Abstract
This study examined the usefulness of airborne light detection and ranging (LiDAR) data for estimating the individual tree health condition in Japanese mountain cherry (Cerasus jamasakura) in Yoshinoyama, Nara Prefecture, Japan. LiDAR variables that represented the ratio of lasers hitting tree components were calculated and their effectiveness was examined by relating them to the results of conventional field-based visual tree health assessments based on ordination, correlation analyses, and generalized linear models. The results showed that many of the LiDAR variables had significant correlations with the variables derived from visually evaluated tree health condition. In particular, the proportion of “only” returns, which represents the ratio of the lasers reflected from the crown surfaces, was the most effective for estimating total health condition in relation to the crown density, one of the key health indicators for representing physical properties. The individuals with large estimation errors had smaller crowns than the individuals with small errors, suggesting that sufficiently large crown sizes are important for more accurate estimations of the tree health condition using airborne LiDAR data.
- Airborne Laser Scanning
- Detrended Correspondence Analysis
- Hemispherical Photography
- Single Tree Level
- Tree Crown Density
- Tree Health Assessment
INTRODUCTION
Tree health condition is one of the most important parameters for monitoring and proper management of trees and forests which provide various ecosystem services (Ferretti 1997; Nielsen et al. 2014). Evaluation procedures of the tree crown health based on visual observations have been developed at the individual tree level, which are provided by the United States Department of Agriculture (Schomaker et al. 2007). Many studies have also evaluated the condition of crowns using visual procedures (Long et al. 1997; Schaberg et al. 2006; Scowcroft et al. 2007; Pontius and Hallett 2014). Similarly, in Japan a visual tree health evaluation method based on several indicators, including the tree form, crown dieback, and crown density, is used widely (Hasegawa et al. 1984; Kozawa and Kobayashi 1999; Imanishi et al. 2011; Iida et al. 2013). These visual methods are often subjective and differences between the evaluations made by observers have been reported (Mizoue and Dobbertin 2003; Nakajima et al. 2011). Passive optical remote sensing data from satellites or aerial images can provide an objective evaluation of tree condition for a large area by calculating vegetation indices, such as the normalized difference vegetation index, although the estimation accuracies are lower than visual assessments (Michez et al. 2016). The passive optical sensors cannot penetrate canopies and therefore fail to acquire precise 3-D (three-dimensional) information, and many of the vegetation indices become saturated at high levels of leaf densities (Chen and Cihlar 1996; Sasaki et al. 2016b).
Light detection and ranging (LiDAR) is an active remote sensing technology which can directly acquire 3-D structure information (Lefsky et al. 2002; Griffith et al. 2015; Sasaki et al. 2016a). The LiDAR system emits numerous laser pulses (tens of thousands of pulses per second) and records one or more discrete returns or a continuous waveform for each laser emission (Lim et al., 2003). LiDAR is less susceptible to weather conditions compared to optical sensors and it is effective for measuring objects with complex 3-D structures such as trees and forests (Wulder et al. 2008). Information about the objects hit by the lasers is converted into a large number of point cloud data, which all have x-, y-, and z-coordinate values, as well as attributes such as the intensity and return number. At the forest stand level, the canopy coverage and porosity have been estimated with high accuracies using variables based on the laser canopy hits and laser penetration rate derived from LiDAR data (Sasaki et al. 2008; Hopkinson and Chasmer 2009; Korhonen et al. 2011; Sasaki et al. 2016b). Thus, physical properties of trees such as their crown density and porosity may be estimated by applying these estimation methods at the individual tree level. Furthermore, the health condition of a tree can be estimated if it correlates with these physical indicators. Although many past studies identified individual trees in forests and estimated tree characteristics using LiDAR data (Popescu et al. 2003; Holmgren et al. 2008; Ørka et al. 2009), relatively few studies monitored trees outside forests (Schell et al. 2015), and none have associated the LiDAR variables with field-based visual tree health evaluations.
In this study, we aimed to verify the effectiveness of using airborne LiDAR data for estimating the individual tree health condition of Japanese mountain cherry (Cerasus jamasakura) planted in Yoshinoyama, Japan. In particular, we focused on the crown density, which is a visual health indicator, and used the crown porosity data obtained from hemispherical photographs to support the objectivity of the visual health assessment.
METHODS
Study Site
Yoshinoyama (34° 22′ N, 135° 52′ E), located in Nara Prefecture, Japan, is known for its cultural landscape containing flowering cherry trees. It has been designated as a UNESCO World Heritage Site because of the Sacred Sites and Pilgrimage Routes in the Kii Mountain Range. It has an elevation range of ca. 200 to 850 m above sea level, where cherry trees (mostly Japanese mountain cherry [Cerasus jamasakura]) have traditionally been planted for about 1,300 years (Imanishi et al. 2016). At present, the total planted area is about 50 hectares (0.5 km2). The four main planted areas are known as Shimosenbon, Nakasenbon, Kamisenbon, and Okusenbon, in order from the lowest elevation (Figure 1). Recently, local residents have become concerned about declined tree health condition, which is inferred to be caused by multiple factors, including soil water retention ability determined by the topography, infection of fungus, and barkfeeding damage by animals (Imanishi et al. 2012).
Field Survey
We selected 11 clusters comprising a total of 324 cherry trees (Figure 1). We aimed to include trees in various health conditions throughout the entire area of Yoshinoyama. In August 2011, each tree was visually classified by two experts according to one of four health ranks (1 = good to 4 = poor) taking into account the following eight indicators: tree vigor, tree form, branch growth, crown dieback, crown density, leaf shape and size, leaf color, and bark condition, in accordance with the Japanese traditional method (Shibata 2007; Table 1).
For each tree, the stem location was measured using differential GPS and total station; the crown projection area was depicted using total station by finding the apexes of branches and connecting them. Each crown projection area on the ground was transformed into a polygon with ArcGIS and its area was calculated. We excluded five trees with excessively small crowns that did not produce any meshes when the polygons were converted into 0.5-m meshes in the later processing step, as well as 31 cherry trees other than C. jamasakura, and we used 288 individual trees in the following analyses (35 in Shimosenbon, 163 in Nakasenbon, 60 in Kamisenbon, and 30 in Okusenbon). The means and standard deviations for height and diameter at breast height (dbh) in the selected 288 individuals were 11.4 m ± 2.8 m and 36.8 cm ± 16.7 cm, respectively. The means and standard deviations for height and dbh in the excluded five small individuals were 5.7 m ± 0.9 m and 8.9 cm ± 3.5 cm, respectively.
Among the 288 C. jamasakura trees, we took hemispherical photographs of 96 individuals (15 in Shimosenbon, 40 in Nakasenbon, 27 in Kamisenbon, and 14 in Okusenbon) and calculated the individual-level crown porosities. The trees with relatively upright stems and evenly distributed branches around the trunks were selected to minimize noises when comparing the field and the LiDAR data. Because the hemispherical photograph taken under the crown includes the trunk considerably, we took two photographs on opposite sides of the trunk for each tree and integrated each half not blocked by the trunk in the following analyses (Figure 2). The photographs were taken using a fish-eye lens (Sigma Circular Fisheye 4.5-mm F2.8 EX DC) attached to a digital camera (Canon EOS 60D), which was leveled on a tripod 1.3 m above the ground, under overcast sky conditions or immediately after sunset in August 2011.
LiDAR Data
The airborne LiDAR data were collected over the study area using a RIEGL LMS-Q560 sensor (Riegl Laser Measurement Systems GmbH, Horn, Austria) mounted on a helicopter platform on August 4, 2011, which was during the leaf-on season. This system projects near-infrared laser beams (1,550 nm) and records the full waveform of the reflection. The pulse frequency was 240 kHz and the scanning angle was ±30°. The flying height was 300 m above ground-level and the beam divergence was 0.5 mrad, thereby yielding a ground footprint measuring approximately 0.15 m in diameter. The flight speed was around 80 km/h−1. A back-and-forth flight pattern was used to survey the entire area. The full-waveform data from the entire area were converted into discrete points by detecting the local amplitude maxima of the waveforms (by Nakanihon Air Service Co. Ltd., Japan). All of the points obtained had x-, y-, and z-coordinate values and any of following attributes: first, intermediate, last, or only returns. First, intermediate, and last refer to the order in which the projected laser hit the tree components while passing through the crown. If all of the energy from a projected laser was returned at the same time, it was recorded as an only return.
Data Processing
Calculation of Crown Porosity
The individual-level crown porosity was calculated based on the hemispherical photographs of each of the 96 trees using CanopOn 2 software (Takenaka 2009), which binarized all the pixels in a hemispherical photograph to classify them as either sky or vegetation, where the hemispherical photograph was divided into 11 annulus rings split at zenith angles of 8.6°, 16.0°, 24.3°, 32.4°, 40.9°, 49.9°, 57.8°, 65.0°, 73.2°, 81.7°, and 90.0°. Thus, 11 gap fraction values were calculated according to the integrated annuli from 0–8.6° to 0–90.0°. For each of the two photographs per tree, we selected the zenith angle to completely cover the crown on the opposite side from the trunk (Figure 2), with the exception of abnormally long branches and overlaps with adjacent trees. The porosity value was calculated by integrating two hemicycles per tree, which we defined as the crown porosity.
LiDAR Data Processing
A 0.5-m mesh digital elevation model (DEM) was created by building a triangulated irregular network according to the public manual (Geospatial Information Authority of Japan 2006). For all the LiDAR points, we derived their height above the ground using the DEM data, and the points with height ≥ 1.3 m were classified as “vegetation” points and the remaining as “ground” points. The numbers of the “vegetation” points and “ground” points were exported with 0.5-m meshes. The numbers of the points with each return type, i.e., first, intermediate, last, and only returns, were also exported with 0.5-m meshes.
The numbers of points in each class and each return type within each crown polygon were calculated. The meshes on the polygon border line were dealt with as the meshes in the polygon if the center of the mesh was included in the polygon. Thus, the mean and standard deviation of the point density inside the crown were 167.0 and 81.4 points/m2 or m−2, respectively (n = 288). We calculated ten LiDAR variables based on previous studies (Morsdorf et al. 2006; Sasaki et al. 2008; Korhonen et al. 2011; Sasaki et al. 2016b) for each tree crown (Table 2), all of which represented the proportion of laser hits within the crowns using the information obtained from the point classes (vegetation or ground) and attributes (first, intermediate, last, or only).
Ordination and Correlation Analysis
For the 288 C. jamasakura trees, each 4-rank health indicator (Table 1) was converted into binary value: for example, if the “Tree vigor (Tvg)” for a tree is rank 3, the rank value is expressed by four binary values; Tvg1: 0, Tvg2: 0, Tvg3: 1, Tvg4: 0. This conversion was conducted for the data of the two experts respectively. Using the produced binary data, detrended correspondence analysis (DCA) was conducted using the vegan package in R version 3.4.1. (R Development Core Team). DCA is an improved method of correspondence analysis (CA) and corrects “arch effect” inherent to CA when using unimodal data (Hill and Gauch 1980). By this method, the health indicators for each rank were ordinated along the axes. The scores of the first axis by the two experts were averaged for each tree and used in the following analyses (hereinafter, referred to as “DCA axis 1 score”). Relationships between the DCA axis 1 score and the ten LiDAR variables (Table 2) were investigated using Pearson’s correlation coefficient.
For the 96 individuals with hemispherical photographs, relationships between the crown porosities in the photographs and the ten LiDAR variables were assessed using Pearson’s correlation coefficients. Correlation between the DCA axis 1 score and the crown porosity was also calculated.
Estimating the Tree Health Condition and Crown Porosity Using LiDAR Variables
We verified whether the DCA axis 1 score could be estimated using the LiDAR variables. A generalized linear model (GLM) was produced using a Gaussian distribution as an error distribution. To eliminate multicollinearity, if the correlation coefficients among the explanatory variables exceeded 0.5, that with a lower correlation with the response variable was eliminated in order. The conceivable regression models were arranged in the ascending order of Akaike’s information criterion (AIC), and the model with the lowest AIC was selected as the best model. The effectiveness of the models was verified using the root mean squared error (RMSE) and the Pearson’s correlation between the estimated values and the measured values.
The crown porosity in the hemispherical photographs was also estimated by the GLM using the LiDAR variables. The binomial distribution and logit were used as an error distribution and a link function, respectively. The relationships between the estimated and measured values were verified using Pearson’s correlation coefficient.
We analyzed influences of tree size on estimation accuracies by comparing the crown areas with the differences between estimated and measured values.
RESULTS
Results of Field Survey and Ordination
The result of DCA ordination was shown in Figure 3. All of the health indicators were arranged from 1 (good) to 4 (poor) in the ascending order of the DCA axis 1 score.
The median, maximum, and minimum values of the DCA axis 1 score for the trees were −0.12, 3.03, and −1.70, respectively (Figure 4). For the crown porosity values obtained from the hemispherical photographs, the median, maximum, and minimum values were 0.174, 0.533, and 0.082, respectively (Figure 4).
Correlation Analyses
The correlation between the DCA axis 1 score and the crown porosity was 0.704 (p < 0.001). Table 3 shows the correlations of the LiDAR variables with DCA axis 1 score and with crown porosity. Many of the LiDAR variables had significant coefficients with each variable. VOVFVO had the strongest correlations with both DCA axis 1 score and crown porosity (−0.600 and −0.617, respectively; p < 0.001).
Estimating Health Condition and Crown Porosity Using LiDAR Variables
Table 4 shows selected explanatory variables and the results obtained by the GLM. VOVFVO and VFO were selected to estimate the DCA axis 1 score. We only used VOVFVO for the GLM because it had the lowest AIC value. The RMSE value was 0.803, and the correlation coefficients between the estimated and measured values was 0.600 (p < 0.001).
VOVFVO and VLLO were selected to estimate the crown porosity. We only used VOVFVO for the GLM because it had the lowest AIC value. The RMSE value was 0.060 and the correlation coefficient between the estimated and measured values was 0.643 (p < 0.001).
Relationships Between the Estimation Accuracy and Crown Area
We analyzed whether the difference between estimated and measured values of the DCA axis 1 score and the crown porosity were influenced by the crown area (Figure 5). The individuals with larger crowns tended to have smaller errors. For the DCA axis 1 score, the individuals with crown areas larger than 100 m2 had errors smaller than 1. For the crown porosity, this tendency was similar, but one individual with a crown area of 90.3 m2 had by far the largest error (0.265). This tree, which had the largest crown porosity, was the outlier of the data (0.533, Figure 4). Another outlier individual with the second largest crown porosity (0.509, Figure 4) had the second largest error (0.187, Figure 5).
DISCUSSION
Possibility of Using Crown Density/Porosity to Estimate Tree Health
The DCA result indicated that all the health indicators showed similar arrangement along axis 1 (Figure 3), thereby suggesting that estimating some health indicators with high accuracy could facilitate the estimation of the total tree health condition. In addition, the DCA axis 1 score appeared to indicate the total health condition of the trees. The correlation coefficient between the DCA axis 1 score and the crown porosity was high (r = 0.704, p < 0.001), indicating objectivity of the traditional visual assessment made by experts. This result suggested that it is possible to estimate the total health condition by estimating the crown density/porosity using LiDAR data.
Relationships Between LiDAR Variables and Indicators Obtained in the Field
Of the LiDAR variables, VOVFVO, representing the ratio of “only” returns among the “vegetation” returns, had the strongest negative correlations with both the DCA axis 1 score and crown porosity (Table 3). It was presumed because healthy trees have closed crowns, thereby increasing the likelihood of lasers being reflected from the crown surfaces. According to a previous study, VOVFVO was effective for estimating the canopy porosity in a forest with simple structure, similar-aged trees, less canopy surface roughness, and few understory trees (Sasaki et al. 2016b). In this study site, only cherry trees have been planted, and there are very few other trees and shrubs under the crowns of the cherry trees, resulting in a simple woodland structure. In addition, we did not need to consider the gaps between tree crowns and the differences of forms and heights among adjacent trees because we calculated the porosity within the crown at the individual tree level. These conditions seemed to lead to results similar to a forest with a simple structure.
On the other hand, VALL, the complement of which has been verified as an effective variable for estimating the leaf area and canopy porosity in various types of forests (Riaño et al. 2004; Sasaki et al. 2008; Richardson et al. 2009; Sasaki et al. 2016b), had smaller correlation coefficients than VOVFVO with both of the two variables (Table 3), although they were significant. VALL counts all the points that hit the forest vegetation components above 1.3 m, including the sub-canopy and shrub layers, so it will increase the stability when measuring complex structures in forests (Sasaki et al. 2016b). However, in this study, it was hardly necessary to consider the vegetation components under the crown, which appeared to increase the noise caused by the lasers hitting the trunk.
GLM Results and Estimation Accuracy
The GLM estimates of the DCA axis 1 score and the crown porosity demonstrated the effectiveness of VOVFVO because for both of the two variables, the AIC values were the lowest, and there were significant correlations between the estimated and measured values when using only VOVFVO as an explanatory variable (Table 4).
According to the estimate of the DCA axis 1 score, the individuals with large crown areas (> 100 m2) had small errors (< 1) between the estimated and measured values (Figure 5). For the crown porosity, this tendency was similar except for the outlier individuals (Figure 5). In previous LiDAR studies, small trees were more difficult to measure compared to large trees when monitoring forest canopies (Zimble et al. 2003; Maltamo et al. 2004). The results of the present study are consistent with those of previous studies, thereby suggesting that airborne LiDAR is effective for individual tree health estimation, but a large crown size is necessary to obtain more accurate estimates.
CONCLUSION
The results of this study suggest that the total tree health condition can be estimated by the crown density or porosity because the DCA axis 1 score from the visually assessed health indicators had a strong correlation with the crown porosity determined from hemispherical photographs. Furthermore, the LiDAR variables had significant relationships with the tree health indicator; in particular, VOVFVO, which represents the proportion of the lasers reflected from the crown surface, seemed to be the most effective. These results may be applied to monitoring trees at the individual level, including trees other than C. jamasakura, because the LiDAR variables used in this study represent physical characteristics. This method has challenges if applied to densely forested areas with complex structures where individual tree delineation is difficult, but it will be acceptable for the diagnosis of trees in urban areas, such as parks or roadsides.
Acknowledgments
We acknowledge the Yoshino town office for providing the LiDAR data and Yoshinoyama-hoshou-kai for their cooperation with field observations.
- © 2019, International Society of Arboriculture. All rights reserved.