Extracting LiDAR indices to characterise multilayered forest structure using mixture distribution functions

Discrete Light Detection and Ranging (LiDAR) data is used to stratify a multilayered eucalyptus forest and characterise the structure of the vertical profile. We present a methodology that may prove useful for a very broad range of forest management applications, particularly for timber inventory evaluation and forest growth modelling. In this study, we use LiDAR data to stratify a multilayered eucalyptus forest and characterise the structure of specific vegetation layers for forest hydrology research, as vegetation dynamics influence a catchment's streamflow yield. A forest stand's crown height, density, depth, and closure, influence aerodynamic properties of the forest structure and the amount of transpiring leaf area, which in turn determine evapotranspiration rates. We present a methodology that produces canopy profile indices of understorey and overstorey vegetation using mixture models with a wide range of theoretical distribution functions. Mixture models provide a mechanism to summarise complex canopy attributes into a short list of parameters that can be empirically analysed against stand characteristics.Few studies have explored theoretical distribution functions to represent the vertical profile of vegetation structure in LiDAR data. All prior studies have focused on a Weibull distribution function, which is unimodal. In a complex native forest ecosystem, the form of the distribution of LiDAR points may be highly variable between forest types and age classes. We compared 44 probability distributions within a two component mixture model to determine the most suitable bimodal distributions for representing LiDAR density estimates of Mountain Ash forests in south-eastern Australia. An elimination procedure identified eleven candidate distributions for representing the eucalyptus component of the mixture model.We demonstrate the methodology on a sample of plots to predict overstorey stand volumes and basal area, and understorey basal area of 18-, 37-, and 70-year old Mountain Ash forest with variable density classes. The 70-year old forest has been subjected to a range of treatments including: thinning of the eucalyptus layer with two distinct retention rates, removal of the understorey, and clear felling of patches that have 37 year old regenerating forest. We demonstrate that the methodology has clear potential, as observed versus predicted values of eucalyptus basal area and stand volume were highly correlated, with bootstrap based r² ranging from 0.61 to 0.89 and 0.67 to 0.88 respectively. Non-eucalyptus basal area r² ranged from 0.5 to 0.91.     

Quantifying forest above ground carbon content using LiDAR remote sensing

The UNFCCC and interest in the source of the missing terrestrial carbon sink are prompting research and development into methods for carbon accounting in forest ecosystems. Here we present a canopy height quantile-based approach for quantifying above ground carbon content (AGCC) in a temperate deciduous woodland, by means of a discrete-return, small-footprint airborne LiDAR. Fieldwork was conducted in Monks Wood National Nature Reserve UK to estimate the AGCC of five stands from forest mensuration and allometric relations. In parallel, a digital canopy height model (DCHM) and a digital terrain model (DTM) were derived from elevation measurements obtained by means of an Optech Airborne Laser Terrain Mapper 1210. A quantile-based approach was adopted to select a representative statistic of height distributions per plot. A forestry yield model was selected as a basis to estimate stemwood volume per plot from these heights metrics. Agreement of r=0.74 at the plot level was achieved between ground-based AGCC estimates and those derived from the DCHM. Using a 20×20 m grids superposed to the DCHM, the AGCC was estimated at the stand level and at the woodland level. At the stand level, the agreement between the plot data upscaled in proportion to area and the LiDAR estimates was r=0.85. At the woodland level, LiDAR estimates were nearly 24% lower than those from the upscaled plot data. This suggests that field-based approaches alone may not be adequate for carbon accounting in heterogeneous forests. Conversely, the LiDAR 20×20 m grid approach has an enhanced capability of monitoring the natural variability of AGCC across the woodland.     

Forest canopy recovery from the 1938 hurricane and subsequent salvage damage measured with airborne LiDAR

