Calibration of area based diameter distribution with individual tree based diameter estimates using airborne laser scanning

Diameter distribution is essential for calculating stem volume and timber assortments of forest stands. A new method was proposed in this study to improve the estimation of stem volume and timber assortments, by means of combining the Area-based approach (ABA) and individual tree detection (ITD), the two main approaches to deriving forest attributes from airborne laser scanning (ALS) data. Two methods, replacement, and histogram matching were employed to calibrate ABA-derived diameter distributions with ITD-derived diameter estimates at plot level. The results showed that more accurate estimates were obtained when calibrations were applied. In view of the highest accuracy between ABA and ITD, calibrated diameter distributions decreased its relative RMSE of the estimated entire growing stock, saw log and pulpwood fractions by 2.81%, 3.05% and 7.73% points at best, respectively. Calibration improved pulpwood fraction significantly, which contributed to the negligible bias of the estimated entire growing stock.     

Estimation of aboveground biomass of Robinia pseudoacacia forest in the Yellow River Delta based on UAV and Backpack LiDAR point clouds

Forest plantations are an important source of terrestrial carbon sequestration. The forest of Robinia pseudoacacia in the Yellow River Delta (YRD) is the largest artificial ecological protection forest in China. However, more than half of the forest has appeared different degrees of dieback and even death since the 1990s. Timely and accurate estimation of the forest aboveground biomass (AGB) is a basis for studying the carbon cycle of forests. Light Detecting and Ranging (LiDAR) has been proved to be one of the most powerful methods for forest biomass estimation. However, because of an irregular and overlapping shape of the broadleaved forest canopy in a growing season, it is difficult to segment individual trees and estimate the tree biomass from airborne LiDAR data. In this study, a new method was proposed to solve this problem of individual tree detection in the Robinia pseudoacacia forest based on a combination of the Unmanned Aerial Vehicle-Light Detecting and Ranging (UAV-LiDAR) with the Backpack-LiDAR. The proposed method mainly consists of following steps: (i) at a plot level, trees in the UAV-LiDAR data were detected by seed points obtained by an individual tree segmentation (ITS) method from the Backpack-LiDAR data; (ii) height and diameter at breast height (DBH) of an individual tree would be extracted from UAV and Backpack LiDAR data, respectively; (iii) the individual tree AGB would be calculated through an allometric equation and the forest AGB at the plot level was accumulated; and (iv) the plot-level forest AGB was taken as a dependent variable, and various metrics extracted from UAV-LiDAR point cloud data as independent variables to estimate forest AGB distribution in the study area by using both multiple linear regression (MLR) and random forest (RF) models. The results demonstrate that: (1) the seed points extracted from Backpack-LiDAR could significantly improve the overall accuracy of individual tree detection (F = 0.99), and thus increase the forest AGB estimation accuracy; (2) compared with MLR model, the RF model led to a higher estimation accuracy (p < 0.05); and (3) LiDAR intensity information selected by both MLR and RF models and laser penetration rate (LP) played an important role in estimating healthy forest AGB.     

机载LiDAR数据估算样地和单木尺度森林地上生物量

利用机载激光雷达点云数据,结合大量实测单木结构信息,分别从样地和单木尺度估算了森林地上生物量AGB。首先,利用局部最大值单木提取算法提取了每个样地内的单木结构参数,并针对样地和单木尺度分别计算了一组激光雷达变量。然后,利用激光雷达变量和地上生物量及其两者的对数形式,从样地和单木尺度分别构建了估算模型。最后,针对两种尺度估算过程中存在的不确定性进行了详细讨论。结果表明:(1)样地和单木尺度模型估算的森林地上生物量与地面实测值都具有明显的相关性,且对数模型估算效果要优于非对数模型;(2)样地尺度模型估算效果(R2=0.84,rRMSE=0.23)明显优于单木尺度模型(R2=0.61,rRMSE=0.46);(3)按树木类型分别进行估算可以提高单木地上生物量的估算精度;(4)不论是样地还是单木尺度地上生物量估算都存在一定的不确定性,与样地尺度相比,单木尺度估算过程的不确定性更大,这种不确定性主要来自单木识别过程。     

Evaluation of machine learning models for rice dry biomass estimation and mapping using quad-source optical imagery

ABSTRACT

Several machine learning regression models have been advanced for the estimation of crop biophysical parameters with optical satellite imagery. However, literature on the comparative performances of such models is still limited in range and scope, especially under multiple data sources, despite the potential of multi-source imagery to improving crop monitoring in cloudy areas. To fill in this knowledge gap, this study explored the synergistic use of Landsat-8, Sentinel-2A, China’s environment and disaster monitoring and forecasting satellites (HJ-1 A and B) and Gaofen-1 (GF-1) data to evaluate four machine learning regression models that include Random Forest (RF), Support Vector Machine (SVM), k-Nearest Neighbor (k-NN), and Gradient Boosting Decision Tree (GBDT), for rice dry biomass estimation and mapping. Taking a major rice cultivation area in southeast China as case study during the 2016 and 2017 growing seasons, a cross-calibrated time series of the Enhanced Vegetation Index (EVI) was obtained from the quad-source optical imagery and on which the aforementioned models were applied, respectively. Results indicate that in the before rice heading scenario, the most accurate dry biomass estimates were obtained by the GBDT model (R² of 0.82 and RMSE of 191.8 g/m²) followed by the RF model (R² of 0.79 and RMSE of 197.8 g/m²). After heading, the k-NN model performed best (R² of 0.43 and RMSE of 452.1 g/m²) followed by the RF model (R² of 0.42 and RMSE of 464.7 g/m²). Whist the k-NN model performed least in the before heading scenario, SVM performed least in the after heading scenario. These findings may suggest that machine learning regression models based on an ensemble of decision trees (RF and GBDT) are more suitable for the estimation of rice dry biomass, at least with optical satellite imagery. Studies that would extend the evaluation of these machine learning models, to other parameters like leaf area index, and to microwave imagery, are hereby recommended.     

Spectral behavior of Persian oak under compound stress of water deficit and dust storm

Persian oak (Quercus Brantii Lindl.) which is the most widely distributed tree in the Zagros Mountain forests is affected by western dust storms, mostly originating in Iraq, and harsh water stress as well. The objective of this research is to analyze the spectral behavior of Persian oak under water and dust stress scenarios, aiming to pave the way for modeling the stresses of drought and dust storms on oak trees using remote sensing images. Experiments were carried out on 54 two-year old oak tree seedlings, using a portable wind tunnel in greenhouse conditions. Water stress was induced on seedlings by means of changes in irrigation practices, i.e. well-watered (100 % field capacity), medium water deficit condition (40 % field capacity), and severe water deficit condition (20 % field capacity) treatments. Dust stress is also investigated by using three different dust particle concentrations, i.e. 350, 750 and 1500 (μg/m³). The spectrometry experiments were carried out at leaf and canopy levels in dark room by Fieldspec-3-ASD spectrometer. Spectral analysis was conducted using four procedures: (i) narrow-band spectral indices analysis, (ii) geometric indicators extraction from absorption features, (iii) Partial Least Squares Regression (PLSR), and SVM classifier. Results show that water stress could be modeled much better using PLSR statistic (R² = 0.87, RMSE = 0.12), narrow-band indices analysis (R²_cv = 0.75, RMSE_cv = 0.17), and continuum removal (R² = 0.71, RMSE = 0.20), respectively. For dust stress, PLSR (R² = 0.83, RMSE = 0.14) and narrow-band indices (R² _cv = 0.7, RMSE _cv = 0.30) showed the best results, respectively. SVM could successfully separate stressed and not-stressed samples and also the stress types at both leaf and canopy levels, but it could not distinguish the different levels of stresses.     

Inter-comparison of remote sensing platforms for height estimation of mango and avocado tree crowns

