Next Article in Journal
A Simple Harmonic Model for FAPAR Temporal Dynamics in the Wetlands of the Volga-Akhtuba Floodplain
Next Article in Special Issue
Dynamics of Fractional Vegetation Coverage and Its Relationship with Climate and Human Activities in Inner Mongolia, China
Previous Article in Journal
Seasonal Separation of African Savanna Components Using Worldview-2 Imagery: A Comparison of Pixel- and Object-Based Approaches and Selected Classification Algorithms
Previous Article in Special Issue
Weak Environmental Controls of Tropical Forest Canopy Height in the Guiana Shield
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Remote Sensing
Volume 8
Issue 9
10.3390/rs8090766
remotesensing-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating Ladder Fuels: A New Approach Combining Field Photography with LiDAR

by
Heather A. Kramer
1,2,*,
Brandon M. Collins
3,
Frank K. Lake
4,
Marek K. Jakubowski
5,
Scott L. Stephens
1 and
Maggi Kelly
1
1
Ecosystem Sciences Division, Department of Environmental Science, Policy, and Management, 130 Mulford Hall, University of California, Berkeley, CA 94720, USA
2
SILVIS Lab, Department of Forest and Wildlife Ecology, University of Wisconsin–Madison, 1630 Linden Drive, Madison, WI 53706, USA
3
University of California Center for Fire Research and Outreach, College of Natural Resources, 130 Mulford Hall, University of California, Berkeley, CA 94720, USA
4
USDA Forest Service, Pacific Southwest Research Station, Orleans Ranger Station, Hwy 96/1 Ishi-Pishi Road, Orleans, CA 95556, USA
5
Marek Jakubowski Consulting, 560 Boulevard PI NE, Atlanta, GA 30308, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2016, 8(9), 766; https://0-doi-org.brum.beds.ac.uk/10.3390/rs8090766
Submission received: 28 April 2016 / Accepted: 12 September 2016 / Published: 17 September 2016
(This article belongs to the Special Issue Remote Sensing of Vegetation Structure and Dynamics)
Browse Figures
Versions Notes

Abstract

:
Forests historically associated with frequent fire have changed dramatically due to fire suppression and past harvesting over the last century. The buildup of ladder fuels, which carry fire from the surface of the forest floor to tree crowns, is one of the critical changes, and it has contributed to uncharacteristically large and severe fires. The abundance of ladder fuels makes it difficult to return these forests to their natural fire regime or to meet management objectives. Despite the importance of ladder fuels, methods for quantifying them are limited and imprecise. LiDAR (Light Detection and Ranging), a form of active remote sensing, is able to estimate many aspects of forest structure across a landscape. This study investigates a new method for quantifying ladder fuel in the field (using photographs with a calibration banner) and remotely (using LiDAR data). We apply these new techniques in the Klamath Mountains of Northern California to predict ladder fuel levels across the study area. Our results demonstrate a new utility of LiDAR data to identify fire hazard and areas in need of fuels reduction.

Graphical Abstract

1. Introduction

Remote sensing has been key to many areas of environmental science, and wildland fire is no exception [1,2,3,4]. Wildland fires consume homes, endanger lives, and consume a large number of taxpayer dollars every year [5,6]. During the past century, many forests in the Western United States have undergone unprecedented change in forest structure and fire behavior due to forest and fire management practices, including resource extraction and fire suppression [7]. One dramatic change in Western forests that historically experienced frequent, low-to moderate-severity fire has been the infilling of ladder fuels that help facilitate the movement of fire from the forest floor to tree canopies [7,8].
Imagery acquired from the Landsat program has informed the study of wildland fire severity [9,10]. Analyses of these Landsat-based fire severity assessments have shown that wildfires over the last 30 years have produced an uncharacteristically large proportions of high severity effects (where over 90% of trees are killed), which have been increasing [10,11]. These altered contemporary fire patterns are due to a number of factors, but considerable increases in surface and ladder fuels play an important role [12]. Ladder fuels allow fire to transition into overstory tree crowns by providing greater vertical fuel continuity, and fire burning ladder fuels also preheats canopy fuels that have not yet ignited [13]. Additionally, dense ladder fuels can make suppression more difficult and increase wildland firefighter exposure to hazardous conditions, especially when fire behavior shifts unexpectedly, by inhibiting escape to safety zones [14]. In addition to the biophysical influences of ladder fuels on fire behavior and fire management, they can also reduce habitat quality through decreasing accessibility and foraging efficiency for wildlife and tribal subsistence gathering [15].
Despite the importance of ladder fuels to fire behavior, effects, and firefighter safety, they have not been directly quantified except in a few cases [13,16,17], all of which are sampled on the ground, with no explicit connection to remote sensing. Because ladder fuels are below the canopy, passive remote sensing platforms are not able to capture their composition, except in very open forests. In most cases, fire models utilize a surrogate for ladder fuels, which is comprised of a combination of canopy base height (CBH) and fuel model (sometimes with an adjustment of fire behavior) [18]. CBH is the height above which there is enough fuel per unit volume to carry the fire upward. The fuel density necessary to carry fire is most commonly 0.012 kg·m−3 [19], but other thresholds have been suggested and used [20,21,22,23].
In addition to a variable fuel density threshold, estimations for CBH are commonly based on allometric equations, moving them further from a consistent direct measurement. Although active remote sensing platforms have been used to estimate CBH in a number of studies, the allometric derivation makes CBH an inherently difficult measure to estimate from remote sensing [24,25,26]. The allometric equations upon which CBH is derived use a combination of the following inputs: species, diameter at breast height (DBH), tree height, crown length, height to live crown base, crown ratio, and crown width [20,27,28,29,30]. While these allometric equations are likely appropriate for plantations with uniform conditions [25,31,32], they are unlikely to maintain accuracy in more natural stands because they often ignore small-diameter trees and shrubs that are critical components of ladder fuels [18,33,34]. This was also the case in Greece, where no satisfactory equation for CBH in Aleppo pine (Pinus halepensis Mill.) forests could be derived from basic stand measurements [35].
Alternatives to the standard methods to measure CBH include those given by Wilson and Baker [29], which estimated the CBH as the midpoint of minimum and average height to live crown base in multilayered stands; and Cruz, Alexander, and Wakimoto [36], which characterized the vertical gap in fuels as an alternative to CBH, and called their new metric the Fuel Strata Gap. Due to CBH’s arbitrary threshold and insensitivity to tree density, an alternate method was developed to estimate the torching index (open windspeed at which fire is expected to move into tree crowns under a given set of conditions) in the Fire and Fuels Extension of the Forest Vegetation Simulator [34]. Despite these inaccuracies, most wildland fire behavior and effects models still rely heavily on CBH [37].
Menning and Stephens [13] recognized the need for a more robust measure of ladder fuels, and developed a categorical system for identifying ladder fuels in forested systems as one of four classes consisting of all combinations of: (1) clumped/unclumped low aerial fuel; and (2) the presence/absence of a 2 m or greater vertical gap in fuel. Kramer et al. [38] related the ground plot measurements from Menning and Stephens [13] to aerial light detection and ranging (LiDAR), showing that LiDAR is able to differentiate between areas with dense and sparse ladder fuels. However, because the categorical, ground-based measures only span four structural distributions to describe every configuration of ladder fuels, relationships could not be built to more specifically estimate these fuels.
Methods to estimate the vertical profile of understory cover may also be adopted to quantify ladder fuels. Vegetation cover boards have long been used in wildlife studies to accurately capture the proportion of cover over different height bands within a few meters of the ground [39,40,41]. The cover estimates derived from these methods have been used to model wildlife habitat and hiding cover [42,43,44,45]. We used a similar approach to estimate ladder fuel cover, which includes shrub and herbaceous cover, but also small trees and lower tree limbs. We use the term ladder fuel cover to convey our intention to estimate the connectivity from the surface fuel layer to the canopy fuel layer. Although there is overlap in methodology with some previous efforts to estimate understory cover, our use of the term ladder fuel cover draws the distinction between our estimates and the more commonly used approach for estimating understory cover based on areal extent within a given plot.
Aerial LiDAR is an active remote sensing technology that is able to penetrate the upper canopy, allowing accurate estimations of forest attributes [46,47] including canopy structure [48,49], canopy bulk density [50], and leaf area index [51,52,53]. Airborne LiDAR presents a distinct advantage over terrestrial and space borne platforms for characterizing landscapes because it is able to cover a wide area, while still maintaining high pulse coverage. LiDAR has been applied to wildlife habitat [54,55] and wildland fire behavior modeling, providing fine-scale input across the landscape for terrain, canopy cover, tree height, crown bulk density, and fuel model [26,56,57,58] with varying degrees of accuracy. However, even though we can now model many inputs to fire behavior models, ladder fuels remain underrepresented. Because LiDAR is able to penetrate the canopy, it presents an excellent data source from which to derive additional information about these ladder fuels. LiDAR has been used to calculate metrics similar to ladder fuels, such as understory vegetation cover [59] and shrub biomass [60], although both of these studies necessitated intense sampling efforts, and neither captured the density and vertical continuity of ladder fuels as a whole.
Here we developed a ground-based methodology for predicting ladder fuel cover that provides more precise, quantitative estimates than the categorical system derived by Menning and Stephens, yet does not require the extensive field data collection efforts required by others [13,59,60]. We then investigate whether LiDAR can predict ladder fuel cover. Specifically, our main research questions for this project are: (1) Does the photo-banner methodology adequately sample ladder fuel cover within each plot? (2) How robustly can aerial LiDAR predict ladder fuel cover across a range of conditions?