The structure of a forest canopy often reflects its disturbance history. Such signatures of past disturbances or legacies can influence how the ecosystem functions across broad spatio-temporal scales. The 1938 hurricane and ensuing salvage operations which swept through New England represent the most recent large, infrequent disturbance (LID) in this region. Though devastating (downing ∼ 70% of the timber at Harvard Forest), the disturbance was not indiscriminate; it left behind a heterogeneous landscape comprised of different levels of canopy damage. We analyzed large-footprint LiDAR, from the Prospect Hill tract at Harvard Forest in central Massachusetts, to assess whether damage to the forest structure from the hurricane and subsequent timber extraction could be discerned after ∼ 65 years. Differences in LiDAR-derived measures of canopy height and vertical diversity were a function of the degree of damage from the 1938 hurricane and the predominant tree species which is, in part, a function of land use history. Higher levels of damage corresponded to slightly shorter canopies with a less even vertical distribution of return from the ground to the top. In addition, differences in canopy topography as revealed by spatial autocorrelation of canopy top heights were found among the damage classes. Less disturbed stands were characterized by lower levels of local autocorrelation for canopy height and higher levels of vertical diversity of LiDAR returns. These differences in canopy structure reveal that the forest tract has not completely recovered from the 1938 LID and salvage regime, which may have implications on arboreal and understory habitat and other ecosystem functions.     

Strengths and limitations of assessing forest density and spatial configuration with aerial LiDAR

Changes in the structural state of forests of the semi-arid U.S.A., such as an increase in tree density, are widely believed to be leading to an ecological crisis, but accurate methods of quantifying forest density and configuration are lacking at landscape scales. An individual tree canopy (ITC) method based on aerial LiDAR has been developed to assess forest structure by estimating the density and spatial configuration of trees in four different height classes. The method has been tested against field measured forest inventory data from two geographically distinct forests with independent LiDAR acquisitions. The results show two distinct patterns: accurate, unbiased density estimates for trees taller than 20 m, and underestimation of density in trees less than 20 m tall. The underestimation of smaller trees is suggested to be a limitation of LiDAR remote sensing. Ecological applications of the method are demonstrated through landscape metrics analysis of density and configuration rasters.     

Fusion of airborne LiDAR and satellite multispectral data for the estimation of timber volume in the Southern Alps

Remote sensing can be considered a key instrument for studies related to forests and their dynamics. At present, the increasing availability of multisensor acquisitions over the same areas, offers the possibility to combine data from different sensors (e.g., optical, RADAR, LiDAR). This paper presents an analysis on the fusion of airborne LiDAR and satellite multispectral data (IRS 1C LISS III), for the prediction of forest stem volume at plot level in a complex mountain area (Province of Trento, Southern Italian Alps), characterized by different tree species, complex morphology (i.e. altitude ranges from 65 m to 3700 m above sea level), and a range of different climates (from the sub-Mediterranean to Alpine type). 799 sample plots were randomly distributed over the 3000 km² of the forested areas of the Trento Province. From each plot, a set of variables were extracted from both LiDAR and multispectral data. A regression analysis was carried out considering two data sources (LiDAR and multispectral) and their combination, and dividing the plot areas into groups according to their species composition, altitude and slope. Experimental results show that the combination of LiDAR and IRS 1C LISS III data, for the estimation of stem volume, is effective in all the experiments considered. The best developed models comprise variables extracted from both of these data sources. The RMSE% on an independent validation set for the stem volume estimation models ranges between 17.2% and 26.5%, considering macro sets of tree species (deciduous, evergreen and mixed), between 17.5% and 29.0%, considering dominant species plots, and between 15.5% and 21.3% considering altitude and slope sets.     

Development of a simulation model to predict LiDAR interception in forested environments