To support the adoption of precision agricultural practices in horticultural tree crops, prior research has investigated the relationship between crop vigour (height, canopy density, health) as measured by remote sensing technologies, to fruit quality, yield and pruning requirements. However, few studies have compared the accuracy of different remote sensing technologies for the estimation of tree height. In this study, we evaluated the accuracy, flexibility, aerial coverage and limitations of five techniques to measure the height of two types of horticultural tree crops, mango and avocado trees. Canopy height estimates from Terrestrial Laser Scanning (TLS) were used as a reference dataset against height estimates from Airborne Laser Scanning (ALS) data, WorldView-3 (WV-3) stereo imagery, Unmanned Aerial Vehicle (UAV) based RGB and multi-spectral imagery, and field measurements. Overall, imagery obtained from the UAV platform were found to provide tree height measurement comparable to that from the TLS (R² = 0.89, RMSE = 0.19 m and rRMSE = 5.37 % for mango trees; R² = 0.81, RMSE = 0.42 m and rRMSE = 4.75 % for avocado trees), although coverage area is limited to 1–10 km² due to battery life and line-of-sight flight regulations. The ALS data also achieved reasonable accuracy for both mango and avocado trees (R² = 0.67, RMSE = 0.24 m and rRMSE = 7.39 % for mango trees; R² = 0.63, RMSE = 0.43 m and rRMSE = 5.04 % for avocado trees), providing both optimal point density and flight altitude, and therefore offers an effective platform for large areas (10 km²–100 km²). However, cost and availability of ALS data is a consideration. WV-3 stereo imagery produced the lowest accuracies for both tree crops (R² = 0.50, RMSE = 0.84 m and rRMSE = 32.64 % for mango trees; R² = 0.45, RMSE = 0.74 m and rRMSE = 8.51 % for avocado trees) when compared to other remote sensing platforms, but may still present a viable option due to cost and commercial availability when large area coverage is required. This research provides industries and growers with valuable information on how to select the most appropriate approach and the optimal parameters for each remote sensing platform to assess canopy height for mango and avocado trees.     

Individual tree crown approach for predicting site index in boreal forests using airborne laser scanning and hyperspectral data

Site productivity is essential information for sustainable forest management and site index (SI) is the most common quantitative measure of it. The SI is usually determined for individual tree species based on tree height and the age of the 100 largest trees per hectare according to stem diameter. The present study aimed to demonstrate and validate a methodology for the determination of SI using remotely sensed data, in particular fused airborne laser scanning (ALS) and airborne hyperspectral data in a forest site in Norway. The applied approach was based on individual tree crown (ITC) delineation: tree species, tree height, diameter at breast height (DBH), and age were modelled and predicted at ITC level using 10-fold cross validation. Four dominant ITCs per 400 m² plot were selected as input to predict SI at plot level for Norway spruce (Picea abies (L.) Karst.) and Scots pine (Pinus sylvestris L.). We applied an experimental setup with different subsets of dominant ITCs with different combinations of attributes (predicted or field-derived) for SI predictions. The results revealed that the selection of the dominant ITCs based on the largest DBH independent of tree species, predicted the SI with similar accuracy as ITCs matched with field-derived dominant trees (RMSE: 27.6% vs 23.3%). The SI accuracies were at the same level when dominant species were determined from the remotely sensed or field data (RMSE: 27.6% vs 27.8%). However, when the predicted tree age was used the SI accuracy decreased compared to field-derived age (RMSE: 27.6% vs 7.6%). In general, SI was overpredicted for both tree species in the mature forest, while there was an underprediction in the young forest. In conclusion, the proposed approach for SI determination based on ITC delineation and a combination of ALS and hyperspectral data is an efficient and stable procedure, which has the potential to predict SI in forest areas at various spatial scales and additionally to improve existing SI maps in Norway.     

Sentinel-2卫星落叶松林龄信息反演

林龄结构信息能够有效反映区域森林群落不同生长阶段的固碳能力,对于评估森林生态系统的健康状况具有重要意义。本研究以中国温带典型优势树种落叶松林为研究对象,分别选择其芽萌动期、展叶期和落叶期时段的Sentinel-2影像,采用多元线性回归（MLR）、随机森林（RF）、支持向量机回归（SVR）、前馈反向传播神经网络（BP）以及多元自适应回归样条（MARS）等5种方法依次构建落叶松林龄反演模型。通过相关性分析首先确定最佳遥感反演物候期,并在此基础上根据相关性差异筛选出5个最优特征变量用于模型反演,分别为冠层含水量（CWC）,归一化水体指数（NDWI）,叶面积指数（LAI）,光合有效辐射吸收率（FAPAR）和植被覆盖度（FVC）。研究结果表明,展叶期为落叶松林最佳遥感反演物候期。除植被衰减指数（PSRI）以及落叶期的NDVI、RVI外,落叶松林龄与各指标之间均呈负相关关系,其中与冠层含水量（CWC）的相关性最高,pearson相关系数达到-0.74（p<0.01）。此外,不同模型反演结果表明,随机森林模型（RF）为最佳落叶松林龄估测模型,其平均决定系数R²和平均均方根误差RMSE分别为0.89和2.91 a;多元线性回归模型（MLR）的林龄估测结果最差,其平均决定系数R²和平均均方根误差RMSE仅为0.57和5.69 a,非线性模型能更好的解释林龄与建模变量之间的关系。     

用地基激光雷达提取单木结构参数——以白皮松为例

以白皮松(Pinus bungeana Zucc)为研究对象,针对地基激光雷达TLS扫描的3维点云数据在单株木垂直方向的分布特征,提出了一种基于体元化方法的树干覆盖度变化检测方法,获取单木枝下高;然后根据获取的枝下高引入2维凸包算法获取垂直方向分层树冠轮廓,并计算树冠体积和冠幅;同时获取的单木参数还有胸径与树高。结果表明:单木枝下高的估测精度较高,R2与RMSE分别为0.97 m和0.21 m;胸径估测结果的R2与RMSE分别为0.79 cm和1.07 cm;采用逐步线性回归方法建立单木树冠体积与其他单木参数的相关关系,模型变量包括冠幅、叶子填充树冠长度和胸径,样本数为20,模型的R2与RMSE分别是0.967 m3和2.64 m3。本文方法能较准确地估测枝下高,TLS数据具有对树冠结构3维建模的潜力。     

Mapping forest aboveground biomass in the reforested Buffelsdraai landfill site using texture combinations computed from SPOT-6 pan-sharpened imagery

Developing models for estimating aboveground biomass (AGB) in naturally growing forests is critical for climate change modelling. AGB models developed using satellite imagery varies with study area, depending on the complexity of vegetation and landscape structure, which affects the upwelling radiance. We assessed the potential of SPOT-6 imagery in predicting AGB of trees planted at different time periods, using image texture combinations. Image texture variables were computed from the SPOT6 pan-sharpened image data, which is characterised by a 1.5 m spatial resolution. In addition, we incorporated the minimal variance technique to select the optimum window sizes that best captures AGB variation in our study area. The results showed that image texture was able to detect AGB for both mature and young trees, however, models detecting mature trees were more superior, with accuracies of R² = 0.70 and 0.25 for 2009–2011 and 2011–2013 plantation phases, respectively. In addition, our results showed that the three band texture ratios yielded the highest accuracy (R² = 0.88 and RMSE = 54.54 kg m⁻²) compared to two texture (R² = 0.85 and RMSE = 60.65 kg m⁻²) and single texture band combinations (R² = 0.64 and RMSE = 94.13 kg m⁻²). A frequency analysis was also run to determine which bands appeared more frequently in the selected texture band models. The frequency analysis revealed that both the red and green bands appeared more frequently on the selected texture band variables, indicating that they were more sensitive to the variation of AGB in our study area. The results showed high variation in AGB within the Buffelsdraai reforestation site, especially due to varying tree plantation phases as well as topography. In essence, the study demonstrated the possibility of image texture combinations computed from the SPOT-6 image in estimating AGB.     

Classifying individual tree species under leaf-off and leaf-on conditions using airborne lidar

In this paper, a methodology for individual tree-based species classification using high sampling density and small footprint lidar data is clarified, corrected and improved. For this purpose, a well-defined directed graph (digraph) is introduced and it plays a fundamental role in the approach. It is argued that there exists one and only one such unique digraph that describes all four pure events and resulting disjoint sets of laser points associated with a single tree in data from a two-return lidar system. However, the digraph is extendable so that it fits an n-return lidar system (n > 2) with higher logical resolution. Furthermore, a mathematical notation for different types of groupings of the laser points is defined, and a new terminology for various types of individual tree-based concepts defined by the digraph is proposed. A novel calibration technique for estimating individual tree heights is evaluated. The approach replaces the unreliable maximum single laser point height of each tree with a more reliable prediction based on shape characteristics of a marginal height distribution of the whole first-return point cloud of each tree. The result shows a reduction of the RMSE of the tree heights of about 20% (stddev = 1.1 m reduced to stddev = 0.92 m). The method improves the species classification accuracy markedly, but it could also be used for reducing the sampling density at the time of data acquisition. Using the calibrated tree heights, a scale-invariant rescaled space for the universal set of points for each tree is defined, in which all individual tree-based geometric measurements are conducted. With the corrected and improved classification methodology the total accuracy raises from 60% to 64% for classifying three leaf-off individual tree deciduous species (N = 200 each) in West Virginia, USA: oaks (Quercus spp.), red maple (Acer rubrum), and yellow poplar (Liriodendron tuliperifera).     

Estimating urban above ground biomass with multi-scale LiDAR

Background

Urban trees have long been valued for providing ecosystem services (mitigation of the “heat island” effect, suppression of air pollution, etc.); more recently the potential of urban forests to store significant above ground biomass (AGB) has also be recognised. However, urban areas pose particular challenges when assessing AGB due to plasticity of tree form, high species diversity as well as heterogeneous and complex land cover. Remote sensing, in particular light detection and ranging (LiDAR), provide a unique opportunity to assess urban AGB by directly measuring tree structure. In this study, terrestrial LiDAR measurements were used to derive new allometry for the London Borough of Camden, that incorporates the wide range of tree structures typical of an urban setting. Using a wall-to-wall airborne LiDAR dataset, individual trees were then identified across the Borough with a new individual tree detection (ITD) method. The new allometry was subsequently applied to the identified trees, generating a Borough-wide estimate of AGB.

Results

Camden has an estimated median AGB density of 51.6 Mg ha^–1 where maximum AGB density is found in pockets of woodland; terrestrial LiDAR-derived AGB estimates suggest these areas are comparable to temperate and tropical forest. Multiple linear regression of terrestrial LiDAR-derived maximum height and projected crown area explained 93% of variance in tree volume, highlighting the utility of these metrics to characterise diverse tree structure. Locally derived allometry provided accurate estimates of tree volume whereas a Borough-wide allometry tended to overestimate AGB in woodland areas. The new ITD method successfully identified individual trees; however, AGB was underestimated by ≤?25% when compared to terrestrial LiDAR, owing to the inability of ITD to resolve crown overlap. A Monte Carlo uncertainty analysis identified assigning wood density values as the largest source of uncertainty when estimating AGB.

Conclusion

Over the coming century global populations are predicted to become increasingly urbanised, leading to an unprecedented expansion of urban land cover. Urban areas will become more important as carbon sinks and effective tools to assess carbon densities in these areas are therefore required. Using multi-scale LiDAR presents an opportunity to achieve this, providing a spatially explicit map of urban forest structure and AGB.

     

Estimating aboveground and organ biomass of plant canopies across the entire season of rice growth with terrestrial laser scanning

Non-destructive and accurate estimation of crop biomass is crucial for the quantitative diagnosis of growth status and timely prediction of grain yield. As an active remote sensing technique, terrestrial laser scanning (TLS) has become increasingly available in crop monitoring for its advantages in recording structural properties. Some researchers have attempted to use TLS data in the estimation of crop aboveground biomass, but only for part of the growing season. Previous studies rarely investigated the estimation of biomass for individual organs, such as the panicles in rice canopies, which led to the poor understanding of TLS technology in monitoring biomass partitioning among organs. The objective of this study was to investigate the potential of TLS in estimating the biomass for individual organs and aboveground biomass of rice and to examine the feasibility of developing universal models for the entire growing season. The field plots experiments were conducted in 2017 and 2018 and involved different nitrogen (N) rates, planting techniques and rice varieties. Three regression approaches, stepwise multiple linear regression (SMLR), random forest regression (RF) and linear mixed-effects (LME) modeling, were evaluated in estimating biomass with extensive TLS and biomass data collected at multiple phenological stages of rice growth across the entire season. The models were calibrated with the 2017 dataset and validated independently with the 2018 dataset.The results demonstrated that growth stage in LME modeling was selected as the most significant random effect on rice growth among the three candidates, which were rice variety, growth stage and planting technique. The LME models grouped by growth stage exhibited higher validation accuracies for all biomass variables over the entire season to varying degrees than SMLR models and RF models. The most pronounced improvement with a LME model was obtained for panicle biomass, with an increase of 0.74 in R² (LME: R² = 0.90, SMLR: R² = 0.16) and a decrease of 1.15 t/ha in RMSE (LME: RMSE =0.79 t/ha, SMLR: RMSE =2.94 t/ha). Compared to SMLR and RF, LME modeling yielded similar estimation accuracies of aboveground biomass for pre-heading stages, but significantly higher accuracies for post-heading stages (LME: R² = 0.63, RMSE =2.27 t/ha; SMLR: R² = 0.42, RMSE =2.42 t/ha; RF: R² = 0.57, RMSE =2.80 t/ha). These findings implied that SMLR was only suitable for the estimation of biomass at pre-heading stages and LME modeling performed remarkably well across all growth stages, especially for post-heading. The results suggest coupling TLS with LME modeling is a promising approach to monitoring rice biomass at post-heading stages at high accuracy and to overcoming the saturation of canopy reflectance signals encountered in optical remote sensing. It also has great potential in the monitoring of other crops in cloud-cover conditions and the instantaneous prediction of grain yield any time before harvest.     