2. Materials and Methods

2.1. Study Area

The Klamath River Basin is located in Northern California, and includes diverse land ownership and management. The 2622 ha (6480 ac) study area is located between Orleans, CA and Happy Camp, CA, centered at 41.5° North, 123.5° West (Figure 1). While LiDAR was collected over a 14,323 ha (35,394 ac) area, four focal areas are utilized by the overarching Western Klamath Restoration Partnership project to concentrate sampling effort in areas of maximum project interest for fuels assessment, and limiting the study area to 2622 ha. The four focal areas represent local forest diversity around geographically distant communities of homes. Topography is steep and complex, with elevation in our study area ranging between 160 and 520 m. Canopy cover ranged between 0 and 100%, but was often quite dense, averaging 90% across the study area (derived from the proportion of LiDAR first returns over breast height, or 1.37 m).
The Klamath has a Mediterranean climate: hot, dry summers, and cool, wet winters, although there is substantial variability in climate across the region [8,61]. Normal daily maximum and minimum temperatures range between 11 °C and 35 °C, and between 1 °C and 11 °C, respectively [8]. The majority of precipitation falls during the winter months, with an average of 143 cm every year [8]. The Klamath includes the highest diversity of conifer forests in North America [62]. Our study area contains areas of forest, shrub, and meadow, with the dominant tree species being tanoak (Notholithocarpus densiflorus) and Douglas-fir (Psedotsuga menziesii). It is also common to find pacific madrone (Arbutus menziesii) and canyon live oak (Quercus chrysolepis). The study area also includes rarer abundance of golden chinquapin (Chrysolepis chrysophylla), California black oak (Quercus kelloggi), California bay (Umbellularia californica), pacific dogwood (Cornus nuttallii), big-leaf maple (Acer macrophyllum), red alder (Alnus rubra), ponderosa pine (Pinus ponderosa), and sugar pine (Pinus lambertiana). Many of these tree species, such as California bay and multiple species of oak, are important sources for tribal gatherers, both historically and in the present [63]. Coniferous species, as well as a high diversity of oaks and other hardwoods present in this system, provide tribal community uses and important habitat for many species [8].
Fire has long been an important part of the Klamath, with many ignitions coming from lightning strikes [8] as well as the indigenous people [64]. Before the 1900s, fires were ignited and managed by indigenous tribes, ranchers, and foresters, maintaining a frequent, mixed severity fire regime [8,64,65]. In 1911, fires were discouraged by the Weeks Act, and by the 1920s, fires were actively suppressed in accessible areas, with remote fires remaining unchecked until 1945, with the advent of enhanced firefighting equipment [8]. Since the implementation of fire suppression, fire frequency decreased dramatically, resulting in a sharp decrease in annual area burned [8,61,66].
The change in fire frequency also led to a change in forest structure, with oaks becoming overtopped, and meadows getting encroached upon by conifers. This resulted in a decrease in forest complexity, with denser, more homogenous forests, and smaller forest openings [61,67].

2.2. Field Data

Field plots were established and measured during the summer of 2015. Plot locations were chosen based on a stratified random sample, which sorted stands in the combined four focal areas by: (1) management history—based on the CalVeg database [68]; (2) insolation; and (3) quadratic mean diameter of the stand. This yielded 10 different strata combinations, or stand types, which are shown in Figure 2.
More common stand types were assigned more plots, but rare types still contained at least three plots. Plot locations were limited to areas with <45% slope, between 25 m and 200 m from the nearest road, at least 25 m from all stand boundaries (where forest type changes), and with at least 100 m between plots. These limits were established to ensure the safety of the field crew and reduce travel time to plots to maximize plots sampled. Sixty plots were established and measured during the 2015 field season.
Plot centers were recorded with a Trimble Geo 7X GPS receiver with a Tempest antenna. One thousand points per plot were collected, followed by differential correction with a PDOP filter of 4. This resulted in plot center coordinates with an average horizontal precision of 87 cm. While the positional accuracies of the differential GPS and of the LiDAR data are high (sub-meter), their combined error creates an inconsistent misalignment between the two datasets. Because the accuracy of our predictive models relies heavily on the ground data lining up with LiDAR, we tested and adjusted plot center coordinates using mapped tall tree crowns. The distance and azimuth from plot center were measured to at least five tall, dominant conifers whose crowns reached above the main canopy and therefore would be apparent in the LiDAR canopy height model (CHM). Many of these identification trees fall outside the plot boundary. The distances to trees were measured using a TruPulse 360 laser rangefinder with an average error of ±30 cm, thus the local accuracy of the stem map’s geometry, as well as the accuracy of the GPS and LiDAR measurements, all fall below 1 m, and should coincide reasonably well. The stem maps—plot centers and locations of non-leaning trees—were incorporated into a GIS database based on the plot centers’ GPS positions, and then compared to the CHM generated by the LiDAR point cloud. All tree characteristics recorded in the field, including a unique ID number, tree height, species, DBH, and lean, were added to the GIS to help with the CHM comparison. The CHM was generated at 25 cm pixel resolution to correct for possible GPS-LiDAR misalignment.
In general, the tree positions in the CHM were close to the tree positions derived from the GPS; however, in a number of instances, the tree configuration apparent in the CHM was clearly different than the ground data. Plot locations were only shifted if an improved location fit was apparent. The combination of individual plot center and trees associated with that plot were moved manually and as a unit. Only x-y shift was applied (no stretching or warping) until the new fit between GPS trees and the CHM was optimal. We cross-checked the alignment of the shifted plot using tree heights as measured on the ground vs. the local maxima of the LiDAR-derived CHM. Table A1 shows the horizontal precision reported for each plot, as well as the distance of plot shift. Of the 60 plots, 45 were shifted an average of 3.3 m. The average shift distance for all plots, including those that were unshifted, was 2.5 m. Figure 3 shows an example case of plot shifting. In this example, some trees and part of the unshifted plot fall on the road, whereas trees in the shifted plot line up with the LiDAR-derived canopy height model. Each plot center was intentionally located at least 25 m from roads, showing that the plot center as measured by the GPS does not line up with the LiDAR data. Had the unshifted plot been used, model accuracy would have been reduced due to the inclusion of road in the LiDAR point cloud.
Plots used a 16.93 m radius to produce a plot area of 900 m2. While a suite of measurements was taken at each plot, the only measures used by this study were a set of four photographs. A banner, measuring 4 m tall × 0.5 m wide with 1 m vertical gradations, was placed at plot center and a photograph was taken from plot edge at each cardinal direction towards the center with a T70h Android Quadcore Rugged Tablet. The banner was similar to a vegetation cover board used in many wildlife studies [39,40]. The photos were taken from eye-level, so camera height varied between 1.5 and 2 m. Figure 4A shows an example of one of these photos in an area with low ladder fuel cover.
Ladder fuel cover in each 1 m vertical segment of each photo was separately analyzed by four different technicians. Although we attempted photo processing automation similar to that used by Bennett, Judd, and Adams [69], we found that we were not able to accurately assess photos where large foreground objects occluded the camera lens, such as that shown in Figure 4B. We therefore utilized a standard ocular estimation technique for estimating canopy cover as a proxy for estimating ladder fuel cover from each photograph. Cover estimates were binned into five classes between 0 and 100%, shown in Figure 5. One cover class was assigned to each 1 m vertical segment of banner that best approximated ladder fuel cover as shown by the photo for that 1 m tall × 0.5 m wide section of banner.
To maintain consistency between estimates, photo-analysis periods were limited to 45 min. At the start of each session, technicians had to pass a cover estimation calibration test—scoring five correct estimates in a row. Correct classifications for analysts’ tests were determined by an accurate assignment of a sample cover image, (one of the 30 shown in Figure 5), to a class designation (the six classes shown in Figure 5). To assist technicians with their estimates, they were provided with a printout of Figure 5, which served as a reference for their analysis.
To maintain a focus on ladder fuel fuels and eliminate the complications of perspective (a leaf near the camera lens sometimes covered an entire 1 m section of the banner), a 10% rule was established, whereby, if the banner area covered by a single object was more than 10% of a 1 m banner section, that area was removed from analysis. This addressed potential error introduced by tree boles in front of the banner, as well as leaves and smaller branches that were close to the camera lens and thus took up a larger than appropriate area on the banner. One example of this rule is demonstrated by Figure 4B (the raw photo) and Figure 4D (with unanalyzed area greyed out), where the banner is partially obscured by tanoak leaves close to the camera lens and the 10% rule was used. In this case, all large tanoak leaves were eliminated from analysis (shown in Figure 4D) because each leaf obscured over 10% of a banner 1 m segment. Ultimately, Figure 4A was analyzed as having 5% cover (the lowest class) in all height bins; Figure 4B was analyzed as having 20% cover from 0–1 m, 33% cover from 1–2 m, too little banner to determine cover from 2–3 m, and 27% cover from 3–4 m.
Because four independent technicians estimated the percent cover in each 1 m section of each photo, estimates often were not in perfect agreement. In cases where estimates differed by more than two consecutive bins (15% of samples), photos were re-assessed by a pair of technicians, who added a fifth measure to the set of estimates. Estimates for each 1 m segment of each photo were then averaged and used for further analysis.

2.3. LiDAR Data and Processing

