Sub-footprint analysis to uncover tree height variation using ICESat/GLAS

Detailed forest height data are an indispensable prerequisite for many forestry and earth science applications. Existing research of using Geoscience Laser Altimeter System (GLAS) data mainly focuses on deriving average or maximum tree heights within a GLAS footprint, i.e. an ellipse with a diameter of 65 m. However, in most forests, it is likely that the tree heights within such ellipse are heterogeneous. Therefore, it is desired to uncover detailed tree height variation within a GLAS footprint. To the best of our knowledge, no such methods have been reported as of now. In this study, we aim to characterize tree heights’ variation within a GLAS footprint as different layers, each of which corresponds to trees with similar heights. As such, we developed a new method that embraces two steps: first, a refined Levenberg–Marquardt (LM) algorithm is proposed to decompose raw GLAS waveform into multiple Gaussian signals, within which it is hypothesized that each vegetation signal corresponds to a particular tree height layer. Second, for each layer, three parameters were first defined: Canopy Top Height (CTH), Crown Length (CL), and Cover Proportion (CP). Then we extracted the three parameters from each Gaussian signal through a defined model. In order to test our developed method, we set up a study site in Ejina, China where the dominant specie is Populus euphratica. Both simulated and field tree height data were adopted. With regard to the simulation data, results presented a very high agreement for the three predefined parameters between our results and simulation data. When our methods were applied to the field data, the respective R² become 0.78 (CTH), CL (R² = 0.76), CP (R² = 0.74). Overall, our studies revealed that large footprint GLAS waveform data have the potentials for obtaining detailed forest height variation.     

Above ground biomass estimation in an African tropical forest with lidar and hyperspectral data

The estimation of above ground biomass in forests is critical for carbon cycle modeling and climate change mitigation programs. Small footprint lidar provides accurate biomass estimates, but its application in tropical forests has been limited, particularly in Africa. Hyperspectral data record canopy spectral information that is potentially related to forest biomass. To assess lidar ability to retrieve biomass in an African forest and the usefulness of including hyperspectral information, we modeled biomass using small footprint lidar metrics as well as airborne hyperspectral bands and derived vegetation indexes. Partial Least Square Regression (PLSR) was adopted to cope with multiple inputs and multicollinearity issues; the Variable of Importance in the Projection was calculated to evaluate importance of individual predictors for biomass. Our findings showed that the integration of hyperspectral bands (R² = 0.70) improved the model based on lidar alone (R² = 0.64), this encouraging result call for additional research to clarify the possible role of hyperspectral data in tropical regions. Replacing the hyperspectral bands with vegetation indexes resulted in a smaller improvement (R² = 0.67). Hyperspectral bands had limited predictive power (R² = 0.36) when used alone. This analysis proves the efficiency of using PLSR with small-footprint lidar and high resolution hyperspectral data in tropical forests for biomass estimation. Results also suggest that high quality ground truth data is crucial for lidar-based AGB estimates in tropical African forests, especially if airborne lidar is used as an intermediate step of upscaling field-measured AGB to a larger area.     

PTrees: A point-based approach to forest tree extraction from lidar data

This paper introduces PTrees, a multi-scale dynamic point cloud segmentation dedicated to forest tree extraction from lidar point clouds. The method process the point data using the raw elevation values (Z) and compute height (H = Z − ground elevation) during post-processing using an innovative procedure allowing to preserve the geometry of crown points. Multiple segmentations are done at different scales. Segmentation criteria are then applied to dynamically select the best set of apices from the tree segments extracted at the various scales. The selected set of apices is then used to generate a final segmentation. PTrees has been tested in 3 different forest types, allowing to detect 82% of the trees with under 10% of false detection rate. Future development will integrate crown profile estimation during the segmentation process in order to both maximize the detection of suppressed trees and minimize false detections.     

3-D voxel-based solid modeling of a broad-leaved tree for accurate volume estimation using portable scanning lidar

