Assessment of the CLASlite forest monitoring system in detecting disturbance from selective logging in the

Detecting and monitoring forest disturbance from selective logging is necessary to develop effective strategies and polices that conserve tropical forests and mitigate climate change. We assessed the potential of using the remote sensing tool, CLASlite forest monitoring system, to detect disturbance from timber harvesting in four community forests (ejidos) of the Selva Maya on the Yucatan Peninsula, Mexico. Selective logging impacts (e.g. felling gaps, skid trails, logging roads and log landings) were mapped using GPS in the 2014 annual cutting areas (ACAs) of each ejido. We processed and analyzed two pre-harvest Landsat images (2001 and 2013) and one post-harvest image (November 2014) with the CLASlite system, producing maps of degraded, deforested and unlogged areas in each ACA. Based on reference points of disturbed (felling and skidding), deforested (log landings and roads) and unlogged areas in each ACA, we applied accuracy assessments which showed very low overall accuracies (<19.1%). Selective logging impacts, mainly from log landings and new logging road construction, were detected in only one ejido which had the highest logging intensity (7 m3 ha–1).


Introduction
Tropical forests are a major global carbon stock and play a key role in climate change mitigation (Putz et al. 2008).However, their extension is being reduced by deforestation, and degradation processes are also impacting these ecosystems, resulting in greater carbon emissions (69%) than from forest cover loss (Baccini et al. 2017).Forest degradation is defined as the reduction in the capacity of the forest to provide products and environmental services due to changes in biomass, structure and composition without implying a change in land cover; it is caused by anthropogenic activities such as selective logging, firewood collection, fire and grazing (Griscom et al. 2009;Sasaki and Putz 2009;Simula and Mansur 2011).In Latin America, selective logging implies the removal of high-value timber from the forest usually present in low densities (1 to 20 trees per hectare: Putz et al. 2008) and is identified as the main cause of forest degradation (Pearson et al. 2017).Thus, detecting and monitoring potential degradation from timber harvesting is vital to developing effective strategies that conserve tropical forests and mitigate climate change (Asner 2009;Asner et al. 2010;Asner 2014;Reimer et al. 2015).
Unfortunately, forest degradation is not as easy to detect and monitor in comparison to deforestation; it is costlier and more complex using remote sensing techniques, and results can often be subjective due to succession dynamics and assumed baselines for non-degraded forest (Putz and Redford 2010;Thompson et al. 2013).Efforts have been pursued to detect and measure anthropogenic and natural forest disturbances applying remote sensing of satellite imagery in order to develop viable methods of monitoring the phenomenon of degradation in the tropics (Asner 2009;Griscom et al. 2009;Thompson et al. 2013).Landsat images have been successfully applied to detect forest disturbance from selective logging in the tropics (Broadbent et al. 2006;Zu et al. 2012;Anwar and Stein 2012;Tritsch et al. 2016), providing a potentially valuable image source for tropical forest managers and conservation practitioners.
The objective of this study is to evaluate the Carnegie Landsat Analysis System-lite (CLA-Slite), a software package developed to identify deforestation and forest degradation from satellite imagery, especially Landsat.CLASlite is a semi-automated system that integrates remote sensing research and algorithms into a user-friendly environment, providing a tool that natural resource managers, governments institutions and conservation organizations can use for forest monitoring and environmental planning (http://claslite.ciw.edu).Although CLASlite has been applied to detect deforestation and forest degradation from logging and other anthropogenic activities in the tropics (Asner et al. 2010;Allnutt et al. 2013;Bryan et al. 2013;Chicas et al. 2016), its effectiveness in detecting disturbances and monitoring potential degradation from selective logging in Mexico has not been evaluated.We tested the accuracy of the automated and user-friendly CLASlite system to detect and monitor selective logging disturbances in the Selva Maya, Mexico as a potential tool for local forest managers and conservation practitioners with limited or no remote sensing experience.

Study area and mapping
We assessed the capability of CLASlite in detecting forest disturbance from selective logging in four community forestry ejidos (Caobas, Felipe Carrillo Puerto, Noh-Bec and Petcacab) located in the state of Quintana Roo, southeast Mexico (Fig. 1).The ejidos are situated in flat karst topography and at low elevations (10 to 150 m a.s.l.).Climate in the region is hot and sub-humid with an annual average precipitation of 1200 mm, with a marked dry season (<60 mm a month) from November to April (Gutiérrez-Granados et al. 2011).The predominant upland forests in the region are between 15 and 25 m in height and are semi-deciduous, loosing around 25 to 50% of its foliage during the dry season (Flores and Espejel 2004;Ellis and Porter-Bolland 2008).
The study ejidos are all involved with selective logging for commercial purposes and apply a polycyclic silvicultural system with a 25-year cutting cycle guided by a required forest management plan (FMP) (Ellis et al. 2015(Ellis et al. , 2017)).Forest management is conducted within annual cutting areas (ACAs) planned out in the FMP.Each ACA can range from 200 to 1000 ha, depending on the size and scale of the forest management operation.Harvestable trees are previously marked by logging crews and then felled by chainsaw and extracted (skidded) using an articulated forestry tractor (skidder), or in the case of Caobas, a modified agricultural tractor.Skidded logs are deposited in log-landings, which are accessed by logging roads to transport timber to neighboring sawmills or markets.
Mapping was conducted in the 2014 ACAs from March to November 2015 after logging activities were concluded.Delimited ACAs and sub-blocks were georeferenced from maps in FMPs or from shapefiles provided by forest technicians.Two 100 ha sub-blocks were randomly selected for field sampling and were mapped during fieldwork using Garmin GPS Map 60csx.In each sub-block, all felled trees (stumps) or felling gaps, skid trails, log landings, and logging roads were georeferenced.

CLASlite analysis
Landsat scenes (path 19, row 47) were selected and downloaded in GEOTIFF format from the United States Geological Survey (USGS) Earth Explorer website (https://earthexplorer.usgs.gov).We used Landsat images due to their free and easy access, keeping in mind that CLASlite could be applied by local forest managers and conservation NGOs with limited budgets and remote sensing capacities.We selected two pre-harvest scenes, a Landsat 7 ETM with minimal cloud cover (<10%) from March 2001 (dry season) and a Landsat 8 OLI with under 20% cloud cover from July 2013 (wet season), and one post-harvest scene, a Landsat 8 OLI from November 2014 (dry season).Very few Landsat 8 OLI scenes with less than 20% cloud cover were available, indicating the difficulty in obtaining quality images of the study area.Nevertheless, CLASlite has been reported as being effective in removing cloud cover and other atmospheric artefacts common in the humid tropics (Reimer et al. 2015).
Forest change detection of deforested and degraded areas was conducted between the pre and post-harvest input imagery using CLASlite 3.3.The CLASlite system automatically performs radiometric and atmospheric correction and produces fractional cover maps of photosynthetic vegetation (PV), non-photosynthetic vegetation (NPV) and bare soil (S) from an available library of spectral signatures and applying Automated Monte Carlo Unmixing algorithm (AutoMCU), also known as Spectral Mixture Analysis (Asner et al. 2005;Asner 2009; Romero-Sanchez and Ponce-Hernandez 2017).CLASlite combines fractional cover maps of both dates to produce results of deforested, degraded and areas with unchanged forest cover.For detection of forest loss (deforestation) and degradation (areas of persistent forest disturbance), CLASlite applies expressions derived from extensive field testing and validation (Asner et al. 2005;Asner et al. 2009).For example, in their default settings, forested pixels have a PV ≥ 80%, deforested pixels are those where PV decreases to 60% or less, S increases 15% or more and NPV increases to 20% or more, while degraded pixels represent those where PV decreased from 80%, but remain greater than 60%, S increases 15% or more and NPV increases by 20% or more.

Accuracy assessment
Based on validation points obtained during field mapping of ACAs, we conducted accuracy assessments of the classification outputs produced by CLASlite 3.3: 1) undisturbed forest cover, 2) deforested and 3) disturbed or degraded forest areas.The validation points extracted and used for accuracy assessments included: 1) 100 felling gap points randomly selected from all georeferenced gaps (stumps) using a random number generator; 2) 25 skid trail points located in randomly selected trails and at random distances from each skid trail entrance, 3) 25 logging road points taken every 500 m on the road during the ACA mapping process, 4) the center point of all log landings mapped in each ACAs (2 to 9) and 5) 25 randomly selected un-logged points taken during field mapping at a random direction and distance (between 50 to 80 m) from each skid trail validation point.We then combined our reference points of felling and skidding into a single "disturbance" class and log landings and logging roads into a single "deforestation" class based on their type of impact on the forest biomass (i.e.damage to forest biomass and small canopy gaps or complete removal of trees and large canopy openings).
For accuracy assessment we applied the method described by Olofsson et al. (2013Olofsson et al. ( , 2014) using the accurAssess tool for QGIS described and available from Mas et al. (2014).The adjusted confusion matrices and producer accuracies were used to determine the accuracy of CLASlite in classifying disturbance, deforestation and unlogged forest in the four annual cutting areas sampled.Confusion matrices generated by the Olofsson method are adjusted to weigh the number of reference points and proportions of each category on the map (Olofsson et al. 2013(Olofsson et al. , 2014;;Mas et al. 2014).

Results
Fig. 2 shows the post-harvest image (2014), pre (2001) and post-harvest (2014) fractional cover maps and CLASlite deforestation and disturbance results for the four ejido ACAs, also indicating validation points of disturbance and deforestation impacts used for accuracy assessments.The ejidos we sampled contain among the largest ACAs (>800 ha) of the community forestry ejidos in the region, and Noh Bec and Petcacab are among the most intensively harvested; yet, there were no indications that even in these larger forest management operations CLASlite could accurately detect selective logging impacts on forest biomass.Even the most discernable harvest impacts, logging roads and log landings, were not sufficiently detected by CLASlite.
Overall accuracies of CLASlite results were very low: 13.9% in Petcacab and Caobas, 14.5% in Felipe Carrillo Puerto, and 19.1% in Noh Bec, the ejido with the most visible impacts of selective logging in the post-harvest image.Table 1 shows the adjusted confusion matrices for the 2001-2014 forest change maps in each ACA (i.e.unchanged forest cover or un-logged, deforested or disturbed forest).As expected, unlogged areas were the most accurately predicted (14 to 16%) with producer accuracies over 99%.However, over 68% of the observed disturbed areas were mis-classified as unlogged pixels in all ejidos.Producer accuracies for classifying deforested areas from building logging roads and log landings and disturbed areas from felling and skidding were practically null in all ejidos (<0.3%) except for Noh Bec where producer accuracies were 1% for disturbed and 3% for deforested areas, showing very little effectiveness of CLASlite in detecting selective logging forest impacts in the Selva Maya of the Yucatan Peninsula.

Discussion
Our results show that even with high logging intensities and in large-scale community forest management operations in the Selva Maya, CLASlite could not accurately detect forest disturbance from timber harvesting using Landsat images (30 m spatial resolution).In contrast, the CLAS algorithms worked well in detecting forest degradation from logging in other countries of the tropics (for example, Belize, Madagascar, Malaysia and Brazil).However, in Madagascar (Allnutt et al. 2013) and Belize (Chicas et al. 2016), the authors did not distinguish selective logging from other disturbance sources, including slash and burn cultivation and road construction.In other studies where degradation caused from selective logging was identified and mapped, logging intensities were much higher, such as in Malaysia (50 to 150 m 3 ha -1 : Bryan et al. 2013) and Brazil (23 m 3 ha -1 : Asner et al. 2005;Asner 2009); whereas in our study site, the ejidos with the highest logging intensities, Nohbec and Petcacab, were only 7 and 5 m 3 ha -1 respectively.Thus, the use of CLASlite by forest managers and conservation organizations to monitor degradation by selective logging on the Yucatan Peninsula would not be feasible in this region with very low logging intensities (Ellis et al. 2019).
Even though, CLASlite could conceivable detect forest disturbance at resolution of 0.1 ha (1000 m 2 ) (Asner et al. 2010;Asner 2014), forest impacts from felling gaps (100 to 500 m 2 ) or skid trails (4 or 5 m wide) in the Selva Maya were not large enough to be detected by CLASlite with Landsat imagery.Even larger deforestation impacts, logging yards (approximately 1000 m 2 ) and recently constructed logging roads, were very poorly detected by CLASlite.Besides very low logging intensities in the Selva Maya, ejidos tend to utilize and maintain old logging roads, minimizing forest impacts from new roads.Moreover, they conduct improved management practices such as directional felling and use of modified agricultural tractor that also reduce biomass impacts (Ellis et al. 2019).The rapid recovery of photosynthetic vegetation after logging in the region (Martínez and García 2007; Toledo 2009) may further complicate monitoring selective logging disturbance using Landsat images (Broadbent et al. 2006) with the CLASlite monitoring system.Impacts from selective logging of up to 400 to 600 m 2 have been detected with Landsat images in several remote sensing studies (Asner et al. 2005;Tritsch et al. 2016;Vogelmann et al. 2017) and higher accuracies in detecting forest disturbance from selective logging in the same study area have been obtained using vegetation indices such as NDVI and Tasseled Cap.However, the automated CLASlite system oriented towards forest managers and conservation practitioners with little remote sensing experience is not sufficiently accurate for detecting selective logging disturbance.
The only study that applied CLASlite and Landsat images to measure forest degradation on the Yucatan Peninsula was conducted in the state of Yucatan in agricultural areas with slash and burn cultivation where the results of fractional and forest cover maps were used to model biomass  -Hernandez (2017) focuses on evaluating temporal biomass changes as a proxy for degradation, evaluating for the most part disturbances from forest clearings typically larger than 1000 m 2 .These clearings are cultivated for 2 to 3 years and then left to fallow for different periods until the cycle repeats, requiring denser time series to detect forest biomass dynamics and potential degradation.In comparison, we evaluate disturbance from selective logging which occurs at a much smaller scale, longer rotations and does not involve a change in forest cover.We therefore recommend CLASlite as tool to detect and monitor forest impacts greater than 900 m 2 resulting from other disturbances, for example, slash and burn agriculture, pasture establishment or fire.

SilvaFig. 2 .
Fig. 2. Images used and results of the analysis using the CLASlite 3.3 software.

Table 1 .
Adjusted confusion matrices of CLASlite 3.3 results for each ejido ACA.Unlogged changes at regional scale (30 by 40 km) (Romero-Sanchez and Ponce-Hernandez 2017).The study of Romero-Sanchez and Ponce