Next Article in Journal
Wood Density and Ring Width in Quercus rotundifolia Trees in Southern Portugal
Previous Article in Journal
Effect of Heavy Machine Traffic on Soil CO2 Concentration and Efflux in a Pinus koraiensis Thinning Stand
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Regional Variation in Forest Canopy Height and Implications for Koala (Phascolarctos cinereus) Habitat Mapping and Forest Management

by
Dave L. Mitchell
1,2,*,
Mariela Soto-Berelov
2 and
Simon D. Jones
2
1
Australian Koala Foundation, GPO Box 2659, Brisbane, QLD 4001, Australia
2
School of Science, RMIT University, Melbourne, VIC 3000, Australia
*
Author to whom correspondence should be addressed.
Forests 2021, 12(11), 1494; https://doi.org/10.3390/f12111494
Submission received: 4 October 2021 / Revised: 23 October 2021 / Accepted: 24 October 2021 / Published: 29 October 2021
(This article belongs to the Section Forest Ecology and Management)

Abstract

:
Previous research has shown that the Koala (Phascolarctos cinereus) prefers larger trees, potentially making this a key factor influencing koala habitat quality. Generally, tree height is considered at regional scales which may overlook variation at patch or local scales. In this study, we aimed to derive a set of parameters to assist in classifying koala habitat in terms of tree height, which can then be used as an overlay for existing habitat maps. To determine canopy height variation within a specific forest community across a broad area in eastern Australia, we used freely available Airborne Laser Scanning (ALS) data and adopted a straightforward approach by extracting maximum-height ALS returns within a total of 288 30 m × 30 m “virtual” ALS plots. Our findings show that while maximum tree heights generally fall within published regional-scale parameters (mean height 33.2 m), they vary significantly between subregions (mean height 28.8–39.0 m), within subregions (e.g., mean height 21.3–29.4 m), and at local scales, the tree heights vary in response to previous land-use (mean height 28.0–34.2 m). A canopy height dataset useful for habitat management needs to recognise and incorporate these variations. To examine how this information might be synthesised into a usable map, we used a wall-to-wall canopy height map derived from ALS to investigate spatial and nonspatial clustering techniques that capture canopy height variability at both intra-subregional (100s of hectares) and local (60 hectare) scales. We found that nonspatial K-medians clustering with three or four height classes is suited to intra-subregional extents because it allows for simultaneous assessment and comparison of multiple forest community polygons. Spatially constrained clustering algorithms are suited to individual polygons, and we recommend the use of the Redcap algorithm because it delineates contiguous height classes recognisable on a map. For habitat management, an overlay combining these height classification approaches as separate attributes would provide the greatest utility at a range of scales. In addition to koala habitat management, canopy height maps could also assist in managing other fauna; identifying forest disturbance, regenerating forest, and old-growth forest; and identifying errors in existing forest maps.

1. Introduction

1.1. Koala Habitat

Wildlife managers are required to balance competing demands, firstly, to provide solutions to ensure the sustainability of wildlife populations, and secondly, to mitigate potential conflicts arising from human resource needs. Resolution of these conflicts requires adequate knowledge about how wildlife uses habitats, this information is collected using, for example, presence/absence surveys, dietary studies, and telemetry. In the case of the Koala (Phascolarctos cinereus, henceforth koala), use of habitat is primarily influenced by the density of preferred food trees, e.g., [1]. This knowledge is transferred to vegetation maps to generate a habitat map which ranks vegetation communities in terms of their relative importance to the koala [2]. These maps are derived from regional-scale vegetation maps covering ~100,000 hectares and might not capture habitat variability at site or landscape scales relevant to the utility needs of both the koala and habitat management stakeholders. We refer to “stakeholders” as a wide group including government agencies, developers, individuals, and community groups; all of whom have an interest in making decisions about koala habitat.
Larger trees have been identified as a factor influencing koala habitat quality by many studies, using methods ranging from; correlations between koala use and tree size at the plot scale; to mixed multi-scale models with additional environmental factors (Table 1). Mixed models have provided varying results, but no studies published have shown that the koala prefers smaller trees.
The preference for larger trees has been attributed to several factors; (a) there may be a threshold tree size that guides the selection of foraging behaviour [4]; (b) larger trees have larger canopies providing increased food resources [7,9] and shelter [5,13]; and (c) tree health and vigour, as determined by environmental gradients [13], in more challenging environments trees are likely to be smaller, with less nutritional quality and greater levels of deterrent foliar chemicals [16], and site competition also affects tree size [17]. In addition, other factors have been recognised, for example, soil fertility influences habitat occupancy at the landscape scale [3,18] and surrogate measures of fertility (slope and elevation) have been incorporated into some koala habitat models (e.g., [19]). A study by Callaghan et al. [3] suggested that, where detailed structural mapping is available, ‘an overlay map could indicate areas where koala habitat quality is enhanced by the predominance of large, older-growth trees.’ Identifying these areas might therefore provide useful information for habitat management at a range of scales.

1.2. The Use of LiDAR (Light Detection and Ranging) for Canopy Height Information

Canopy height, defined as the statistical aggregation of individual tree heights over a nominated area [20] is one of the most important forest metrics [21] It can be a categorical or continuous variable, and is scale-independent whether used at plot, landscape, regional or even global scales [22] Canopy heights are traditionally measured from the ground using a clinometer and measuring tape, which can be a slow and laborious technique. In dense forests, observing the exact canopy top is difficult because of visual obstruction, and tree heights can be underestimated by up to 2 m [23] Since 2000, remote sensing of forests using LiDAR (Light Detection and Ranging) has developed apace, and LiDAR can be used to assess forest structure.
LiDAR data collection methods include terrestrial, unmanned aerial vehicle (UAV), airborne (ALS), and, more recently, satellite platforms. LiDAR offers several useful products: digital terrain models, tree canopy or height models, tree mapping, and plant projective cover maps [24]; and provides opportunities to explore wildlife interactions with their habitat in the third dimension, e.g., canopy height and forest structure [25,26]. Typically, ALS data is pre-classified into discrete-return classes, for example, ground, buildings, roads, shrubs, and trees. ALS has been used to develop canopy height models, and, with a sufficient density of within-canopy returns, internal structure of individual trees can be modelled [27]. ALS provides broad-scale wall-to-wall coverage at individual tree resolution and is freely available over large areas of the koala’s geographic range in eastern Australia [28].
ALS research generally focusses on two different methods to obtain canopy height, firstly, area-based methods, and secondly, individual tree detection (ITD) [27] which focusses on the delineation of individual tree crowns, or individual tree vertical profiles [29]. Both methods have different uses, for example, area-based methods provide maximum canopy height while ITD can be used to obtain average canopy height. ITD methods begin with the assumption that local elevation maxima determine the top of the tree, and then, for example, use pixel-based region-growing methods to identify “watersheds” or “valley bottoms” to determine the base of the canopy [27] or Delaunay triangulation [30]. Stopping the algorithm when the base of the tree crown is sufficiently identified can be a complicated, multi-step process [31]. Over broader regional and landscape areas, ITD, because it can detect smaller trees, reduces errors of omission, and helps maintain consistent analysis between structurally variable forests [32].
Area-based approaches, e.g., within a nominated plot or grid cell, are computationally less intensive and provide similar results as ITD for certain metrics, including the height of the largest trees [32]. Area-based approaches include Canopy Height Models (CHMs), a raster-based product with the height value determined by ALS maximum return height within the grid cell [33]. Other definitions of a dominant-height metric have been formulated, e.g., the 50 highest trees/ha and the 95th percentile of all returns [22]; the predominant height based on means of local maxima [34]; the 99th percentile [35]; modal height, mean height, and 95th percentile height [36]; the arithmetic mean of returns above the 80th percentile [35]; and the median height of all returns [37]. A recent development is “cloth simulation”, where a virtual cloth is draped over a fine-scale CHM to extract large trees across different scales by varying cloth tension [33].
Virtual ALS plots have been proposed as a proxy for field-based data collection and offer efficient and cost-effective data collection over large or inaccessible areas with provision of information, including canopy height, equivalent to the accuracy of field-based plots [38]. Furthermore, ALS plots can be matched to field-based plots when necessary, for example, when data calibration is required [39]. With efficient sampling design, virtual plots can be used to analyse forest structure in wall-to-wall ALS coverage [38].
In Australia, forestry research using LiDAR began in the early 2000s, and is commonly used by forestry agencies, for example, in biomass assessment of plantations and harvest planning [40] as well as by land management agencies [41]. However, most ITD and CHM methods have been developed in northern-hemisphere boreal forests with Pinus species and are not readily applicable to the more structurally-heterogenous Eucalyptus species forests found in Australia. These forests typically have sparse, highly clumped canopies and erectophile leaf angle distributions, allowing ALS pulses to penetrate deeper into the canopy; therefore, they may underestimate dominant or maximum height [42]. Conversely, denser canopies in some tree species affect the accuracy of the underlying terrain model from which canopy heights are obtained [43]. Considering these two factors, point densities above 0.5 pts/m2 provide forest metrics, including height, that did not differ significantly from those obtained with a density of 10 pts/m2 [43]. ALS metadata typically specifies a target pulse density of 1 pt/m2. At broader scales, ALS has been used to calibrate space-borne LiDAR data to produce vegetation community structural classifications in Queensland [36], and across Australia [44].
Koalas prefer vegetation communities with a high proportion of food trees [2,3], and they also prefer habitats with larger trees (Table 1). However, no previous research has attempted to combine these factors in a habitat management map. Given extensive ALS coverage is available in the southeast Queensland and the New South Wales north coast regions where koalas are widely distributed, our aim is to evaluate the utility of ALS in developing high resolution canopy height models which can complement existing low resolution regional-scale habitat maps and further assist wildlife managers to better manage koalas. To do this, we ask the following research questions: (1) do published herbarium descriptions adequately capture height variability at regional, sub-regional, and intra-subregional scales, and; (2) how can canopy height variability best be incorporated into a koala habitat map that can be useful to wildlife managers across different spatial scales.

2. Materials and Methods

Using herbarium descriptions as a baseline, we first examine height variability at different scales, and then, how this information could be incorporated into existing koala habitat management maps to provide greater utility at a range of scales. We make the following hierarchical scale definitions:
  • Regional: The whole study region encompassing southeast Queensland and the adjoining Tweed Local Government Area in New South Wales (2,500,000 hectares),
  • Subregional Group: Any identified area which is geographically separate from other areas within the study region (30,000–230,000 ha),
  • Intra-subregional Group: Any identified area which is geographically separate from other areas within a subregional group (400–1600 ha), or our subregional study area,
  • Polygon Group: For map production, we used contiguous groups of 25 m × 25 m cells derived from a vector polygon within intra-subregional groups (<100 ha).

2.1. Variation in Forest Canopy Height

We chose a floristically variable vegetation community with tree species known to be used by the koala [2,45] within the Southeast Queensland (SEQ) Interim Bioregion [46], which covers SEQ and part of the New South Wales north coast. In SEQ, the community is mapped as Regional Ecosystem (RE) 12.11.3 [47]. Its floristic composition is described by the Regional Ecosystem Descriptions Database [48] as mainly composed of Eucalyptus siderophloia (Grey Ironbark) and E. propinqua (Small-fruited Grey Gum), with other patchily distributed species including E. microcorys (Tallowwood), Lophostemon confertus (Brush Box), Corymbia intermedia (Pink Bloodwood), E. biturbinata (a Grey Gum), E. acmenoides (White Mahogany), and E. tereticornis (Queensland Blue Gum). Vegetation communities similar to RE 12.11.3 extend into New South Wales [49,50].
RE 12.11.3 forest structure is described in a Technical Description [51], and varies from Tall Open Forest to Woodland, with overstorey heights ranging between 15–37 m (average 28.1 m). The Description uses 66 Herbarium CORVEG plots, data collected at these plots includes floristic composition and crown cover from all forest structural layers at relatively undisturbed and representative plots, CORVEG data was collected between 1991–2011, 30 of these plots were collected within 10 years of ALS acquisition which began in 2009. We aimed to use all 38 plots within the ALS footprint but discarded five sites which had been disturbed by development post-CORVEG data collection.
RE 12.11.3 in southeast Queensland and the associated NSW community covers 76,390 ha. Our study was limited to those areas containing ALS data (49,680 ha) (Figure 1).
We developed a workflow (Figure 2) providing a multiscale approach to assess the forest height variability within the vegetation community (henceforth referred to as RE 12.11.3) at the regional, subregional, and intra-subregional scales.
Figure 1. We examined a forest community used by koalas (Regional Ecosystem 12.11.3) occurring across Southeast Queensland and the New South Wales north coast (Tweed). Within this study area, six distinct subregions (dotted black ellipses) of RE 12.11.3 can be recognised by spatial proximity or separation by rivers or state borders. We selected the Koala Coast to examine the intra-subregion canopy height variability (Figure 3) and an additional site on the Sunshine Coast to examine the height variability resulting from historical forest disturbances (blue dot; see Figure 4). ALS coverage is shown (pale yellow background); vegetation community RE 12.11.3 is covered by ALS (dark green, with similar areas in NSW), and RE 12.11.3 not covered by ALS (light green). Random plots (205) within ALS-covered RE 12.11.3 (white circles) and CORVEG plots (33) within ALS-covered RE 12.11.3 (black circles) are shown.
Figure 1. We examined a forest community used by koalas (Regional Ecosystem 12.11.3) occurring across Southeast Queensland and the New South Wales north coast (Tweed). Within this study area, six distinct subregions (dotted black ellipses) of RE 12.11.3 can be recognised by spatial proximity or separation by rivers or state borders. We selected the Koala Coast to examine the intra-subregion canopy height variability (Figure 3) and an additional site on the Sunshine Coast to examine the height variability resulting from historical forest disturbances (blue dot; see Figure 4). ALS coverage is shown (pale yellow background); vegetation community RE 12.11.3 is covered by ALS (dark green, with similar areas in NSW), and RE 12.11.3 not covered by ALS (light green). Random plots (205) within ALS-covered RE 12.11.3 (white circles) and CORVEG plots (33) within ALS-covered RE 12.11.3 (black circles) are shown.
Forests 12 01494 g001

2.1.1. Regional Scale