Airborne scanning LiDAR systems are used to predict a range of forest attributes. However, the accuracy with which this can be achieved is highly dependent on the sensor configuration and the structural characteristics of the forest examined. As a result, there is a need to understand laser light interactions with forest canopies so that LiDAR sensor configurations can be optimised to assess particular forest types. Such optimisation will not only ensure the targeted forest attributes can be accurately and consistently quantified, but may also minimise the cost of data acquisition and indicate when a survey configuration will not deliver information needs.In this paper, we detail the development and application of a model to simulate laser interactions within forested environments. The developed model, known as the LiDAR Interception and Tree Environment (LITE) model, utilises a range of structural configurations to simulate trees with variable heights, crown dimensions and foliage clumping. We developed and validated the LITE model using field data obtained from three forested sites covering a range of structural classes. Model simulations were then compared to coincident airborne LiDAR data collected over the same sites. Results indicate that the LITE model can be used to produce comparable estimates of maximum height of trees within plots (differences < 2.42 m), mean heights of first return data (differences < 2.27 m), and canopy height percentiles (r² = 0.94, p < 0.001) when compared to airborne LiDAR data. In addition, the distribution of airborne LiDAR hits through the canopy profile was closely matched by model predictions across the range of sites. Importantly, this demonstrates that the structural differences between forest stands can be characterised by LITE. Models that are capable of interpreting the response of small-footprint LiDAR waveforms can facilitate algorithm development, the generation of corrections for actual LiDAR data, and the optimisation of sensor configurations for differing forest types, benefiting a range of experimental and commercial LiDAR applications. As a result, we also performed a scenario analysis to demonstrate how differences in forest structure, terrain, and sensor configuration can influence the interception of LiDAR beams.     

Assessing the utility of airborne hyperspectral and LiDAR data for species distribution mapping in the coastal Pacific Northwest, Canada

To effectively manage forested ecosystems an accurate characterization of species distribution is required. In this study we assess the utility of hyperspectral Airborne Imaging Spectrometer for Applications (AISA) imagery and small footprint discrete return Light Detection and Ranging (LiDAR) data for mapping 11 tree species in and around the Gulf Islands National Park Reserve, in coastal South-western Canada. Using hyperspectral imagery yielded producer's and user's accuracies for most species ranging from > 52-95.4 and > 63-87.8%, respectively. For species dominated by definable growth stages, pixel-level fusion of hyperspectral imagery with LiDAR-derived height and volumetric canopy profile data increased both producer's (+ 5.1-11.6%) and user's (+ 8.4-18.8%) accuracies. McNemar's tests confirmed that improvements in overall accuracies associated with the inclusion of LiDAR-derived structural information were statistically significant (p < 0.05). This methodology establishes a specific framework for mapping key species with greater detail and accuracy then is possible using conventional approaches (i.e., aerial photograph interpretation), or either technology on its own. Furthermore, in the study area, acquisition and processing costs were lower than a conventional aerial photograph interpretation campaign, making hyperspectral/LiDAR fusion a viable replacement technology.     

Fusion of LiDAR and imagery for estimating forest canopy fuels

Due to increased fuel loading as a result of fire suppression, land managers in the American west are in need of precise information about the fuels they manage, including canopy fuels. Canopy fuel metrics such as canopy height (CH), canopy base height (CBH), canopy bulk density (CBD) and available canopy fuel (ACF) are specific inputs for wildfire behavior models such as FARSITE and emission models such as FOFEM. With finer spatial resolution data, accurate quantification of these metrics with detailed spatial heterogeneity can be accomplished. Light Detection and Ranging (LiDAR) and color near-infrared imagery are active and passive systems, respectively, that have been utilized for measuring a range of forest structure characteristics at high resolution. The objective of this research was to determine which remote sensing dataset can estimate canopy fuels more accurately and whether a fusion of these datasets produces more accurate estimates. Regression models were developed for ponderosa pine (Pinus ponderosa) stand representative of eastern Washington State using field data collected in the Ahtanum State Forest and metrics derived from LiDAR and imagery. Strong relationships were found with LiDAR alone and LiDAR was found to increase canopy fuel accuracy compared to imagery. Fusing LiDAR with imagery and/or LiDAR intensity led to small increases in estimation accuracy over LiDAR alone. By improving the ability to estimate canopy fuels at higher spatial resolutions, spatially explicit fuel layers can be created and used in wildfire behavior and smoke emission models leading to more accurate estimations of crown fire risk and smoke related emissions.     