We developed a method to produce a 3-D voxel-based solid model of a tree based on portable scanning lidar data for accurate estimation of the volume of the woody material. First, we obtained lidar measurements with a high laser pulse density from several measurement positions around the target, a Japanese zelkova tree. Next, we converted lidar-derived point-cloud data for the target into voxels. The voxel size was 0.5 cm × 0.5 cm × 0.5 cm. Then, we used differences in the spatial distribution of voxels to separate the stem and large branches (diameter > 1 cm) from small branches (diameter ⩽ 1 cm). We classified the voxels into sets corresponding to the stem and to each large branch and then interpolated voxels to fill out their surfaces and their interiors. We then merged the stem and large branches with the small branches. The resultant solid model of the entire tree was composed of consecutive voxels that filled the outer surface and the interior of the stem and large branches, and a cloud of voxels equivalent to small branches that were discretely scattered in mainly the upper part of the target. Using this model, we estimated the woody material volume by counting the number of voxels in each part and multiplying the number of voxels by the unit voxel volume (0.13 cm³). The percentage error of the volume of the stem and part of a large branch was 0.5%. The estimation error of a certain part of the small branches was 34.0%.     

Assessing plantation canopy condition from airborne imagery using spectral mixture analysis and fractional abundances

Pine plantations in Australia are subject to a range of abiotic and biotic damaging agents that affect tree health and productivity. In order to optimise management decisions, plantation managers require regular intelligence relating to the status and trends in the health and condition of trees within individual compartments. Remote sensing technology offers an alternative to traditional ground-based assessment of these plantations. Automated estimation of foliar crown health, especially in degraded crowns, can be difficult due to mixed pixels when there is low or fragmented vegetation cover. In this study we apply a linear spectral unmixing approach to high spatial resolution (50 cm) multispectral imagery to quantify the fractional abundances of the key image endmembers: sunlit canopy, shadow, and soil. A number of Pinus radiata tree crown attributes were modelled using multiple linear regression and endmember fraction images. We found high levels of significance (r² = 0.80) for the overall crown colour and colour of the crown leader (r² = 0.79) in tree crowns affected by the fungal pathogen Sphaeropsis sapinea, which produces both needle necrosis and chlorosis. Results for stands associated with defoliation and chlorosis through infestation by the aphid Essigella californica were lower with an r² = 0.33 for crown transparency and r² = 0.31 for proportion of crown affected. Similar analysis of data from a nitrogen deficient site produced an outcome somewhat in between the other two damaging agents. Overall the sunlit canopy image fraction has been the most important variable used in the modelling of forest condition for all damaging agents.     

Detecting forest damage after a low-severity fire using remote sensing at multiple scales

Remote sensing technologies are an ideal platform to examine the extent and impact of fire on the landscape. In this study we assess that capacity of the RapidEye constellation and Landsat (Thematic Mapper and Operational Land Imager to map fine-scale burn attributes for a small, low severity prescribed fire in a dry Western Canadian forest. Estimates of burn severity from field data were collated into a simple burn index and correlated with a selected suite of common spectral vegetation indices. Burn severity classes were then derived to map fire impacts and estimate consumed woody surface fuels (diameter ≥2.6 cm). All correlations between the simple burn index and vegetation indices produced significant results (p < 0.01), but varied substantially in their overall accuracy. Although the Landsat Soil Adjusted Vegetation Index provided the best regression fit (R² = 0.56), results suggested that RapidEye provided much more spatially detailed estimates of tree damage (Soil Adjusted Vegetation Index, R² = 0.51). Consumption estimates of woody surface fuels ranged from 3.38 ± 1.03 Mg ha⁻¹ to 11.73 ± 1.84 Mg ha⁻¹, across four derived severity classes with uncertainties likely a result of changing foliage moisture between the before and after fire images. While not containing spectral information in the short wave infrared, the spatial variability provided by the RapidEye imagery has potential for mapping and monitoring fine scale forest attributes, as well as the potential to resolve fire damage at the individual tree level.     

A new 500-m resolution map of canopy height for Amazon forest using spaceborne LiDAR and cloud-free MODIS imagery

In the present study, we aimed to map canopy heights in the Brazilian Amazon mainly on the basis of spaceborne LiDAR and cloud-free MODIS imagery with a new method (the Self-Organizing Relationships method) for spatial modeling of the LiDAR footprint. To evaluate the general versatility, we compared the created canopy height map with two different canopy height estimates on the basis of our original field study plots (799 plots located in eight study sites) and a previously developed canopy height map. The compared canopy height estimates were obtained by: (1) a stem diameter at breast height (D) – tree height (H) relationship specific to each site on the basis of our original field study, (2) a previously developed D–H model involving environmental and structural factors as explanatory variables (Feldpausch et al., 2011), and (3) a previously developed canopy height map derived from the spaceborne LiDAR data with different spatial modeling method and explanatory variables (Simard et al., 2011). As a result, our canopy height map successfully detected a spatial distribution pattern in canopy height estimates based on our original field study data (r = 0.845, p = 8.31 × 10⁻³) though our canopy height map showed a poor correlation (r = 0.563, p = 0.146) with the canopy height estimate based on a previously developed model by Feldpausch et al. (2011). We also confirmed that the created canopy height map showed a similar pattern with the previously developed canopy height map by Simard et al. (2011). It was concluded that the use of the spaceborne LiDAR data provides a sufficient accuracy in estimating the canopy height at regional scale.     

Quantifying post-fire fallen trees using multi-temporal lidar

Massive tree-felling due to root damage is a common fire effect on burnt areas in Scandinavia, but has so far not been analyzed in detail. Here we explore if pre- and post-fire lidar data can be used to estimate the proportion of fallen trees. The study was carried out within a large (14,000 ha) area in central Sweden burnt in August 2014, where we had access to airborne lidar data from both 2011 and 2015. Three data-sets of predictor variables were tested: POST (post-fire lidar metrics), DIF (difference between post- and pre-fire lidar metrics) and combination of those two (POST_DIF). Fractional logistic regression was used to predict the proportion of fallen trees. Training data consisted of 61 plots, where the number of fallen and standing trees was calculated both in the field and with interpretation of drone images. The accuracy of the best model was tested based on 100 randomly selected validation plots with a size of 25 × 25 m.Our results showed that multi-temporal lidar together with field-collected training data can be used for quantifying post-fire tree felling over large areas. Several height-, density- and intensity metrics correlated with the proportion of fallen trees. The best model combined metrics from both datasets (POST_DIF), resulting in a RMSE of 0.11. Results were slightly poorer in the validation plots with RMSE of 0.18 using pixel size of 12.5 m and RMSE of 0.15 using pixel size of 6.25 m. Our model performed least well for stands that had been exposed to high-intensity crown fire. This was likely due to the low amount of echoes from the standing black tree skeletons. Wall-to-wall maps produced with this model can be used for landscape level analysis of fire effects and to explore the relationship between fallen trees and forest structure, soil type, fire intensity or topography.     

Object-based class modelling for multi-scale riparian forest habitat mapping

Object-based class modelling allows for mapping complex, hierarchical habitat systems. The riparian zone, including forests, represents such a complex ecosystem. Forests within riparian zones are biologically high productive and characterized by a rich biodiversity; thus considered of high community interest with an imperative to be protected and regularly monitored. Satellite earth observation (EO) provides tools for capturing the current state of forest habitats such as forest composition including intermixture of non-native tree species. Here we present a semi-automated object based image analysis (OBIA) approach for the mapping of riparian forests by applying class modelling of habitats based on the European Nature Information System (EUNIS) habitat classifications and the European Habitats Directive (HabDir) Annex 1. A very high resolution (VHR) WorldView-2 satellite image provided the required spatial and spectral details for a multi-scale image segmentation and rule-base composition to generate a six-level hierarchical representation of riparian forest habitats. Thereby habitats were hierarchically represented within an image object hierarchy as forest stands, stands of homogenous tree species and single trees represented by sunlit tree crowns. 522 EUNIS level 3 (EUNIS-3) habitat patches with a mean patch size (MPS) of 12,349.64 m² were modelled from 938 forest stand patches (MPS = 6868.20 m²) and 43,742 tree stand patches (MPS = 140.79 m²). The delineation quality of the modelled EUNIS-3 habitats (focal level) was quantitatively assessed to an expert-based visual interpretation showing a mean deviation of 11.71%.     

Above ground biomass and tree species richness estimation with airborne lidar in tropical Ghana forests

Estimates of forest aboveground biomass are fundamental for carbon monitoring and accounting; delivering information at very high spatial resolution is especially valuable for local management, conservation and selective logging purposes. In tropical areas, hosting large biomass and biodiversity resources which are often threatened by unsustainable anthropogenic pressures, frequent forest resources monitoring is needed. Lidar is a powerful tool to estimate aboveground biomass at fine resolution; however its application in tropical forests has been limited, with high variability in the accuracy of results. Lidar pulses scan the forest vertical profile, and can provide structure information which is also linked to biodiversity. In the last decade the remote sensing of biodiversity has received great attention, but few studies focused on the use of lidar for assessing tree species richness in tropical forests.This research aims at estimating aboveground biomass and tree species richness using discrete return airborne lidar in Ghana forests. We tested an advanced statistical technique, Multivariate Adaptive Regression Splines (MARS), which does not require assumptions on data distribution or on the relationships between variables, being suitable for studying ecological variables.We compared the MARS regression results with those obtained by multilinear regression and found that both algorithms were effective, but MARS provided higher accuracy either for biomass (R² = 0.72) and species richness (R² = 0.64). We also noted strong correlation between biodiversity and biomass field values. Even if the forest areas under analysis are limited in extent and represent peculiar ecosystems, the preliminary indications produced by our study suggest that instrument such as lidar, specifically useful for pinpointing forest structure, can also be exploited as a support for tree species richness assessment.     

Spectral resampling based on user-defined inter-band correlation filter: C3 and C4 grass species classification

In this paper, a user-defined inter-band correlation filter function was used to resample hyperspectral data and thereby mitigate the problem of multicollinearity in classification analysis. The proposed resampling technique convolves the spectral dependence information between a chosen band-centre and its shorter and longer wavelength neighbours. Weighting threshold of inter-band correlation (WTC, Pearson's r) was calculated, whereby r = 1 at the band-centre. Various WTC (r = 0.99, r = 0.95 and r = 0.90) were assessed, and bands with coefficients beyond a chosen threshold were assigned r = 0. The resultant data were used in the random forest analysis to classify in situ C₃ and C₄ grass canopy reflectance. The respective WTC datasets yielded improved classification accuracies (kappa = 0.82, 0.79 and 0.76) with less correlated wavebands when compared to resampled Hyperion bands (kappa = 0.76). Overall, the results obtained from this study suggested that resampling of hyperspectral data should account for the spectral dependence information to improve overall classification accuracy as well as reducing the problem of multicollinearity.     

Assessing the utility WorldView-2 imagery for tree species mapping in South African subtropical humid forest and the conservation implications: Dukuduku forest patch as case study

Indigenous forest biome in South Africa is highly fragmented into patches of various sizes (most patches < 1 km²). The utilization of timber and non-timber resources by poor rural communities living around protected forest patches produce subtle changes in the forest canopy which can be hardly detected on a timely manner using traditional field surveys. The aims of this study were to assess: (i) the utility of very high resolution (VHR) remote sensing imagery (WorldView-2, 0.5–2 m spatial resolution) for mapping tree species and canopy gaps in one of the protected subtropical coastal forests in South Africa (the Dukuduku forest patch (ca.3200 ha) located in the province of KwaZulu-Natal) and (ii) the implications of the map products to forest conservation. Three dominant canopy tree species namely, Albizia adianthifolia, Strychnos spp. and Acacia spp., and canopy gap types including bushes (grass/shrubby), bare soil and burnt patches were accurately mapped (overall accuracy = 89.3 ± 2.1%) using WorldView-2 image and support vector machine classifier. The maps revealed subtle forest disturbances such as bush encroachment and edge effects resulting from forest fragmentation by roads and a power-line. In two stakeholders’ workshops organised to assess the implications of the map products to conservation, participants generally agreed amongst others implications that the VHR maps provide valuable information that could be used for implementing and monitoring the effects of rehabilitation measures. The use of VHR imagery is recommended for timely inventorying and monitoring of the small and fragile patches of subtropical forests in Southern Africa.     

Assessing the potential of multi-seasonal WorldView-2 imagery for mapping West African agroforestry tree species

High resolution satellite systems enable efficient and detailed mapping of tree cover, with high potential to support both natural resource monitoring and ecological research. This study investigates the capability of multi-seasonal WorldView-2 imagery to map five dominant tree species at the individual tree crown level in a parkland landscape in central Burkina Faso. The Random Forest algorithm is used for object based tree species classification and for assessing the relative importance of WorldView-2 predictors. The classification accuracies from using wet season, dry season and multi-seasonal datasets are compared to gain insights about the optimal timing for image acquisition. The multi-seasonal dataset produced the most accurate classifications, with an overall accuracy (OA) of 83.4%. For classifications based on single date imagery, the dry season (OA = 78.4%) proved to be more suitable than the wet season (OA = 68.1%). The predictors that contributed most to the classification success were based on the red edge band and visible wavelengths, in particular green and yellow. It was therefore concluded that WorldView-2, with its unique band configuration, represents a suitable data source for tree species mapping in West African parklands. These results are particularly promising when considering the recently launched WorldView-3, which provides data both at higher spatial and spectral resolution, including shortwave infrared bands.     

Camera derived vegetation greenness index as proxy for gross primary production in a low Arctic wetland area

The Arctic is experiencing disproportionate warming relative to the global average, and the Arctic ecosystems are as a result undergoing considerable changes. Continued monitoring of ecosystem productivity and phenology across temporal and spatial scales is a central part of assessing the magnitude of these changes. This study investigates the ability to use automatic digital camera images (DCIs) as proxy data for gross primary production (GPP) in a complex low Arctic wetland site. Vegetation greenness computed from DCIs was found to correlate significantly (R² = 0.62, p < 0.001) with a normalized difference vegetation index (NDVI) product derived from the WorldView-2 satellite. An object-based classification based on a bi-temporal image composite was used to classify the study area into heath, copse, fen, and bedrock. Temporal evolution of vegetation greenness was evaluated and modeled with double sigmoid functions for each plant community. GPP at light saturation modeled from eddy covariance (EC) flux measurements were found to correlate significantly with vegetation greenness for all plant communities in the studied year (i.e., 2010), and the highest correlation was found between modeled fen greenness and GPP (R² = 0.85, p < 0.001). Finally, greenness computed within modeled EC footprints were used to evaluate the influence of individual plant communities on the flux measurements. The study concludes that digital cameras may be used as a cost-effective proxy for potential GPP in remote Arctic regions.     

Retrieval of effective leaf area index (LAIe) and leaf area density (LAD) profile at individual tree level using high density multi-return airborne LiDAR

As an important canopy structure indicator, leaf area index (LAI) proved to be of considerable implications for forest ecosystem and ecological studies, and efficient techniques for accurate LAI acquisitions have long been highlighted. Airborne light detection and ranging (LiDAR), often termed as airborne laser scanning (ALS), once was extensively investigated for this task but showed limited performance due to its low sampling density. Now, ALS systems exhibit more competing capacities such as high density and multi-return sampling, and hence, people began to ask the questions like—“can ALS now work better on the task of LAI prediction?” As a re-examination, this study investigated the feasibility of LAI retrievals at the individual tree level based on high density and multi-return ALS, by directly considering the vertical distributions of laser points lying within each tree crown instead of by proposing feature variables such as quantiles involving laser point distribution modes at the plot level. The examination was operated in the case of four tree species (i.e. Picea abies, Pinus sylvestris, Populus tremula and Quercus robur) in a mixed forest, with their LAI-related reference data collected by using static terrestrial laser scanning (TLS). In light of the differences between ALS- and TLS-based LAI characterizations, the methods of voxelization of 3D scattered laser points, effective LAI (LAIe) that does not distinguish branches from canopies and unified cumulative LAI (ucLAI) that is often used to characterize the vertical profiles of crown leaf area densities (LADs) was used; then, the relationships between the ALS- and TLS-derived LAIes were determined, and so did ucLAIs. Tests indicated that the tree-level LAIes for the four tree species can be estimated based on the used airborne LiDAR (R² = 0.07, 0.26, 0.43 and 0.21, respectively) and their ucLAIs can also be derived. Overall, this study has validated the usage of the contemporary high density multi-return airborne LiDARs for LAIe and LAD profile retrievals at the individual tree level, and the contribution are of high potential for advancing forest ecosystem modeling and ecological understanding.     

Mapping land cover gradients through analysis of hyper-temporal NDVI imagery

The green cover of the earth exhibits various spatial gradients that represent gradual changes in space of vegetation density and/or in species composition. To date, land cover mapping methods differentiate at best, mapping units with different cover densities and/or species compositions, but typically fail to express such differences as gradients. Present interpretation techniques still make insufficient use of freely available spatial-temporal Earth Observation (EO) data that allow detection of existing land cover gradients. This study explores the use of hyper-temporal NDVI imagery to detect and delineate land cover gradients analyzing the temporal behavior of NDVI values. MODIS-Terra MVC-images (250 m, 16-day) of Crete, Greece, from February 2000 to July 2009 are used. The analysis approach uses an ISODATA unsupervised classification in combination with a Hierarchical Clustering Analysis (HCA). Clustering of class-specific temporal NDVI profiles through HCA resulted in the identification of gradients in landcover vegetation growth patterns. The detected gradients were arranged in a relational diagram, and mapped. Three groups of NDVI-classes were evaluated by correlating their class-specific annual average NDVI values with the field data (tree, shrub, grass, bare soil, stone, litter fraction covers). Multiple regression analysis showed that within each NDVI group, the fraction cover data were linearly related with the NDVI data, while NDVI groups were significantly different with respect to tree cover (adj. R² = 0.96), shrub cover (adj. R² = 0.83), grass cover (adj. R² = 0.71), bare soil (adj. R² = 0.88), stone cover (adj. R² = 0.83) and litter cover (adj. R² = 0.69) fractions. Similarly, the mean Sorenson dissimilarity values were found high and significant at confidence interval of 95% in all pairs of three NDVI groups. The study demonstrates that hyper-temporal NDVI imagery can successfully detect and map land cover gradients. The results may improve land cover assessment and aid in agricultural and ecological studies.     

Texture augmented detection of macrophyte species using decision trees

Image classification using multispectral sensors has shown good performance in detecting macrophytes at the species level. However, species level classification often does not utilize the texture information provided by high resolution images. This study investigated whether image texture provides useful vector(s) for the discrimination of monospecific stands of three floating macrophyte species in Quickbird imagery of the South Nation River. Semivariograms indicated that window sizes of 5 × 5 and 13 × 13 pixels were the most appropriate spatial scales for calculation of the grey level co-occurrence matrix and subsequent texture attributes from the multispectral and panchromatic bands. Of the 214 investigated vectors (13 Haralick texture attributes * 15 bands + 9 spectral bands + 10 transformations/indices), feature selection determined which combination of spectral and textural vectors had the greatest class separability based on the Mann–Whitney U-test and Jefferies–Matusita distance. While multispectral red and near infrared (NIR) performed satisfactorily, the addition of panchromatic-dissimilarity slightly improved class separability and the accuracy of a decision tree classifier (Kappa: red/NIR/panchromatic-dissimilarity – 93.2% versus red/NIR – 90.4%). Class separability improved by incorporating a second texture attribute, but resulted in a decrease in classification accuracy. The results suggest that incorporating image texture may be beneficial for separating stands with high spatial heterogeneity. However, the benefits may be limited and must be weighed against the increased complexity of the classifier.     

Modeling biophysical properties of broad-leaved stands in the hyrcanian forests of Iran using fused airborne laser scanner data and ultraCam-D images

Inventories of mixed broad-leaved forests of Iran mainly rely on terrestrial measurements. Due to rapid changes and disturbances and great complexity of the silvicultural systems of these multilayer forests, frequent repetition of conventional ground-based plot surveys is often cost prohibitive. Airborne laser scanning (ALS) and multispectral data offer an alternative or supplement to conventional inventories in the Hyrcanian forests of Iran. In this study, the capability of a combination of ALS and UltraCam-D data to model stand volume, tree density, and basal area using random forest (RF) algorithm was evaluated. Systematic sampling was applied to collect field plot data on a 150 m × 200 m sampling grid within a 1100 ha study area located at 36°38′- 36°42′N and 54°24′–54°25′E. A total of 308 circular plots (0.1 ha) were measured for calculation of stand volume, tree density, and basal area per hectare. For each plot, a set of variables was extracted from both ALS and multispectral data. The RF algorithm was used for modeling of the biophysical properties using ALS and UltraCam-D data separately and combined. The results showed that combining the ALS data and UltraCam-D images provided a slight increase in prediction accuracy compared to separate modeling. The RMSE as percentage of the mean, the mean difference between observed and predicted values, and standard deviation of the differences using a combination of ALS data and UltraCam-D images in an independent validation at 0.1-ha plot level were 31.7%, 1.1%, and 84 m³ ha⁻¹ for stand volume; 27.2%, 0.86%, and 6.5 m² ha⁻¹ for basal area, and 35.8%, −4.6%, and 77.9 n ha⁻¹ for tree density, respectively. Based on the results, we conclude that fusion of ALS and UltraCam-D data may be useful for modeling of stand volume, basal area, and tree density and thus gain insights into structural characteristics in the complex Hyrcanian forests.     

A sun–crown–sensor model and adapted C-correction logic for topographic correction of high resolution forest imagery

Canopy shadowing mediated by topography is an important source of radiometric distortion on remote sensing images of rugged terrain. Topographic correction based on the sun–canopy–sensor (SCS) model significantly improved over those based on the sun–terrain–sensor (STS) model for surfaces with high forest canopy cover, because the SCS model considers and preserves the geotropic nature of trees. The SCS model accounts for sub-pixel canopy shadowing effects and normalizes the sunlit canopy area within a pixel. However, it does not account for mutual shadowing between neighboring pixels. Pixel-to-pixel shadowing is especially apparent for fine resolution satellite images in which individual tree crowns are resolved. This paper proposes a new topographic correction model: the sun–crown–sensor (SCnS) model based on high-resolution satellite imagery (IKONOS) and high-precision LiDAR digital elevation model. An improvement on the C-correction logic with a radiance partitioning method to address the effects of diffuse irradiance is also introduced (SCnS + C). In addition, we incorporate a weighting variable, based on pixel shadow fraction, on the direct and diffuse radiance portions to enhance the retrieval of at-sensor radiance and reflectance of highly shadowed tree pixels and form another variety of SCnS model (SCnS + W). Model evaluation with IKONOS test data showed that the new SCnS model outperformed the STS and SCS models in quantifying the correlation between terrain-regulated illumination factor and at-sensor radiance. Our adapted C-correction logic based on the sun–crown–sensor geometry and radiance partitioning better represented the general additive effects of diffuse radiation than C parameters derived from the STS or SCS models. The weighting factor W_t also significantly enhanced correction results by reducing within-class standard deviation and balancing the mean pixel radiance between sunlit and shaded slopes. We analyzed these improvements with model comparison on the red and near infrared bands. The advantages of SCnS + C and SCnS + W on both bands are expected to facilitate forest classification and change detection applications.     

Aboveground total and green biomass of dryland shrub derived from terrestrial laser scanning

Sagebrush (Artemisia tridentata), a dominant shrub species in the sagebrush-steppe ecosystem of the western US, is declining from its historical distribution due to feedbacks between climate and land use change, fire, and invasive species. Quantifying aboveground biomass of sagebrush is important for assessing carbon storage and monitoring the presence and distribution of this rapidly changing dryland ecosystem. Models of shrub canopy volume, derived from terrestrial laser scanning (TLS) point clouds, were used to accurately estimate aboveground sagebrush biomass. Ninety-one sagebrush plants were scanned and sampled across three study sites in the Great Basin, USA. Half of the plants were scanned and destructively sampled in the spring (n = 46), while the other half were scanned again in the fall before destructive sampling (n = 45). The latter set of sagebrush plants was scanned during both spring and fall to further test the ability of the TLS to quantify seasonal changes in green biomass. Sagebrush biomass was estimated using both a voxel and a 3-D convex hull approach applied to TLS point cloud data. The 3-D convex hull model estimated total and green biomass more accurately (R² = 0.92 and R² = 0.83, respectively) than the voxel-based method (R² = 0.86 and R² = 0.73, respectively). Seasonal differences in TLS-predicted green biomass were detected at two of the sites (p < 0.001 and p = 0.029), elucidating the amount of ephemeral leaf loss in the face of summer drought. The methods presented herein are directly transferable to other dryland shrubs, and implementation of the convex hull model with similar sagebrush species is straightforward.