Monitoring aboveground forest biomass dynamics over three decades using Landsat time-series and single-date inventory data

Understanding forest biomass dynamics is crucial for carbon and environmental monitoring, especially in the context of climate change. In this study, we propose a robust approach for monitoring aboveground forest biomass (AGB) dynamics by combining Landsat time-series with single-date inventory data. We developed a Random Forest (RF) based kNN model to produce annual maps of AGB from 1988 to 2017 over 7.2 million ha of forests in Victoria, Australia. The model was internally evaluated using a bootstrapping technique. Predictions of AGB and its change were then independently evaluated using multi-temporal Lidar data (2008 and 2016). To understand how natural and anthropogenic processes impact forest AGB, we analysed trends in relation to the history of disturbance and recovery. Specifically, change metrics (e.g., AGB loss and gain, Years to Recovery - Y2R) were calculated at the pixel level to characterise the patterns of AGB change resulting from forest dynamics. The imputation model achieved a RMSE value of 132.9 Mg ha⁻¹ (RMSE% = 46.3%) and R² value of 0.56. Independent assessments of prediction maps in 2008 and 2016 using Lidar-based AGB data achieved relatively high accuracies, with a RMSE of 108.6 Mg ha⁻¹ and 135.9 Mg ha⁻¹ for 2008 and 2016, respectively. Annual validations of AGB maps using un-changed, homogenous Lidar plots suggest that our model is transferable through time (RMSE ranging from 109.65 Mg ha⁻¹ to 112.27 Mg ha⁻¹ and RMSE% ranging from 25.38% to 25.99%). In addition, changes in AGB values associated with forest disturbance and recovery (decrease and increase, respectively) were captured by predicted maps. AGB change metrics indicate that AGB loss and Y2R varied across bioregions and were highly dependent on levels of disturbance severity (i.e., a greater loss and longer recovery time were associated with a higher severity disturbance). On average, high severity fire burnt from 200 Mg ha⁻¹ to 550 Mg ha⁻¹ of AGB and required up to 15 years to recover while clear-fell logging caused a reduction in 250 Mg ha⁻¹ to 600 Mg ha⁻¹ of AGB and required nearly 20 years to recover. In addition, AGB within un-disturbed forests showed statistically significant but monotonic trends, suggesting a mild gradual drop over time across most bioregions. Our methods are designed to support forest managers and researchers in developing forest monitoring systems, especially in developing regions, where only a single date forestry inventory exists.     

Estimation of higher chlorophylla concentrations using field spectral measurement and HJ-1A hyperspectral satellite data in Dianshan Lake,China

Based on in situ water sampling and field spectral measurements in Dianshan Lake, a semi-analytical three-band algorithm was used to estimate Chlorophylla (Chla) content in case II waters. The three bands selected to estimate Chla for high concentrations included 653, 691 and 748 nm. An equation, based on the difference in reciprocal reflectance between 653 and 691 nm, multiplied by reflectance at 748 nm as [R_rs⁻¹(653) − R_rs⁻¹ (691)] R_rs(748), explained 85.57% of variance in Chla concentration with a root mean square error (RMSE) of <6.56 mg/m³. In order to test the utility of this model with satellite data, HJ-1A Hyperspectral Imager (HSI) data were analyzed using comparable wavelengths selected from the in situ data [B₆₇⁻¹(656) − B₈₀⁻¹(716)] B₈₇(753). This model accounted for 84.3% of Chla variation, estimating Chla concentrations with an RMSE of <4.23 mg/m³. The results illustrate that, based on the determined wavelengths, the spectrum-based model can achieve a high estimation accuracy and can be applied to hyperspectral satellite imagery especially for higher Chla concentration waters.     

Assessing the sensitivity of site index models developed using bi-temporal airborne laser scanning data to different top height estimates and grid cell sizes

Site productivity and forest growth are critical inputs into projecting wood volume and biomass accumulation over time. Site productivity, which is determined most commonly using site index models is also the primary criterion to consider many forest management decisions. Most of the previous research utilizing the remote sensing data for assessment of site index with forest height are based on the existing site index models developed with traditional dendrometric methods. However, these traditional methods are both time-consuming and expensive. This study demonstrates how bi-temporal airborne laser scanning (ALS) data collected within the 8-year period can be used for the development of site index models for Scots pine. The accuracy of ALS-derived models was assessed by comparison to the reference site index model developed based on data from stem analysis of 174 felled Scots pine trees. We evaluated the effect of different height metrics and grid cell size on the trajectory of site index models developed from ALS-derived measurements. Four methods of estimating top height from ALS point clouds were evaluated: 95th, 99th and 100th percentiles of point clouds and an individual tree detection approach (ITD). The models were created for a range of grid cell sizes: 10 × 10 m, 30 × 30 m, and 50 × 50 m. The results indicate that bitemporal ALS data could substitute traditional methods that have been applied to date for stand growth modelling. It was found that top height increment can be estimated by using both ITD approach and the 100th percentile of point cloud giving an appropriate top height (TH) increment estimation. Observed growth curves of reference trees agreed best with the trajectories that were obtained based on TH calculated using ITD method (R² = 0.892) and 100th percentile (R² = 0.797). In case of TH obtained from 99th and 95th percentiles only weak correlation was found: R² = 0.358 and R² = 0.213, accordingly. The height growth models developed with 95th and 99th percentiles of point cloud were not compatible with the reference model. We also found that grid cell size did not affect the model height growth trajectories. Irrespective of the grid cell size, the obtained model trajectories for the given method of TH estimation are nearly identical for cells 10 × 10, 30 × 30 and 50 × 50 m.     

Integration of WorldView-2 and airborne LiDAR data for tree species level carbon stock mapping in Kayar Khola watershed,Nepal

Integration of WorldView-2 satellite image with small footprint airborne LiDAR data for estimation of tree carbon at species level has been investigated in tropical forests of Nepal. This research aims to quantify and map carbon stock for dominant tree species in Chitwan district of central Nepal. Object based image analysis and supervised nearest neighbor classification methods were deployed for tree canopy retrieval and species level classification respectively. Initially, six dominant tree species (Shorea robusta, Schima wallichii, Lagerstroemia parviflora, Terminalia tomentosa, Mallotus philippinensis and Semecarpus anacardium) were able to be identified and mapped through image classification. The result showed a 76% accuracy of segmentation and 1970.99 as best average separability. Tree canopy height model (CHM) was extracted based on LiDAR’s first and last return from an entire study area. On average, a significant correlation coefficient (r) between canopy projection area (CPA) and carbon; height and carbon; and CPA and height were obtained as 0.73, 0.76 and 0.63, respectively for correctly detected trees. Carbon stock model validation results showed regression models being able to explain up to 94%, 78%, 76%, 84% and 78% of variations in carbon estimation for the following tree species: S. robusta, L. parviflora, T. tomentosa, S. wallichii and others (combination of rest tree species).     

Global, 30-m resolution continuous fields of tree cover: Landsat-based rescaling of MODIS vegetation continuous fields with lidar-based estimates of error

Abstract