Quantum Spatial collected aerial LiDAR between 23 May and 26 May 2015. A Leica ALS70 laser system mounted on a Cessna was used to collect LiDAR points utilizing a scan angle of ±15° from nadir. A pulse rate of 172 kHz, with unlimited returns per pulse (though often no more than 5) was used, averaging eight pulses per m2. IPAS TC v.3.1, ALS Post Processing Software v.2.75, Waypoint Inertial Explorer v.8.5, Leica Cloudpro v.1.2.1, and the TerraSolid software suite v.14 and v.15 were used to calculate point positions, classify points, and test spatial accuracy of the point cloud by Quantum Spatial [71]. The vendor completed real time kinematic (RTK) and post-processed kinematic (PPK) surveys and reported that average vertical and horizontal accuracy were 5.7 cm and 0.9 cm, respectively, based on the mean divergence of points from ground survey point coordinates (60 ground survey points were collected and compared to measured LiDAR ground points across the study area).
LiDAR metrics were extracted for each 900 m2 plot using LasTools [72]. Metrics included basic point statistics (six metrics), percentile heights (eight metrics), and relative percent cover values (40 metrics), yielding a total of 54 LiDAR-derived metrics (see Appendix Table A2 for a detailed breakdown). Relative percent cover values were calculated for a range of 20 strata each for first returns and all returns. A visual representation of this suite of values is shown in Figure 6. Rasters of each LiDAR-derived metric were also generated at a spatial resolution of 30 m to be used to make predictions across the landscape.

2.4. Statistical Analysis

Our workflow through plot and LiDAR processing is summarized in Figure 7. Average ladder fuel cover was calculated for each 1 m height band in each photo by averaging the estimates of the four photo-interpretation technicians. To determine whether plots were sufficiently sampled, the R software program [73] was used to calculate the standard error of ladder fuel cover among the four photos for a given plot, for each 1 m height stratum. This was done for a random sample of two-, three-, and four-photo samples to examine whether repeated sampling per plot decreased standard error and increased the reliability of this measure.
For all subsequent analyses, the four estimates for each 1 m height band (one from each cardinal direction) were averaged to produce a single estimate of ladder cover for each 1 m height band in each plot. These data were used to examine the distribution of ladder fuel cover in each 1 m height band across the study area. Ladder fuel cover estimates between 1 and 4 m were then averaged for each plot to produce a dependent variable, upon which to build a regression model using LiDAR. Cover below 1 m was left out of this measure since fuels in this range are commonly defined as surface fuels [74]. LiDAR is also less precise at differentiating ground points from low points within 1 m of the ground, and these points are often removed from analysis, though many different thresholds have been used [75,76,77,78].
Because LiDAR-derived independent variables were highly correlated, an iterative model building approach was taken. (1) The best regression models using one to four independent variables were chosen using the leaps package in R [79]; (2) Each model was evaluated for variable collinearity (correlation >0.6 in either Pearson or Spearman correlations was used to indicate collinear variables) and variable significance (p < 0.05). Only models where all independent variables were significant and none were collinear with one another were kept; (3) The leaps package was used to find the best non-collinear independent variable to add to each model; (4) The process was repeated for all potential models until independent variables were no longer all significant; (5) The model with the best adjusted R-squared value was chosen from these models; (6) This model was tested for heteroscedasticity (where the variability across the range of a regression is not consistent across the range of values) using the Breusch–Pagan test [80]. R-squared, root mean squared error (RMSE), and a 10-fold cross-validation error using the CVTools package in R [81] were calculated to evaluate the fit of the final model. We repeated this model building for shifted, as well as unshifted plots to determine the difference in model fit between the two methods for this study.
After a model was created, LiDAR-derived rasters of each independent variable used by the model were generated at 30 m resolution. This spatial resolution was chosen to match the plot size of 900 m2. A new 30 m raster of predicted ladder fuel cover was then generated for the entire study area.

3. Results

3.1. Assessment of Photo-Banner

To determine whether each plot had been sufficiently sampled, we calculated the standard error between two-, three-, and four-photo estimates of ladder fuel cover at the same height within a single plot. Table 1 shows the distribution of these standard error measures of increasing sampling density across all plots for each height band.
The variability of samples in the lowest height class (0–1 m) did not change much with increased sampling. With the lowest sampling density of only two photos per plot, standard error was highest for the highest height class, at 9.94. However, for all but the height class within 1 m of the ground, standard error decreased with every subsequent photo sample added. Since we do not have more than four photos per plot, it is difficult to speculate where standard errors level off. However, based on the decrease in standard error of less than 1% from three to four photos in all height strata above 1 m, it appears that four photos may be sufficient. Average standard error for plots over 1 m, with four photo samples per plot, was only 6.98%. Considering that the estimates of ladder fuel cover were derived from six classes (shown in Figure 5), with most classes representing a range of 20%, the low standard error was encouraging. These low values suggest reasonable consistency between measurements, even between only two samples per plot.
Average ladder fuel cover is summarized in Figure 8 for each height class, as well as a combination of height classes from 1–4 m. Average 1–4 m ladder fuel cover was used in subsequent analyses to build a predictive model from LiDAR-derived metrics. The majority of plots had 1–4 m ladder fuel cover between 21% and 45%, with a maximum value of 91%, and minimum value of 7%, showing a good representation of sparse ladder fuels, and a weak, though still present, representation of dense ladder fuels.

3.2. Linking LiDAR to Ground-Based Measures

Using the iterative model-building approach described, we built a multiple regression using the shifted and unshifted plot centers. Shifted plot centers produced the model shown in Equation (1), where independent variables included the percentage of all LiDAR returns between 1 and 8 m (COV1_8), the standard deviation of LiDAR point heights above two meters (STD), and the percentage of first return LiDAR points between 8 and 16 m (FCOV8_16). Table 2 shows the correlation matrix of each of these variables, as well as the percentage of all LiDAR returns between 1 and 4 m (COV1_4), and 4 and 8 m (COV4_8).
Ladder fuel cover = 20.41 + 0.873 × (COV1_8) − 1.73 × (STD) − 0.189 × (FCOV8_16)
For the above model, all independent variables, as well as the model as a whole, were significant at p < 0.05, with an R-squared value of 0.73, root mean squared error of 9.92, and 10-fold cross-validation error of 10.57. Figure 9 shows a scatterplot of the predicted versus measured values of ladder fuel cover, using Equation (1).
A model was also created using data from the unshifted plot locations. Model fit was nearly identical, with significance at p < 0.05, an R-squared value of 0.73, and 10-fold cross-validation error of 11.02. While percentage of all LiDAR returns between 1 and 8 m still emerged as the primary model driver, other variables composed the rest of the model: height of the 90th percentile of LiDAR returns above 2 m, percentage of LiDAR first returns over 8 m, and percentage of LiDAR returns over breast height (1.37 m).
Equation (1) was used to predict ladder fuel cover across the entire study area, and is shown in Figure 10. Predicted ladder fuel cover ranged from −8% to 99% across the four focal areas, with an average of 36.6% and standard deviation of 15.9%. Although some areas had predicted ladder fuel cover under 0%, these only represented 0.12% of the total study area. While most areas have moderate to low ladder fuel cover, a few areas have pockets with very high ladder fuel cover. The top 10% of ladder fuel cover in the study area encompassed predicted densities between 58% and 100%, and were concentrated in focal areas C and D. Table 3 displays the average ladder fuel cover for each focal area, and the proportion of each focal area with exceptionally dense ladder fuel cover (in the top 10% of predicted values for the study area). Average ladder fuel cover ranged between 31% and 40% across the four focal areas, but area D had, by far, the greatest proportion of dense ladder fuel cover, having over four times as much as area A where ladder fuel cover was over the 90th percentile. Since area D has the highest human population density, it is especially critical to identify these pockets of dense ladder fuel cover for fire planning and mitigation.

4. Discussion

We developed an improved ground-based method for quantifying ladder fuels using an ocular estimate of ladder fuel cover as a surrogate. LiDAR-derived variables included in our model were the percentage of all points between 1 and 8 m, the standard deviation of point heights above two meters, and the percentage of first return points between 8 and 16 m. The primary model driver was the percentage of all points between 1 and 8 m, a good representation of understory and midstory vegetation cover. This variable was the single best predictor of ladder fuel cover in models derived for both shifted and unshifted plots. It is curious that the most logical indicator of our field-derived ladder fuel cover, percentage of points between 1 and 4 m, was not a driver of the model. While the correlation between the LiDAR-derived 1–4 m and 4–8 m relative cover was not strong (see Table 2), the correlation of both of these variables with 1–8 m relative cover is high, indicating that this could better capture both the understory and midstory. This midstory may influence vegetative ladder fuels through light availability and the amount of competition with other plants. Light penetration can also influence the number of LiDAR pulses that reach the lower stratum with enough energy to register a return from understory vegetation. Less important independent variables are more contentious, since many LiDAR-derived variables are highly correlated (and one can become a better predictor than another with just a small change in plot data). However, standard deviation of point heights is a good indication of the overall forest structure, and has been used to predict forest structure [82]. Percentage of first return points is commonly used to estimate cover [60,75], so a stratum between 8 and 16 m is a good representation of upper canopy cover, since all but a few trees were under 16 m in height. While this is not a direct representation of ladder fuel cover, upper canopy cover could be important for light availability and LiDAR pulse penetration for the same reasons described above. Shifting plot centers to better match the LiDAR point cloud did not improve model accuracy, likely due to the large size of the plots and raster cells, and the small shift in most plots. However, we emphasize the importance plot shifting could have for smaller plots and still believe it to be a good check for larger plot shifts or those that could negatively influence the model.
This method is a more quantitative approach to ladder fuels measurement than the five-class system developed by Menning et al. [13], and is also a promising method for change detection. Because our estimation relied on calibrated field photographs, our methods could be repeated in plots over time or after treatment or fire. The resulting series of photos taken in field plots could be utilized to record a quantitative change in ladder fuel cover, but could also serve as a visual indication of change, with a consistent and apparent measure of scale. In our study, area ladder fuel cover was fairly consistent within most plots, with low standard error between estimates. However, to minimize standard error in all plots in our study area, our results suggested collecting at least three samples per plot to bring the average standard error in plots (for 1 m bands over 1 m) below 8% and bring the maximum standard error near 20% (the size of most of our ladder fuel cover bins).
Using the relative ladder fuel cover derived from this ground-based measure, we built a model that predicted this value from the LiDAR point cloud with reasonable precision. This landscape-scale prediction also identified areas where there are pockets of uncharacteristically dense ladder fuels, indicating dense ladder fuels. Although it is out of the scope of this study to further investigate these areas, they indicate locations where efforts to mitigate ladder fuel hazard could be prioritized. While our measurement is not a direct input for current fire models, we believe that it holds much potential for landscape monitoring, fuels management, cultural resource identification and protection, and the evaluation of wildland fire behavior. This method is robust, quantitative, repeatable, fairly inexpensive (10 person-minutes in the field plus 16 in the office per plot (assuming four independent assessments per photo)), and is based on direct observation (i.e., does not rely on allometry).