Modeling forest songbird species richness using LiDAR-derived measures of forest structure

Conservation of biodiversity requires information at many spatial scales in order to detect and preserve habitat for many species, often simultaneously. Vegetation structure information is particularly important for avian habitat models and has largely been unavailable for large areas at the desired resolution. Airborne LiDAR, with its combination of relatively broad coverage and fine resolution provides existing new opportunities to map vegetation structure and hence avian habitat. Our goal was to model the richness of forest songbirds using forest structure information obtained from LiDAR data. In deciduous forests of southern Wisconsin, USA, we used discrete-return airborne LiDAR to derive forest structure metrics related to the height and density of vegetation returns, as well as composite variables that captured major forest structural elements. We conducted point counts to determine total forest songbird richness and the richness of foraging, nesting, and forest edge-related habitat guilds. A suite of 35 LiDAR variables were used to model bird species richness using best-subsets regression and we used hierarchical partitioning analysis to quantify the explanatory power of each variable in the multivariate models. Songbird species richness was correlated most strongly with LiDAR variables related to canopy and midstory height and midstory density (R² = 0.204, p < 0.001). Richness of species that nest in the midstory was best explained by canopy height variables (R² = 0.197, p < 0.001). Species that forage on the ground responded to mean canopy height and the height of the lower canopy (R² = 0.149, p < 0.005) while aerial foragers had higher richness where the canopy was tall and dense and the midstory more sparse (R² = 0.216, p < 0.001). Richness of edge-preferring species was greater where there were fewer vegetation returns but higher density in the understory (R² = 0.153, p < 0.005). Forest interior specialists responded positively to a tall canopy, developed midstory, and a higher proportion of vegetation returns (R² = 0.195, p < 0.001). LiDAR forest structure metrics explained between 15 and 20% of the variability in richness within deciduous forest songbird communities. This variability was associated with vertical structure alone and shows how LiDAR can provide a source of complementary predictive data that can be incorporated in models of wildlife habitat associations across broad geographical extents.     

Characterizing forest ecological structure using pulse types and heights of airborne laser scanning

Characterizing forest structure is an important part of any comprehensive biodiversity assessment. However, current methods for measuring structural complexity require a laborious process that involves many logistically expensive point based measurements. An automated or semi-automated method would be ideal. In this study, the utility of airborne laser scanning (LiDAR; Light Detection and Ranging) for characterizing the ecological structure of a forest landscape is examined. The innovation of this paper is to use different laser pulse return properties from a full waveform LiDAR to characterize forest ecological structure. First the LiDAR dataset is stratified into four vertical layers: ground, low vegetation (0-1 m from the ground), medium vegetation (1-5 m from the ground) and high vegetation (> 5 m). Subsequently the “Type” of LiDAR return is analysed: Type 1 (singular returns); Type 2 (first of many returns); Type 3 (intermediate returns); and Type 4 (last of many returns). A forest characterization scheme derived from LiDAR point clouds is proposed. A validation of the scheme is then presented using a network of field sites that recorded commonly used metrics of biodiversity. The proposed forest characterization categories allow for quantification of gaps (above bare ground, low vegetation and medium vegetation), canopy cover and its vertical density as well as the presence of various canopy strata (low, medium and high). Regression analysis showed that LiDAR derived variables were good predictors of field recorded variables (R² = 0.82, P < 0.05 between LiDAR derived presence of low vegetation and field derived LAI for low vegetation). The proposed scheme clearly shows the potential of full waveform LiDAR to provide information on the complexity of habitat structure.     

Evaluating uncertainty in mapping forest carbon with airborne LiDAR

Airborne LiDAR is increasingly used to map carbon stocks in tropical forests, but our understanding of mapping errors is constrained by the spatial resolution (i.e., plot size) used to calibrate LiDAR with field data (typically 0.1-0.36 ha). Reported LiDAR errors range from 17 to 40 Mg C ha^− 1, but should be lower at coarser resolutions because relative errors are expected to scale with (plot area)^-1/2. We tested this prediction empirically using a 50-ha plot with mapped trees, allowing an assessment of LiDAR prediction errors at multiple spatial resolutions. We found that errors scaled approximately as expected, declining by 38% (compared to 40% predicted from theory) from 0.36- to 1-ha resolution. We further reduced errors at all spatial resolutions by accounting for tree crowns that are bisected by plot edges (not typically done in forestry), and collectively show that airborne LiDAR can map carbon stocks with 10% error at 1-ha resolution — a level comparable to the use of field plots alone.     

Simulated impact of sample plot size and co-registration error on the accuracy and uncertainty of LiDAR-derived estimates of forest stand biomass

Regression has been widely applied in Light Detection And Ranging (LiDAR) remote sensing to spatially extend predictions of total aboveground biomass (TAGB) and other biophysical properties over large forested areas. Sample (field) plot size has long been considered a key sampling design parameter and focal point for optimization in forest surveys, because of its impact on sampling effort and the estimation accuracy of forest inventory attributes. In this study, we demonstrate how plot size and co-registration error interact to influence the estimation of LiDAR canopy height and density metrics, regression model coefficients, and the prediction accuracy of least-squares estimators of TAGB. We made use of simulated forest canopies and synthetic LiDAR point clouds, so that we could maintain strict control over the spatial scale and complexity of forest scenes, as well as the magnitude and type of planimetric error inherent in ground-reference and LiDAR datasets. Our results showed that predictions of TAGB improved markedly as plot size increased from 314 (10 m radius) to 1964 m² (25 m radius). The co-registration error (spatial overlap) between ground-reference and LiDAR samples negatively impacted the estimation of LiDAR metrics, regression model fit, and the prediction accuracy of TAGB. We found that larger plots maintained a higher degree of spatial overlap between ground-reference and LiDAR datasets for any given GPS error, and were therefore more resilient to the ill effects of co-registration error compared to small plots. The impact of co-registration error was more pronounced in tall, spatially heterogeneous stands than short, homogeneous stands. We identify and briefly discuss three possible ways that LiDAR data could be used to optimize plot size, sample selection, and the deployment of GPS resources in forest biomass surveys.     

Nearest neighbor imputation of species-level, plot-scale forest structure attributes from LiDAR data

Meaningful relationships between forest structure attributes measured in representative field plots on the ground and remotely sensed data measured comprehensively across the same forested landscape facilitate the production of maps of forest attributes such as basal area (BA) and tree density (TD). Because imputation methods can efficiently predict multiple response variables simultaneously, they may be usefully applied to map several structural attributes at the species-level. We compared several approaches for imputing the response variables BA and TD, aggregated at the plot-scale and species-level, from topographic and canopy structure predictor variables derived from discrete-return airborne LiDAR data. The predictor and response variables were associated using imputation techniques based on normalized and unnormalized Euclidean distance, Mahalanobis distance, Independent Component Analysis (ICA), Canonical Correlation Analysis (aka Most Similar Neighbor, or MSN), Canonical Correspondence Analysis (aka Gradient Nearest Neighbor, or GNN), and Random Forest (RF). To compare and evaluate these approaches, we computed a scaled Root Mean Square Distance (RMSD) between observed and imputed plot-level BA and TD for 11 conifer species sampled in north-central Idaho. We found that RF produced the best results overall, especially after reducing the number of response variables to the most important species in each plot with regard to BA and TD. We concluded that RF was the most robust and flexible among the imputation methods we tested. We also concluded that canopy structure and topographic metrics derived from LiDAR surveys can be very useful for species-level imputation.     