Within the study area, we nominated six sub-regions (Figure 1). The Sunshine Coast and Western subregions have 38% and 29% ALS coverage, respectively, but the other subregions are entirely covered. Within the region covered by ALS, in addition to the 33 CORVEG plots, we used ArcMap [52] to generate 205 random 30 m × 30 m “virtual” plots within RE 12.11.3, which were at least 100 m apart.
The technical description for RE 12.11.3 provides the average canopy height for the RE by estimating the heights of at least three trees representing the median canopy height in each of the 66 [53] representative 20 m × 50 m plots [54]. Owing to the absence of individual plot orientation, we could not replicate this procedure, so we used 30 m × 30 m plots centred on the CORVEG plots and extracted the maximum canopy height in each of these plots using ALS.

2.1.2. Subregional-Scale

Within the wider region, we nominated six spatially separate subregional groups (Figure 1) and ensured that each subregion had at least 30 of the 205 virtual plots to enable an assessment of the cross-region height variation and compare it against the CORVEG plot heights.

2.1.3. Intra-Subregional-Scale

To examine the variability at finer scales, we selected the geographically distinct Koala Coast subregion (Figure 1), which also has the highest density of koalas in the wider region. Within the Koala Coast, we nominated seven distinct intra-subregional groups based on spatial separation, with two close groups demarcated by differences in elevation (Figure 3).
For the Koala Coast intra-subregional groups, we generated a CHM covering all RE 12.11.3 polygons and populated an aligned rectangular 25 m × 25 m polygon grid with CHM height values.
To further investigate the intra-subregional variation, and to examine the influence of disturbance on the forest height, we used data from a 60-hectare study site in the Sunshine Coast subregion. This focus study site encompassed two tenures with different land use histories—a private property “Quinlans” and Mapleton National Park. Quinlans was a dairy farm up until the 1960s, with the vegetation mostly cleared; the subsequent regrowth vegetation is now classified as Remnant Vegetation [47]. Mapleton NP was declared a forest reserve in 1923, and its vegetation is also mapped as Remnant. Figure 4 shows the vegetation cover in 1940; in some areas of Quinlans, the vegetation was almost completely removed, and the remaining forest areas were extensively thinned to an open woodland. A 1958 aerial photo (not shown) indicates that tree regrowth was actively occurring in Quinlans at that time. Areas of possible disturbance in Mapleton NP 1940 are also shown and might be because of previous logging, fire, and/or storm damage or topographic influences (dryer northerly aspect).
Within this site, we used 50 30 m × 30 m virtual plots spaced on a 100-m map grid (Map Grid of Australia MGA94 Zone 56). In May 2019, we collected data, including tree species, diameter at breast height (DBH), and MGA coordinates, for every tree over 10-cm DBH at 44 sites. Data collection was aimed at obtaining detailed knowledge of both tree species variability and structural variability across the site and to provide a context to assist the interpretation of the ALS data.

2.1.4. Height Variability—Extracting Maximum Canopy Height from ALS

We obtained all ALS data from ELVIS (Elevation Information System) [28] at https://elevation.fsdf.org.au/, accessed on 24 June 2021. This site collects elevation data from many Australian government agencies and includes classified ALS point clouds. Our regional study area encompassed multiple ALS acquisition projects (Table 2); all available metadata showed that the minimum spacing requirement of 0.5 returns/m2, appropriate for extracting forest metrics, such as maximum canopy height, was met [43].
For all virtual plots, we clipped each ALS tile to the 30 m × 30 m plot footprints using the lasclip tool [55] in ArcMap [52] and obtained the maximum heights using lasheight [55]. We then manually extracted the maximum height from the Arc Map’s layer display while visually checking all plots for ALS misclassifications (e.g., birds) and correcting when necessary.

2.1.5. Statistical Analysis

We first examined the datasets for sample size, normality, and variance (Table 3). The sample sizes varied considerably; the majority were not distributed normally, and uneven variances occurred at all the scales. For two-group comparisons, we used Welch’s t-test, and for more than two groups, we determined that ANOVA with Games-Howell post hoc tests (with bootstrapping) were suitable for analysis, since these tests do not assume normality, equal variances, or sample sizes [56].

2.2. Map Production—Canopy Height Model (CHM)

This part of our research examined the implications of the regional variations in tree height in the context of providing a forest height map that could be used for habitat or forest management at the regional, subregional, and intra-subregional scales. We considered using individual tree detection (ITD), but since we only required the maximum canopy height, we decided that a canopy height model (CHM) would suffice for our purposes, with the added advantage of this approach being much easier to implement over broad areas. For the Koala Coast subregion and the Sunshine Coast study site, we extracted all ALS data overlapping RE 12.11.3 polygons and generated a 25-m resolution CHM using the canopymodel batch function in Fusion software [57] to extract the maximum height for each 25 m × 25 m cell. We again visually checked all tiles and removed misclassifications (birds, telecommunications towers, and a high-voltage powerline). We generated a 25 m × 25 m “fishnet” aligned with the CHM and extracted the CHM height values to the fishnet for analysis. We chose this cell size (0.0625 hectares) as a compromise between the Regional Ecosystem minimum mapping unit (0.5–1 ha) in Southeast Queensland [54] and some smaller units in the Tweed dataset (e.g., 0.025 ha). This cell size also provides sufficient resolution for decision-making at the habitat management and development scales.

Producing a Classified Height Map