4.1. Immediate Implications for Managers

Managers can use our methods for ground-based ladder fuel characterization to establish a baseline surrogate for ladder fuels across their management area and return to these plots to remeasure ladder fuel cover after change, similar to photographic monitoring used for evaluating fuels treatment and wildland fire effects [83]. Although our photographic methods require some analysis in the office, they are efficient to implement in the field, and provide a visual reference for each plot for further analysis or evaluation of change over time. Mechanical and/or prescribed fire fuel reduction treatments are planned throughout this study area, likely including several areas that overlap our plots. Plot remeasurement is planned post-treatment to assist in evaluating treatment effectiveness. Not only is this a quantitative measure of change, but also a visual indication, since plot photos, with consistent scale, are byproducts of plot measurement. This will allow for robust quantification of change in ladder fuel cover at the plot scale, and if LiDAR data collection is repeated, at the landscape scale as well.
Knowledge of the locations of areas with especially dense ladder fuels (such as those with ladder fuel cover in the 90th percentile for the study area) can be invaluable for prioritizing fuel reduction to facilitate community safety and fire preparedness. If areas of dense ladder fuels are near communities or evacuation routes (especially along routes that represent the only evacuation corridor for residents), these areas of dense ladder fuel cover could be targeted for fuel reduction or flagged as potentially hazardous for crown fire if a wildland fire were to occur in the area. Making these high-risk areas known to residents and fire managers could improve safety on large fires that threaten these communities in the wildland urban interface. Information about dense pockets of ladder fuels can also be useful to firefighters on the ground, for risk assessment of crown fire, ease of cutting fireline, and accessibility to and mobility through these areas, especially if they are to be traversed to get to a safety zone. Our results directly inform land management, including treatment prioritization at the landscape scale within the Western Klamath Restoration Partnership’s Geographic Zones of Agreement [84].
Our work can augment work that addresses other aspects of wildland fire risk and modeling that can be predicted by LiDAR. Aerial LiDAR is an excellent predictor of topography (thus elevation, slope, and aspect), tree height, and canopy cover, all used extensively by wildland fire models [1]. Canopy bulk density has also been predicted by LiDAR [50,58,60,85]. Surface fuel model has been predicted in some areas [56,57,86], but still presents a challenge in denser forest types on steep, topographically complex terrain [26]. CBH has also been estimated from LiDAR [25,87], completing the suite of landscape metrics required by many wildland fire behavior models. However, because neither fuel model nor CBH is a direct measure of ladder fuels, we encourage managers and modelers to consider incorporating a metric that is better able to account for ladder fuels, both for land management and fire modeling. While ladder fuel cover is also not a perfect surrogate, it represents an alternative method of measuring ladder fuels that focuses on the fuels, rather than tree crowns, to estimate ladder fuels.
For wildlife that require understory structural complexity such as fishers, these data could also be valuable [88]. From a cultural perspective, which is especially important in the tribal communities surrounding the study area, these results will be used to mitigate potential impacts from fire near valuable cultural resources, both historic and present [64,89]. These resources could include archaeological, heritage, and tribal sacred sites, as well as legacy oaks and pines that provide a traditional food source, but are also critical habitat elements [90,91]. Thinning and prescribed fire treatments to enhance these cultural resources that are at elevated risk of intense wildland fire will be facilitated by our predictive map of ladder fuel cover across the landscape [92]. Because areas with lower densities of surface and ladder fuels facilitate enhanced accessibility, these areas are preferred by traditional gatherers, as suggested by Hummel and Lake [89]. Our work to predict ladder fuel cover could be used to identify especially open understories that may be valuable for identifying gathering locations for (non-timber) forest products and associated ecological services for the community, tribes and public [93,94].
Beyond the immediate implications of this work, this photographic technique, as well as the relationship we built between ground measures and the LiDAR point cloud, represents a new method to quantify an important fuel source that has been so elusive to previous fire behavior models that an allometric proxy was implemented in the place of a quantitative measure based on fuels actually present. Although Kane et al. [95] were not able to find a significant relationship between standard forest structure metrics and fire severity, we encourage further investigation into whether this estimation of ladder fuel cover has any predictive power over fire behavior and effects. It will be particularly useful to test whether these methods for estimating ladder fuel cover can be used to explain observed fire severity in areas where pre-fire LiDAR data are available.

4.2. Study Limitations

While photo interpretation methods were robust, we did not show that the estimate of ladder fuel cover relates to any specific measure on the ground, such as ladder fuel bulk density. Therefore, this model was only useful for predicting the relative ladder fuel cover, and should not be used to predict specific quantities of fuel. This method did not account for vertical continuity of fuels (such as a tree crown reaching all the way to the ground and providing a continuous “ladder” for fire to climb). Vertical continuity was addressed by Menning and Stephens [13], but remains a difficult parameter to incorporate, both for sampling methods and sampling density, and was out of the scope of this study. Photo interpretation methods also contained numerous inherent sources of error, including: (1) slope distorting the apparent height of vegetation; (2) objects in the foreground appearing disproportionately large; (3) using standard focal length settings, which could alter relative object sizes in the photo [96]; and (4) photo interpreter error. We also noticed a moderate degree of inaccuracy, averaging 2.5 m (up to 15 m for some plots), between the LiDAR and the ground-based measures of tree locations. While we attempted to mitigate this error by shifting plot center GPS coordinates based on the mapped height of tall trees in relation to plot center, there are likely still inaccuracies between the plot measured on the ground and the LiDAR point cloud extracted from each plot perimeter. Furthermore, only a single study area was used, introducing the potential for locational bias in our analysis. Because our study area was a single sample of the landscape, results from this study should not be applied to other areas without first collecting ground reference data, testing the relationship, and adjusting where necessary to ensure accurate model estimates.
This research also emphasizes the need for better accuracy assessment of plot center coordinates for LiDAR-based studies. LiDAR vertical accuracy was reported to be 5.7 cm. Despite our use of a high-accuracy GPS with antenna, collecting at least 1000 estimates per plot, using differential correction, and receiving an average horizontal precision of 87 cm, plot centers were still mis-aligned with the LiDAR data by an average of over 2.5 m, and up to at least 15 m. We urge those that work with LiDAR to assess how well recorded plot centers line up with the LiDAR point cloud and take measures to better match the plot to the LiDAR before undertaking further analyses. This is especially important when plots are located in heterogeneous landscapes or near roads, and when plot sizes are small (so a small shift in plot center equates to a dramatically different area of the LiDAR point cloud). While some LiDAR-derived metrics are not as sensitive to shifts in plot center location [97], many likely are, and plot accuracy should be better accounted for in all cases.

4.3. Future Research

This research demonstrated the importance of matching plot data to the LiDAR point cloud. We urge others to pursue further research into the best practices of measuring and addressing error when working with the LiDAR point cloud.
We are eager to further explore the ability of LiDAR to predict ladder fuel cover, and plan to work more extensively with land managers and community members in this area to discuss ways to best utilize the products from these analyses. We urge others to experiment with our methods, not only as a way to estimate ladder fuels, which are important to the field of wildland fire and tribal cultural resources [92], but also as a monitoring method for plot vegetation and fuel loading, as well as for the enhancement of access to and the quality of cultural resources.
We encourage others to test this metric against fire behavior to determine its influence on the likelihood of crown fire, and hope that fire behavior models will incorporate some version of this metric to further refine future wildland fire behavior and effects models.

5. Conclusions

Ladder fuels are an important piece of fire behavior and effects modeling. However, because they are difficult to quantify in the field, ladder fuels are indirectly measured through allometric equations that estimate canopy base height, leaving out the shrubs and small trees that are a large component of ladder fuels in many forests. Our work presents a new ground-based method for estimating ladder fuels and establishes a strong predictive link with aerial LiDAR. We were able to predict our ground-based estimates of ladder fuel from aerial LiDAR with an R-squared value of 0.73. We encourage others to explore this relationship in other areas and use our estimates of landscape scale ladder fuel cover to reduce wildfire risk to communities in our study area through targeted treatment. We also demonstrate that plot center accuracy must be independently tested, and urge LiDAR users to critically compare their plot locations to the LiDAR point cloud.

Acknowledgments