The effects of rectification and Global Positioning System errors on satellite image-based estimates of forest area

Satellite image-based maps of forest attributes are of considerable interest and are used for multiple purposes such as international reporting by countries that have no national forest inventory and small area estimation for all countries. Construction of the maps typically entails, in part, rectifying the satellite images to a geographic coordinate system, observing ground plots whose coordinates are obtained from Global Positioning System (GPS) receivers that are calibrated to the same geographic coordinate system, and then matching ground plots to image pixels containing the centers of the ground plots. Errors in rectification and GPS coordinates cause observations of ground attributes to be associated with spectral values of incorrect pixels which, in turn, introduces classification errors into the resulting maps. The most important finding of the study is that for common magnitudes of rectification and GPS errors, as many as half the ground plots may be assigned to incorrect pixels. The effects on areal estimates obtained by aggregating class predictions for individual pixels are deviation of the estimates from their true values, erroneous confidence intervals, and incorrect inferences. Results are reported in detail for both probability-based (design-based) and model-based approaches to inference for proportion forest area using maps constructed from Landsat imagery, forest inventory plot observations and a logistic regression model.     

Integrating airborne multispectral imagery and airborne LiDAR data for enhanced lithological mapping in vegetated terrain

Practical and financial constraints associated with traditional field-based lithological mapping are often responsible for the generation of maps with insufficient detail and inaccurately located contacts. In arid areas with well exposed rocks and soils, high-resolution multi- and hyperspectral imagery is a valuable mapping aid as lithological units can be readily discriminated and mapped by automatically matching image pixel spectra to a set of reference spectra. However, the use of spectral imagery in all but the most barren terrain is problematic because just small amounts of vegetation cover can obscure or mask the spectra of underlying geological substrates. The use of ancillary information may help to improve lithological discrimination, especially where geobotanical relationships are absent or where distinct lithologies exhibit inherent spectral similarity. This study assesses the efficacy of airborne multispectral imagery for detailed lithological mapping in a vegetated section of the Troodos ophiolite (Cyprus), and investigates whether the mapping performance can be enhanced through the integration of LiDAR-derived topographic data. In each case, a number of algorithms involving different combinations of input variables and classification routine were employed to maximise the mapping performance. Despite the potential problems posed by vegetation cover, geobotanical associations aided the generation of a lithological map - with a satisfactory overall accuracy of 65.5% and Kappa of 0.54 - using only spectral information. Moreover, owing to the correlation between topography and lithology in the study area, the integration of LiDAR-derived topographic variables led to significant improvements of up to 22.5% in the overall mapping accuracy compared to spectral-only approaches. The improvements were found to be considerably greater for algorithms involving classification with an artificial neural network (the Kohonen Self-Organizing Map) than the parametric Maximum Likelihood Classifier. The results of this study demonstrate the enhanced capability of data integration for detailed lithological mapping in areas where spectral discrimination is complicated by the presence of vegetation or inherent spectral similarities.     

Multivariate forest structure modelling and mapping using high resolution airborne imagery and topographic information