We developed a global, 30-m resolution dataset of percent tree cover by rescaling the 250-m MOderate-resolution Imaging Spectroradiometer (MODIS) Vegetation Continuous Fields (VCF) Tree Cover layer using circa- 2000 and 2005 Landsat images, incorporating the MODIS Cropland Layer to improve accuracy in agricultural areas. Resulting Landsat-based estimates maintained consistency with the MODIS VCF in both epochs (RMSE =8.6% in 2000 and 11.9% in 2005), but showed improved accuracy in agricultural areas and increased discrimination of small forest patches. Against lidar measurements, the Landsat-based estimates exhibited accuracy slightly less than that of the MODIS VCF (RMSE=16.8% for MODIS-based vs. 17.4% for Landsat-based estimates), but RMSE of Landsat estimates was 3.3 percentage points lower than that of the MODIS data in an agricultural region. The Landsat data retained the saturation artifact of the MODIS VCF at greater than or equal to 80% tree cover but showed greater potential for removal of errors through calibration to lidar, with post-calibration RMSE of 9.4% compared to 13.5% in MODIS estimates. Provided for free download at the Global Land Cover Facility (GLCF) website (www.landcover.org), the 30-m resolution GLCF tree cover dataset is the highest-resolution multi-temporal depiction of Earth's tree cover available to the Earth science community.     

Non-linear partial least square regression increases the estimation accuracy of grass nitrogen and phosphorus using in situ hyperspectral and environmental data

Grass nitrogen (N) and phosphorus (P) concentrations are direct indicators of rangeland quality and provide imperative information for sound management of wildlife and livestock. It is challenging to estimate grass N and P concentrations using remote sensing in the savanna ecosystems. These areas are diverse and heterogeneous in soil and plant moisture, soil nutrients, grazing pressures, and human activities. The objective of the study is to test the performance of non-linear partial least squares regression (PLSR) for predicting grass N and P concentrations through integrating in situ hyperspectral remote sensing and environmental variables (climatic, edaphic and topographic). Data were collected along a land use gradient in the greater Kruger National Park region. The data consisted of: (i) in situ-measured hyperspectral spectra, (ii) environmental variables and measured grass N and P concentrations. The hyperspectral variables included published starch, N and protein spectral absorption features, red edge position, narrow-band indices such as simple ratio (SR) and normalized difference vegetation index (NDVI). The results of the non-linear PLSR were compared to those of conventional linear PLSR. Using non-linear PLSR, integrating in situ hyperspectral and environmental variables yielded the highest grass N and P estimation accuracy (R² = 0.81, root mean square error (RMSE) = 0.08, and R² = 0.80, RMSE = 0.03, respectively) as compared to using remote sensing variables only, and conventional PLSR. The study demonstrates the importance of an integrated modeling approach for estimating grass quality which is a crucial effort towards effective management and planning of protected and communal savanna ecosystems.     

Terrestrial laser scanning to estimate plot-level forest canopy fuel properties

This paper evaluates the potential of a terrestrial laser scanner (TLS) to characterize forest canopy fuel characteristics at plot level. Several canopy properties, namely canopy height, canopy cover, canopy base height and fuel strata gap were estimated. Different approaches were tested to avoid the effect of canopy shadowing on canopy height estimation caused by deployment of the TLS below the canopy. Estimation of canopy height using a grid approach provided a coefficient of determination of R² = 0.81 and an RMSE of 2.47 m. A similar RMSE was obtained using the 99th percentile of the height distribution of the highest points, representing the 1% of the data, although the coefficient of determination was lower (R² = 0.70). Canopy cover (CC) was estimated as a function of the occupied cells of a grid superimposed upon the TLS point clouds. It was found that CC estimates were dependent on the cell size selected, with 3 cm being the optimum resolution for this study. The effect of the zenith view angle on CC estimates was also analyzed. A simple method was developed to estimate canopy base height from the vegetation vertical profiles derived from an occupied/non-occupied voxels approach. Canopy base height was estimated with an RMSE of 3.09 m and an R² = 0.86. Terrestrial laser scanning also provides a unique opportunity to estimate the fuel strata gap (FSG), which has not been previously derived from remotely sensed data. The FSG was also derived from the vegetation vertical profile with an RMSE of 1.53 m and an R² = 0.87.