We thank the Western Klamath Restoration Partnership for sharing their LiDAR dataset and assisting with plot data collection, with specific thanks to Kenny Sauve, Jill Beckmann (Karuk Tribe, GIS analysts), and the Mid Klamath Watershed Council plot crew (Chris Root, Skip Lowrey, Eric Nelson, Shan Davis, and Tony Dennis) for their hard work establishing and collecting plot data. We acknowledge Deer Creek GIS for helping get this project started. We also thank Alexis Bernal and Aseem Singh for their assistance with photo interpretation. Much thanks to Martin Isenburg and LAStools, as well as the LASmoons grant, for use of the software program, LAStools. Joan Canfield, and Philip Riggan, and anonymous reviewers contributed comments that improved this paper. This work is supported by the Western Klamath Restoration Partnership, with funding through the US Forest Service, Karuk Tribe and Mid Klamath Watershed Council. This material is also based upon work supported by the National Science Foundation Graduate Research Fellowship Program under Grant No. DGE 1106400.

Author Contributions

Heather A. Kramer, Frank K. Lake, Marek K. Jakubowski, Brandon M. Collins, Maggi Kelly, and Scott L. Stephens conceived and designed the experiments; Heather A. Kramer and Frank K. Lake performed the experiments; Heather A. Kramer, Marek K. Jakubowski, and Brandon M. Collins analyzed the data; Frank K. Lake, Marek K. Jakubowski, Maggi Kelly, and Scott L. Stephens contributed analysis tools; and Heather A. Kramer wrote the paper.

Conflicts of Interest

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

Abbreviations

The following abbreviations are used in this manuscript:
LiDARLight Detection and Ranging
CBHCanopy base height
CHMCanopy height model
FIAForest Inventory and Analysis
RMSERoot mean squared error
COV1_8Percentage of all points between 1 and 8 m
STDStandard deviation of point heights above 2 m
FCOV8_16Percentage of first return points between 8 and 16 m
COV1_4Percentage of all LiDAR returns between 1 and 4 m
COV4_8Percentage of all LiDAR returns between 4 and 8 m
DBHDiameter at breast height

Appendix A

Table
Table A1. Reported plot center precision and distance of plot shift when stem-mapped trees were compared to the LiDAR-derived canopy height model are presented for each plot.
Table A1. Reported plot center precision and distance of plot shift when stem-mapped trees were compared to the LiDAR-derived canopy height model are presented for each plot.
Plot IDHorizontal Precision (m)Distance Shifted (m)
D0010.73.2
D1040.20
D1050.60.9
D0111.10
D1160.81.6
D1170.82.9
D1180.60
D1242.25.1
D1260.80
D1330.21.9
D1361.326.4
D01513.0
D15110
D0180.92.9
D20410
D2340.80
D2351.32.7
D2500.93.1
D2550.92.6
D0270.71.0
D0392.52.5
D0400.70.6
D0050.90
B0021.210.8
B02013.3
B0210.30
B0220.91.7
B0240.90
B0260.80.8
B0290.92.1
B0341.13.3
B0040.71.9
B0410.91.5
B0420.81.9
B0450.93.2
B0460.76.6
B0490.80
B0500.92.8
C0330.80
C1100.80
C1140.82.8
C1310.73.1
C13211.3
C13911.9
C0360.81.2
C1200.72.0
C13011.7
C1370.22.9
C1520.71.4
A0100.60.8
A0170.94.0
A2070.72.4
A2160.83.0
A0230.65.8
A2360.61.6
A23818.2
A04312.8
A0471.10
A0070.92.8
A0090.90
Table
Table A2. Details each LiDAR-derived variable.
Table A2. Details each LiDAR-derived variable.
Variable GroupIndividual Variable Description
Basic point statisticsMaximum point height
Average of point height above 2 m
Quadratic average of point height above 2 m
Skewness of point heights above 2 m
Standard deviation of point heights above 2 m
Kurtosis of point heights above 2 m
Cover by vertical strata
calculated for first and all returns
= # points in strata _
 # points in and below strata		Percentage of points between 1 and 2 m
Percentage of points between 2 and 3 m
Percentage of points between 3 and 4 m
Percentage of points between 2 and 4 m
Percentage of points between 1 and 4 m
Percentage of points between 1 and 8 m
Percentage of points between 2 and 8 m
Percentage of points between 4 and 8 m
Percentage of points between 2 and 16 m
Percentage of points between 4 and 16 m
Percentage of points between 8 and 16 m
Percentage of points between 2 and 32 m
Percentage of points between 4 and 32 m
Percentage of points between 8 and 32 m
Percentage of points between 16 and 32 m
Percentage of points over breast height (1.37 m)
Percentage of points over 2 m
Percentage of points over 8 m
Percentage of points over 16 m
Percentage of points over 32 m
Percentile HeightsHeight of 5th percentile of points above 2 m
Height of 10th percentile of points above 2 m
Height of 25th percentile of points above 2 m
Height of 50th percentile of points above 2 m
Height of 75th percentile of points above 2 m
Height of 90th percentile of points above 2 m
Height of 95th percentile of points above 2 m
Height of 99th percentile of points above 2 m