For those involved in koala habitat or forest management decision-making, classified maps are a valuable tool that enable a feature of interest to be interpreted, compared, and evaluated in a spatial context. We now address our second research aim, “How can canopy height variability best be incorporated into a map?” Our approach was to use statistical clustering methods. Spatial clusters can be defined as a grouping of similar objects in a geographic space that are more homogenous than those from different clusters and that occupy contiguous regions [58,59] a cluster analysis uses analytical pattern recognition to group similar data points [60]. There are many different clustering algorithms, and a selection of the most suitable method requires some consideration of the data structure, e.g., data density and arrangement, and purpose, i.e., whether we sought out global or more localised cluster identification. Our dataset had distinctive characteristics: (a) height was the only variable, (b) the sample size (i.e., CHM) was 100%, (c) the data was in the form of unconnected polygons, each composed of connected groups of cells, and (d) the polygon groups had different sizes and shapes, i.e., they may be compact, linear, or complex. Unconnected polygon groups required a nonspatial global approach, whereas individual polygon groups (i.e., with contiguous cells) were also suited to spatially constrained algorithms, since they met the contiguity requirement. Our goal was to evaluate several clustering algorithms and recommend those that (a) displayed canopy height information in a user-friendly map format useful for habitat management, (b) were reproducible using readily available software, and (c) were amenable to a high degree of automation. For map production, we used two polygon groups drawn from the Koala Coast and our Sunshine Coast study site.
For these areas, we firstly derived maps showing interquartile height ranges (i.e., four classes with equal numbers of data points) and visually identified six apparent clusters of higher canopy heights. Cluster maps were generated using algorithms implemented in Geoda [61], chosen because (a) it is freely available and regularly updated software with a graphical user interface (https://geodacenter.github.io/download.html, accessed on 1 October 2021), (b) the data formats are compatible with most GIS software, and (c) a wide selection of clustering algorithms is available. We chose eight candidate algorithms suited to our cellular data and used default values for each algorithm. For each algorithm we produced maps with two, three, and four classes, and to assess the map utility and interpretability, a three-member expert panel subjectively evaluated each output using the criteria shown in Table 4.

3. Results

3.1. Regional Variation

The RE 12.11.3 technical description [53] uses canopy heights obtained from an average of three to five trees at each CORVEG plot [54] and summarises the data from 66 plots as the mean, maximum, and minimum heights (Table 5). We were not able to replicate this procedure, so we recorded the maximum height in each of our 33-plot ALS subsamples of the CORVEG plots, which we used to compare with all the other ALS virtual plots. Since we measured the maximum height at each CORVEG plot (Table 4), the ALS plots were substantially higher (maximum height 51.4 m) than the technical description (maximum average height 37 m), and the average maximum ALS height (33.2 m) was higher than the published mean height (28 m) [53].
These results indicate that the CORVEG plots were not significantly different in height to the random plots (p = 0.057) and so, for RE 12.11.3, the Herbarium technical description adequately captures the height variations across the wider Southeast Queensland region and the adjoining Tweed Coast in New South Wales.

3.2. Subregional Variation

The 205 randomly selected virtual plots across the Southeast Queensland Bioregion were assigned to six identified subregional groups to, firstly, examine the variations between the subregions and, secondly, to identify those subregions outside the CORVEG height ranges. The average maximum ALS canopy heights measured at the CORVEG plots were substantially below the average of all the subregional group heights, except the Koala Coast and Western. The Sunshine Coast and Tweed subregional groups had higher canopy heights with a larger range of heights; in contrast, the Western group had the lowest canopy. The Koala Coast canopies were also lower, with the smallest variations in height range (Figure 5).
The Koala Coast had substantially lower mean maximum heights (31.32 m) than the average of the other near-coastal regional groups and a lower standard deviation in maximum height (3.89 m) than the average of these groups (6.5 m). The Koala Coast and Western subregions had lower average canopy heights than CORVEG; all the other subregions had higher average heights than CORVEG. The height distribution was positively skewed, particularly for the Gold Coast and Western Regions (Table 6).
The differences between CORVEG and the subregional plots were not significant for the Sunshine Coast and Tweed subregions, and between the subregions, the Koala Coast and Western subregions were not significantly different to each other, but both were significantly different to the other four subregions. These four near-coastal subregions (Brisbane, Gold Coast, Sunshine Coast, and Tweed) were not significantly different from each other (Table 7).

3.3. Intra-Subregional Variation

3.3.1. Intra-Subregional Variation (Koala Coast)

We previously identified (Figure 3) seven intra-subregional groups within the Koala Coast subregional group, and for this part of our analysis, we used a 100% sample rate, i.e., all 25 m × 25 m cells within each group. Figure 6 shows similar within-group variations, except for Group 6, where there was a greater range of height values.
Table 8 shows height variations between these groups; compared to CORVEG, these groups showed a more compact canopy (standard error, standard deviation, and variance were all lower), and the CORVEG distribution was strongly positively skewed. The statistical comparisons are shown in Table 9.

3.3.2. Intra-Subregional Variation (Sunshine Coast Study Site)

The effects of forest disturbances are readily apparent at the Sunshine Coast intra-subregional study site, which encompasses two land tenures with different land use histories. The canopy heights fall within the CORVEG height ranges but are substantially less than Sunshine Coast virtual plots, and the canopy height in Quinlans averages 6.2 m below the adjoining Mapleton National Park (Figure 7 and Table 10) but still falls within the CORVEG parameters (e.g., range and standard deviation from the mean). Pairwise differences between CORVEG and the entire study site, and CORVEG and Mapleton National Park, were not significantly different, but all other pairwise comparisons were significantly different (Table 11).
While it may technically fall within the definition of Remnant Forest (Neldner et al. 2019; Table 4), compared to CORVEG and all the other Sunshine Coast datasets, Quinlans is lower in height, and we can conclude that Quinlans has a regenerating forest. The difference in canopy heights between the Mapleton National Park and Quinlans is illustrated in Figure 8 using the ALS profile.

3.4. Clustering and Sample Map Production

The Sunshine Coast study site and two polygon groups in the Koala Coast were selected to evaluate clustering algorithms. We chose these because; (a) we have detailed disturbance history and field data for the study site; (b); the two Koala Coast groups represent higher and lower canopy cohorts; and (c) the three polygon groups are representative of compact, linear, and complex shapes which might, potentially, influence clustering algorithms in different ways. The three groups are shown in Figure 9, coloured by the interquartile ranges (IQR) for each group of cells, with median heights for each IQR shown. By visual inspection we identified six distinct height canopy areas, marked A–E (higher canopy), and an ecotone identified during fieldwork marked F (moderate canopy height) in Figure 9.
For 90 output cluster maps, an expert panel of three visually evaluated each map to assess how well the algorithm identified the six nominated features for individual polygons, results are summarised in Table 12. Elbow plots (see Appendix A) show the class variance explained by each algorithm, indicate the appropriate number of clusters, and algorithm efficiency. The complete evaluation process is detailed in Appendix A.
Three-class cluster maps with the highest evaluation scores are shown in Figure 10. Nonspatial K-means and hierarchical clustering successfully identified five of six features, but cluster areas varied substantially between algorithms for Feature D. Feature F, the forest ecotone community, was only partially captured by K-medians.
Four-class cluster maps are shown in Figure 11. K-medians again ranked highest within the nonspatial algorithms, but several features have less well-defined clusters than three-class maps. Redcap [62] had the highest score from the expert panel, with tightly defined contiguous clusters for all features. Redcap extended the area of feature A, and effectively shows the height difference between Mapleton National Park and regenerating forest in Quinlans (see Figure 4). Redcap also accurately captures the extent of the forest ecotone in Quinlans (Feature F). However, compared to Figure 9 and other cluster maps (Figure 10), Redcap overestimated the extent of one feature (D) and underestimated another (E).
The panel evaluation for multiple polygon cluster maps produced using all Koala Coast subregional data is summarised in Table 13.
K-means and K-medians were ranked equally highly for three classes (Figure 12), but K-medians was ranked higher for four classes. Features D and E were identified, but because features B and C were not high compared to D and E, were assigned to lower canopy height classes. In addition, areas of High Value Regrowth Forest [63] were identified by both algorithms.

4. Discussion

Current koala habitat mapping uses regional-scale vegetation maps incorporating a high degree of floristic variability, which may limit their application at the local scale (1–300 hectares), where tree species variability, i.e., food resource availability, becomes more apparent. Since koalas prefer larger trees [13,16], incorporating this variable into a map overlay might be useful for habitat management at this scale [3]. Our study was aimed at how best to incorporate larger trees into the current habitat maps using heights obtained from ALS. To derive the parameters for an overlay map, we used published information as a baseline and found that, at the regional scale, forest canopy heights are within published parameters [53] but, at the subregional and intra-subregional scales, the height variability is generally outside these parameters. Additionally, the variability is both significantly different between subregions and within subregions. Two height parameters demonstrated the effects of previous land use. In regenerating forests, not only is the canopy height lower, but the canopy has a lower standard deviation in height, i.e., the canopy is “flatter”, with shade-tolerant and shade-intolerant species in active competition [17]. The range in the standard deviation therefore has potential for identifying similar areas of regenerating forests.
To incorporate ALS canopy height information into a habitat map overlay, we examined several clustering algorithms and found that nonspatial K-medians best characterised the nominated high-canopy areas for multiple polygons at the intra-subregional scale. This approach might be useful for height comparisons between different areas for broader management purposes but produced clusters suffering from a lack of contiguity, limiting their use at finer scales. At the site scale, i.e., individual habitat map polygons, the spatially constrained hierarchical clustering algorithm Redcap [61,62] was able to: (a) provide greater map utility by identifying contiguous clusters; (b) separate historically disturbed, lower canopy height forests from relatively undisturbed higher canopy forests, and (c) capture a forest ecotone with a higher proportion of preferred koala food resources. Redcap further showed its utility in another polygon by capturing an area independently mapped as Mature Regrowth Forest [63]. We used algorithms requiring no user input so that our approach could be automated, allowing any number of forest map polygons to be clustered at any desired scale. The number of algorithms we used, and the consistent results between many of them, confirmed that most of our visually nominated features were, indeed, real clusters of similar heights and that our algorithm recommendations were valid. This multifaceted approach was recommended by Ketchen and Shook [64], who warned against adopting any single clustering algorithm. However, we noted that two of our nominated features were not recognised as being statistically significant by any algorithm, which shows that arbitrary classifications of canopy height (a) may not be meaningful and (b) may not be repeatable.
The previous use of multivariate modelling that included ALS-derived canopy height as a variable has shown that this metric features prominently within these analyses, for example, by using multivariate K-means clustering to portray forest structures, one study found that most landscape indices were rendered redundant by the median and standard deviation of the canopy height [65], i.e., ALS provides directly measured predictor variables and improves the accuracy of maps compared to inferred variables. In Alberta, Canada, a cluster analysis with ALS variables (canopy height, cover, and height variation) was used to identify areas of high conservation values within mapped forest polygons at the regional scale [66]. In another example, a Northern Spotted Owl (Strix accidentalis caurina) nesting habitat in Oregon, USA, mapped using traditional approaches (aerial photo interpretation, satellite imagery, and environmental predictor variables), was substantially improved with the addition of the ALS-derived large tree density and average stand height [67], and the importance of height was confirmed by another study in Northern California [68]. The maximum canopy height was found to be the best predictor of Spotted Owl habitats in Northern California (Strix occidentalis) [69]. In The Netherlands, a study found that, for species distribution models, the ALS-derived canopy height was a better predictor of bird species presence/absence than a suite of land cover variables [70] and had advantages at multiple scales where species habitat requirements are either not known or vary according to particular needs (e.g., feeding, breeding, and dispersal) [71].
We have recommended the use of Redcap with four classes for individual polygons, but to further maximise the overlay map utility at broader scales, for example, in conservation planning, we recommend that a height overlay also includes height classes derived from K-medians. For at-risk species dependent on large hollow-bearing trees, such as the greater glider (Petauroides volans) [72], this multiscale approach might be more appropriate. Koala researchers have previously stratified field sites based on low-resolution vegetation community types and/or geology maps [18]; sampling strategies using height clusters could reduce field efforts by effectively sampling areas of mature forest, regrowth, or ecotones within these broader map units. A similar approach using LiDAR has been used to increase the efficiency of survey plot sampling designs in pine plantations [73].
In Queensland, the distinction between Remnant and Regrowth vegetation is important, because, e.g., Remnant Vegetation [54] is subject to stricter legislation governing land use compared to Regrowth or other undefined forest; therefore, it is necessary to accurately differentiate between these forest types. We identified potential mapping errors within our sample forest communities by identifying areas below the minimum height required for classification as Remnant Vegetation [54]; the application of our methods could highlight these mapping errors so they can be corrected. During our research, we noted gaps and overlaps between the Remnant and Mature Regrowth datasets, and because we successfully captured mapped areas of Mature Regrowth [63], our methods could also be used to rectify these mapping errors.
We acknowledge that our study was limited in scope, with the Sunshine Coast and Western areas having only approximately 30% LiDAR coverage. As time goes on, we can expect greater LiDAR coverage, not only from ALS but from newer LiDAR platforms, e.g., the recently commissioned Global Ecosystem Dynamics Investigation (GEDI) satellite [74]. Our study only examined one regional vegetation community and assessed three community polygons with clustering, so our results may not be applicable to all forest communities. However, by choosing a floristically variable community with a likely corresponding variability in canopy height, and with disturbance information for one polygon and partial disturbance information for another, we regard our study as a valuable contribution to this topic.

5. Conclusions

In this study we showed that, because forest canopy heights vary at the subregional and intra-subregional scales, height overlay maps need to incorporate this variability. A regional framework provides an overall picture of the forest canopy height, but at the subregional scale and intra-regional scale relevant to habitat management, the overlay can incorporate attributes derived from the canopy height at these scales. This approach can be implemented over broad areas to produce a multiscale map that can be used to capture an important koala habitat attribute, i.e., that koalas prefer larger trees at any given location.
For koalas, which constantly face development pressure, habitat maps depicting areas of large trees could provide valuable assistance in identifying valuable habitats for protection. Our clustering approach could provide unambiguous and independent tools to assess development at both the site and local government authority scales, reducing the reliance of development approval authorities on information supplied by development proponents. Other potential uses include the identification of candidate old-growth forests, mapping areas of larger hollow-bearing trees required for fauna, conservation reserve planning, and improvement in the spatial accuracy of forest community mapping.

Author Contributions

Conceptualization, D.L.M., S.D.J. and M.S.-B.; methodology, D.L.M., S.D.J. and M.S.-B.; formal analysis, D.L.M.; investigation, D.L.M.; data curation, D.L.M.; writing—original draft preparation, D.L.M.; writing—review and editing, D.L.M., M.S.-B. and S.D.J.; supervision, S.D.J. and M.S.-B.; project administration, S.D.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Australian Koala Foundation, GPO Box 2659, Brisbane, Queensland 4001, Australia, www.savethekoala.com (accessed on 10 October 2021).

Acknowledgments

We particularly thank Deborah Tabart OAM (AKF Chief Executive Officer) and the AKF Board for their ongoing dedication and support. We would also like to acknowledge the international volunteers who assisted with sometimes arduous field work (Ben Masters, Elisa Chapuis, Katie O’Connell, Yulan Lu). Special thanks to Derek Johnson (QTree Vegetation Assessment) and David Staines (AKF) for fieldwork assistance and cluster map evaluation; and Tim Ryan (Queensland Herbarium) for advice on using REDD, CORVEG, and the Technical Descriptions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Appendix A.1. Introduction

We regarded simple height classification using arbitrary height classes unsatisfactory for two reasons: (a) class definitions are arbitrary and (b) may not identify actual areas of interest. Clustering algorithms are capable of separating heights in a statistically robust and repeatable manner and require little or no user input. With the aim of producing classified canopy height maps suitable for Koala habitat management and other forest management purposes, we firstly used clustering algorithms implemented in Geoda software (https://geodacenter.github.io/cheatsheet.html, accessed on 24 June 2021) to produce three separate cluster maps for individual polygon groups representing compact, linear, and complex shapes. Then, for the two of these polygons within the Koala Coast, we ran those algorithms capable of clustering multiple polygons to (a) enable comparisons between any polygon groups and (b) compare these with maps generated with individual polygon groups.
We evaluated the maps in three ways. Firstly, we constructed elbow plots, which (a) indicate the appropriate number of clusters and, (b) by displaying the variance explained by the number of clusters for each algorithm, provide comparisons between algorithms. Secondly, an expert panel composed of one of the authors (DM) and two analysts scored each map visually using three simple criteria. Then, the maps were assessed by the magnitude of separation between the highest canopy class and the next-highest class. The three methods were individually ranked, and finally, a total rank score was used to choose the best map.
The cluster maps, all with two, three, or four classes, are shown below (Appendix A.2), together with the panel score for each map (see Appendix A.5 for the score allocation).

Appendix A.2. Individual Polygon Clustering

Our aim was to evaluate which Geoda algorithm best detects higher canopy areas (features A–F) in each of the three polygons depicted in Figure 12, reproduced below (Figure A1). The algorithms for the individual polygon groups used only the height data within that group.
Figure A1. Compact (a), linear (b), and complex (c) polygon groups, showing the interquartile ranges (IQR) and interquartile median canopy heights within each group. We used this map to assess and compare the clustering algorithm outputs which follow. Spatial patterns are evident, for, e.g., patches of the highest canopy (A–E, light blue) and an area of moderate height (F, dark green), surrounded by a canopy with a slightly lower height (light green). The height appears more randomly distributed in other areas, except for lower canopy height areas (dark blue). All maps use the same scale.
Figure A1. Compact (a), linear (b), and complex (c) polygon groups, showing the interquartile ranges (IQR) and interquartile median canopy heights within each group. We used this map to assess and compare the clustering algorithm outputs which follow. Spatial patterns are evident, for, e.g., patches of the highest canopy (A–E, light blue) and an area of moderate height (F, dark green), surrounded by a canopy with a slightly lower height (light green). The height appears more randomly distributed in other areas, except for lower canopy height areas (dark blue). All maps use the same scale.
Forests 12 01494 g0a1
Figure A2. Nonspatial algorithms with 2 classes, using K-means (ac), K-medians (df) and nonspatial hierarchical clustering (gi). Letters A–F indicate areas of higher canopy. features A–E are captured. K-medians best captures the areal extent of these features and best captures F, an unmapped area of the koala food tree Eucalyptus tereticornis (d). Panel scores: K-means = 16, K-medians = 17, and Hierarchical = 12.
Figure A2. Nonspatial algorithms with 2 classes, using K-means (ac), K-medians (df) and nonspatial hierarchical clustering (gi). Letters A–F indicate areas of higher canopy. features A–E are captured. K-medians best captures the areal extent of these features and best captures F, an unmapped area of the koala food tree Eucalyptus tereticornis (d). Panel scores: K-means = 16, K-medians = 17, and Hierarchical = 12.
Forests 12 01494 g0a2
Figure A3. Spatially constrained K-means (ac) and K-medians (df) with 2 classes. These algorithms use spatial weights (x, y), which, in all cases, outweigh those of the object (h = height). Letters A–F indicate areas of higher canopy, only one of the designated features (E) is adequately captured using these algorithms (f). Panel scores: K-means = 6 and K-medians = 8.
Figure A3. Spatially constrained K-means (ac) and K-medians (df) with 2 classes. These algorithms use spatial weights (x, y), which, in all cases, outweigh those of the object (h = height). Letters A–F indicate areas of higher canopy, only one of the designated features (E) is adequately captured using these algorithms (f). Panel scores: K-means = 6 and K-medians = 8.
Forests 12 01494 g0a3
Figure A4. Spatially constrained hierarchical (SCHC, (ac), Skater (df), and Redcap (gi)) clustering maps with 2 classes. Letters A–F indicate areas of higher canopy, features A and B are adequately captured by these algorithms (a,b,d,e,g,h), but the areal extents of D and E (e,f,i) are too large. C and F are not recognised. Redcap separates an area of regrowth forest on land previously used for farming (g, dark blue) adjoining a relatively undisturbed national park (light blue), which encompasses Feature A. Panel scores: Hierarchical clustering = 9, Skater = 9, and Redcap = 9.
Figure A4. Spatially constrained hierarchical (SCHC, (ac), Skater (df), and Redcap (gi)) clustering maps with 2 classes. Letters A–F indicate areas of higher canopy, features A and B are adequately captured by these algorithms (a,b,d,e,g,h), but the areal extents of D and E (e,f,i) are too large. C and F are not recognised. Redcap separates an area of regrowth forest on land previously used for farming (g, dark blue) adjoining a relatively undisturbed national park (light blue), which encompasses Feature A. Panel scores: Hierarchical clustering = 9, Skater = 9, and Redcap = 9.
Forests 12 01494 g0a4
Figure A5. Nonspatial algorithm maps (K-means, ac; K-medians, df; nonspatial hierarchical, gi) with 3 classes. Letters A–F indicate areas of higher canopy. K-medians (df) best captures features A, D, C, all algorithms capture B and C. No algorithm adequately captures F, the unmapped area of Eucalyptus tereticornis. Panel scores: K-means = 12, K-medians = 16, and Hierarchical = 11.
Figure A5. Nonspatial algorithm maps (K-means, ac; K-medians, df; nonspatial hierarchical, gi) with 3 classes. Letters A–F indicate areas of higher canopy. K-medians (df) best captures features A, D, C, all algorithms capture B and C. No algorithm adequately captures F, the unmapped area of Eucalyptus tereticornis. Panel scores: K-means = 12, K-medians = 16, and Hierarchical = 11.
Forests 12 01494 g0a5
Figure A6. Spatially constrained K-means (ac) and K-medians (df) cluster maps with 3 classes. Letters A–F indicate areas of higher canopy. Again, these algorithms use spatial weights (x, y) which outweigh those of the object (h = height). The only features adequately captured by these algorithms are A (a,d) and E (c,f). Panel scores: K-means = 9 and K-medians = 9.
Figure A6. Spatially constrained K-means (ac) and K-medians (df) cluster maps with 3 classes. Letters A–F indicate areas of higher canopy. Again, these algorithms use spatial weights (x, y) which outweigh those of the object (h = height). The only features adequately captured by these algorithms are A (a,d) and E (c,f). Panel scores: K-means = 9 and K-medians = 9.
Forests 12 01494 g0a6
Figure A7. Spatially constrained hierarchical (SCHC, (ac), Skater (df), and Redcap (gi)) clustering maps with 3 classes. Letters A–F indicate areas of higher canopy, features A and B are adequately captured by all algorithms, but the areal extent of E is too small (c,i) or too large (f). C and D are not recognised. Redcap (g) separates an area of regrowth forest on land previously used for farming (dark blue) adjoining a relatively undisturbed national park (light blue), which encompasses Feature A and recognises F, the unmapped area of Eucalyptus tereticornis. Panel scores: Hierarchical clustering = 7, Skater = 9, and Redcap = 9.
Figure A7. Spatially constrained hierarchical (SCHC, (ac), Skater (df), and Redcap (gi)) clustering maps with 3 classes. Letters A–F indicate areas of higher canopy, features A and B are adequately captured by all algorithms, but the areal extent of E is too small (c,i) or too large (f). C and D are not recognised. Redcap (g) separates an area of regrowth forest on land previously used for farming (dark blue) adjoining a relatively undisturbed national park (light blue), which encompasses Feature A and recognises F, the unmapped area of Eucalyptus tereticornis. Panel scores: Hierarchical clustering = 7, Skater = 9, and Redcap = 9.
Forests 12 01494 g0a7
Figure A8. Nonspatial algorithm cluster maps, K-means (ac), K-medians (df) and hierarchical (gi) with 4 classes. Letters A–F indicate areas of higher canopy, K-medians (df) best captures all the features, but D and F are not contiguous. B and C are captured in all maps, with little difference between algorithms. Panel scores: K-means = 11, K-medians = 15, and Hierarchical = 12.
Figure A8. Nonspatial algorithm cluster maps, K-means (ac), K-medians (df) and hierarchical (gi) with 4 classes. Letters A–F indicate areas of higher canopy, K-medians (df) best captures all the features, but D and F are not contiguous. B and C are captured in all maps, with little difference between algorithms. Panel scores: K-means = 11, K-medians = 15, and Hierarchical = 12.
Forests 12 01494 g0a8
Figure A9. Spatially constrained K-means (ac) and K-medians (df) cluster maps with 4 classes. Letters A–F indicate areas of higher canopy. Again, these algorithms use spatial weights (x, y) that far outweigh those of the object (h = height). As with 3 classes for these algorithms, the only features adequately captured are A (a,d) and E (c,f). Panel scores: K-means = 7 and K-medians = 7.
Figure A9. Spatially constrained K-means (ac) and K-medians (df) cluster maps with 4 classes. Letters A–F indicate areas of higher canopy. Again, these algorithms use spatial weights (x, y) that far outweigh those of the object (h = height). As with 3 classes for these algorithms, the only features adequately captured are A (a,d) and E (c,f). Panel scores: K-means = 7 and K-medians = 7.
Forests 12 01494 g0a9
Figure A10. Spatially constrained hierarchical (SCHC, (ac), Skater (df), and Redcap (gi)) clustering maps with 4 classes. Letters A–F indicate areas of higher canopy. Redcap (g) separates an area of regrowth forest on land previously used for farming (A, dark blue) adjoining a relatively undisturbed national park (light blue). Redcap also recognises F, the unmapped area of Eucalyptus tereticornis (dark green), which is slightly lower in height than the national park canopy but higher than the surrounding forest (dark blue). All algorithms portray C (b,e,h) as slightly lower than B, but the areal extent of D is exaggerated by Skater (i). Skater is the algorithm that best captures the areal extent of E (i). Panel scores: Hierarchical clustering = 14, Skater = 15, and Redcap = 16.
Figure A10. Spatially constrained hierarchical (SCHC, (ac), Skater (df), and Redcap (gi)) clustering maps with 4 classes. Letters A–F indicate areas of higher canopy. Redcap (g) separates an area of regrowth forest on land previously used for farming (A, dark blue) adjoining a relatively undisturbed national park (light blue). Redcap also recognises F, the unmapped area of Eucalyptus tereticornis (dark green), which is slightly lower in height than the national park canopy but higher than the surrounding forest (dark blue). All algorithms portray C (b,e,h) as slightly lower than B, but the areal extent of D is exaggerated by Skater (i). Skater is the algorithm that best captures the areal extent of E (i). Panel scores: Hierarchical clustering = 14, Skater = 15, and Redcap = 16.
Forests 12 01494 g0a10

Appendix A.3. Multiple Polygon Clustering

Non-contiguous polygons can only be clustered by nonspatial algorithms, i.e., K-means, K-medians, and hierarchical clustering. For the Koala Coast area, we selected two polygon groups, one from intra-subregional Group 6 (212 cells, 13.3 ha) and one from Group 7 (1357 cells, 84.8 ha), i.e., a polygon from each of the high- and low-canopy cohorts. We again produced a canopy height map classified by interquartile height ranges (IQR) as a reference point but with height data from 6007 cells from all seven Koala Coast intra-subregional groups (Figure A11). Some clustering is apparent in the maps—in particular, in Map (b), two areas of higher canopy (D and E), and a section of lower canopy in the southeast. In Map (a), the higher canopy areas (B and C) are differentiated from the low-canopy height (dark blue, IQR 1). Figure A11 thus identifies the same features as Figure A1 but provides a subregional context, i.e., whilst B and C are higher than the remainder of their polygon group, they are significantly lower than D and E and similar areas across the subregion.
Figure A11. Selected polygon groups from Koala Coast intra-subregional Groups 6 (a, low-canopy cohort) and 7 (b, high-canopy cohort). Letters B–E indicate areas of higher canopy These maps are similar to Figure A1, with features B–E recognisable, but, because they encompass the entire canopy height range within the Koala Coast, B and C (a) are in middle IQRs (interquartile ranges), D and E (b) are still in the upper quartile. Lower canopy heights (b, dark blue) adjoining B and C are pronounced. Cell colours represent the IQRs of all height data within the Koala Coast, i.e., intra-subregional Groups 1–7. IQR median heights are also shown.
Figure A11. Selected polygon groups from Koala Coast intra-subregional Groups 6 (a, low-canopy cohort) and 7 (b, high-canopy cohort). Letters B–E indicate areas of higher canopy These maps are similar to Figure A1, with features B–E recognisable, but, because they encompass the entire canopy height range within the Koala Coast, B and C (a) are in middle IQRs (interquartile ranges), D and E (b) are still in the upper quartile. Lower canopy heights (b, dark blue) adjoining B and C are pronounced. Cell colours represent the IQRs of all height data within the Koala Coast, i.e., intra-subregional Groups 1–7. IQR median heights are also shown.
Forests 12 01494 g0a11
Figure A12. Nonspatial algorithm (K-means, (a,b), K-medians (c,d), hierarchical (e,f)) cluster maps with 2 classes. Letters B–E indicate areas of higher canopy. Features B and C (a,c,e) are captured, but the areal extent of D and E (b,d,f) is not as well-defined as in Figure A11 Panel scores: K-means = 6, K-medians = 6, and Hierarchical = 7.
Figure A12. Nonspatial algorithm (K-means, (a,b), K-medians (c,d), hierarchical (e,f)) cluster maps with 2 classes. Letters B–E indicate areas of higher canopy. Features B and C (a,c,e) are captured, but the areal extent of D and E (b,d,f) is not as well-defined as in Figure A11 Panel scores: K-means = 6, K-medians = 6, and Hierarchical = 7.
Forests 12 01494 g0a12
Figure A13. Nonspatial algorithm (K-means, a,b; K-medians c,d; hierarchical, e,f) cluster maps with 3 classes. Letters B–E indicate areas of higher canopy. Features B and C are captured by all algorithms, but D is best captured by K-means (b). Panel scores: K-means = 10, K-medians = 10, and Hierarchical = 8.
Figure A13. Nonspatial algorithm (K-means, a,b; K-medians c,d; hierarchical, e,f) cluster maps with 3 classes. Letters B–E indicate areas of higher canopy. Features B and C are captured by all algorithms, but D is best captured by K-means (b). Panel scores: K-means = 10, K-medians = 10, and Hierarchical = 8.
Forests 12 01494 g0a13
Figure A14. Nonspatial algorithm (K-means, a,b; K-medians c,d; hierarchical, e,f) cluster maps with 4 classes. Letters B–E indicate areas of higher canopy. Features B and C are captured by all algorithms, but D and E are best captured by K-medians (c,d). Panel scores: K-means = 6, K-medians = 10, and Hierarchical = 8.
Figure A14. Nonspatial algorithm (K-means, a,b; K-medians c,d; hierarchical, e,f) cluster maps with 4 classes. Letters B–E indicate areas of higher canopy. Features B and C are captured by all algorithms, but D and E are best captured by K-medians (c,d). Panel scores: K-means = 6, K-medians = 10, and Hierarchical = 8.
Forests 12 01494 g0a14

Appendix A.4. Elbow Plots

Our assumption was that two, three, or four cluster class maps would have the greatest utility in terms of user practicality, but we generated maps with 2–10 clusters to allow use of the elbow method (Figure A15, Figure A16 and Figure A17) to estimate the optimal number of clusters [64]. The plots indicate that the elbow position is not clearly demarcated but occurs between two and four clusters; after four clusters, additional clusters explain less and less variance. Nonspatial methods show the greatest consistency between algorithms, and differences between spatial algorithms may be due to different polygon group shapes, i.e., compact, linear, and complex.
Figure A15. Plot of the percentage of variance explained as a function of the number of clusters generated by eight algorithms for the Sunshine Coast study site. This elbow plot indicates that 2–4 classes is the optimum, e.g., for nonspatial algorithms, 3 classes explain 82.7% (K-means), 79.2% (K-medians), and 80.9% (Hierarchical clustering) of the variance, and 4 classes explain 90.3% (K-means), 89.2% (K-medians), and 90.0% (Hierarchical clustering) of the variance. In comparison, spatially constrained algorithms (sc) explain a maximum of 56% (3 classes) and 63% (4 classes) of the variance.
Figure A15. Plot of the percentage of variance explained as a function of the number of clusters generated by eight algorithms for the Sunshine Coast study site. This elbow plot indicates that 2–4 classes is the optimum, e.g., for nonspatial algorithms, 3 classes explain 82.7% (K-means), 79.2% (K-medians), and 80.9% (Hierarchical clustering) of the variance, and 4 classes explain 90.3% (K-means), 89.2% (K-medians), and 90.0% (Hierarchical clustering) of the variance. In comparison, spatially constrained algorithms (sc) explain a maximum of 56% (3 classes) and 63% (4 classes) of the variance.
Forests 12 01494 g0a15
Figure A16. Plot of percentage of the variance explained as a function of the number of clusters generated by eight algorithms for the Koala Coast intra-subregional Group 6 polygon group. This elbow plot indicates that 2 to 3 classes is the optimum, e.g., for nonspatial algorithms, 2 classes explain 86.9% (K-means), 86.3% (K-medians), and 84.3% (Hierarchical clustering) of the variance. For spatial algorithms (sc), there is no clear elbow, with the Redcap algorithm explaining a maximum of 33.2% of the variances for 3 classes.
Figure A16. Plot of percentage of the variance explained as a function of the number of clusters generated by eight algorithms for the Koala Coast intra-subregional Group 6 polygon group. This elbow plot indicates that 2 to 3 classes is the optimum, e.g., for nonspatial algorithms, 2 classes explain 86.9% (K-means), 86.3% (K-medians), and 84.3% (Hierarchical clustering) of the variance. For spatial algorithms (sc), there is no clear elbow, with the Redcap algorithm explaining a maximum of 33.2% of the variances for 3 classes.
Forests 12 01494 g0a16
Figure A17. Plot of percentage of the variance explained as a function of the number of clusters generated by eight algorithms for the Koala Coast intra-subregional Group 7 polygon group. This elbow plot indicates that 2–4 classes is the optimum, e.g., for nonspatial algorithms, 3 classes explain 78.1% (K-means), 76.3% (K-medians), and 76.6% (Hierarchical clustering) of the variance, and 4 classes explain 86.1% (K-means), 85.4% (K-medians), and 84.3% (Hierarchical clustering) of the variance. For spatial algorithms (sc), Redcap explains a maximum of 32.1% (3 classes) and 38.3% (4 classes) of the variance.
Figure A17. Plot of percentage of the variance explained as a function of the number of clusters generated by eight algorithms for the Koala Coast intra-subregional Group 7 polygon group. This elbow plot indicates that 2–4 classes is the optimum, e.g., for nonspatial algorithms, 3 classes explain 78.1% (K-means), 76.3% (K-medians), and 76.6% (Hierarchical clustering) of the variance, and 4 classes explain 86.1% (K-means), 85.4% (K-medians), and 84.3% (Hierarchical clustering) of the variance. For spatial algorithms (sc), Redcap explains a maximum of 32.1% (3 classes) and 38.3% (4 classes) of the variance.
Forests 12 01494 g0a17

Appendix A.5. Visual Evaluation

To evaluate the utility of each algorithm, we used three criteria: (a) Does the algorithm identify the six features shown in Figure A1? (b) Is the spatial extent of each feature similar? (c) Is the feature internally contiguous, i.e., are all the cells connected? The maps (Figure A2, Figure A3, Figure A4, Figure A5, Figure A6, Figure A7, Figure A8, Figure A9, Figure A10, Figure A12, Figure A13 and Figure A14) were evaluated by an expert panel composed of one of the authors (DM), an ecologist specialising in forest mapping (DJ), and a GIS analyst (DS), i.e., potential height cluster map users. To provide a “hard” rule and reduce the subjectivity, we scored these criteria using a yes/no binary (Table A1).
Table A1. Procedure used to rate each cluster map. Binary scoring (yes or no) for each identified feature is designed to discourage subjectivity. Letters P, A, and C refer to codes in Table 12 and Table 13.
Table A1. Procedure used to rate each cluster map. Binary scoring (yes or no) for each identified feature is designed to discourage subjectivity. Letters P, A, and C refer to codes in Table 12 and Table 13.
CriterionEvaluationScore
Is the feature identified?yes1
P = presentno0
Is the areal extent similar?Within ±20% of the IQR * map1
A = areaMore than 25% different to the IQR map0
Is the feature internally contiguous?all cells joined (edges and/or corners)1
C = contiguityno better than the IQR map0
* IQR Interquartile range (Figure A1 and Figure A11).
The panellists first independently allocated scores, their disagreements were highlighted, and the panellists then re-evaluated their scores. These scores are shown with each of the above figures. The results for all algorithms are shown in Table 12; for nonspatial clustering, K-medians was assigned the highest rank for the two, three, and four classes. In comparison, the spatially constrained algorithms (Hierarchical, Skater, Redcap) performed poorly for the two and three classes but slightly better than the nonspatial algorithms for the four class. No single algorithm was judged to capture all the features, and the panellists had the most disagreement with feature D. The spatially constrained K-means and K-medians performed particularly poorly, identifying only one feature.
For two or more non-contiguous polygon groups, only the nonspatial algorithms (K-means, K-medians, and hierarchical clustering) are suitable clustering methods. For two sample polygon groups (Figure A11) selected from the highest (Group 7) and lowest (Group 6) canopy height cohorts, we used the entire Koala Coast canopy height model to evaluate these algorithms. Comparing multiple polygon groups within a subregion allows for the generation of multiscale clustering maps. Algorithms with two classes did not capture any features. K-medians captured features D and E for three and four classes and was assigned the highest rank (Table 13).

Appendix A.6. Cluster Height Differentiation

The mapped highest-height clusters may be significantly higher than the next-highest height, but these differences need to have practical significance, for example, a 4-m spread would be apparent to field observers and is well outside the error range of the LiDAR sensor. The spread between classes is directly influenced by the number of height classes, for example, for K-means, the average spread varies from 8.13 m (two classes) to 6.42 m (three classes) and 5.48 m (four classes) (Table A2). For the spatial algorithms, there is a trade-off between class separability and location, for example, spatial K-means and K-medians and Redcap have well-defined map clusters, but at the expense of height differentiation, which is less than 4 m.
Table A2. Height differences (metres) between the highest cluster median height and next-lowest cluster median height. Spatially constrained K-means and K-medians have the smallest spread between these height classes, while the other spatial algorithms have a larger spread but are still less than the nonspatial methods. n.a. = algorithm not applicable.
Table A2. Height differences (metres) between the highest cluster median height and next-lowest cluster median height. Spatially constrained K-means and K-medians have the smallest spread between these height classes, while the other spatial algorithms have a larger spread but are still less than the nonspatial methods. n.a. = algorithm not applicable.
AlgorithmClassesFigure A1Figure A1Figure A1Figure A11AverageRank
Map (a)Map (b)Map (c)Maps (a), (b)
K-means (ns)28.2910.145.868.258.132
K-medians (ns)26.288.984.685.426.345
Hierarchical Clustering (ns)210.6610.165.968.448.811
K-means (sc)27.504.991.64n.a.4.717
K-medians (sc)24.405.472.55n.a.4.148
Hierarchical Clustering (s)27.599.114.79n.a.7.163
Skater (s)27.488.552.98n.a.6.346
Redcap (s)27.128.554.77n.a.6.824
K-means (ns)37.617.354.955.766.423
K-medians (ns)35.926.223.543.934.906
Hierarchical Clustering (ns)38.698.365.894.666.902
K-means (sc)37.432.142.33n.a.3.978
K-medians (sc)37.623.811.25n.a.4.237
Hierarchical Clustering (s)36.737.247.10n.a.7.021
Skater (s)37.026.681.60n.a.5.105
Redcap (s)32.346.757.05n.a.5.384
K-means (ns)46.625.564.764.975.481
K-medians (ns)46.075.633.903.304.723
Hierarchical Clustering (ns)46.275.634.724.665.322
K-means (sc)45.562.762.61n.a.3.648
K-medians (sc)45.194.821.91n.a.3.976
Hierarchical Clustering (s)45.281.745.24n.a.4.095
Skater (s)46.362.834.91n.a.4.704
Redcap (s)43.252.585.19n.a.3.677

References

  1. McAlpine, C.A.; Bowen, M.E.; Callaghan, J.G.; Lunney, D.; Rhodes, J.R.; Mitchell, D.L.; Pullar, D.V.; Possingham, H.P. Testing Alternative Models for the Conservation of Koalas in Fragmented Rural-Urban Landscapes. Austral Ecol. 2006, 31, 529–544. [Google Scholar] [CrossRef]
  2. Department of Environment and Science (DES). Spatial Modelling for Koalas in South East Queensland v1.0; Biodiversity Assessment Team, Queensland Herbarium, Science and Technology Division, Department of Environment and Science: Brisbane, QLD, Australia, 2020. Available online: https://environment.des.qld.gov.au/wildlife/animals/living-with/koalas/mapping/koalamaps#toc-2 (accessed on 3 February 2021).
  3. Callaghan, J.; McAlpine, C.; Thompson, J.; Mitchell, D.; Bowen, M.; Rhodes, J.; De Jong, C.; Sternberg, R.; Scott, A. Ranking and Mapping Koala Habitat Quality for Conservation Planning on the Basis of Indirect Evidence of Tree-Species Use: A Case Study of Noosa Shire, South-Eastern Queensland. Wildl. Res. 2011, 38, 89–102. [Google Scholar] [CrossRef]
  4. Ellis, W.A.H.; Melzer, A.; Bercovitch, F.B. Spatiotemporal Dynamics of Habitat Use by Koalas: The Checkerboard Model. Behav. Ecol. Sociobiol. 2009, 63, 1181–1188. [Google Scholar] [CrossRef]
  5. Crowther, M.S.; Lunney, D.; Lemon, J.; Stalenberg, E.; Wheeler, R.; Madani, G.; Ross, K.A.; Ellis, M. Climate-Mediated Habitat Selection in an Arboreal Folivore. Ecography 2014, 37, 336–343. [Google Scholar] [CrossRef]
  6. McAlpine, C.A.; Rhodes, J.R.; Bowen, M.E.; Lunney, D.; Callaghan, J.G.; Mitchell, D.L.; Possingham, H.P. Can Multiscale Models of Species’ Distribution Be Generalized from Region to Region? A Case Study of the Koala. J. Appl. Ecol. 2008, 45, 558–567. [Google Scholar] [CrossRef]
  7. Moore, B.D.; Lawler, I.R.; Wallis, I.R.; Beale, C.M.; Foley, W.J. Palatability Mapping: A Koala’s Eye View of Spatial Variation in Habitat Quality. Ecology 2010, 91, 3165–3176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Dargan, J.R.; Moriyama, M.; Mella, V.; Lunney, D.; Crowther, M.S. The Challenge for Koala Conservation on Private Land: Koala Habitat Use Varies with Season on a Fragmented Rural Landscape. Anim. Conserv. 2019, 22, 543–555. [Google Scholar] [CrossRef]
  9. Hindell, M.A.; Lee, A.K. Habitat Use and Tree Preferences of Koalas in a Mixed Eucalypt Forest. Wildl. Res. 1987, 14, 349–360. [Google Scholar] [CrossRef]
  10. Santamaria, F.; Keatley, M.; Schlagloth, R. Does Size Matter: Tree Use by Translocated Koalas. Vic. Nat. 2005, 122, 4–13. [Google Scholar]
  11. Thompson, J. The Comparative Ecology and Population Dynamics of Koalas in the Koala Coast Region of South-East Queensland. Ph.D. Thesis, School of Integrative Biology, University of Queensland, Brisbane, QLD, Australia, 2006. [Google Scholar] [CrossRef] [Green Version]
  12. Januchowski, S.R.; McAlpine, C.A.; Callaghan, J.G.; Griffin, C.B.; Bowen, M.; Mitchell, D.; Lunney, D. Identifying Multiscale Habitat Factors Influencing Koala (Phascolarctos cinereus) Occurrence and Management in Ballarat, Victoria, Australia. Ecol. Manag. Restor. 2008, 9, 134–142. [Google Scholar] [CrossRef]
  13. Smith, A.G.; Mcalpine, C.; Rhodes, J.R.; Seabrook, L.; Baxter, G.; Lunney, D.; Bradley, A. At What Spatial Scales Does Resource Selection Vary? A Case Study of Koalas in a Semi-Arid Region. Austral Ecol. 2013, 38, 230–240. [Google Scholar] [CrossRef]
  14. Moore, B.; Foley, W. Tree Use by Koalas in a Chemically Complex Landscape. Nature 2005, 435, 488–490. [Google Scholar] [CrossRef] [PubMed]
  15. Matthews, A.; Lunney, D.; Gresser, S.; Maitz, W. Tree Use by Koalas (Phascolarctos cinereus) after Fire in Remnant Coastal Forest. Wildl. Res. 2007, 34, 84–93. [Google Scholar] [CrossRef]
  16. Moore, B.D.; Foley, W.J.; Wallis, I.R.; Cowling, A.; Handasyde, K.A. Eucalyptus Foliar Chemistry Explains Selective Feeding by Koalas. Biol. Lett. 2005, 1, 64–67. [Google Scholar] [CrossRef] [Green Version]
  17. Florence, R.G. Ecology and Silviculture of Eucalypt Forests, 1st ed.; CSIRO: Collingwood, VIC, Australia, 1996. [CrossRef]
  18. Phillips, S.; Callaghan, J.G. Tree Species Preferences of Koalas (Phascolarctos cinereus) in the Campbelltown Area South-West of Sydney, New South Wales. Wildl. Res. 2000, 27, 509–516. [Google Scholar] [CrossRef]
  19. Gutteridge Haskins & Davey (GHD). South East Queensland Koala Habitat Assessment and Mapping Project; GHD: Brisbane, QLD, Australia, 2009; Available online: https://qldgov.softlinkhosting.com.au/liberty/opac/search.do?queryTerm=GHDReport%3ASouthEastQueenslandkoalaHabitatAssessmentandMappingProjectandmode=ADVANCEDand=undefinedandmodeRadio=KEYWORDandoperator=ANDandactiveMenuItem=false (accessed on 24 June 2021).
  20. Parker, G.G.; Brown, M.J. Forest Canopy Stratification—Is It Useful? Am. Nat. 2000, 155, 473–484. [Google Scholar] [CrossRef]
  21. Jones, S.; Haywood, A.; Suárez, L.; Wilkes, P.; Woodgate, W.; Soto-Berelov, M.; Mellor, A.; Axelsson, C. Literature Review for Determining Optimal Data Primitives for Characterising Australian Woody Vegetation and Scalable for Landscape-Level Woody Vegetation Feature Generation, 2013. CRCSI Project Report. Available online: http://goo.gl/knst81 (accessed on 3 February 2021).
  22. Suarez, L.; Jones, S.D.; Haywood, A.; Wilkes, P.; Woodgate, W.; Soto-Berelov, M.; Mellor, A. Woody Vegetation Landscape Feature Generation from Multispectral and LiDAR Data (A CRCSI 2.07 Woody Attribution Paper). In Proceedings of the International Geoscience and Remote Sensing Symposium (IGARSS), Melbourne, VIC, Australia, 21–26 July 2013; pp. 262–265. [Google Scholar] [CrossRef]
  23. Coops, N.C.; Wulder, M.A.; Culvenor, D.S.; St-Onge, B. Comparison of Forest Attributes Extracted from Fine Spatial Resolution Multispectral and Lidar Data. Can. J. Remote Sens. 2004, 30, 855–866. [Google Scholar] [CrossRef]
  24. Arroyo, L.A.; Johansen, K.; Armston, J.; Phinn, S. Integration of LiDAR and QuickBird Imagery for Mapping Riparian Biophysical Parameters and Land Cover Types in Australian Tropical Savannas. For. Ecol. Manag. 2010, 259, 598–606. [Google Scholar] [CrossRef]
  25. Vogeler, J.C.; Cohen, W.B. A Review of the Role of Active Remote Sensing and Data Fusion for Characterizing Forest in Wildlife Habitat Models. Rev. Teledetección 2016, 45, 1–14. [Google Scholar] [CrossRef] [Green Version]
  26. Carrasco, L.; Giam, X.; Papȩs, M.; Sheldon, K.S. Metrics of Lidar-Derived 3D Vegetation Structure Reveal Contrasting Effects of Horizontal and Vertical Forest Heterogeneity on Bird Species Richness. Remote Sens. 2019, 11, 743. [Google Scholar] [CrossRef] [Green Version]
  27. Wallace, L.; Lucieer, A.; Watson, C.W. Evaluating Tree Detection and Segmentation Routines on Very High Resolution UAV LiDAR Data. IEEE Trans. Geosci. Remote Sens. 2014, 52, 7619–7628. [Google Scholar] [CrossRef]
  28. Geosciences Australia. ELVIS—Elevation and Depth—Foundation Spatial Data Version: 0.7.0. 2020. Available online: https://elevation.fsdf.org.au/ (accessed on 24 June 2021).
  29. Coops, N.C.; Hilker, T.; Wulder, M.A.; St-Onge, B.; Newnham, G.; Siggins, A.; Trofymow, J.A. Estimating Canopy Structure of Douglas-Fir Forest Stands from Discrete-Return LiDAR. Trees–Struct. Funct. 2007, 21, 295–310. [Google Scholar] [CrossRef] [Green Version]
  30. Alexander, C. Delineating Tree Crowns from Airborne Laser Scanning Point Cloud Data Using Delaunay Triangulation. Int. J. Remote Sens. 2015, 30, 3843–3848. [Google Scholar] [CrossRef]
  31. Zhen, Z.; Quackenbush, L.J.; Zhang, L. Impact of Tree-Oriented Growth Order in Marker-Controlled Region Growing for Individual Tree Crown Delineation Using Airborne Laser Scanner (ALS) Data. Remote Sens. 2014, 6, 555–579. [Google Scholar] [CrossRef] [Green Version]
  32. Jeronimo, S.M.; Kane, V.R.; Churchill, D.J.; Mcgaughey, R.J.; Franklin, J.F. Applying LiDAR Individual Tree Detection to Management of Structurally Diverse Forest Landscapes. J. For. 2018, 116, 336–346. [Google Scholar] [CrossRef] [Green Version]
  33. Zhang, W.; Cai, S.; Liang, X.; Shao, J.; Hu, R.; Yu, S.; Yan, G. Cloth Simulation-Based Construction of Pit-Free Canopy Height Models from Airborne LiDAR Data. For. Ecosyst. 2020, 7, 1–14. [Google Scholar] [CrossRef] [Green Version]
  34. Lovell, J.L.; Jupp, D.L.; Culvenor, D.S.; Coops, N.C. Using Airborne and Ground-Based Ranging Lidar to Measure Canopy Structure in Australian Forests. Can. J. Remote Sens. 2003, 29, 607–622. [Google Scholar] [CrossRef]
  35. Jenkins, R.B. Airborne Laser Scanning for Vegetation Structure Quantification in a South East Australian Scrubby Forest-Woodland. Austral Ecol. 2012, 37, 44–55. [Google Scholar] [CrossRef]
  36. Scarth, P.; Armston, J.; Goodwin, N. If You Climb Up A Tree, You Must Climb Down The Same Tree. But How High Was It? In Proceedings of the 15th Australasian Remote Sensing and Photogrammetry Conference, Alice Springs, NT, Australia, 13–19 September 2010; Available online: https://www.researchgate.net/publication/309211791_If_You_Climb_Up_A_Tree_You_Must_Climb_Down_The_Same_Tree_But_How_High_Was_It (accessed on 24 June 2021).
  37. Liao, Z.; Van Dijk, A.I.; He, B.; Larraondo, P.R.; Scarth, P.F. Woody vegetation cover, height and biomass at 25-m resolution across Australia derived from multiple site, airborne and satellite observations. Int. J. Appl. Earth Obs. Geoinf. 2020, 93, 102209, ISSN 0303-2434. [Google Scholar] [CrossRef]
  38. Wulder, M.A.; White, J.C.; Bater, C.W.; Coops, N.C.; Hopkinson, C.; Chen, G.; Wulder, M.A. Lidar Plots—A New Large-Area Data Collection Option: Context, Concepts, and Case Study. Can. J. Remote Sens. 2014, 38, 600–618. [Google Scholar] [CrossRef]
  39. Wilkes, P.; Jones, S.D.; Suarez, L.; Haywood, A.; Soto-Berelov, M.; Mellor, A.; Axelsson, C.; Woodgate, W. Deriving Metrics of Vertical Structure at the Plot Level for Use in Regional Characterisation of S.E. Australian Forests. CEUR Workshop Proc. 2012, 1328. Available online: https://researchrepository.rmit.edu.au/esploro/outputs/conferenceProceeding/Deriving-metrics-of-vertical-structure-at-the-plot-level-for-use-in-regional-characterisation-of-SE-Australian-forests/9921862129601341 (accessed on 24 June 2021).
  40. Saremi, H.; Kumar, L.; Turner, R.; Stone, C. Airborne LiDAR Derived Canopy Height Model Reveals a Significant Difference in Radiata Pine (Pinus Radiata, D. Don) Heights Based on Slope and Aspect of Sites. Trees-Struct. Funct. 2014, 28, 733–744. [Google Scholar] [CrossRef]
  41. Nguyen, T.H.; Jones, S.D.; Soto-berelov, M.; Haywood, A.; Hislop, S. Monitoring Aboveground Forest Biomass Dynamics over Three Decades Using Landsat Time-Series and Single-Date Inventory Data. Int. J. Appl. Earth Obs. Geoinf. 2020, 84, 101952. [Google Scholar] [CrossRef]
  42. Lee, A.C.; Lucas, R.M. A LiDAR-Derived Canopy Density Model for Tree Stem and Crown Mapping in Australian Forests. Remote Sens. Environ. 2007, 111, 493–518. [Google Scholar] [CrossRef]
  43. Wilkes, P.; Jones, S.D.; Suarez, L.; Haywood, A.; Woodgate, W.; Soto-berelov, M.; Mellor, A.; Skidmore, A.K. Understanding the Effects of ALS Pulse Density for Metric Retrieval across Diverse Forest Types. Photogramm. Eng. Remote Sens. 2015, 81, 625–635. [Google Scholar] [CrossRef]
  44. Scarth, P.; Armston, J.; Lucas, R.; Bunting, P. A Structural Classification of Australian Vegetation Using ICESat/GLAS, ALOS PALSAR, and Landsat Sensor Data. Remote Sens. 2019, 11, 147. [Google Scholar] [CrossRef] [Green Version]
  45. Office of Environment and Heritage (OEH). A Review of Koala Tree Use across New South Wales; Office of Environment and Heritage: Sydney, NSW, Australia, 2018. Available online: https://www.environment.nsw.gov.au/research-and-publications/publications-search/a-review-of-koala-tree-use-across-new-south-wales (accessed on 20 November 2020).
  46. Department of Agriculture, Water and the Environment (DAWE). Interim Biogeographic Regionalisation for Australia (IBRA). 2012. Available online: https://www.environment.gov.au/land/nrs/science/ibra#ibra (accessed on 24 February 2021).
  47. Department of Environment and Science (DES). Biodiversity Status of Remnant Regional Ecosystems-Queensland, Version 11. State of Queenxland, Department of Environment and Science. 2018. Available online: http://qldspatial.information.qld.gov.au/catalogue/custom/detail.page?fid=%7B8FDF54D2-654C-4822-8295-1D8E8E772373%7D (accessed on 21 April 2020).
  48. Department of Environment and Science (DES). Regional Ecosystem Description Database (REDD). Version 11. State of Queenxland, Department of Environment and Science. 2018. Available online: https://www.qld.gov.au/__data/assets/excel_doc/0021/91551/redd-v11-1-2019.xlsx (accessed on 21 April 2020).
  49. Kingston, M.; Turnbill, J.; Hall, P. Tweed Vegetation Management Strategy 2004; Ecograph: Limpinwood, NSW, Australia. Available online: https://www.google.com/Url?Sa=T&Rct=J&Q=&Esrc=S&Source=Web&Cd=&Cad=Rja&Uact=8&Ved=2ahukewiawuwyv6_Xahxy4nmbhw6qdrwqfjaaegqibhad&Url=Https%3a%2f%2fwww.Tweed.Nsw.Gov.Au%2fdownload.Aspx%3fpath%3d~%2fdocuments%2fenvironment%2fflora%2520and%2520fauna%2ftsc00219_Tweed_Vegetation_Management_Strategy_Volume_3_Of_3_Appendices.Pdf&Usg=Aovvaw3y4xkb16tpotzq6ljkzxxg (accessed on 24 June 2021).
  50. Fisher, M.; Body, M.; Gill, J. The Vegetation of the Coffs Harbour City Council LGA. A Report to Coffs Harbour City Council. Vegetation Community Profiles-Volume 2. 1996. Available online: https://www.coffsharbour.nsw.gov.au/environment/Plants-and-Animals/Pages/Class5VegetationMapforCoffsHarbour.aspx (accessed on 24 June 2021).
  51. Ryan, T.S. (Ed.) Technical Descriptions of Regional Ecosystems of Southeast Queensland; Queensland Herbarium, Science and Technology Division, Department of Environment and Science: Brisbane, QLD, Australia, 2012. Available online: https://www.publications.qld.gov.au/dataset/re-technical-descriptions/resource/1771f391-44b9-4d25-8315-f92033610a9b (accessed on 24 June 2021).
  52. ESRI. ArcMap (Version 10.7.1); Environmental Systems Research Institute (ESRI): Redlands, CA, USA, 2019; Available online: https://www.esri.com/en-us/arcgis/products/arcgis-desktop/overview (accessed on 12 May 2020).
  53. Ryan, T.S. (Ed.) Draft Technical Description of 12.11.3; Queensland Herbarium, Science and Technology Division, Department of Environment and Science: Brisbane, QLD, Australia, 2018.
  54. Neldner, V.J.; Wilson, B.A.; Dilleward, H.A.; Ryan, T.S.; Butler, D.W.; McDonald, W.J.F.; Addicot, E.P.; Appelman, C.N. Methodology for Surveying and Mapping Regional Ecosystems and Vegetation Communities in Queensland. Version 5.1; Queensland Herbarium, Queensland Department of Environment and Science, Brisbane, Australia. 2019. Available online: https://www.publications.qld.gov.au/dataset/redd/resource/6dee78ab-c12c-4692-9842-b7257c2511e4 (accessed on 3 March 2020).
  55. Isenburg, M. ArcGIS LAStools Toolbox for LiDAR Processing; Rapidlasso GmbH, Gilching, Germany. 2013. Available online: https://rapidlasso.com/2013/07/02/arcgis-lastools-toolbox-for-lidar-processing/ (accessed on 24 June 2021).
  56. Ruxton, G.D.; Beauchamp, G. Time for some a priori thinking about post hoc testing. Behav. Ecol. 2008, 19, 690–693. [Google Scholar] [CrossRef] [Green Version]
  57. McGaughey, R.J. FUSION/LDV: Software for LIDAR Data Analysis and Visualization, FUSION Version 4.0. US Dep. Agric. 2020, 170. Available online: http://scholar.google.com/scholar?hl=en&btnG=Search&q=intitle:FUSION+/+LDV+:+Software+for+LIDAR+Data+Analysis+and+Visualization#0 (accessed on 20 February 2021).
  58. Assunção, R.M.; Neves, M.C.; Câmara, G.; Da Costa Freitas, C. Efficient Regionalization Techniques for Socio-Economic Geographical Units Using Minimum Spanning Trees. Int. J. Geogr. Inf. Sci. 2006, 20, 797–811. [Google Scholar] [CrossRef]
  59. Waller, L.A. Detection of Clustering in Spatial Data. In The SAGE Handbook of Spatial Analysis; Sage Publications: London, UK, 2012; pp. 298–320. Available online: https://uk.sagepub.com/en-gb/eur/the-sage-handbook-of-spatial-analysis/book227940 (accessed on 24 June 2021). [CrossRef]
  60. Jacquez, G.M. Spatial Cluster Analysis. Handb. Geogr. Inf. Sci. 2008, 395–416. [Google Scholar] [CrossRef]
  61. Anselin, L.; Syabri, I.; Kho, K. GeoDa: An Introduction to Spatial Data Analysis. Geogr. Anal. 2006, 38, 5–22. [Google Scholar] [CrossRef]
  62. Guo, D. Regionalization with Dynamically Constrained Agglomerative Clustering and Partitioning (REDCAP). Int. J. Geogr. Inf. Sci. 2008, 22, 801–823. [Google Scholar] [CrossRef]
  63. Department of Environment and Science (DES). Environmental Protection Act 1994—High Value Regrowth. Department of Environment and Science (DES). 2012. Available online: http://qldspatial.information.qld.gov.au/catalogue/custom/search.page (accessed on 3 February 2021).
  64. Ketchen, D.J.; Shook, C.J. The Application of Cluster Analysis in Strategic Management Research: An Analysis and Critique. Strateg. Manag. J. 1996, 17, 441–458. [Google Scholar] [CrossRef]
  65. Pascual, C.; García-Abril, L.G.; García-Montero, L.; Martín-Fernández, S.; Cohen, W.B. Object-Based Semi-Automatic Ap-proach for Forest Structure Characterization Using Lidar Data in Heterogeneous Pinus Sylvestris Stands. For. Ecol. Manag. 2008, 255, 3677–3685. [Google Scholar] [CrossRef] [Green Version]
  66. Guo, X.; Coops, N.C.; Tompalski, P.; Nielsen, S.E.; Bater, C.W.; John Stadt, J. Regional Mapping of Vegetation Structure for Biodiversity Monitoring Using Airborne Lidar Data. Ecol. Inform. 2017, 38, 50–61. [Google Scholar] [CrossRef]
  67. Ackers, S.H.; Davis, R.J.; Olsen, K.A.; Dugger, K.M. The Evolution of Mapping Habitat for Northern Spotted Owls (Strix Occidentalis Caurina): A Comparison of Photo-Interpreted, Landsat-Based, and Lidar-Based Habitat Maps. Remote Sens. Environ. 2015, 156, 361–373. [Google Scholar] [CrossRef]
  68. Hagar, J.C.; Yost, A.; Haggerty, P.K. Incorporating LiDAR Metrics into a Structure-Based Habitat Model for a Canopy-Dwelling Species. Remote Sens. Environ. 2020, 236, 111499. [Google Scholar] [CrossRef]
  69. North, M.P.; Kane, J.T.; Kane, V.R.; Asner, G.P.; Berigan, W.; Churchill, D.J.; Conway, S.; Gutiérrez, R.; Jeronimo, S.; Keane, J.; et al. Cover of Tall Trees Best Predicts California Spotted Owl Habitat. For. Ecol. Manag. 2017, 405, 166–178. [Google Scholar] [CrossRef]
  70. Ficetola, G.F.; Bonardi, A.; Mücher, C.A.; Gilissen, N.L.; Padoa-Schioppa, E. How Many Predictors in Species Distribution Models at the Landscape Scale? Land Use versus LiDAR-Derived Canopy Height. Int. J. Geogr. Inf. Sci. 2014, 28, 1723–1739. [Google Scholar] [CrossRef]
  71. Seavy, N.E.; Viers, J.H.; Wood, J.K.; Seavy, N.E.; Viers, J.H.; Julian, K.W. Riparian Bird Response to Vegetation Structure: A Multiscale Analysis Using LiDAR Measurements of Canopy Height. Ecol. Appl. 2016, 19, 1848–1857. [Google Scholar] [CrossRef] [Green Version]
  72. Eyre, T.J. Regional Habitat Selection of Large Gliding Possums at Forest Stand and Landscape Scales in Southern Queensland, Australia. I. Greater Glider (Petauroides volans). For. Ecol. Manag. 2006, 235, 270–282. [Google Scholar] [CrossRef]
  73. Melville, G.; Stone, C.; Turner, R. Application of LiDAR Data to Maximise the Efficiency of Inventory Plots in Softwood Plantations. N. Z. J. For. Sci. 2015, 45, 9. [Google Scholar] [CrossRef] [Green Version]
  74. Dubayah, R.; Blair, J.B.; Goetz, S.; Fatoyinbo, L.; Hansen, M.; Healey, S.; Hofton, M.; Hurtt, G.; Kellner, J.; Luthcke, S.; et al. The Global Ecosystem Dynamics Investigation: High-Resolution Laser Ranging of the Earth’s Forests and Topography. Sci. Remote Sens. 2020, 1, 100002. [Google Scholar] [CrossRef]
Figure 2. Workflow implemented to examine the forest height variability at the regional, subregional, and intra-subregional scales for the nominated forest community (Regional Ecosystem 12.11.3). We assessed Herbarium field plots (CORVEG) to examine how well the published information captures the height variabilities within RE 12.11.3 at these scales.
Figure 2. Workflow implemented to examine the forest height variability at the regional, subregional, and intra-subregional scales for the nominated forest community (Regional Ecosystem 12.11.3). We assessed Herbarium field plots (CORVEG) to examine how well the published information captures the height variabilities within RE 12.11.3 at these scales.
Forests 12 01494 g002
Figure 3. RE 12.11.3 Koala Coast intra-subregional groups, numbered 1–7. We used spatial separation to define most groups and elevation to separate Groups 6 and 7.
Figure 3. RE 12.11.3 Koala Coast intra-subregional groups, numbered 1–7. We used spatial separation to define most groups and elevation to separate Groups 6 and 7.
Forests 12 01494 g003
Figure 4. Map 1: Location of the study site in the Sunshine Coast subregion (see, also, Figure 1). Map 2: 2006 aerial photo (Photo: QAP6183058.jpg, Source: https://qimagery.information.qld.gov.au/, accessed on 27 June 2018) showing the forest cover as near-contemporaneous with the ALS acquisition in 2009. Map 3: Focus study area (thick black outline, rainforest excluded from study) in a 1940 aerial photo (Photo: 1940jul08_MAP648_76044_2011.tif, Source: https://qimagery.information.qld.gov.au/, accessed on 19 August 2018). Thinner black lines outline the tenure boundaries. Labels show disturbance types from the interpretations of the 1940 photo. Areas in Mapleton National Park marked “Possibly Disturbed Forest” might reflect the topographic position (dryer northerly aspect) and not disturbance per se. Yellow squares in Maps 2 and 3 show locations of 50 ALS virtual plots.
Figure 4. Map 1: Location of the study site in the Sunshine Coast subregion (see, also, Figure 1). Map 2: 2006 aerial photo (Photo: QAP6183058.jpg, Source: https://qimagery.information.qld.gov.au/, accessed on 27 June 2018) showing the forest cover as near-contemporaneous with the ALS acquisition in 2009. Map 3: Focus study area (thick black outline, rainforest excluded from study) in a 1940 aerial photo (Photo: 1940jul08_MAP648_76044_2011.tif, Source: https://qimagery.information.qld.gov.au/, accessed on 19 August 2018). Thinner black lines outline the tenure boundaries. Labels show disturbance types from the interpretations of the 1940 photo. Areas in Mapleton National Park marked “Possibly Disturbed Forest” might reflect the topographic position (dryer northerly aspect) and not disturbance per se. Yellow squares in Maps 2 and 3 show locations of 50 ALS virtual plots.
Forests 12 01494 g004
Figure 5. ALS virtual plot (VP) maximum tree heights for 33 Qld Herbarium CORVEG plots and 205 VPs from identified subregional areas in the Southeast Queensland Bioregion (SEQ).
Figure 5. ALS virtual plot (VP) maximum tree heights for 33 Qld Herbarium CORVEG plots and 205 VPs from identified subregional areas in the Southeast Queensland Bioregion (SEQ).
Forests 12 01494 g005
Figure 6. Canopy height variations within the seven Koala Coast intra-subregional groups. Groups 2, 4, 5, and 7 are generally taller than Groups 1, 3, and 6. Heights below 15 m are typically near the polygon group edges and can be attributed to mapping errors or edge disturbances.
Figure 6. Canopy height variations within the seven Koala Coast intra-subregional groups. Groups 2, 4, 5, and 7 are generally taller than Groups 1, 3, and 6. Heights below 15 m are typically near the polygon group edges and can be attributed to mapping errors or edge disturbances.
Forests 12 01494 g006
Figure 7. Canopy height variations for the Sunshine Coast intra-subregional study area. The study site, with 50 virtual plots, shows similar distribution to CORVEG as a whole, but Quinlans (Q) is substantially lower than CORVEG, and all the datasets are lower than the 32 Sunshine Coast random virtual plots.
Figure 7. Canopy height variations for the Sunshine Coast intra-subregional study area. The study site, with 50 virtual plots, shows similar distribution to CORVEG as a whole, but Quinlans (Q) is substantially lower than CORVEG, and all the datasets are lower than the 32 Sunshine Coast random virtual plots.
Forests 12 01494 g007
Figure 8. ALS profile (southwest to northeast) across the Sunshine Coast local study site illustrating the forest canopy height differences between the regenerating forest (Quinlans) and mature/old growth forest (Mapleton National Park). Rainforest area is excluded from our research. Height information from Table 8. Transect length 1300 m, vertical exaggeration 7×.
Figure 8. ALS profile (southwest to northeast) across the Sunshine Coast local study site illustrating the forest canopy height differences between the regenerating forest (Quinlans) and mature/old growth forest (Mapleton National Park). Rainforest area is excluded from our research. Height information from Table 8. Transect length 1300 m, vertical exaggeration 7×.
Forests 12 01494 g008
Figure 9. Polygon groups representing compact (a), linear (b), and complex (c) shapes, with interquartile ranges (IQR) and interquartile median height. We used these maps to assess clustering algorithms. Higher-canopy spatial patterns are evident, for e.g., patches of higher canopy (A–E, light blue), and a canopy area of moderate height (F, dark green, some light green and light blue), surrounded by canopy with slightly lower height (light green). Heights appear more randomly distributed in other areas (light green, dark green) except for areas of lowest canopy height (dark blue). All maps use the same scale.
Figure 9. Polygon groups representing compact (a), linear (b), and complex (c) shapes, with interquartile ranges (IQR) and interquartile median height. We used these maps to assess clustering algorithms. Higher-canopy spatial patterns are evident, for e.g., patches of higher canopy (A–E, light blue), and a canopy area of moderate height (F, dark green, some light green and light blue), surrounded by canopy with slightly lower height (light green). Heights appear more randomly distributed in other areas (light green, dark green) except for areas of lowest canopy height (dark blue). All maps use the same scale.
Forests 12 01494 g009
Figure 10. Nonspatial K-medians (ac) and nonspatial hierarchical cluster maps (df) for individual polygons with 3 height classes. Letters A–F indicate areas of higher canopy as shown in Figure 9. Five of six features were identified, and cluster areas are similar except for feature D (c,f). Feature F was partially captured by K-medians (d).
Figure 10. Nonspatial K-medians (ac) and nonspatial hierarchical cluster maps (df) for individual polygons with 3 height classes. Letters A–F indicate areas of higher canopy as shown in Figure 9. Five of six features were identified, and cluster areas are similar except for feature D (c,f). Feature F was partially captured by K-medians (d).
Forests 12 01494 g010
Figure 11. Nonspatial K-medians (ac), nonspatial hierarchical (df) and Redcap [64] spatial hierarchical (gi) cluster maps for individual polygons (refer Figure 9) with 4 height classes. Letters A–F indicate areas of higher canopy as shown in Figure 9. All six features are identified, but features A, D, and F are not contiguous. In comparison, the spatially constrained algorithm Redcap (gi) identifies all the features as contiguous but assigns C, D, and F to a lower height class. Compared to Figure 9, Redcap overestimates D and underestimates E; A is also larger, because it also captures height differences reflecting the land use history of Quinlans compared to Mapleton National Park. Redcap also accurately captures the extent of an ecotone identified during fieldwork (F); this ecotone has a higher proportion of a preferred koala food tree than the surrounding area (dark blue), and which is absent from the national park.
Figure 11. Nonspatial K-medians (ac), nonspatial hierarchical (df) and Redcap [64] spatial hierarchical (gi) cluster maps for individual polygons (refer Figure 9) with 4 height classes. Letters A–F indicate areas of higher canopy as shown in Figure 9. All six features are identified, but features A, D, and F are not contiguous. In comparison, the spatially constrained algorithm Redcap (gi) identifies all the features as contiguous but assigns C, D, and F to a lower height class. Compared to Figure 9, Redcap overestimates D and underestimates E; A is also larger, because it also captures height differences reflecting the land use history of Quinlans compared to Mapleton National Park. Redcap also accurately captures the extent of an ecotone identified during fieldwork (F); this ecotone has a higher proportion of a preferred koala food tree than the surrounding area (dark blue), and which is absent from the national park.
Forests 12 01494 g011
Figure 12. Nonspatial K-medians cluster maps for the Group 6 (low canopy) polygon (a,c) and Group 7 (high canopy) polygon (b,d) using all the height data from Koala Coast subregional Groups 1–7 (Figure 3). Letters B–E indicate areas of higher canopy as shown in Figure 9. Features D and E are clearly identified but B and C are not, because they belong to the second- and third-highest height classes across the wider Koala Coast subregion. B and C are still significantly higher than the lowest height class (a,c, dark blue). Dark blue areas in (a,c) accurately capture areas of Mature Regrowth vegetation [63], outlined in black.
Figure 12. Nonspatial K-medians cluster maps for the Group 6 (low canopy) polygon (a,c) and Group 7 (high canopy) polygon (b,d) using all the height data from Koala Coast subregional Groups 1–7 (Figure 3). Letters B–E indicate areas of higher canopy as shown in Figure 9. Features D and E are clearly identified but B and C are not, because they belong to the second- and third-highest height classes across the wider Koala Coast subregion. B and C are still significantly higher than the lowest height class (a,c, dark blue). Dark blue areas in (a,c) accurately capture areas of Mature Regrowth vegetation [63], outlined in black.
Forests 12 01494 g012
Table 1. Previous studies investigating use of larger trees by the koala. DBH = diameter at breast height.
Table 1. Previous studies investigating use of larger trees by the koala. DBH = diameter at breast height.
MethodScaleFindingReference
Use vs. Availability by DBH classplotpositive correlation with food tree, negative with non-food trees[3]
tree girth and koala usekoala radio trackingpositive correlation with food tree size[4]
generalised linear mixed modelkoala radio trackingtaller non-eucalypts on hot days[5]
mixed-effects logistic model in 3 regions with 9 habitat variablestreetree size importance varies regionally. Noosa: Ranked 1 of 9. Port Stephens: 7 of 9. Ballarat: 6 of 9[6]
koalas observed in trees, generalised additive models, generalised additive mixed modelssite (7.6 ha)increasing DBH highly positively significant[7]
generalised linear mixed effects model, 3 seasonal surveystransectpositively associated with DBH, negative with tree height[8]
Use vs. Availability by DBH classtransectpositively associated with large trees[9]
2-way ANOVAkoala radio trackingpositively associated with larger trees within plots[10]
DBH, t-testplotpositively associated with larger trees within plots[11]
crown size, t-testplotpositively associated with crown size within plots[11]
Logistic regression and hierarchical partitioning of 12 variablesplotlarger DBH had strongest independent effect, but negative[12]
9 variables, mixed-effects resource selection probability functionsplot/landscapetree species and tree height most important factors at both scales[13]
Observed and expected tree visitation based on DBHtreekoalas use trees significantly larger than expected[14]
one-way ANOVAkoala radio trackinglarge trees selected disproportionately to availability[15]
Table 2. Available metadata for ALS data used in this paper. Point density was not included, but most datasets include ground point spacing. Following Reference [43], we regard all datasets to be suitable for our study. All data was downloaded from ELVIS (Elevation Information System) at https://elevation.fsdf.org.au/, accessed on 24 June 2021.
Table 2. Available metadata for ALS data used in this paper. Point density was not included, but most datasets include ground point spacing. Following Reference [43], we regard all datasets to be suitable for our study. All data was downloaded from ELVIS (Elevation Information System) at https://elevation.fsdf.org.au/, accessed on 24 June 2021.
Project and YearHeight (m)Swath Width (m)Ground pt Spacing (m)Pulse Rate (kHz)Sensor
Lismore 2017334019591.62 ALS80
Murwillumbah 2013212010241.51 ALS50
Tweed Heads 2013212010241.51 ALS50
Fraser Coast 2009150010842145.3ALS60
Inland Towns Stage 4 2012 1.74125–150Gemini/Pegasus
Inland Towns Stage 6 2014 300varies
Noosa-Gympie 2015 300Q1560
Sunshine Coast LGA 200810009241.74125–150Optech ALTM
Sunshine Coast LGA 201470011555.08300Trimble Harrier 68i INS
Logan LGA 2009 0.7126.2ALS50-II
Logan LGA 201314009921100ALS50-II
Moreton Bay LGA 2014 7383.2144ALS60
Scenic Rim 201130–1160 2.03 Riegl LMS-Q56
Southeast Queensland 2009 0.7126.2ALS50-II
Southeast Queensland 20142060240 ALS50, SN069, Riegl 680i
Toowoomba LGA 2010120010001 Optech ALTM Gemini
Table 3. Virtual ALS plot datasets are characterised by uneven sample sizes, non-normal data distribution, and unequal variances. We used both the Kolmogorov–Smirnov (with Lilliefors Significance Correction) and Anderson–Darling normality tests. Non-normal datasets are shown in blue, and highly unequal variances are shown in orange. * denotes maximum significance level reported by software package.
Table 3. Virtual ALS plot datasets are characterised by uneven sample sizes, non-normal data distribution, and unequal variances. We used both the Kolmogorov–Smirnov (with Lilliefors Significance Correction) and Anderson–Darling normality tests. Non-normal datasets are shown in blue, and highly unequal variances are shown in orange. * denotes maximum significance level reported by software package.
Regional GroupsSamplesKolmogorov–Smirnov Anderson–Darling
StatisticSig.Normality ANormality pVariance
CORVEG virtual plots330.1090.200 *0.4760.2342.2
Regional virtual plots (VPs)2050.0710.021.0640.0055.6
Subregional Groups
Brisbane VPs350.0720.200 *0.1760.9133.4
Gold Coast VPs310.1480.081.1000.0030.0
Koala Coast VPs300.1110.200 *0.4650.2415.1
Sunshine Coast VPs320.1210.200 *0.4410.2750.9
Tweed VPs400.1270.10.6960.0746.4
Western VPs370.1710.011.2490.0051.9
Intra-subregional Groups—Koala Coast
Koala Coast Group 1 VPs380.1140.200 *0.3350.4913.4
Koala Coast Group 2 VPs2410.0320.200 *0.4710.2413.2
Koala Coast Group 3 VPs6250.082<0.0018.9600.0025.9
Koala Coast Group 4 VPs10610.041<0.0012.6450.0013.9
Koala Coast Group 5 VPs16670.062<0.0019.2970.0023.6
Koala Coast Group 6 VPs3090.09<0.0012.7160.0037.2
Koala Coast Group 7 VPs20660.0220.021.4310.0015.3
Intra-subregional Groups—Sunshine Coast
Sunshine Coast VPs320.1210.200 *0.4410.2750.9
Study site all VPs500.1310.030.7530.0522.8
Mapleton National Park VPs220.1340.200 *0.3970.3421.6
Quinlans VPs280.090.200 *0.2950.577.9
Table 4. Procedure used by an expert panel to rate each cluster map. Binary scoring (yes or no) for each identified feature was designed to discourage any assessment subjectivity.
Table 4. Procedure used by an expert panel to rate each cluster map. Binary scoring (yes or no) for each identified feature was designed to discourage any assessment subjectivity.
CriterionEvaluationScore
Is the feature identified?yes1
no0
Is the areal extent similar?Within ±20% of the IQR * map1
More than 25% different to the IQR map0
Is the feature internally contiguous?all cells joined (edges and/or corners)1
no better than the IQR map0
* IQR Interquartile range.
Table 5. RE 12.11.3 mean, maximum, and minimum canopy heights from 66 CORVEG plots [53] compared with maximum heights obtained from the 33 CORVEG plots available to ALS, and these 33 plots were compared to the 205 ALS virtual plots within the Southeast Queensland Bioregion (SEQ), randomly sampled from six identified subregional areas within SEQ. Height differences were not significantly different (Welch’s t-test, df = 46.66, F = 1.205, p = 0.057). All measurements are in metres. Max- = maximum, Min. = minimum, Std = standard, Std Dev. = standard deviation.
Table 5. RE 12.11.3 mean, maximum, and minimum canopy heights from 66 CORVEG plots [53] compared with maximum heights obtained from the 33 CORVEG plots available to ALS, and these 33 plots were compared to the 205 ALS virtual plots within the Southeast Queensland Bioregion (SEQ), randomly sampled from six identified subregional areas within SEQ. Height differences were not significantly different (Welch’s t-test, df = 46.66, F = 1.205, p = 0.057). All measurements are in metres. Max- = maximum, Min. = minimum, Std = standard, Std Dev. = standard deviation.
n PlotsMean height Max. HeightMin. Height
6628 3715
Mean max. heightStd ErrorMax. max. heightMin. max. heightMedian heightRangeVarianceStd Dev.
3333.151.1351.3622.3632.3629.0042.246.50
20537.031.0357.0218.4736.2127.1835.095.92
Table 6. Variations in the canopy height between CORVEG and six geographically distinct subregional groups within the study area (Figure 1).
Table 6. Variations in the canopy height between CORVEG and six geographically distinct subregional groups within the study area (Figure 1).
MeasureCORVEGBrisbaneGold CoastKoala CoastSunshine CoastTweedWestern
Mean33.1536.9736.9731.3239.9839.5528.82
Standard Error1.130.980.980.711.261.081.18
Maximum51.3652.9750.8641.1257.0251.6454.99
Minimum22.3625.7929.0522.2129.9528.0118.47
Median32.3636.2136.1231.05539.0138.0526.95
Range29.027.221.818.927.123.636.5
Standard Deviation6.505.785.473.897.136.817.21
Variance42.2433.3929.9715.1050.8946.3551.94
Skewness0.870.381.210.430.420.211.59
Count33353130324037
Table 7. Differences in the maximum heights between the CORVEG virtual ALS plots and subregional virtual plots; significant p-values (ANOVA with post hoc tests, df = 231, F = 16.065) are shown in bold. The CORVEG plots were significantly different to the Sunshine Coast and Tweed subregional groups, and the Koala Coast and Western groups were significantly different to all the other subregional groups.
Table 7. Differences in the maximum heights between the CORVEG virtual ALS plots and subregional virtual plots; significant p-values (ANOVA with post hoc tests, df = 231, F = 16.065) are shown in bold. The CORVEG plots were significantly different to the Sunshine Coast and Tweed subregional groups, and the Koala Coast and Western groups were significantly different to all the other subregional groups.
Subregional Group(C)BGCKCSCTW
(CORVEG)10.1580.1610.8150.0030.0020.128
Brisbane 11.000<0.0010.4980.570<0.001
Gold Coast 1<0.0010.5000.573<0.001
Koala Coast 1<0.001<0.0010.543
Sunshine Coast 11.000<0.001
Tweed 1<0.001
Western 1
Table 8. Variations in the canopy heights for seven geographically distinct intra-subregional groups within the Koala Coast regional group. Compared to CORVEG, several groups have negatively skewed distributions and lower variances.
Table 8. Variations in the canopy heights for seven geographically distinct intra-subregional groups within the Koala Coast regional group. Compared to CORVEG, several groups have negatively skewed distributions and lower variances.
MeasureCORVEGGroup 1Group 2Group 3Group 4Group 5Group 6Group 7All Groups
Mean33.1523.4327.1323.7528.1328.7221.3329.4127.86
Standard Error1.130.590.230.200.110.120.350.090.06
Maximum51.3632.3536.9637.8946.8543.8033.2544.2646.85
Minimum22.3615.5415.561.2014.457.825.416.621.20
Median32.3623.6027.0524.1228.0129.2922.0229.4528.18
Range29.0016.8121.4036.6932.4035.9827.8437.6445.65
Standard Deviation6.503.663.635.083.734.866.103.914.96
Variance42.2413.4113.1825.8513.9023.6237.2315.2624.59
Skewness0.870.18−0.07−0.970.27−0.49−0.500.12−0.65
Count33382416251061166730920666607
Table 9. Differences in the maximum heights between the CORVEG and Koala Coast intra-subregional virtual plots; significant p-values (ANOVA with post hoc tests, df 1 = 7, df2 = 6032, F = 231.483) are shown in bold. We also compared the CORVEG plots to the pooled group plots (Welch’s t-test, df = 6038, F = 4.434). Except for Group 7, all the groups were significantly different to CORVEG, and Group 1 was not significantly different to Groups 3 and 6. n.a. = not applicable.
Table 9. Differences in the maximum heights between the CORVEG and Koala Coast intra-subregional virtual plots; significant p-values (ANOVA with post hoc tests, df 1 = 7, df2 = 6032, F = 231.483) are shown in bold. We also compared the CORVEG plots to the pooled group plots (Welch’s t-test, df = 6038, F = 4.434). Except for Group 7, all the groups were significantly different to CORVEG, and Group 1 was not significantly different to Groups 3 and 6. n.a. = not applicable.
Koala Coast Group(CORVEG)Group 1Group 2Group 3Group 4Group 5Group 6Group 7All Groups
(CORVEG)1<0.001<0.001<0.0010.0020.010<0.0010.0430.002
Group 1 1<0.0010.999<0.001<0.0010.063<0.001n.a.
Group 2 1<0.0010.004<0.001<0.001<0.001n.a.
Group 3 1<0.001<0.001<0.001<0.001n.a.
Group 4 10.009<0.001<0.001n.a.
Group 5 1<0.001<0.001n.a.
Group 6 1<0.001n.a.
Group 7 1n.a.
All Groups 1
Table 10. Regional CORVEG virtual plots (VPs) and subregional Sunshine Coast VP canopy heights compared to the Sunshine Coast intra-subregional study site. Quinlans (Q) has a lower canopy characterised by a lower standard deviation and low variance.
Table 10. Regional CORVEG virtual plots (VPs) and subregional Sunshine Coast VP canopy heights compared to the Sunshine Coast intra-subregional study site. Quinlans (Q) has a lower canopy characterised by a lower standard deviation and low variance.
MeasureCORVEGSC VPsStudy Site VPsMNP VPsQ VPs
Mean33.1539.9830.7134.1928.00
Standard Error1.131.260.670.990.53
Maximum51.3657.0245.2145.2134.97
Minimum22.3629.9523.2925.3123.29
Median32.3639.0129.6234.7528.09
Range29.0027.0721.9219.9011.68
Standard Deviation6.507.134.784.652.82
Variance42.2450.8922.8221.657.95
Skewness0.870.420.780.250.43
Count3332512228
Table 11. Differences in the maximum heights between the CORVEG and Sunshine Coast intra-subregional study sites; significant p-values (Welch’s t-test) are shown in bold.
Table 11. Differences in the maximum heights between the CORVEG and Sunshine Coast intra-subregional study sites; significant p-values (Welch’s t-test) are shown in bold.
Pairwise CombinationFdfp
CORVEG-Study site2.16354.060.088
CORVEG-Mapleton NP2.18952.680.503
CORVEG-Quinlans10.69245.130.001
Sunshine Coast VPs-Study site9.81648.520.001
Sunshine Coast VPs-Mapleton NP7.76851.900.001
Sunshine Coast VPs-Quinlans26.36541.530.001
Mapleton NP-Quinlans3.14332.740.001
Table 12. Visual evaluation by a panel composed of one author (DM), a GIS analyst (DS) and an ecologist (DJ); for maps produced for individual polygons by clustering algorithms. Maps were ranked according to how well six nominated features (A–F in Figure 10) were captured by each algorithm, each with 2, 3, and 4 classes. Panellists first independently allocated scores, disagreements were highlighted, and the panellists re-evaluated their scores. Key: P—feature identified (1), not identified (0); A—feature has similar area (1), smaller or larger (0); C—feature contiguous (1), not contiguous (0), ‘s’ prefix denotes spatially constrained algorithm. 1 1 1 = panellists agree that feature recognised by algorithm using all criteria- 0 0 0 feature not recognised using any criteria- Dark grey cells indicate panellists agree with individual criteria, 0/1 one panellist disagrees. For all cells, there was 86.3% agreement by all panellists. Panellists most often disagreed on whether feature D was captured. Nonspatial K-medians (all classes) has the highest overall score, followed by nonspatial K-means. For 4 classes, Redcap ranked slightly higher than nonspatial K-medians.
Table 12. Visual evaluation by a panel composed of one author (DM), a GIS analyst (DS) and an ecologist (DJ); for maps produced for individual polygons by clustering algorithms. Maps were ranked according to how well six nominated features (A–F in Figure 10) were captured by each algorithm, each with 2, 3, and 4 classes. Panellists first independently allocated scores, disagreements were highlighted, and the panellists re-evaluated their scores. Key: P—feature identified (1), not identified (0); A—feature has similar area (1), smaller or larger (0); C—feature contiguous (1), not contiguous (0), ‘s’ prefix denotes spatially constrained algorithm. 1 1 1 = panellists agree that feature recognised by algorithm using all criteria- 0 0 0 feature not recognised using any criteria- Dark grey cells indicate panellists agree with individual criteria, 0/1 one panellist disagrees. For all cells, there was 86.3% agreement by all panellists. Panellists most often disagreed on whether feature D was captured. Nonspatial K-medians (all classes) has the highest overall score, followed by nonspatial K-means. For 4 classes, Redcap ranked slightly higher than nonspatial K-medians.
Feature A B C D E F
AlgorithmClassesPACPACPACPACPACPACSumRank
K-means2111111111101111110162
K-medians2111111111101111111171
Hier. Clust.2101111111101101000123
s K-means210110110100000000066
s K-medians210110100010110100085
s Hier. Clust.210111100010110100094
s Skater210111100010110100094
s Redcap210111100010110100094
K-means3101111111100111000122
K-medians3111111111111111100161
Hier. Clust.3101111111100101000113
s K-means310110100010111100094
s K-medians310110100010111100094
s Hier. Clust.310111100000010100075
s Skater310111100010110100094
s Redcap310111100000010111094
K-means4101111111000101100115
K-medians4111111111100111101152
Hier. Clust.4110111111100101100124
s K-means410100000010111100076
s K-medians410100000010111100076
s Hier. Clust.4101111111101101101143
s Skater4101111111101111101152
s Redcap4111111111101101111161
Disagreements(Sum)2422 3 132113834 2
Table 13. Visual evaluation of the clustering algorithms using all height data within the Koala Coast subregional group (6007 cells), scored using the same methods as Table 12. Maps were ranked according to how well four nominated features (B–E in Figure A11) were captured by each algorithm, each with 2, 3, and 4 classes. For all cells, there was 84.3% agreement by all panellists, and the panellists disagreed most often on whether the feature area (column A) was captured. K-means and K-medians are ranked equally for 3 classes, but K-medians is ranked higher for 4 classes.
Table 13. Visual evaluation of the clustering algorithms using all height data within the Koala Coast subregional group (6007 cells), scored using the same methods as Table 12. Maps were ranked according to how well four nominated features (B–E in Figure A11) were captured by each algorithm, each with 2, 3, and 4 classes. For all cells, there was 84.3% agreement by all panellists, and the panellists disagreed most often on whether the feature area (column A) was captured. K-means and K-medians are ranked equally for 3 classes, but K-medians is ranked higher for 4 classes.
AlgorithmFeature B C D E SumRank
ClassesPACPACPACPAC
K-means210010010110162
K-medians210010010110162
Hier. Clust.211010010110171
K-means3110110111111101
K-medians3110110111111101
Hier. Clust.311011010110182
K-means410010010011163
K-medians4110110111111101
Hier. Clust.411011010011182
Disagreements(Sum) 3 14 2 151
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mitchell, D.L.; Soto-Berelov, M.; Jones, S.D. Regional Variation in Forest Canopy Height and Implications for Koala (Phascolarctos cinereus) Habitat Mapping and Forest Management. Forests 2021, 12, 1494. https://doi.org/10.3390/f12111494

AMA Style

Mitchell DL, Soto-Berelov M, Jones SD. Regional Variation in Forest Canopy Height and Implications for Koala (Phascolarctos cinereus) Habitat Mapping and Forest Management. Forests. 2021; 12(11):1494. https://doi.org/10.3390/f12111494

Chicago/Turabian Style

Mitchell, Dave L., Mariela Soto-Berelov, and Simon D. Jones. 2021. "Regional Variation in Forest Canopy Height and Implications for Koala (Phascolarctos cinereus) Habitat Mapping and Forest Management" Forests 12, no. 11: 1494. https://doi.org/10.3390/f12111494

APA Style

Mitchell, D. L., Soto-Berelov, M., & Jones, S. D. (2021). Regional Variation in Forest Canopy Height and Implications for Koala (Phascolarctos cinereus) Habitat Mapping and Forest Management. Forests, 12(11), 1494. https://doi.org/10.3390/f12111494

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop