Creating Landscape-Scale Site Index Maps for the Southeastern US Is Possible with Airborne LiDAR and Landsat Imagery
Abstract
:1. Introduction
- How effective is ALS data collected via disparate projects coupled with Landsat-based forest age products in predicting site index for pine plantations?
- What is the distribution and regional spatial pattern of site productivity?
2. Methods
2.1. Study Area
2.2. Remote Sensing Datasets
- ALS data: We considered many ALS data acquisition projects that: (1) were available in the public domain as of 2011; (2) were from the southeastern states of the USA; (3) had recorded multiple returns (i.e., two or more) for each pulse; and (4) contained at least one FIA plot with pine plantation, in their area of coverage (see below). This resulted in 38 projects being finally selected, spread across eight southeastern states (Virginia, North Carolina, South Carolina, Georgia, Florida, Alabama, Mississippi and Texas), as shown in Figure 1. We used georeferenced point-cloud data from these projects. Their acquisition dates ranged from 2005 to 2012, and 31 of these projects used instruments capable of recording up to four returns.
- Landsat-based prediction of plantation age: The vegetation change tracker (VCT) algorithm is a fairly well-established and extensively tested technique that effectively leverages multitemporal Landsat data (both Landsat 5 TM and Landsat 7 ETM+; 1984–2011) to detect disturbances to forest stands (such as clearcutting) [33]. Recently, additional steps that partitioned such disturbances into stand-clearing ones (i.e., clearcuts) and non-stand-clearing ones (such as thinning, insect damage and ice damage) were proposed as part of an enhanced vegetation change tracker algorithm (henceforth, eVCT) [30]. Identifying and demarcating the stand-clearing disturbance helps delineate the stand, and the number of years since disturbance is a good proxy for stand age. All pixels in a particular delineated stand are given a single value of year of disturbance (and hence, age).
- Landsat-based identification of planted pine stands: We used a map of planted pine forest stands in the southeastern USA generated by a combination of spectral and multitemporal Landsat data [31]. The following were the prominent datasets used in this classification: (1) National Land Cover Database 2011, a Landsat-based land cover map; (2) summer and winter Normalized Difference Vegetation Index (NDVI) (and their difference) from Landsat; (3) the 30 m resolution global forest change dataset by Hansen et al. [34] (2000–2013); and (4) the Vegetation Change Tracker dataset (VCT; 1985–2011 [33]). A decision tree model was developed using these datasets as well as a combination of plot-level forest type data and Google Earth imagery for training and verification. From the generated map, it was estimated that about 28% of the total forest area of southeastern USA is covered by pine plantations.
2.3. Canopy Height Model
- The plots had to be single-condition. “Condition” is a term used by the FIA to refer to major land-use types on the plots [35]. A consequence of FIA’s plot center randomizing strategy is that some plots can straddle multiple land-use types (e.g., a plot may be 70% pine stand, 20% hardwoods stand, and 10% water). Hence, we filtered out such non-homogeneous plots.
- The plot was marked to have clear evidence of artificial regeneration (i.e., it is a plantation).
- The pine species basal area of the plot (percent) was ≥95% of total basal area on the plot.
- ALS data within two years from the year of visit and measurement of the FIA plots are available. That is, there is a maximum of ±2 years difference between the years of ALS acquisition and of FIA field measurement.
- Ground classification, using the lasground tool of lastools (http://rapidlasso.com/lastools/).
- The understory points were identified, using a threshold of 3.0 m. This threshold was determined after inspecting 12 FIA plots (of the 211) for which understory height information was available from the FIA.
- ALS distributional metrics: The standard percentiles of height above ground (, etc.) were extracted from the point cloud (using MATLAB v7.14, 2012), after excluding understory points.
- The dominant height of all trees on our field plot (subplot 1 of the FIA) was calculated as the average height of the five tallest trees on that subplot. We used the FIA’s HT variable, in the TREE table.
2.4. Generation of a Sample Site Index Map from ALS and Landsat Data
2.5. Accuracy Assessment of Site Index Predictions
- Dominant height: This would be a normal distribution of mean 0.0 and standard deviation () equal to the root mean square error of the OLS-based linear model.
- Age, predicted by eVCT: For this, we used the following independent ground-truth sources:
- (a)
- Appomattox-Buckingham State Forest plot set: This is a set of 23 loblolly pine forest plots located in stands of known age in central Virginia, USA [41]. We have good estimates of the age of these stands, based on quality-checked planting records. However, these plots are relatively few in number and their coverage is geographically limited.
- (b)
- FIA plot set: This is a set of 75 FIA plantation pine plots for which ages (year of planting) were available. The accuracy of these is variable, as the availability of planting records and reliable tree core data is limited. The lack of agreement between the FIA field estimates of age and eVCT is assumed to be due to erroneous eVCT ages. Note that this may not always be true, especially given the uncertainty in FIA age estimates.
The error distribution in the eVCT-based age was then estimated by randomly sampling among these plots, and getting the difference between the ages.
2.6. HU-Based Spatial Aggregation
3. Results
3.1. Canopy Height Model
3.2. Accuracy Assessment of Site Index Map
- The Appomattox-Buckingham State Forest plot set: The bias of calculated site index was estimated to be −1.0 m, the RMSE estimate was 3.0 m, and the relative RMSE was 15.1%.
- FIA plot set: In this case, the estimated bias was −0.28 m, the RMSE was 3.8 m, and the relative RMSE was 19.7%.
4. Discussion
4.1. Efficacy of Site Index Maps Based on Disparate ALS Projects Data and Landsat
4.2. Distribution and Regional Patterns of Productivity
4.3. Contrasting Current and Historical Site Index Levels
4.4. Forest Productivity Maps: Other Utilization Avenues
5. Conclusions
Author Contributions
Funding
Conflicts of Interest
References
- Coops, N.C. Characterizing forest growth and productivity using remotely sensed data. Curr. For. Rep. 2015, 1, 195–205. [Google Scholar] [CrossRef]
- Splechtn, B.E. Height growth and site index models for Pacific silver fir in southwestern British Columbia. J. Ecosyst. Manag. 2001, 1, 1–14. [Google Scholar]
- Dyck, B. Precision forestry—The path to increased profitability. In Proceedings of the 2nd International Precision Forestry Symposium, Seattle, WA, USA, 15–17 June 2003; pp. 3–8. [Google Scholar]
- Turner, D.P.; Koerper, G.J.; Harmon, M.E.; Lee, J.J. A carbon budget for forests of the conterminous United States. Ecol. Appl. 1995, 5, 421–436. [Google Scholar] [CrossRef]
- Han, F.X.; Plodinec, M.J.; Su, Y.; Monts, D.L.; Li, Z. Terrestrial carbon pools in southeast and south-central United States. Clim. Chang. 2007, 84, 191–202. [Google Scholar] [CrossRef]
- USFS. Future of America’s Forest and Rangelands: Forest Service 2010 Resources Planning Act Assessment; General Technical Report GTR-WO-87; USDA Forest Service: Washington, DC, USA, 2012.
- Skovsgaard, J.P.; Vanclay, J.K. Forest site productivity: A review of the evolution of dendrometric concepts for even-aged stands. Forestry 2008, 81, 13–31. [Google Scholar] [CrossRef]
- Sellers, P.J.; Shuttleworth, W.J.; Dorman, J.L.; Dalcher, A.; Roberts, J.M. Calibrating the simple biosphere model for Amazonian tropical forest using field and remote sensing data. Part I: Average calibration with field data. J. Appl. Meteorol. 1989, 28, 727–759. [Google Scholar] [CrossRef]
- Launay, M.; Guerif, M. Assimilating remote sensing data into a crop model to improve predictive performance for spatial applications. Agric. Ecosyst. Environ. 2005, 111, 321–339. [Google Scholar] [CrossRef]
- Avery, T.E.; Burkhart, H.E. Forest Measurements; Waveland Press: Long Grove, IL, USA, 2015. [Google Scholar]
- Zhao, D.; Kane, M.; Teskey, R.; Fox, T.R.; Albaugh, T.J.; Allen, H.L.; Rubilar, R. Maximum response of loblolly pine plantations to silvicultural management in the southern United States. For. Ecol. Manag. 2016, 375, 105–111. [Google Scholar] [CrossRef] [Green Version]
- Wear, D.; Greis, J. The Southern Forest Futures Project: Summary Report; Technical Report SRS-GTR-168; USDA-Forest Service, Southern Research Station: Asheville, NC, USA, 2012.
- Sabatia, C.O.; Burkhart, H.E. Predicting site index of plantation loblolly pine from biophysical variables. For. Ecol. Manag. 2014, 326, 142–156. [Google Scholar] [CrossRef]
- Bryars, C.; Maier, C.; Zhao, D.; Kane, M.; Borders, B.; Will, R.; Teskey, R. Fixed physiological parameters in the 3-PG model produced accurate estimates of loblolly pine growth on sites in different geographic regions. For. Ecol. Manag. 2013, 289, 501–514. [Google Scholar] [CrossRef]
- Carlson, C.A.; Fox, T.R.; Creighton, J.; Dougherty, P.M.; Johnson, J.R. Nine-year growth responses to planting density manipulation and repeated early fertilization in a loblolly pine stand in the Virginia Piedmont. South. J. Appl. For. 2009, 33, 109–114. [Google Scholar]
- Thomas, R.Q.; Brooks, E.B.; Jersild, A.L.; Ward, E.J.; Wynne, R.H.; Albaugh, T.J.; Dinon-Aldridge, H.; Burkhart, H.E.; Domec, J.C.; Fox, T.R. Leveraging 35 years of Pinus taeda research in the southeastern US to constrain forest carbon cycle predictions: Regional data assimilation using ecosystem experiments. Biogeosciences 2017, 14, 3525. [Google Scholar] [CrossRef]
- Jokela, E.J.; Dougherty, P.M.; Martin, T.A. Production dynamics of intensively managed loblolly pine stands in the southern United States: A synthesis of seven long-term experiments. For. Ecol. Manag. 2004, 192, 117–130. [Google Scholar] [CrossRef]
- Miller, J.H.; Zutter, B.R.; Zedaker, S.M.; Edwards, M.B.; Newbold, R.A. Growth and yield relative to competition for loblolly pine plantations to midrotation: A southeastern United States regional study. South. J. Appl. For. 2003, 27, 237–252. [Google Scholar]
- Gatziolis, D. LiDAR-derived site index in the US Pacihic Northwest: Challenges and opportunities. In Proceedings of the ISPRS Workshop on Laser Scanning 2007 and SilviLaser 2007, Espoo, Finland, 12–14 September 2007. [Google Scholar]
- Wulder, M.A.; White, J.C.; Stinson, G.; Hilker, T.; Kurz, W.A.; Coops, N.C.; St-Onge, B.; Trofymow, J.T. Implications of differing input data sources and approaches upon forest carbon stock estimation. Environ. Monit. Assess. 2010, 166, 543–561. [Google Scholar] [CrossRef] [PubMed]
- Tompalski, P.; Coops, N.C.; White, J.C.; Wulder, M.A.; Pickell, P.D. Estimating forest site productivity using airborne laser scanning data and Landsat time series. Can. J. Remote Sens. 2015, 41, 232–245. [Google Scholar] [CrossRef]
- Kandare, K.; Ørka, H.O.; Dalponte, M.; Næsset, E.; Gobakken, T. Individual tree crown approach for predicting site index in boreal forests using airborne laser scanning and hyperspectral data. Int. J. Appl. Earth Obs. Geoinf. 2017, 60, 72–82. [Google Scholar] [CrossRef]
- Packalén, P.; Mehtätalo, L.; Maltamo, M. ALS-based estimation of plot volume and site index in a eucalyptus plantation with a nonlinear mixed-effect model that accounts for the clone effect. Ann. For. Sci. 2011, 68, 1085. [Google Scholar] [CrossRef]
- Watt, M.S.; Dash, J.P.; Bhandari, S.; Watt, P. Comparing parametric and non-parametric methods of predicting Site Index for radiata pine using combinations of data derived from environmental surfaces, satellite imagery and airborne laser scanning. For. Ecol. Manag. 2015, 357, 1–9. [Google Scholar] [CrossRef]
- Socha, J.; Pierzchalski, M.; Bałazy, R.; Ciesielski, M. Modelling top height growth and site index using repeated laser scanning data. For. Ecol. Manag. 2017, 406, 307–317. [Google Scholar] [CrossRef]
- Ham, S.A.; Mikhailova, E.A.; Gering, L.R.; Post, C.J.; Bridges, W.C.; Cox, S.K. Temporal analysis of field, SSURGO, and LiDAR-derived site indices in the southeastern United States. Soil Sci. 2013, 178, 325–334. [Google Scholar] [CrossRef]
- Gopalakrishnan, R.; Thomas, V.A.; Coulston, J.W.; Wynne, R.H. Prediction of canopy heights over a large region using heterogeneous lidar datasets: Efficacy and challenges. Remote Sens. 2015, 7, 11036–11060. [Google Scholar] [CrossRef]
- Sexton, J.O.; Bax, T.; Siqueira, P.; Swenson, J.J.; Hensley, S. A comparison of lidar, radar, and field measurements of canopy height in pine and hardwood forests of southeastern North America. For. Ecol. Manag. 2009, 257, 1136–1147. [Google Scholar] [CrossRef]
- Popescu, S.C.; Wynne, R.H. Seeing the trees in the forest: Using lidar and multispectral data fusion with local filtering and variable window size for estimating tree height. Photogramm. Eng. Remote Sens. 2004, 70, 589–604. [Google Scholar] [CrossRef]
- Kauffman, J.S.; Prisley, S.P. Automated estimation of forest stand age using Vegetation Change Tracker and machine learning. Math. Comput. For. Nat. Resour. Sci. 2016, 8, 4. [Google Scholar]
- Fagan, M.E.; Morton, D.C.; Cook, B.D.; Masek, J.; Zhao, F.; Nelson, R.F.; Huang, C. Mapping pine plantations in the southeastern US using structural, spectral, and temporal remote sensing data. Remote Sens. Environ. 2018, 216, 415–426. [Google Scholar] [CrossRef]
- Sugarbaker, L.J.; Constance, E.W.; Heidemann, H.K.; Jason, A.L.; Lukas Vicki, S.D.; Stoker, J.M. The 3D Elevation Program Initiative: A Call for Action; US Geological Survey: Reston, VA, USA, 2014.
- Huang, C.; Goward, S.N.; Masek, J.G.; Thomas, N.; Zhu, Z.; Vogelmann, J.E. An automated approach for reconstructing recent forest disturbance history using dense Landsat time series stacks. Remote Sens. Environ. 2010, 114, 183–198. [Google Scholar] [CrossRef]
- Hansen, M.C.; Potapov, P.V.; Moore, R.; Hancher, M.; Turubanova, S.A.; Tyukavina, A.; Thau, D.; Stehman, S.V.; Goetz, S.J.; Loveland, T.R. High-resolution global maps of 21st-century forest cover change. Science 2013, 342, 850–853. [Google Scholar] [CrossRef] [PubMed]
- Bechtold, W.A.; Patterson, P.L. The Enhanced Forest Inventory and Analysis Program: National Sampling Design and Estimation Procedures; US Department of Agriculture Forest Service, Southern Research Station: Asheville, NC, USA, 2005.
- McRoberts, R.E. The effects of rectification and Global Positioning System errors on satellite image-based estimates of forest area. Remote Sens. Environ. 2010, 114, 1710–1717. [Google Scholar] [CrossRef]
- White, J.C.; Wulder, M.A.; Vastaranta, M.; Coops, N.C.; Pitt, D.; Woods, M. The utility of image-based point clouds for forest inventory: A comparison with airborne laser scanning. Forests 2013, 4, 518–536. [Google Scholar] [CrossRef]
- Schultz, R.P. Loblolly Pine: The Ecology and Culture of Loblolly Pine (Pinus taeda L.); Department of Agriculture, Forest Service: Washington, DC, USA, 1997.
- Diéguez-Aranda, U.; Burkhart, H.E.; Amateis, R.L. Dynamic site model for loblolly pine (Pinus taeda L.) plantations in the United States. For. Sci. 2006, 52, 262–272. [Google Scholar]
- Lobell, D.B.; Asner, G.P.; Ortiz-Monasterio, J.I.; Benning, T.L. Remote sensing of regional crop production in the Yaqui Valley, Mexico: Estimates and uncertainties. Agric. Ecosyst. Environ. 2003, 94, 205–220. [Google Scholar] [CrossRef]
- Quirino, V.F. Evaluating the Potential for Estimating Age of Even-Aged Loblolly Pine Stands Using Active and Passive Remote Sensing Data. Ph.D. Dissertation, Virginia Tech., Blacksburg, VA, USA, 2014. [Google Scholar]
- Seaber, P.R.; Kapinos, F.P.; Knapp, G.L. Hydrologic Unit Maps; US Government Printing Office: Denver, CO, USA, 1987.
- Burkhart, H.E.; Brooks, E.B.; Dinon-Aldridge, H.; Sabatia, C.O.; Gyawali, N.; Wynne, R.H.; Thomas, V.A. Regional Simulations of Loblolly Pine Productivity with CO2 Enrichment and Changing Climate Scenarios. For. Sci. 2018, 64, 349–357. [Google Scholar] [CrossRef]
- Omernik, J.M. Ecoregions of the conterminous United States. Ann. Assoc. Am. Geogr. 1987, 77, 118–125. [Google Scholar] [CrossRef]
- Belote, R.T.; Prisley, S.; Jones, R.H.; Fitzpatrick, M.; de Beurs, K. Forest productivity and tree diversity relationships depend on ecological context within mid-Atlantic and Appalachian forests (USA). For. Ecol. Manag. 2011, 261, 1315–1324. [Google Scholar] [CrossRef]
- Noormets, A.; Gavazzi, M.J.; McNulty, S.G.; Domec, J.C.; Sun, G.E.; King, J.S.; Chen, J. Response of carbon fluxes to drought in a coastal plain loblolly pine forest. Glob. Chang. Biol. 2010, 16, 272–287. [Google Scholar] [CrossRef]
- NRCS. Land Resource Regions and Major Land Resource Areas of the United States, the Caribbean, and the Pacific Basin; US Department of Agriculture Handbook: Washington, DC, USA, 2006.
- Carmean, W.H.; Hahn, J.T.; Jacobs, R.D. Site Index Curves for Forest Tree Species in the Eastern United States; General Technical Report NC-128; US Department of Agriculture, Forest Service, North Central Forest Experiment Station: St. Paul, MN, USA, 1989.
- Noordermeer, L.; Bollandsås, O.M.; Gobakken, T.; Næsset, E. Direct and indirect site index determination for Norway spruce and Scots pine using bitemporal airborne laser scanner data. For. Ecol. Manag. 2018, 428, 104–114. [Google Scholar] [CrossRef]
- Næsset, E. Predicting forest stand characteristics with airborne scanning laser using a practical two-stage procedure and field data. Remote Sens. Environ. 2002, 80, 88–99. [Google Scholar] [CrossRef]
- Riaño, 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]
- Hopkinson, C.; Chasmer, L.; Lim, K.; Treitz, P.; Creed, I. Towards a universal lidar canopy height indicator. Can. J. Remote Sens. 2006, 32, 139–152. [Google Scholar] [CrossRef]
- Fox, T.R.; Jokela, E.J.; Allen, H.L. The development of pine plantation silviculture in the southern United States. J. For. 2007, 105, 337–347. [Google Scholar]
- Fox, T.R.; Kyle, K.H.; Andrews, L.J.; Aust, W.M.; Burger, J.A.; Hansen, G.H. Long-term effects of drainage, bedding, and fertilization on growth of loblolly pine (Pinus taeda L.) in the Coastal Plain of Virginia. South. J. Appl. For. 2005, 29, 205–214. [Google Scholar]
- Albaugh, T.J.; Fox, T.R.; Allen, H.L.; Rubilar, R.A. Juvenile southern pine response to fertilization is influenced by soil drainage and texture. Forests 2015, 6, 2799–2819. [Google Scholar] [CrossRef]
- NRCS. National Soil Survey Handbook (Title 430-VI); U.S. Department of Agriculture: Washington, DC, USA, 2007.
- Norby, R.J.; Zak, D.R. Ecological lessons from free-air CO2 enrichment (FACE) experiments. Annu. Rev. Ecol. Evol. Syst. 2011, 42, 181. [Google Scholar] [CrossRef]
- Murthy, R.; Dougherty, P.M.; Zarnoch, S.J.; Allen, H.L. Effects of carbon dioxide, fertilization, and irrigation on photosynthetic capacity of loblolly pine trees. Tree Physiol. 1996, 16, 537–546. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Thomas, R.Q.; Canham, C.D.; Weathers, K.C.; Goodale, C.L. Increased tree carbon storage in response to nitrogen deposition in the US. Nat. Geosci. 2010, 3, 13. [Google Scholar] [CrossRef]
- Subedi, S.; Fox, T.R.; Wynne, R.H. Determination of fertility rating (FR) in the 3-PG model for loblolly pine plantations in the southeastern United States based on site index. Forests 2015, 6, 3002–3027. [Google Scholar] [CrossRef]
© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
Gopalakrishnan, R.; Kauffman, J.S.; Fagan, M.E.; Coulston, J.W.; Thomas, V.A.; Wynne, R.H.; Fox, T.R.; Quirino, V.F. Creating Landscape-Scale Site Index Maps for the Southeastern US Is Possible with Airborne LiDAR and Landsat Imagery. Forests 2019, 10, 234. https://doi.org/10.3390/f10030234
Gopalakrishnan R, Kauffman JS, Fagan ME, Coulston JW, Thomas VA, Wynne RH, Fox TR, Quirino VF. Creating Landscape-Scale Site Index Maps for the Southeastern US Is Possible with Airborne LiDAR and Landsat Imagery. Forests. 2019; 10(3):234. https://doi.org/10.3390/f10030234
Chicago/Turabian StyleGopalakrishnan, Ranjith, Jobriath S. Kauffman, Matthew E. Fagan, John W. Coulston, Valerie A. Thomas, Randolph H. Wynne, Thomas R. Fox, and Valquiria F. Quirino. 2019. "Creating Landscape-Scale Site Index Maps for the Southeastern US Is Possible with Airborne LiDAR and Landsat Imagery" Forests 10, no. 3: 234. https://doi.org/10.3390/f10030234
APA StyleGopalakrishnan, R., Kauffman, J. S., Fagan, M. E., Coulston, J. W., Thomas, V. A., Wynne, R. H., Fox, T. R., & Quirino, V. F. (2019). Creating Landscape-Scale Site Index Maps for the Southeastern US Is Possible with Airborne LiDAR and Landsat Imagery. Forests, 10(3), 234. https://doi.org/10.3390/f10030234