References

  1. Lefsky, M.A.; Cohen, W.B.; Parker, G.G.; Harding, D.J. Lidar remote sensing for ecosystem studies. Bioscience 2002, 52, 19–30. [Google Scholar] [CrossRef]
  2. Dubayah, R.O.; Drake, J.B. Lidar remote sensing for forestry. J. For. 2000, 98, 44–46. [Google Scholar]
  3. Arroyo, L.A.; Pascual, C.; Manzanera, J.A. Fire models and methods to map fuel types: The role of remote sensing. For. Ecol. Manag. 2008, 256, 1239–1252. [Google Scholar] [CrossRef] [Green Version]
  4. Hyyppä, J.; Hyyppä, H.; Inkinen, M.; Engdahl, M.; Linko, S.; Zhu, Y.-H. Accuracy comparison of various remote sensing data sources in the retrieval of forest stand attributes. For. Ecol. Manag. 2000, 128, 109–120. [Google Scholar] [CrossRef]
  5. Husari, S.; Nichols, T.H.; Sugihara, N.G.; Stephens, S.L. Fire and Fuel Management; University of California Press: Berkeley, CA, USA, 2006; pp. 444–465. [Google Scholar]
  6. US Forest Service. The Rising Cost of Fire Operations: Effects on the Forest Service’s Non-Fire Work; US Department of Agriculture: Washington, DC, USA, 2015; p. 16.
  7. Hessburg, P.F.; Agee, J.K.; Franklin, J.F. Dry forests and wildland fires of the inland Northwest USA: Contrasting the landscape ecology of the pre-settlement and modern eras. For. Ecol. Manag. 2005, 211, 117–139. [Google Scholar] [CrossRef]
  8. Skinner, C.N.; Taylor, A.H.; Agee, J.K. Klamath Mountains Bioregion. In Fire in California Ecosystems; Sugihara, N.S., Van Wagtendonk, J.W., Fites-Kaufman, J., Shaffer, K.E., Thode, A., Eds.; University of California Press: Berkeley, CA, USA, 2006; pp. 170–194. [Google Scholar]
  9. Miller, J.; Thode, A. Quantifying burn severity in a heterogeneous landscape with a relative version of the delta Normalized Burn Ratio (dNBR). Remote Sens. Environ. 2007, 109, 66–80. [Google Scholar] [CrossRef]
  10. Miller, J.; Knapp, E.; Key, C.; Skinner, C.; Isbell, C.; Creasy, R.; Sherlock, J. Calibration and validation of the relative differenced Normalized Burn Ratio (RdNBR) to three measures of fire severity in the Sierra Nevada and Klamath Mountains, California, USA. Remote Sens. Environ. 2009, 113, 645–656. [Google Scholar] [CrossRef]
  11. Miller, J.; Safford, H. Trends in wildfire severity: 1984 to 2010 in the Sierra Nevada, Modoc Plateau, and Southern Cascades, California, USA. Fire Ecol. 2012, 8, 41–57. [Google Scholar] [CrossRef]
  12. Agee, J.K.; Skinner, C.N. Basic principles of forest fuel reduction treatments. For. Ecol. Manag. 2005, 211, 83–96. [Google Scholar] [CrossRef]
  13. Menning, K.M.; Stephens, S.L. Fire climbing in the forest: A semiqualitative, semiquantitative approach to assessing ladder fuel hazards. West. J. Appl. For. 2007, 22, 88–93. [Google Scholar]
  14. Hirsch, K.; Martell, D. A review of initial attack fire crew productivity and effectiveness. Int. J. Wildland Fire 1996, 6, 199–215. [Google Scholar] [CrossRef]
  15. Norgaard, K.M. The Politics of fire and the social impacts of fire exclusion on the Klamath. Humboldt J. Soc. Relat. 2014, 36, 77–101. [Google Scholar]
  16. Prichard, S.J.; Sandberg, D.V.; Ottmar, R.D.; Eberhardt, E.; Andreu, A.; Eagle, P.; Swedin, K. Fuel Characteristic Classification System Version 3.0: Technical Documentation. USDA Forest Service; General Technical Report PNW-GTR-887; Pacific Northwest Research Station: Portland, OR, USA, 2013. [Google Scholar]
  17. Wright, C.S.; Ottmar, R.; Vihnanek, R. Stereo Photo Series for Quantifying Natural Fuels. Volume VIII: Hardwood, Pitch Pine, and Red Spruce/Basam Fire Types in the Northeastern United States; Seventh Symposium on Fire and Forest Meteorology, 2007; US National Wildfire Coordinating Group, National Interagency Fire Center: Boise, ID, USA, 2007; p. 91.
  18. Scott, J.H.; Reinhardt, E.D. Assessing Crown Fire Potential by Linking Models of Surface and Crown Fire Behavior; Research Paper RMRS-RP-29; USDA Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2001.
  19. Reinhardt, E.; Lutes, D.; Scott, J. FuelCalc: A method for estimating fuel characteristics. In Proceedings of the Fuels Management-How to Measure Success: Conference Proceedings, Portland, OR, USA, 28–30 March 2006; USDA Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2006. [Google Scholar]
  20. Sando, R.W.; Wick, C.H. A Method of Evaluating Crown Fuels in Forest Stands; RP-NC-84; USDA Forest Service: Ogden, UT, USA, 1972.
  21. Mitsopoulos, I.D.; Dimitrakopoulos, A.P. Canopy fuel characteristics and potential crown fire behavior in Aleppo pine (Pinus halepensis Mill.) forests. Ann. For. Sci. 2007, 64, 287–299. [Google Scholar] [CrossRef]
  22. Ottmar, R.D.; Vihnanek, R.E.; Wright, C.S. Stereo Photo Series for Quantifying Natural Fuels, Volume 1: Mixed-Conifer with Mortality, Western Juniper, Sagebrush, and Grassland Types in the Interior Pacific Northwest; US National Wildfire Coordinating Group, National Interagency Fire Center: Boise, ID, USA, 1998; p. 73.
  23. Fernandes, P.M. Combining forest structure data and fuel modelling to classify fire hazard in Portugal. Ann. For. Sci. 2009, 66, 1–9. [Google Scholar] [CrossRef]
  24. Erdody, T.L.; Moskal, L.M. Fusion of LiDAR and imagery for estimating forest canopy fuels. Remote Sens. Environ. 2010, 114, 725–737. [Google Scholar] [CrossRef]
  25. Andersen, H.-E.; McGaughey, R.J.; Reutebuch, S.E. Estimating forest canopy fuel parameters using LIDAR data. Remote Sens. Environ. 2005, 94, 441–449. [Google Scholar] [CrossRef]
  26. Jakubowski, M.K.; Guo, Q.; Collins, B.; Stephens, S.; Kelly, M. Predicting surface fuel models and fuel metrics using lidar and CIR imagery in a dense, mountainous forest. Photogramm. Eng. Remote Sens. 2013, 79, 37–49. [Google Scholar] [CrossRef]
  27. Brown, J.K. Weight and Density of Crowns of Rocky Mountain Conifers; Research Paper INT-197; USDA Forest Service, Intermountain Forest and Range Experiment Station: Ogden, UT, USA, 1978.
  28. Brown, J.K.; Johnston, C.M. Debris Prediction System; Fuel Science RWU 2104; USDA Forest Service, Intermountain Forest and Range Experiment Station: Missoula, MT, USA, 1976.
  29. Wilson, J.; Baker, P. Mitigating fire risk to late-successional forest reserves on the east slope of the Washington Cascade Range, USA. For. Ecol. Manag. 1998, 110, 59–75. [Google Scholar] [CrossRef]
  30. Cruz, M.G.; Alexander, M.E.; Wakimoto, R.H. Assessing canopy fuel stratum characteristics in crown fire prone fuel types of western North America. Int. J. Wildland Fire 2003, 12, 39–50. [Google Scholar] [CrossRef]
  31. McAlpine, R.; Hobbs, M. Predicting the height to live crown base in plantations of four boreal forest species. Int. J. Wildland Fire 1994, 4, 103–106. [Google Scholar] [CrossRef]
  32. Dean, T.J.; Cao, Q.V.; Roberts, S.D.; Evans, D.L. Measuring heights to crown base and crown median with LiDAR in a mature, even-aged loblolly pine stand. For. Ecol. Manag. 2009, 257, 126–133. [Google Scholar] [CrossRef]
  33. Hall, S.A.; Burke, I.C.; Box, D.O.; Kaufmann, M.R.; Stoker, J.M. Estimating stand structure using discrete-return lidar: An example from low density, fire prone ponderosa pine forests. For. Ecol. Manag. 2005, 208, 189–209. [Google Scholar] [CrossRef]
  34. Rebain, S.A. The Fire and Fuels Extension to the Forest Vegetation Simulator: Updated Model. Documentation; USDA Forest Service Internal Report; Forest Management Service Center: Fort Collins, CO, USA, 2010.
  35. Mitsopoulos, I.; Dimitrakopoulos, A. Estimation of canopy fuel characteristics of Aleppo pine (Pinus halepensis Mill.) forests in Greece based on common stand parameters. Eur. J. For. Res. 2014, 133, 73–79. [Google Scholar] [CrossRef]
  36. Cruz, M.G.; Alexander, M.E.; Wakimoto, R.H. Modeling the likelihood of crown fire occurrence in conifer forest stands. For. Sci. 2004, 50, 640–658. [Google Scholar]
  37. Hall, S.A.; Burke, I.C. Considerations for characterizing fuels as inputs for fire behavior models. For. Ecol. Manag. 2006, 227, 102–114. [Google Scholar] [CrossRef]
  38. Kramer, H.A.; Collins, B.M.; Kelly, M.; Stephens, S.L. Quantifying ladder fuels: A new approach using LiDAR. Forests 2014, 5, 1432–1453. [Google Scholar] [CrossRef]
  39. Nudds, T.D. Quantifying the vegetative structure of wildlife cover. Wildl. Soc. Bull. 1977, 5, 113–117. [Google Scholar]
  40. Jones, R.E. A board to measure cover used by prairie grouse. J. Wildl. Manag. 1968, 32, 28–31. [Google Scholar] [CrossRef]
  41. Griffith, B.; Youtie, B.A. Two devices for estimating foliage density and deer hiding cover. Wildl. Soc. Bull. 1988, 16, 206–210. [Google Scholar]
  42. Brown, J.D.; Benson, T.J.; Bednarz, J.C. Vegetation characteristics of Swainson’s warbler habitat at the white river national wildlife refuge, Arkansas. Wetlands 2009, 29, 586–597. [Google Scholar] [CrossRef]
  43. Duebbert, H.F.; Lokemoen, J.T. Duck nesting in fields of undisturbed grass-legume cover. J. Wildl. Manag. 1976, 40, 39–49. [Google Scholar] [CrossRef]
  44. Huegel, C.N.; Dahlgren, R.B.; Gladfelter, H.L. Bedsite selection by white-tailed deer fawns in Iowa. J. Wildl. Manag. 1986, 50, 474–480. [Google Scholar] [CrossRef]
  45. Hyde, E.J.; Simons, T.R. Sampling plethodontid salamanders: Sources of variability. J. Wildl. Manag. 2001, 65, 624–632. [Google Scholar] [CrossRef]
  46. Lefsky, M.A.; Cohen, W.B.; Acker, S.A.; Parker, G.G.; Spies, T.A.; Harding, D. Lidar remote sensing of the canopy structure and biophysical properties of Douglas-fir Western Hemlock Forests. Remote Sens. Environ. 1999, 70, 339–361. [Google Scholar] [CrossRef]
  47. Kelly, M.; Di Tommaso, S. Mapping forests with Lidar provides flexible, accurate data with many uses. Calif. Agric. 2015, 69, 14–20. [Google Scholar] [CrossRef]
  48. 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 2007, 21, 295–310. [Google Scholar] [CrossRef]
  49. Kane, V.R.; Bakker, J.D.; McGaughey, R.J.; Lutz, J.A.; Gersonde, R.F.; Franklin, J.F. Examining conifer canopy structural complexity across forest ages and elevations with LiDAR data. Can. J. For. Res. 2010, 40, 774–787. [Google Scholar] [CrossRef]
  50. Riaño, D.; Chuvieco, E.; Condés, S.; González-Matesanz, J.; Ustin, S.L. Generation of crown bulk density for Pinus sylvestris L. from lidar. Remote Sens.Environ. 2004, 92, 345–352. [Google Scholar] [CrossRef]
  51. Richardson, J.J. Estimating leaf area index from multiple return, small-footprint aerial lidar at the Washington Park Arboretum. Agric. For. Meteorol. 2008, 149, 1152–1160. [Google Scholar] [CrossRef]
  52. Jensen, J.; Humes, K.; Vierling, L.; Hudak, A. Discrete return lidar-based prediction of leaf area index in two conifer forests. Remote Sens. Environ. 2008, 112, 3947–3957. [Google Scholar] [CrossRef]
  53. Kane, V.R.; North, M.P.; Lutz, J.A.; Churchill, D.J.; Roberts, S.L.; Smith, D.F.; McGaughey, R.J.; Kane, J.T.; Brooks, M.L. Assessing fire effects on forest spatial structure using a fusion of Landsat and airborne LiDAR data in Yosemite National Park. Remote Sens. Environ. 2014, 151, 89–101. [Google Scholar] [CrossRef]
  54. Zhao, F.; Sweitzer, R.A.; Guo, Q.; Kelly, M. Characterizing habitats associated with fisher den structures in the Southern Sierra Nevada, California using discrete return lidar. For. Ecol. Manag. 2012, 280, 112–119. [Google Scholar] [CrossRef]
  55. Garcı, C.; Tempel, D.J.; Kelly, M. LiDAR as a tool to characterize wildlife habitat: California spotted Owl nesting habitat as an example. J. For. 2011, 109, 436–443. [Google Scholar]
  56. Mutlu, M.; Popescu, S.C.; Zhao, K. Sensitivity analysis of fire behavior modeling with LIDAR-derived surface fuel maps. For. Ecol. Manag. 2008, 256, 289–294. [Google Scholar] [CrossRef]
  57. Mutlu, M.; Popescu, S.; Stripling, C.; Spencer, T. Mapping surface fuel models using lidar and multispectral data fusion for fire behavior. Remote Sens. Environ. 2008, 112, 274–285. [Google Scholar] [CrossRef]
  58. Riano, D.; Meier, E.; Allgöwer, B.; Chuvieco, E.; Ustin, S.L. Modeling airborne laser scanning data for the spatial generation of critical forest parameters in fire behavior modeling. Remote Sens. Environ. 2003, 86, 177–186. [Google Scholar] [CrossRef]
  59. Wing, B.M.; Ritchie, M.W.; Boston, K.; Cohen, W.B.; Gitelman, A.; Olsen, M.J. Prediction of understory vegetation cover with airborne lidar in an interior ponderosa pine forest. Remote Sens. Environ. 2012, 124, 730–741. [Google Scholar] [CrossRef]
  60. Skowronski, N.; Clark, K.; Nelson, R.; Hom, J.; Patterson, M. Remotely sensed measurements of forest structure and fuel loads in the Pinelands of New Jersey. Remote Sens. Environ. 2007, 108, 123–129. [Google Scholar] [CrossRef]
  61. Taylor, A.H.; Skinner, C.N. Fire history and landscape dynamics in a late-successional reserve, Klamath Mountains, California, USA. For. Ecol. Manag. 1998, 111, 285–301. [Google Scholar] [CrossRef]
  62. Cheng, S. Forest Service Research Natural Areas in California; General Technical Report PSW-GTR-188; Pacific Southwest Research Station, Forest Service, U.S. Department of Agriculture: Albany, CA, USA, 2004.
  63. Davis, B.; Hendryx, M. Plants and the People: The Ethnobotany of the Karuk Tribe; Siskiyou County Museum: Yreka, CA, USA, 2004; p. 166. [Google Scholar]
  64. Lake, F.K. Historical and cultural fires, tribal management and research issue in Northern California: Trails, fires and tribulations. Interdiscip. Stud. Humanit. 2013, 5, 22. [Google Scholar]
  65. Whittaker, R.H. Vegetation of the Siskiyou mountains, Oregon and California. Ecol. Monogr. 1960, 30, 279–338. [Google Scholar] [CrossRef]
  66. Miller, J.; Skinner, C.; Safford, H.; Knapp, E.E.; Ramirez, C. Trends and causes of severity, size, and number of fires in northwestern California, USA. Ecol. Appl. 2012, 22, 184–203. [Google Scholar] [CrossRef] [PubMed]
  67. Skinner, C.N. Change in spatial characteristics of forest openings in the Klamath Mountains of northwestern California, USA. Landsc. Ecol. 1995, 10, 219–228. [Google Scholar] [CrossRef]
  68. USDA Forest Service—Pacific Southwest Region—Remote Sensing Lab. CalvegTiles_Ecoregions07_4. 2010. Available online: http://www.fs.fed.us/r5/rsl/clearinghouse/gis-download.shtml (accessed on 18 May 2015).
  69. Bennett, L.T.; Judd, T.S.; Adams, M.A. Close-range vertical photography for measuring cover changes in perennial grasslands. J. Range Manag. 2000, 53, 634–641. [Google Scholar] [CrossRef]
  70. Schomaker, M.E.; Zarnoch, S.J.; Bechtold, W.A.; Latelle, D.J.; Burkman, W.G.; Cox, S.M. Crown-Condition Classification: A Guide to Data Collection and Analysis; Department of Agriculture, Forest Service, Southern Research Station: Asheville, NC, USA, 2007; p. 78.
  71. Quantum Spatial. Lower Klamath Watersheds LiDAR & Digital Imagery: Technical Data Report Summary; QSI Environmental: Corvallis, OR, USA, 2015; p. 42. [Google Scholar]
  72. Isenburg, M. LAStools—Efficient Tools for LiDAR Processing. 2011. Available online: http://lastools.org (accessed on 16 October 2015).
  73. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2008. [Google Scholar]
  74. Keane, R.E.; Gray, K.; Bacciu, V.; Leirfallom, S. Spatial scaling of wildland fuels for six forest and rangeland ecosystems of the northern Rocky Mountains, USA. Landsc. Ecol. 2012, 27, 1213–1234. [Google Scholar] [CrossRef]
  75. Kane, V.R.; Lutz, J.A.; Roberts, S.L.; Smith, D.F.; McGaughey, R.J.; Povak, N.A.; Brooks, M.L. Landscape-scale effects of fire severity on mixed-conifer and red fir forest structure in Yosemite National Park. For. Ecol. Manag. 2013, 287, 17–31. [Google Scholar] [CrossRef]
  76. McGaughey, R. FUSION/LDV: Software for LIDAR Data Analysis and Visualization, Version 3.01; USDA Forest Service, Pacific Northwest Research Station, University of Washington: Seattle, WA, USA, 2012.
  77. Hudak, A.T.; Crookston, N.L.; Evans, J.S.; Falkowski, M.J.; Smith, A.M.S.; Gessler, P. Regression modeling and mapping of coniferous forest basal area and tree density from discrete-return LiDAR and multispectral satellite data. Can. J. Remote Sens. 2006, 32, 126–138. [Google Scholar] [CrossRef]
  78. Hudak, A.T.; Strand, E.K.; Vierling, L.A.; Byrne, J.C.; Eitel, J.U.; Martinuzzi, S.; Falkowski, M.J. Quantifying aboveground forest carbon pools and fluxes from repeat LiDAR surveys. Remote Sens. Environ. 2012, 123, 25–40. [Google Scholar] [CrossRef]
  79. Lumley, T.; Miller, A. Leaps: Regression Subset Selection. R package version 2.9. 2009. Available online: http://CRAN.R-project.org/package=leaps (accessed on 12 November 2014).
  80. Breusch, T.S.; Pagan, A.R. A simple test for heteroscedasticity and random coefficient variation. Econom. J. Econom. Soc. 1979, 47, 1287–1294. [Google Scholar] [CrossRef]
  81. Alfons, A. cvTools: Cross-Validation Tools for Regression Models. R Package Version 0.3.2. 2012. Available online: https://cran.r-project.org/web/packages/cvTools/index.html (accessed on 12 November 2014).
  82. Zimble, D.A.; Evans, D.L.; Carlson, G.C.; Parker, R.C.; Grado, S.C.; Gerard, P.D. Characterizing vertical forest structure using small-footprint airborne LiDAR. Remote Sens. Environ. 2003, 87, 171–182. [Google Scholar] [CrossRef]
  83. Vaillant, N.M.; Noonan-Wright, E.; Dailey, S.; Ewell, C.; Reiner, A. Effectiveness and longevity of fuel treatments in coniferous forests across California. JFSP Res. Proj. Rep. 2013, 57, 29. [Google Scholar]
  84. Harling, W.; Tripp, B. Western Klamath Restoration Partnership: A Plan for Restoring Fire Adapted Landscapes; Western Klamath Restoration Partnership; MKWC: Orleans, CA, USA, 2014; p. 57. [Google Scholar]
  85. Skowronski, N.S.; Clark, K.L.; Duveneck, M.; Hom, J. Three-dimensional canopy fuel loading predicted using upward and downward sensing LiDAR systems. Remote Sens. Environ. 2011, 115, 703–714. [Google Scholar] [CrossRef]
  86. Seielstad, C.A.; Queen, L.P. Using airborne laser altimetry to determine fuel models for estimating fire behavior. J. For. 2003, 101, 10–15. [Google Scholar]
  87. Popescu, S.C.; Zhao, K. A voxel-based lidar method for estimating crown base height for deciduous and pine trees. Remote Sens. Environ. 2008, 112, 767–781. [Google Scholar] [CrossRef]
  88. Zielinski, W.J.; Dunk, J.R.; Yaeger, J.S.; LaPlante, D.W. Developing and testing a landscape-scale habitat suitability model for fisher (Martes pennanti) in forests of interior northern California. For. Ecol. Manag. 2010, 260, 1579–1591. [Google Scholar] [CrossRef]
  89. Hummel, S.; Lake, F. Forest site classification for cultural plant harvest by tribal weavers can inform management. J. For. 2015, 113, 30–39. [Google Scholar] [CrossRef]
  90. Anderson, K. Tending the Wild: Native American Knowledge and the Management of California's Natural Resources; University of California Press: Oakland, CA, USA, 2005. [Google Scholar]
  91. Long, J.W.; Quinn-Davidson, L.; Goode, R.W.; Lake, F.K.; Skinner, C.N. Restoring California Black Oak to Support Tribal Values and Wildlife. In Proceedings of the 7th California Oak Symposium: Managing Oak Woodlands in a Dynamic World; Department of Agriculture, Forest Service, Pacific Southwest Research Station: Albany, CA, USA, 2015; pp. 113–122. [Google Scholar]
  92. Lake, F.K.; Long, J.W. Fire and Tribal Cultural Resources; Department of Agriculture, Forest Service, Pacific Southwest Research Station: Albany, CA, USA, 2014; pp. 173–186.
  93. Burger, J.; Gochfeld, M.; Pletnikoff, K.; Snigaroff, R.; Snigaroff, D.; Stamm, T. Ecocultural attributes: Evaluating ecological degradation in terms of ecological goods and services versus subsistence and tribal values. Risk Anal. 2008, 28, 1261–1272. [Google Scholar] [CrossRef] [PubMed]
  94. Barbour, R.J.; Singleton, R.; Maguire, D.A. Evaluating forest product potential as part of planning ecological restoration treatments on forested landscapes. Landsc. Urban Plan. 2007, 80, 237–248. [Google Scholar] [CrossRef]
  95. Kane, V.R.; Cansler, C.A.; Povak, N.A.; Kane, J.T.; McGaughey, R.J.; Lutz, J.A.; Churchill, D.J.; North, M.P. Mixed severity fire effects within the Rim fire: Relative importance of local climate, fire weather, topography, and forest structure. For. Ecol. Manag. 2015, 358, 62–79. [Google Scholar] [CrossRef]
  96. Hoiem, D.; Efros, A.A.; Hebert, M. Putting objects in perspective. Int. J. Comput. Vis. 2008, 80, 3–15. [Google Scholar] [CrossRef]
  97. Kramer, H.A.; Collins, B.M.; Gallagher, C.V.; Keane, J.; Kelly, M.; Stephens, S.L. Accessible LiDAR: Estimating large tree density for habitat identification. 2016; under review. [Google Scholar]
Remotesensing 08 00766 g001 1024
Figure 1. Klamath River study area between Orleans and Happy Camp, CA. Within the 14,323 ha LiDAR footprint, four focal areas, totaling 2622 ha (6480 ac), were chosen to sample the forest.
Figure 1. Klamath River study area between Orleans and Happy Camp, CA. Within the 14,323 ha LiDAR footprint, four focal areas, totaling 2622 ha (6480 ac), were chosen to sample the forest.
Remotesensing 08 00766 g001
Remotesensing 08 00766 g002 1024
Figure 2. A flowchart showing the stratification of plots is shown.
Figure 2. A flowchart showing the stratification of plots is shown.
Remotesensing 08 00766 g002
Remotesensing 08 00766 g003 1024
Figure 3. The original plot location and new, shifted plot location are shown with the locations of mapped tall trees in relation to plot center. Tree heights are shown with the shifted tree locations.
Figure 3. The original plot location and new, shifted plot location are shown with the locations of mapped tall trees in relation to plot center. Tree heights are shown with the shifted tree locations.
Remotesensing 08 00766 g003
Remotesensing 08 00766 g004 1024
Figure 4. Examples of photos used to analyze ladder fuel cover (A,B); and their respective analysis boundaries (C,D). Areas left out of analysis are faded out in (C,D), whereas areas included in analysis are outlined in black. Figure 4A,B shows plots with low ladder fuel cover. Figure 4A’s banner is unobstructed, while the banner in Figure 4B is partially obscured by tanoak leaves close to the camera lens, leading to only part of the banner being used for analysis, shown in Figure 4D.
Figure 4. Examples of photos used to analyze ladder fuel cover (A,B); and their respective analysis boundaries (C,D). Areas left out of analysis are faded out in (C,D), whereas areas included in analysis are outlined in black. Figure 4A,B shows plots with low ladder fuel cover. Figure 4A’s banner is unobstructed, while the banner in Figure 4B is partially obscured by tanoak leaves close to the camera lens, leading to only part of the banner being used for analysis, shown in Figure 4D.
Remotesensing 08 00766 g004
Remotesensing 08 00766 g005 1024
Figure 5. Binned cover estimates, used to approximate ladder fuel cover. A canopy cover image series was adapted from Forest Inventory and Analysis (FIA) protocol [70], and binned into six cover classes. The cover estimates for each image set is shown at the bottom, and the average cover in each bin is shown at the top of each class. This was also the diagram used to assist technicians in assigning cover estimates.
Figure 5. Binned cover estimates, used to approximate ladder fuel cover. A canopy cover image series was adapted from Forest Inventory and Analysis (FIA) protocol [70], and binned into six cover classes. The cover estimates for each image set is shown at the bottom, and the average cover in each bin is shown at the top of each class. This was also the diagram used to assist technicians in assigning cover estimates.
Remotesensing 08 00766 g005
Remotesensing 08 00766 g006 1024
Figure 6. Pictorial representation of the 20 relative percent cover metrics calculated from our LiDAR data. Each vertical bar represents a different height range for which relative percent cover was calculated. For example, the most familiar metric is the percent cover, shown in yellow in the figure above. This metric is calculated as the number of LiDAR hits over breast height (1.37–50 m) divided by the total number of LiDAR hits. These 20 metrics represent only a portion of LiDAR-derived metrics calculated, in addition to basic point statistics and percentile heights. These relative percent cover values were calculated for both first and all LiDAR returns (representing 40 of the 54 metrics calculated).
Figure 6. Pictorial representation of the 20 relative percent cover metrics calculated from our LiDAR data. Each vertical bar represents a different height range for which relative percent cover was calculated. For example, the most familiar metric is the percent cover, shown in yellow in the figure above. This metric is calculated as the number of LiDAR hits over breast height (1.37–50 m) divided by the total number of LiDAR hits. These 20 metrics represent only a portion of LiDAR-derived metrics calculated, in addition to basic point statistics and percentile heights. These relative percent cover values were calculated for both first and all LiDAR returns (representing 40 of the 54 metrics calculated).
Remotesensing 08 00766 g006
Remotesensing 08 00766 g007 1024
Figure 7. Our detailed workflow illustrates how we used plots and LiDAR to answer our key questions.
Figure 7. Our detailed workflow illustrates how we used plots and LiDAR to answer our key questions.
Remotesensing 08 00766 g007
Remotesensing 08 00766 g008 1024
Figure 8. A box-and-whisker plot that displays the distribution of ladder fuel cover within each 1 m class across all plots is shown. We also include the average ladder fuel cover between 1 and 4 m in each plot.
Figure 8. A box-and-whisker plot that displays the distribution of ladder fuel cover within each 1 m class across all plots is shown. We also include the average ladder fuel cover between 1 and 4 m in each plot.
Remotesensing 08 00766 g008
Remotesensing 08 00766 g009 1024
Figure 9. Ladder fuel cover derived from the LiDAR point cloud is compared to ground-based measurements. This model is significant at the p < 0.05 level, with an R-squared value of 0.73, root mean squared error of roughly 10%, and 10-fold cross-validation error of roughly 11%. The ground-based measurements range between 7% and 91% ladder fuel cover.
Figure 9. Ladder fuel cover derived from the LiDAR point cloud is compared to ground-based measurements. This model is significant at the p < 0.05 level, with an R-squared value of 0.73, root mean squared error of roughly 10%, and 10-fold cross-validation error of roughly 11%. The ground-based measurements range between 7% and 91% ladder fuel cover.
Remotesensing 08 00766 g009
Remotesensing 08 00766 g010 1024
Figure 10. Ladder fuel cover, as predicted from the LiDAR point cloud using Equation (1), is displayed across the study area.
Figure 10. Ladder fuel cover, as predicted from the LiDAR point cloud using Equation (1), is displayed across the study area.
Remotesensing 08 00766 g010
Table
Table 1. Average standard error between measures of percent ladder fuel cover within a given plot and height bin are shown. We report the average standard error for two-, three-, and four-photo samples for each 1 m band in each plot.
Table 1. Average standard error between measures of percent ladder fuel cover within a given plot and height bin are shown. We report the average standard error for two-, three-, and four-photo samples for each 1 m band in each plot.
2 Photos3 Photos4 Photos
Height binAverageMaximumAverageMaximumAverageMaximum
0–1 m8.2430.009.0924.738.1019.51
1–2 m8.8441.257.9523.757.7219.46
2–3 m9.3432.507.5421.796.4415.72
3–4 m9.9437.507.6722.426.7915.54
Table
Table 2. Correlation matrix (using Pearson correlation) of variables included in the model, as well as the percentage of all LiDAR returns between 1 and 4 m (COV1_4), and 4 and 8 m (COV4_8).
Table 2. Correlation matrix (using Pearson correlation) of variables included in the model, as well as the percentage of all LiDAR returns between 1 and 4 m (COV1_4), and 4 and 8 m (COV4_8).
Ladder Fuel CoverFCOV8_16COV1_4COV1_8COV4_8STD
Ladder fuel cover1-----
FCOV8_160.361----
COV1_40.690.591---
COV1_80.790.370.881--
COV4_80.740.050.540.861-
STD0.220.340.250.050.251
Table
Table 3. Average ladder fuel cover in the four focal areas, as well as the proportion of each focal area with ladder fuel cover in the top 10% (58%–100% ladder fuel cover) for the entire study area, are shown.
Table 3. Average ladder fuel cover in the four focal areas, as well as the proportion of each focal area with ladder fuel cover in the top 10% (58%–100% ladder fuel cover) for the entire study area, are shown.
Focal AreaAverage Ladder Fuel Cover (%)% Area Above 90th Percentile Ladder Fuel Cover
A31.533.89
B36.276.59
C36.6110.23
D40.9916.02

Share and Cite

MDPI and ACS Style

Kramer, H.A.; Collins, B.M.; Lake, F.K.; Jakubowski, M.K.; Stephens, S.L.; Kelly, M. Estimating Ladder Fuels: A New Approach Combining Field Photography with LiDAR. Remote Sens. 2016, 8, 766. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8090766

AMA Style

Kramer HA, Collins BM, Lake FK, Jakubowski MK, Stephens SL, Kelly M. Estimating Ladder Fuels: A New Approach Combining Field Photography with LiDAR. Remote Sensing. 2016; 8(9):766. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8090766

Chicago/Turabian Style

Kramer, Heather A., Brandon M. Collins, Frank K. Lake, Marek K. Jakubowski, Scott L. Stephens, and Maggi Kelly. 2016. "Estimating Ladder Fuels: A New Approach Combining Field Photography with LiDAR" Remote Sensing 8, no. 9: 766. https://0-doi-org.brum.beds.ac.uk/10.3390/rs8090766

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