Remote sensing has been widely used for modelling and mapping individual forest structural attributes, such as LAI and stem density, however the development and evaluation of methods for simultaneously modelling and mapping multivariate aspects of forest structure are comparatively limited. Multivariate representation of forest structure can be used as a means to infer other environmental attributes such as biodiversity and habitat, which have often been shown to be enhanced in more structurally diverse or complex forests. Image-based modelling of multivariate forest structure is useful in developing an understanding of the associations between different aspects of vertical and horizontal structure and image characteristics. Models can also be applied spatially to all image pixels to produce maps of multivariate forest structure as an alternative to sample-based field assessment. This research used high spatial resolution multispectral airborne imagery to provide spectral, spatial, and object-based information in the development of a model of multivariate forest structure as represented by twenty-four field variables measured in plots within a temperate hardwood forest in southern Quebec, Canada. Redundancy Analysis (RDA) was used to develop a model that explained a statistically significant proportion of the variance of these structural attributes. The resulting model included image variables representing mostly within-crown and within-shadow brightness variance (texture) as well as elevation, taken from a DEM of the study area. It was applied spatially across the entire study area to produce a continuous map of predicted multivariate forest structure. Bootstrapping validation of the model provided an RMSE of 19.9%, while independent field validation of mapped areas of complex and simple structure showed accuracies of 89% and 69%, respectively. Multiscale testing using resampled imagery suggested that the methods could possibly be used with current pan-sharpened, or future sub-metre resolution, multispectral satellite imagery, which would provide much greater spatial coverage and reduced image processing compared to airborne imagery.     

Estimation of above- and below-ground biomass across regions of the boreal forest zone using airborne laser

Regression models relating variables derived from airborne laser scanning (ALS) to above-ground and below-ground biomass were estimated for 1395 sample plots in young and mature coniferous forest located in ten different areas within the boreal forest zone of Norway. The sample plots were measured as part of large-scale operational forest inventories. Four different ALS instruments were used and point density varied from 0.7 to 1.2 m^− 2. One variable related to canopy height and one related to canopy density were used as independent variables in the regressions. The statistical effects of area and age class were assessed by including dummy variables in the models. Tree species composition was treated as continuous variables. The proportion of explained variability was 88% for above- and 85% for below-ground biomass models. For given combinations of ALS-derived variables, the differences between the areas were up to 32% for above-ground biomass and 38% for below-ground biomass. The proportion of spruce had a significant impact on both the estimated models. The proportion of broadleaves had a significant effect on above-ground biomass only, while the effect of age class was significant only in the below-ground biomass model. Because of local effects on the biomass-ALS data relationships, it is indicated by this study that sample plots distributed over the entire area would be needed when using ALS for regional or national biomass monitoring.     

Non-parametric and parametric methods using satellite images for estimating growing stock volume in alpine and Mediterranean forest ecosystems

This paper describes applications of non-parametric and parametric methods for estimating forest growing stock volume using Landsat images on the basis of data measured in the field, integrated with ancillary information. Several k-Nearest Neighbors (k-NN) algorithm configurations were tested in two study areas in Italy belonging to Mediterranean and Alpine ecosystems. Field data were acquired by the regional forest inventory and forest management plans, and satellite images are from Landsat 5 TM and Landsat 7 ETM+. The paper describes the data used, the methodologies adopted and the results achieved in terms of pixel level accuracy of forest growing stock volume estimates. The results show that several factors affect estimation accuracy when using the k-NN method. For the two test areas a total of 3500 different configurations of the k-NN algorithm were systematically tested by changing the number and type of spectral and ancillary input variables, type of multidimensional distance measures, number of nearest neighbors and methods for spectral feature extraction using the leave-one-out (LOO) procedure. The best k-NN configurations were then used for pixel level estimation; the accuracy was estimated with a bootstrapping procedure; and the results were compared to estimates obtained using parametric regression methods implemented on the same data set.

The best k-NN growing stock volume pixel level estimates in the Alpine area have a Root Mean Square Error (RMSE) ranging between 74 and 96 m³ ha^− 1 (respectively, 22% and 28% of the mean measured value) and between 106 and 135 m³ ha^− 1 (respectively, 44% and 63% of the mean measured value) in the Mediterranean area. On the whole, the results cast a promising light on the use of non-parametric techniques for forest attribute estimation and mapping with accuracy high enough to support forest planning activities in such complex landscapes. The results of the LOO analyses also highlight the importance of a local empirical optimization phase of the k-NN procedure before defining the best algorithm configuration. In the tests performed the pixel level accuracy increased, depending on the k-NN configuration, as much as 100%.     

Lithological mapping of the Troodos ophiolite, Cyprus, using airborne LiDAR topographic data

Traditional field-based lithological mapping can be a time-consuming, costly and challenging endeavour when large areas need to be investigated, where terrain is remote and difficult to access and where the geology is highly variable over short distances. Consequently, rock units are often mapped at coarse-scales, resulting in lithological maps that have generalised contacts which in many cases are inaccurately located. Remote sensing data, such as aerial photographs and satellite imagery are commonly incorporated into geological mapping programmes to obtain geological information that is best revealed by overhead perspectives. However, spatial and spectral limitations of the imagery and dense vegetation cover can limit the utility of traditional remote sensing products. The advent of Airborne Light Detection And Ranging (LiDAR) as a remote sensing tool offers the potential to provide a novel solution to these problems because accurate and high-resolution topographic data can be acquired in either forested or non-forested terrain, allowing discrimination of individual rock types that typically have distinct topographic characteristics. This study assesses the efficacy of airborne LiDAR as a tool for detailed lithological mapping in the upper section of the Troodos ophiolite, Cyprus. Morphometric variables (including slope, curvature and surface roughness) were derived from a 4 m digital terrain model in order to quantify the topographic characteristics of four principal lithologies found in the area. An artificial neural network (the Kohonen Self-Organizing Map) was then employed to classify the lithological units based upon these variables. The algorithm presented here was used to generate a detailed lithological map which defines lithological contacts much more accurately than the best existing geological map. In addition, a separate map of classification uncertainty highlights potential follow-up targets for ground-based verification. The results of this study demonstrate the significant potential of airborne LiDAR for lithological discrimination and rapid generation of detailed lithological maps, as a contribution to conventional geological mapping programmes.     

Effects of different sensors, flying altitudes, and pulse repetition frequencies on forest canopy metrics and biophysical stand properties derived from small-footprint airborne laser data

Canopy height distributions were created from small-footprint airborne laser scanner (ALS) data collected over 40 field sample plots with size 1000 m² located in mature conifer forest. ALS data were collected with two different instruments, i.e., the ALTM 1233 and ALTM 3100 laser scanners (Optech Inc.). The ALTM 1233 data were acquired at a flying altitude of 1200 m and a pulse repetition frequency (PRF) of 33 kHz. Three different acquisitions were carried out with ALTM 3100, i.e., (1) a flying altitude of 1100 m and a PRF of 50 kHz, (2) a flying altitude of 1100 m and a PRF of 100 kHz, and (3) a flying altitude of 2000 m and a PRF of 50 kHz. Height percentiles, mean and maximum height values, coefficients of variation of the heights, and canopy density at different height intervals above the ground were derived from the four different ALS datasets and for single + first and last echoes of the ALS data separately. The ALS-derived height- and density variables were assessed in pair-wise comparisons to evaluate the effects of (a) instrument, (b) flying altitude, and (c) PRF. A systematic shift in height values of up to 0.3 m between sensors when the first echoes were compared was demonstrated. Also the density-related variables differed significantly between the two instruments. Comparisons of flying altitudes and PRFs revealed upwards shifted canopy height distributions for the highest flying altitude (2000 m) and the lowest PRF (50 kHz). The distribution of echoes on different echo categories, i.e., single and multiple (first and last) echoes, differed significantly between acquisitions. The proportion of multiple echoes decreased with increasing flying altitude and PRF. Different echo categories have different properties since it is likely that single echoes tend to occur in the densest parts of the tree crowns, i.e., near the apex where the concentration of biological matter is highest and distance to the ground is largest. To assess the influence of instrument, flying altitude, and PRF on biophysical properties derived from ALS data, regression analysis was carried out to relate ALS-derived metrics to mean tree height (h_L) and timber volume (V). Cross validation revealed only minor differences in precision for the different ALS acquisitions, but systematic differences between acquisitions of up to 2.5% for h_L and 10.7% for V were found when comparing data from different acquisitions.