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 r2 ranging from 0.61 to 0.89 and 0.67 to 0.88 respectively. Non-eucalyptus basal area r2 ranged from 0.5 to 0.91.

[1]  E. Næsset.  Estimating timber volume of forest stands using airborne laser scanner data , 1997 .

[2]  G. Golub,et al.  Good Ridge Parameter , 1979 .

[3]  W. Cohen,et al.  Geographic variability in lidar predictions of forest stand structure in the Pacific Northwest , 2005 .

[4]  J. Gove,et al.  A Finite Mixture Model for Characterizing the Diameter Distributions of Mixed-Species Forest Stands , 2002, Forest Science.

[5]  M. Hutchinson.  A new procedure for gridding elevation and stream line data with automatic removal of spurious pits , 1989 .

[6]  R. Vertessy,et al.  Factors determining relations between stand age and catchment water balance in mountain ash forests , 2001 .

[7]  J. Monteith.  Evaporation and environment. , 1965, Symposia of the Society for Experimental Biology.

[8]  L. Zhang,et al.  Fitting irregular diameter distributions of forest stands by Weibull, modified Weibull, and mixture Weibull models , 2006, Journal of Forest Research.

[9]  Geoffrey J. McLachlan,et al.  Finite Mixture Models , 2019, Annual Review of Statistics and Its Application.

[10]  J. Hyyppä,et al.  Estimation of timber volume and stem density based on scanning laser altimetry and expected tree size distribution functions , 2004 .

[11]  Juha Hyyppä,et al.  Identifying and quantifying structural characteristics of heterogeneous boreal forests using laser scanner data , 2005 .

[12]  Scott D. Roberts,et al.  Measuring heights to crown base and crown median with LiDAR in a mature, even-aged loblolly pine stand , 2009 .

[13]  Francesco Sepic,et al.  Automatic detection of dominated vegetation under canopy using Airborne Laser Scanning data , 2008 .

[14]  D. Donoghue,et al.  Remote sensing of species mixtures in conifer plantations using LiDAR height and intensity data , 2007 .

[15]  Alan Brown,et al.  Criteria and Indicators for Sustainable Forest Management , 2001 .

[16]  Åsa Persson,et al.  Identifying species of individual trees using airborne laser scanner , 2004 .

[17]  R. Tibshirani,et al.  Improvements on Cross-Validation: The 632+ Bootstrap Method , 1997 .

[18]  S. P. Lloyd,et al.  Least squares quantization in PCM , 1982, IEEE Trans. Inf. Theory.

[19]  S. Magnussen,et al.  Derivations of stand heights from airborne laser scanner data with canopy-based quantile estimators , 1998 .

[20]  Patrick D. Gerard,et al.  Characterizing vertical forest structure using small-footprint airborne LiDAR , 2003 .

[21]  K. Kraus,et al.  Determination of terrain models in wooded areas with airborne laser scanner data , 1998 .

[22]  I. Burke,et al.  Estimating stand structure using discrete-return lidar: an example from low density, fire prone ponderosa pine forests , 2005 .

[23]  S. Ustin,et al.  Modeling airborne laser scanning data for the spatial generation of critical forest parameters in fire behavior modeling , 2003 .

[24]  Peter Axelsson,et al.  Processing of laser scanner data-algorithms and applications , 1999 .

[25]  J. Gove,et al.  A finite mixture of two Weibull distributions for modeling the diameter distributions of rotated-sigmoid, uneven-aged stands , 2001 .

[26]  Hans C. van Houwelingen,et al.  The Elements of Statistical Learning, Data Mining, Inference, and Prediction. Trevor Hastie, Robert Tibshirani and Jerome Friedman, Springer, New York, 2001. No. of pages: xvi+533. ISBN 0‐387‐95284‐5 , 2004 .

[27]  George Kuczera,et al.  Prediction of water yield reductions following a bushfire in ash-mixed species eucalypt forest , 1987 .

[28]  H. Akaike.  A new look at the statistical model identification , 1974 .

[29]  N. Coops,et al.  Estimating canopy structure of Douglas-fir forest stands from discrete-return LiDAR , 2007, Trees.

[30]  E. Næsset,et al.  Laser scanning of forest resources: the nordic experience , 2004 .

[31]  W. Cohen,et al.  Lidar Remote Sensing of the Canopy Structure and Biophysical Properties of Douglas-Fir Western Hemlock Forests , 1999 .

[32]  A. E. Hoerl,et al.  Ridge Regression: Applications to Nonorthogonal Problems , 1970 .

[33]  Trevor Hastie,et al.  The Elements of Statistical Learning , 2001 .

[34]  E. Næsset.  Determination of mean tree height of forest stands using airborne laser scanner data , 1997 .

[35]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[36]  Robert Tibshirani,et al.  The Elements of Statistical Learning: Data Mining, Inference, and Prediction, 2nd Edition , 2001, Springer Series in Statistics.