BNL-113514-2017-JA
Long-term Post-disturbance Forest Recovery in the Greater Yellowstone Ecosystem Analyzed using Landsat Time Series Stack
Feng R. Zhao, Ran Meng, Chengquan Huang, Maosheng Zhao, Feng A. Zhao, Peng Gong, Le Yu, and Zhiliang Zhu
Accepted for publication in Remote Sensing
June 2016
Environmental & Climate Science Dept. Brookhaven National Laboratory
U.S. Department of Energy USDOE Office of Science (SC), Biological and Environmental Research (BER) (SC-23)
Notice: This manuscript has been authored by employees of Brookhaven Science Associates, LLC under Contract No. DE-SC0012704 with the U.S. Department of Energy. The publisher by accepting the manuscript for publication acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. DISCLAIMER
This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, nor any of their contractors, subcontractors, or their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or any third party’s use or the results of such use of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof or its contractors or subcontractors. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. 1 of 25
1 Article 2 Long-term Post-disturbance Forest Recovery in the 3 Greater Yellowstone Ecosystem Analyzed using Landsat 4 Time Series Stack
5 Feng R. Zhao 1, *, Ran Meng 2 , Chengquan Huang 1, Maosheng Zhao 1, Feng A. Zhao 1, Peng Gong 3, Le 6 Yu 3, and Zhiliang Zhu 4
7 1 Department of Geographical Sciences, University of Maryland, College Park, MD 20742, USA, E-Mails: 8 [email protected] (F.R.Z.) ; [email protected] (C.H.) ; [email protected] (M.Z.) ; [email protected] (F.A.Z.) 9 2 Environmental and Climate Sciences Department, Brookhaven National Laboratory, Bldg. 490A, Upton, NY 10 11973, USA, E-Mail: [email protected] (R.M.) 11 3 Ministry of Education Key Laboratory for Earth System Modeling, Center for Earth System Science, Tsinghua 12 University, Haidian, Beijing 10083, China, E-Mails : [email protected] (P.G.) ; [email protected] 13 (L.Y.) 14 4 U.S. Geological Survey, Reston, VA 20192, USA, E-Mail: [email protected] (Z.Z.) 15 * Correspondence: [email protected]; Tel.: +01-240-485-6468
16 Academic Editor: name 17 Received: date; Accepted: date; Published: date
18 Abstract: Forest recovery from past disturbance is an integral process of ecosystem carbon cycles, and 19 remote sensing provides an effective tool for tracking forest disturbance and recovery over large areas. 20 Although the disturbance products (tracking the conversion from forest to non-forest type) derived 21 using the Landsat Time Series Stack-Vegetation Change Tracker (LTSS-VCT) algorithm have been 22 validated extensively for mapping forest disturbances across the United States, the ability of this 23 approach to characterize long-term post-disturbance recovery (the conversion from non-forest to forest 24 type) has yet to be assessed. In this study, the LTSS-VCT approach was applied to examine long-term 25 (up to 24 years) post-disturbance forest spectral recovery following stand-clearing disturbances (fire and 26 harvests) in the Greater Yellowstone Ecosystems (GYE). Using high spatial resolution images from 27 Google Earth, we validated the detectable forest spectral recovery status mapped by VCT by year 2011. 28 Validation results show that the VCT was able to map long-term post-disturbance forest spectral 29 recovery with overall accuracy of ~80% for different disturbance types and forest types in the GYE. 30 Harvested areas in the GYE have higher percentages of forest spectral recovery than burned areas by 31 year 2011, and National Forests land generally has higher spectral recovery rates compared with 32 National Parks. With the availability of the VCT product nationwide, this approach can also be applied 33 to examine long-term post-disturbance forest spectral recovery in other study regions across the U.S.
34 Keywords: Wildland Fires; Timber Harvest; Detectable Forest Spectral Recovery; 1988 Yellowstone Fires
35 1. Introduction 36 Forests in the Greater Yellowstone Ecosystem (GYE) experience frequent natural (e.g., wildfires, 37 insect and disease outbreaks, snow and wind damage) and anthropogenic (land use changes and timber 38 harvesting) disturbance events [1,2]. Recovery from past disturbance is an integral process in the carbon 39 and nutrient cycles [3,4]. Inclusion of the forest recovery process following disturbance is critical to 40 calculating regional C fluxes and can better inform policy makers on both the importance and uncertainty 41 of disturbances in the regulation of the regional and global carbon cycle [5]. Many ecosystem models 2 of 25
42 assume post-disturbance forest recovery occurs immediately or homogenously across the landscape. This 43 modeling hypothesis, however, has not been supported by ground observations [6]. Therefore, there is 44 urgent need to efficiently examine forest recovery conditions at large scales. 45 Despite the urgent needs in quantifying post-disturbance forest recovery, a universal definition of 46 forest recovery is not available [7]. Terminologies such as “regeneration”, “regrowth”, and “recovery” are 47 often used, sometimes interchangeably, to represent the return of vegetation (including grass, shrub, and 48 forests) following disturbances [7-9]. In the context of forest recovery, it is necessary to clearly define 49 forest, and the recovery of forest. In this study, we defined forest following the definition of the Food and 50 Agriculture Organization (FAO) of the United Nations, i.e., greater than 10% tree cover (FAO 1998). 51 Forest recovery was then defined as the return of forest cover to more than 10% of the area. 52 This definition of forest has allowed remote sensing as a tool for examining forest disturbance and 53 recovery over large areas [10-14], especially for fires [15-18] and logging [19-22]. Many remote sensing- 54 based studies have made contributions to our understanding of vegetation spectral recovery following 55 disturbances [8,23]. In particular, a series of Landsat mission satellites have been imaging the Earth’s 56 surface since 1972, creating a time series of Landsat observations that is highly valuable for tracking land 57 change history for over four decades [24-26]. Although numerous change detection procedures have been 58 developed [27-32], with some having succeeded in mapping forest disturbances and disturbance types 59 over large areas [31,33-40], characterizing long-term forest spectral recovery following specific 60 disturbance types remains challenging. Depending on local environmental conditions and post- 61 disturbance management practices, it often takes years to decades for forests to recover its functions, 62 following a stand-clearing disturbance event [41,42]. It is valuable to determine at what point in time 63 during this process young trees have started to grow back and restore forest ecosystem services. 64 One of the algorithms, the Vegetation Change Tracker (VCT), was designed to both detect forest 65 disturbance and track post-disturbance spectral recovery [43] using annual or biannual Landsat time 66 series stacks (LTSSs) [44]. Although the disturbance products (tracking the conversion from forest to non- 67 forest type) derived using this LTSS-VCT approach have been validated extensively for mapping forest 68 disturbances across the United States [45-48], the ability of this approach to characterize post-disturbance 69 spectral recovery has yet to be assessed. Based on the FAO forest definition, we can approximate the 70 recovery of forest cover to more than 10% of the area with spectral values within pre-determined 71 thresholds. Notably, a spectral recovery is not synonymous with ecological definitions of forest recovery 72 [43]. Spectral recovery signal is also highly dependent on the choice of remote sensing index and the 73 recovery definition. In addition, recovery does not imply that the pixel is occupied by the original species. 74 In the GYE, for example, a post-fire Whitebark Pine forest might grow back as birch forest in the first few 75 decades after the fire [49]. 76 In the GYE, forest recovery trends vary from one forest type to another. This mainly has to do with 77 species-specific regrowth attributes. Lodgepole Pine (Pinus contorta) and Aspen (Populus tremuloides) 78 are fast-growing species that usually occupy lower elevations and may grow back in 5 - 20 years after a 79 stand-clearing disturbance [50,51]. Forests that grow at higher elevations, such as Whitebark Pine (Pinus 80 albicaulis), Engelmann Spruce (Picea engelmannii) and Subalpine Fir (Abies lasiocarpa), usually take 81 longer time to recover [50]. Forest recovery trends in the GYE also differ between various land ownership 82 and disturbance types. The post-harvest forest recovery in national forest land is generally faster than the 83 post-fire forest recovery in both national parks and national forests, likely due to the following reasons. 1) 84 Land productivity in national forests is generally higher than that of the national parks and wilderness 85 areas [52]. Most fires occur at high-elevation low-productivity sites with shorter growing seasons, 86 whereas harvested sites generally have longer growing seasons and higher site productivity. 2) Post- 87 harvest forest management activities helped increase forest recovery rates in the national forests. In 88 addition to other forest management techniques used to aid forest growth, planting is a common practice 89 in national forests, especially for harvested forests after natural regeneration fails [53]. 3) The post- 3 of 25
90 disturbance forest recovery in national parks and wilderness areas is subject to threats from various wild 91 animals, insects and diseases, and national forests usually have higher levels of protection against these 92 factors.
93 The GYE has been a focal point for many post-disturbance vegetation recovery studies [9,54-59], yet 94 most previous studies relied on plot-level data and were limited in their ability to densely sampling in the 95 remote, high-elevation areas of the Yellowstone Caldera. No systematic measurement based assessment 96 of forest recovery following disturbances has been conducted across the whole region in recent decades 97 (1980s to the present). The specific goals of our study, as described in this paper, were as follows: 1) to 98 validate the VCT long-term (up to 24 years) recovery product focused on tracking post-fire and post- 99 harvest forest recovery in GYE, using year 2011 (last year of the VCT product) as a reference year for the 100 recovery; 2) to use the validated recovery product to analyze the recovery rates for major fires and 101 harvests that occurred in the 1980s (to allow adequate time for recovery to occur). Results from this study 102 shed lights to regional post-disturbance forest spectral recovery trends and bring in a time-series 103 perspective in analyzing post-disturbance forest spectral recovery in the GYE. With the availability of 104 VCT product nationally, this approach can be applied to other study regions in the U.S.
105 2. Study Area
106 The 91,758-km2 study area of the GYE region includes Yellowstone National Park (YNP) and Grand 107 Teton National Park in the center and seven surrounding national forests (Bridger-Teton, Caribou- 108 Targhee, Gallatin, Shoshone, Custer, Helena, and Beaverhead-DeerLodge) (Figure 1). 4 of 25
Figure 1. Boundary, ownership, and forest fires and harvests between 1985 and 2011 (Zhao et al., 2015) in the Greater Yellowstone Ecosystem (GYE). NPs stand for National Parks and NFs represent National Forests.
109 The GYE features distinct gradients in elevation, climate and soil. Vegetation distribution highly 110 depends on topographical variations such as changes of elevation and aspects, and the effects of 111 topography are manifested through its associations with temperature and moisture availability [60]. 112 Mean annual temperature varies from 7.6 °C at lower elevations (< 1400 m) to 0.13 °C at higher elevations 113 (> 2300 m) [52]. Precipitation mostly falls as snow and generally increases with elevation; mean annual 114 precipitation ranges from 1368 mm to 2414 mm [60]. The growing season in the GYE varies from less than 115 three months at higher elevations to around six months at lower elevations [1]. A large portion of the 116 national parks (NPs), including the Yellowstone Plateau and surrounding mountain ranges, lies at 117 relatively high elevations. The national forest (NF) lands are mostly at moderate and low elevations on 118 the flanks of the plateau. The soils at higher elevations are largely composed of nutrient poor rhyolites 119 and andesites with low water-holding capacities [60]. The valley bottoms and floodplains contain glacial 120 outwash and alluvium soils that generally feature higher nutrients and water-holding capacities in 121 relative terms [52].
122 Natural forest vegetation in the study area is a mosaic of major coniferous species [61]. Lodgepole 123 Pine is widespread in YNP and dominates about 70% of the forested NP area (approximately 5295 out of 5 of 25
124 7355 km2), followed by forest species such as Whitebark Pine, Subalpine Fir, Engelmann Spruce, Douglas- 125 fir (Pseudotsuga menziesii) and Aspen [60]. Whitebark Pine occupies approximately 15% of the YNP 126 forested area, especially at the higher elevations. Engelmann Spruce and Subalpine Fir often co-exist 127 below the elevation zone of Whitebark Pine, with Douglas-fir dominating the lowest elevations. 128 Lodgepole Pine forests between 2000 m and 2600 m are supported by rhyolite soils, and Douglas-fir can 129 be found up to 2300 m on andesitic soils and in relative warm conditions [1]. Above these elevations in all 130 soil types, forests types such as Subalpine fir, Engelmann spruce, and Whitebark pine dominate. About 131 70% of the GYE forested area is located in national forests (Figure 1), which contain wilderness areas 132 designated by the Wilderness Act of 1964 (~22% of the GYE forested area) and areas managed for timber 133 production (~47.7% of GYE forested area). National parks occupy more than a quarter of the GYE 134 forested area and the remaining forested areas are under other ownership, such as state or private forests.
135 The recent history and composition of disturbance events during the study interval (1984-2011) are 136 shown in Figure 1, which also shows the effects of ownership pattern on disturbances. In the GYE 137 national parks and wilderness area, fire was the most dominant disturbance agent, affecting over 37% of 138 the forested area in the GYE national parks. Active harvest event was a major human-induced 139 disturbance type occurring in the national forests. In particular, the harvested area in the Caribou- 140 Targhee National Forest (west to the National Park) was four times larger than the burned areas during 141 the study period [2].
142 3. Materials and Methods
143 3.1 LTSS assembling 144 LTSSs from both Landsat TM and ETM sensors were assembled for the 8 Landsat World Reference 145 System-2 (WRS-2) path/row locations required to cover the GYE (Figure 1, Landsat scenes for the study 146 region include p37r29, p37r30, p38r28, p38r29, p38r30, p38r31, p39r28, and p39r29). The Landsat images 147 were downloaded from the USGS Global Visualization Viewer (GloVis) at the 30-m resolution. The 148 original images were first converted to surface reflectance using the Landsat Ecosystem Disturbance 149 Adaptive Processing System (LEDAPS) algorithm [62]. Geometrically, no additional correction was 150 performed on these images because they had already been ortho-rectified by the USGS to achieve 151 subpixel geolocation accuracy [44,63]. A detailed description of the procedures involved in assembling 152 LTSSs has been provided in a previous study [44].
153 Each LTSS contained one image per year for the years between 1984 and 2011 that had at least one 154 clear view (cloud-free or nearly cloud-free with cloud cover less than 5%) image acquired during the leaf- 155 on growing season (May to September). If no such image was available in certain years, multiple partly 156 cloudy images acquired during the leaf-on growing season of that year were used to produce a composite 157 image [43,48]. The compositing algorithm identifies and replaces the cloud and shadow contaminated 158 pixels and adjusts the phenological differences among the affected images [44]. The cloud and cloud 159 shadow were identified using an automated masking algorithm developed by Huang et al. (2010). 160 Omission errors for the masking algorithm were around 1% for the cloud class, although the errors were 161 higher for low cloud cover and semiarid environment, potentially leading to higher forest and 162 disturbance mapping errors [64]. If no more than 1 clear-view observation was available in a year at a 163 given pixel location, the pixel with the maximum NDVI value was selected [48,65]. If more than one clear- 6 of 25
164 view observation was available, the clear-view observation that had the highest brightness temperature 165 was selected [48,66]. Here, clear-view observations referred to those that were not contaminated by 166 clouds or shadows and did not have other data quality problems [43]. 167 3.2 Forest disturbance and recovery mapping
168 The LTSSs assembled in Section 3.1 were analyzed using the VCT algorithm to map forest 169 disturbance and recovery. The VCT uses an integrated forest z-score (IFZ) index to track forest changes at 170 each pixel location: