Next Article in Journal
Building Resistance and Resilience: Regeneration Should Not be Left to Chance
Next Article in Special Issue
Mapping Progression and Severity of a Southern Colorado Spruce Beetle Outbreak Using Calibrated Image Composites
Previous Article in Journal
Machine Learning Approaches for Estimating Forest Stand Height Using Plot-Based Observations and Airborne LiDAR Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Detecting and Attributing Drivers of Forest Disturbance in the Colombian Andes Using Landsat Time-Series

by
Paulo J. Murillo-Sandoval
1,2,3,*,
Thomas Hilker
1,†,
Meg A. Krawchuk
4 and
Jamon Van Den Hoek
3
1
Department of Forest Engineering, Resources and Management, Oregon State University, Corvallis, OR 97331, USA
2
Departamento de Topografía, Facultad de Tecnologías, Universidad del Tolima, Ibagué 73000-6299, Colombia
3
Geography and Geospatial Science, College of Earth, Ocean and Atmospheric Sciences, Oregon State University, Corvallis, OR 97331, USA
4
Department of Forest Ecosystems and Society, Oregon State University, Corvallis, OR 97331, USA
*
Author to whom correspondence should be addressed.
Deceased.
Forests 2018, 9(5), 269; https://doi.org/10.3390/f9050269
Submission received: 12 March 2018 / Revised: 23 April 2018 / Accepted: 10 May 2018 / Published: 15 May 2018
(This article belongs to the Special Issue Mapping Forest Health Using Moderate Resolution Satellites)

Abstract

:
The Colombian Andes foothills have seen an expansion of forest disturbance since the 1950s. While understanding the drivers of disturbance is important for quantifying the implications of land use change on regional biodiversity, methods for attributing disturbance to specific drivers of change at a high temporal and spatial resolution are still lacking in the Andes region, in part due to persistent cloud cover. Using 20 years of Landsat images (1996–2015) covering Picachos National Park in the Colombian Andes, we detected sub-annual forest cover disturbances using the Breaks For Additive Season and Trend (BFAST) Monitor algorithm; characterized different types of disturbance using spectral, spatial, and topographic indicators; and attributed causes of forest disturbance such as conversion to pasture, conversion to agriculture, and non-stand replacing disturbance (i.e., thinning) using a Random Forest (RF) classifier. Conversion to pasture has been the main driver of forest disturbance in Picachos, responsible for 11,395 ± 72 ha (17%) of forest cover loss, followed by non-stand replacing disturbance and conversion to agriculture. Disturbance detection had 96% overall agreement with validation data, although we had a high omission error of 21% primarily associated with forest to agriculture conversion. Other change drivers had a much more reliable attribution with forest to pasture conversion or non-stand-replacing disturbance, showing only 1–5% commission and 2–14% omission errors. Our results provide spatially-explicit information on sub-annual disturbances and associated drivers of change that are necessary for evaluating and improving domestic conservation efforts and establishing systematic ecological observations, which is currently absent from Colombia. While effective at revealing forest change dynamics in a geographically remote and socio-politically complex region like Picachos, our approach is highly automated and it can be easily extended to the rest of Andes-Amazon transition belt where low availability of remote sensing data and high cloud cover impede efforts at consistent monitoring of forest cover change dynamics and drivers.

Graphical Abstract

1. Introduction

Colombia is a “megadiverse” country given its high level of endemism and species richness [1]. However, over the last 50 years, the spread of informal settlements into forested regions, internal social conflict, and the growing demand for agricultural land have increased the pressure on Colombia’s forest ecosystems even in remote or protected regions [2]. As a consequence, between 1990 and 2015, Colombia lost more than six million hectares of forest, and more significantly, witnessed a 44% increase in deforestation in 2016 compared to 2015 [3]. Despite the duration and geographic breadth of these changes, there remains limited understanding of the relative influence and timing of different drivers on forest cover change across much of Colombia.
The Colombian Andes foothills were rather inaccessible three decades ago but have seen recent forest disturbance [4,5,6] due to favorable incentives for cattle ranching in the 1980s [7], the enlargement of coca cultivation in the 1990s through 2002 [2], and the more recent expansion of cattle ranching and other agricultural land uses [8,9]. Though National Protected Areas make up 12% of Colombia’s land mass [10], protected areas are also threatened by migration and displacement, economic development, and illicit cultivation of coca (Erythroxylum coca Lam.) [11], with attendant effects on socio-ecological relationships [4,7,12,13,14].
Despite many methods having been recently developed for change detection [15], there remains a lack of studies that identify drivers of disturbances with high spatial and temporal detail in rich biodiversity regions, and the Colombian Andes-Amazon area is no exception. Previous studies have identified drivers of Colombian forest cover change at the sub-national level through comparisons of two or three Landsat images collected five or more years apart [2,10]. However, a bi- or tri-temporal change assessment limits our ability to monitor more temporally complex or spectrally subtle land cover change processes that can only be detected by annual or sub-annual time-series analysis [16]. Other studies have used more regularly collected Moderate Resolution Imaging Spectroradiometer images for measuring the long-term effectiveness of protected areas at limiting deforestation [17,18] or detecting the main drivers of forest cover change at the municipality-level [12]. However, those analyses have incorporated coarse spatial resolution data (>250 m) that are unable to capture small-scale processes of change common in mountains and inaccessible places [19,20]. Because of these limitations, such studies have not been able to gauge the effect of drivers of forest cover change with adequate spatial or temporal detail to support regional forest conservation efforts.
Remote sensing time-series analysis that exploits the full archive of available imagery offers the means of capturing disturbance processes over their entire duration, whether short- or long-term. A diversity of time-series algorithms has been developed to monitor forest disturbance including LandTrendr [21], the vegetation change tracker (VCT) [22], and trend detection approaches based on “best-available-pixel” composites [23,24,25]. However, these approaches produce annual disturbance maps, which may miss short-lived or otherwise sub-annual changes [26] such as post-disturbance recovery or changes in condition (low impact selective logging, small induced fires, and foliage condition) that are better captured using dense time-series [19,27,28]. Further, while such algorithms have been reliable in regions with relatively low cloud cover and a relatively uninterrupted Landsat archive of >30 years (e.g., United States and Canada; [21,23,24,29], these disturbance detection algorithms are broadly untested in regions with regular cloud cover or satellite image scarcity.
The development of algorithms such as the Noise Insensitive Trajectory Algorithm [30], the Detecting Breakpoints and Estimating Segments in Trend [31], and the Breaks For Additive Season and Trend (BFAST)-Monitor [28,32] have improved the capacity for disturbance detection despite persistent cloud cover by building a dense time-series stack of all possible dates of satellite observation. Of these, BFAST-Monitor has had the widest adoption for identifying disturbance events, having been used to model seasonal and trend components in tropical forested regions [33] with a single vegetation index (e.g., NDVI) [34,35], multiple vegetation indices [36], or by taking into account spatial context [37,38].
The detection of forest disturbance events and processes has been highly operationalized at the global scale [39], yet the identification of proximate drivers of forest disturbance remains a challenge. Most driver detection studies have focused on North American forests and used spectral-temporal metrics to map common drivers such as fires, clear-cut harvests, insect outbreaks, and road expansion [21,23,24]. While these approaches have been shown to be effective in data rich contexts, they have not been broadly tested in other regions complicated by fewer satellite observations, persistent cloud cover, or smaller disturbance patches.
In this study, our objective is the detection and attribution of drivers of forest disturbance in an area with abrupt and subtle gradual changes using Landsat time-series (1998–2015). Using the foothills of Picachos National Park as a case of study, we (1) detected forest disturbances using BFAST-Monitor and (2) attributed the disturbance driver using a Random forest classifier and various time-series metrics. We focused on the most common drivers of disturbance, including conversion to pasture, conversion to agriculture, and non-stand replacing disturbance. These drivers were selected given their recognized negative impact on landscape fragmentation, forest structure, and biodiversity [14,40]. Drivers of forest disturbance were characterized using metrics derived from time-series trajectories such as spectral metrics (e.g., average spectral values before and after disturbance event), trend metrics related to magnitude of change, spatial pattern of disturbance (e.g., patch area and perimeter), and topographic metrics (e.g., elevation). Using previous metrics and training data, we identified the most probable driver of disturbance. The results contribute to the adoption of conservation monitoring strategies (e.g., zonification processes) in tropical forested areas that are threatened by unregulated and often undocumented human encroachment, which are difficult to monitor due to regular cloud cover and smaller, although progressive, disturbance patches.

2. Study Area

The study covers 67,284 ha in the foothills of Colombia’s Picachos National Park and includes the basins of the Templado, Platanillo, and Yulo rivers (Figure 1). Picachos is quite rugged with an elevational range of 300–2500 masl, has a temperature range of 5–25 °C, and annual accumulated rainfall varying from 1050–3250 mm [14]. The Park has great hydrological significance, as river headwaters here feed the basins of the Orinoco and Amazon rivers. A majority (62%) of the Picachos’ foothills are part of the Andes-Orinoquia biome characterized by broad-leaf evergreen forest, with additional land covers including pasture and mosaic cropland (shrub, herbaceous cover) [14].
Picachos has one of the largest, relatively intact regions of forest cover within Colombia, connecting the Orinoquia, Amazon, and Andes ecosystems [14], and was designated as a Category II Park by the International Union for Conservation of Nature (IUCN) in 1977 to conserve large-scale ecological processes, species, and ecosystems for current and future generations [41]. Although people are not legally allowed to reside within the Park, more than 200 families and 119 small houses and plots of 50–250 ha were documented in Picachos as of 2016 [14]. The settlement of peasant communities results from historical processes of colonization and forced displacement and has expanded present-day grazing land and subsistence agricultural production [8,14,42,43].

3. Materials and Methods

3.1. Image Processing

Imagery from the Landsat archive was selected given its long record of regular observation and intra-annual observation frequency of nearly eight images per year over Picachos, which supports sub-annual forest disturbance [6]; and its 30 m (0.09 ha) spatial resolution, which is suitable for detecting small disturbance patches common in Picachos (mean =2 ha). We used a time-series of 149 Level-1 Terrain-Corrected (L1T) Landsat 5, 7, and 8 surface reflectance images (path/row 8/58) from 1996 to 2015 with cloud cover less than or equal to 50%; 225 additional images with more than 50% cloud cover were removed prior to analysis. Clouds and cloud shadows, water, and snow were masked using the CFmask cloud-shadow mask product [44]. Cloud confidence and pixel description from pixel Quality Assessment (QA) band’s quality were used to remove medium and high cloud confidence, but also pixels surrounding high cloud confidence.
For each image in the time-series, we calculated NDMI (Normalized Difference Moisture Index; Equation (1)) [45], which combines short-wave infrared (SWIR) and near infrared (NIR) to measure canopy moisture content [46]. For the NDMI calculation, spectral bands B4 and B5 for Landsat 5 and 7 were used, whereas B5 and B6 were employed for the Landsat 8 case.
N D M I = N I R S W I R N I R + S W I R
NDMI has been successfully used to monitor disturbance and regrowth in tropical forests [47], has shown better sensitivity for disturbance detection than greenness indices such as Enhanced Vegetation Index (EVI) and the Normalized Difference Vegetation Index (NDVI) [48], and has provided a similarly high disturbance detection accuracy in tropical forests compared to Tasseled Cap Wetness (TCW) [49] or the Normalized Difference Fraction Index (NDFI) [33,50]. Additionally, NDMI is effective even with the level of intra-annual variability in image availability common in the Andes [33,34]. With NDMI, we produced a binary forest/non-forest mask using a Support Vector Machine algorithm in ENVI 5.1 on a cloud-free image on 22 February 1998. This mask represents all forested pixels at the beginning of the study period and was applied to all subsequent image dates to support the identification of forest disturbances, excluding other land cover dynamics.

3.2. Disturbance Detection Using BFAST-Monitor

Given its proven robustness for monitoring changes across various tropical forests [34,35,38], we selected the BFAST-Monitor algorithm to map disturbances in Picachos. BFAST-Monitor iteratively decomposes time-series data into trend, seasonal, and noise components, and detects significant breakpoints where the difference between the actual and expected signal deviates beyond a 95% significance boundary [32,51].
We applied BFAST-Monitor in one-year sequential periods starting in 1999 using a first-order harmonic curve over the NDMI time-series; the 1996–1998 period was used as a historical baseline. Due to the irregularity of clear-sky Landsat observations, the trend component was omitted to avoid erroneous NDMI values in the monitoring period that may cause the detection of false breakpoints [34]. Although the breakpoint detection approach is robust to noise and other variability in the time-series [32], a threshold based on NDMI change magnitude was applied to avoid the detection of spurious breakpoints due to outliers in the monitoring period. To establish this threshold, we visually compared potential breakpoints with a WorldView-1 true color (0.5 m resolution) image from 12 December 2009, and selected an NDMI change magnitude threshold of −0.01. Therefore, all breakpoints with an NDMI decline lower in magnitude than this threshold were classified as disturbed. Disturbances identified by BFAST-Monitor involve three components: location, change date, and change magnitude. Individual disturbance pixels were grouped into patches based on Rook’s case adjacency; patches with an area less than 1 ha (~9 pixels) were discarded to avoid isolated pixels and the salt-pepper effect. Forest pixels that did not show any disturbance over our study period were considered stable forest (non-disturbed).

3.3. Disturbance Validation Using TimeSync

An accuracy assessment identifies the errors during the classification of disturbance, and the sample data can be used for estimating both accuracy and area along with the uncertainty of these estimates [52]. We used TimeSync [34,53], an interactive tool that supports the visual interpretation of forest cover change across the time-series, to validate the disturbance timing and location. We sampled 614 sites from across Picachos using a stratified sampling procedure following Cochran’s recommendations [54]; twenty-two sites were removed because of excessive cloud cover in either Landsat or high-resolution Google Maps-hosted reference. Within TimeSync and referencing Google Earth and RapidEye high-resolution images, we interpreted annual mean Landsat NDVI and NDMI composite time-series to identify disturbance dates. Of our 592 sample sites, we interpreted 211 disturbed sites and 381 non-disturbed sites. The initial area estimates were derived directly from a disturbed/non-disturbed map (i.e., pixel counting), and the area-adjusted error was measured using a stratified estimator for every class to remove the biased effect of unequal omission and commission errors (see Supplementary Material). For disturbances, a confusion matrix was created by cross-tabulating the BFAST Monitor-derived disturbance information with TimeSync-derived validation data. For drivers of disturbance, a confusion matrix from the Random Forest output was used. For both matrices, total, producer’s, and user’s accuracies were calculated.

3.4. Characterizing Drivers of Change

Disturbance detection offers fundamental information about changes in landscape form and function that result from drivers of change, which differ in their spectral-temporal response function [19,55]. Long-term, high magnitude disturbances are commonly related to land cover conversion while short-term, low magnitude disturbances are more commonly a result of land cover alteration (Figure 2). For example, conversion to agriculture might present an abrupt change followed by an annual cycling of spectral values indicative of vegetation harvest and growth. The spectral trend of non-stand replacing disturbance, on the other hand, has a subtle but prolonged decrease in condition, and conversion to pasture commonly presents an abrupt break followed by a sustained period of relatively low vegetative condition (e.g., Figure 3).
To effectively characterize each driver, we measured patch-level metrics based on spectral indicators, trend, spatial pattern, and topographic condition (Table 1). First, spectral metrics consisted of average and standard deviation of near infrared (NIR), shortwave infrared 1 (SWIR1), shortwave infrared 2 (SWIR2), and NDMI reflectance, before and after the date of disturbance [23,56]. Second, trend metrics were calculated using the mean NDVI trend at each disturbance pixel for each year using ordinary least squares (OLS) regression [26]; we also included the kurtosis and skewness for both the magnitude of change and NDVI trend [24]. Pixel-level trends and magnitudes were then averaged across each patch. Third, for the spatial pattern, we measured the patch area, perimeter [57], shape index [58], and fractal dimension [59]. Finally, for topographic condition, we measured the elevation, slope, and topographic position index defined as the difference between the elevation of a cell and the mean value of its eight surrounding cells [21].

3.5. Attributing Drivers of Change

From the total population of 3450 disturbance patches, we visually interpreted and categorized 576 (17%) patches as follows: conversion to agriculture (50 patches, 146 ha total); non-stand replacing disturbance (147 patches, 361 ha total); and conversion to pasture (379 patches, 1171 ha total). This categorization was carried out using land cover maps based on Corine Land Cover (CLC) data generated by Colombian government agencies for 2002, 2007, and 2012 [8,14], which provide the most reliable reference data for interpreting land cover conversion in Colombia. Additionally, we used high resolution (0.65 m) images hosted by Google Earth (17 January 2003; 31 December 2009) and RapidEye imagery (23 October 2014; 15 December 2014; 5 January 2015) to visually validate drivers in our sampled patches.
Patch-level metrics (Table 1) were used for the training and validation (60% and 40% of sampled patches, respectively) of a Random Forests (RF) classification of drivers for each visually interpreted disturbed patch [21,23,36,60]. With 36 input metrics, we limited the potential for multicollinearity by removing nine metrics with Pearson’s r values greater than 0.70. Drivers attributed to each disturbed patch using majority votes from the RF were compared with the TimeSync validation sample to measure the attribution accuracy. We also calculated the distance to the second most voted class as a measure of classification confidence, calculated as the ratio d2/d1, where d2 corresponds to the second most voted class and d1 is the assigned class [61]. Resulting values are between 0 (when d2 < d1) and 1 (when d2 = d1), and patches with classification confidence values higher than 0.6 were re-labeled as unclassified given the unreliable driver attribution [23].

4. Results

4.1. Disturbance Detection Agreement

The total disturbed forest area within the Picachos foothills from 1998–2015 was 14,213 ± 1095 ha, while stable forest was recorded across 53,076 ± 1094 ha. Our estimate of disturbed forest within Picachos is 14% more than the 12,128 ha of forest cover loss estimated by official reports based on CLC maps during the similar study period of 1999–2012. Using visual interpretation results at 592 validation pixels, we estimated a very high total accuracy (TA) of 96% ± 1 (see pixel-level confusion matrix and full area-weight matrix Tables S1 and S2 in supplementary material). The disturbance class had user’s (UA) and producer’s accuracies (PA) of 96% ± 2 and 79% ± 8, respectively; the stable forest class had UA of 94% ± 2 and PA of 99% ± 0.4. Disturbed pixels had low commission (4%) and moderate omission errors (21%), whereas stable forest pixels were more affected by commission (6%) and less so by omission errors (1%) (Table 2). The higher rate of omission for the disturbance class is a consequence of low image availability (e.g., average of 51 cloud-free Landsat observations in the period 1999–2015), which is much lower compared to other tropical forest disturbance studies that have used BFAST-Monitor [33,47].

4.2. Characterization of Drivers’

Overall, spectral and trend metrics had the most predictive power, followed by topographic and spatial pattern indicators. The mean value post-change NDMI was the most important metric in the RF classification, followed by post-change standard deviation SWIR2, NDVI trend, and mean pre-change NDMI (Figure 4). If these metrics were excluded from the RF model, the accuracy would decrease by 33%, 21%, 15%, and 13%, respectively. Other metrics such as mean post-change NIR (11%), standard deviation of NDVI trend (8%), elevation (7.5%), and median magnitude of NDMI (4.6%) had much lower respective importance.
The most important metrics for characterizing drivers were those related to post-disturbance NDMI and SWIR2, which reflects the variation in post-disturbance spectral trends across drivers. Spatial pattern metrics such as patch area or perimeter displayed little contribution to the classification, likely because the area of a disturbed patch was similar between drivers, and larger patches affected by other drivers agents (e.g., wildfires) are not common in the study area [10].
Drivers had varying degrees of separability across metrics. Considering the six most important metrics, each driver showed considerable variation in median and interquartile ranges (Figure 5). For less important metrics, there was greater overlap in the interquartile range. For example, non-stand replacing and conversion to agriculture are distributed in a wider elevation range than conversion to pasture, which is mostly located at lower elevations. Conversion to pasture was well distinguished from the other two drivers across metrics, while conversion to agriculture and non-stand replacing disturbance showed less spectral separability. Of note, non-stand replacing disturbance had the smallest standard deviation in NDVI, indicating the subtle spectral changes generated by forest thinning over time. Conversion to agriculture exhibited narrowed values, especially for NDVI trend and pre-change NDMI, with the latter indicating stable values before disturbance detection.

4.3. Driver Attribution Agreement

The total agreement of driver attribution at the weight-area level was 98%±0.008 using a 95% level of confidence (Table 3) (patch -level confusion matrix and full area-weight matrix are in Supplementary Tables S3 and S4, respectively). Based on area-weighted methods, conversion to pasture exhibited the lowest omission and commission (1%) errors and was most readily discriminated from non-stand replacing disturbance and conversion to agriculture. Non-stand replacing disturbance had 5% and 12% omission and commission errors, respectively, and the lowest agreement was found in the conversion to agriculture class, which exhibited high omission (30%) and commission (4%) errors. Conversion to agriculture was commonly confused with non-stand replacing disturbance, likely due to spectral similarity in agricultural cycling and subtle short- and long-term changes in the forest (Figure 3). Making a clear distinction between both classes is difficult since agriculture has a relatively low coverage area and is generally surrounded by other land covers, which can introduce a mixed pixel effect, reducing the potential for accurate classification. Moreover, despite non-stand replacing disturbance often preceding deforestation driven by subsistence agriculture [16,36] or pasture expansion [62,63], we were unable to identify causal linkages between specific drivers.

4.4. Driver Dynamics

The annual area of disturbance increased from 1999 through 2007, except for 2005–2006 when the rate of disturbance was relatively low and 2003 when imagery was not available (Figure 6). This increase in disturbed area was mainly the result of pasture expansion for cattle grazing by people living within Picachos. Much of this conversion took place in favorable flat and productive soil along the Templado, Platanillo, and Yulo rivers, and was likely stimulated by economic incentives conferred by the Colombian government to cattle ranchers to expand their production in the 1980s and 1990s [14]. After 2007, the rate of annual disturbance had a more inconsistent trajectory. Of note, the increase in disturbed area in 2007 is not a product of increased Landsat image availability, since other years within our temporal study window had a similar or higher number of images.
Conversion to pasture led to the greatest area of disturbed land (17% of Picachos foothills; 11,425 ± 59 ha), which is slightly more than the CLC-calculated area of 10,827 ha during the similar period of 1999–2012. The second most dominant driver was non-stand replacing disturbance with an area of 1085 ± 76 ha; there is no documented precedent for this driver in CLC assessments. Conversion to agriculture contributed 290 ± 50 ha, which is also in good agreement with the CLC estimate of 380 ha. Finally, 1175 ha, corresponding to 8% of the total disturbed area, could not be confidently attributed and was left as unclassified.
Disturbance was most prominent along the margins of rivers with the Platanillo watershed experiencing the most disturbance; the earliest disturbances took place close to these rivers and have since expanded outwards (Figure 7). Conversion to pasture was found along the Templado, Plantanillo, and Yulo rivers, while conversion to agriculture and unclassified patches were mainly found near Yulo River. As Yulo River forms the southern border of Picachos, disturbance here likely reflects the influence of agricultural expansion and grazing from nearby peasant reserves. On the other hand, non-stand replacing disturbance was more dispersed over the entire study area and even evident at higher elevations. Notably, conversion to pasture and non-stand replacing disturbance are commonly in close proximity to each other, suggesting a gradual land conversion process.

5. Discussion

This study is the most comprehensive analysis of disturbances and drivers in the Andes, and highlights the value of using the full time-series of Landsat observations within an automated disturbance detection and driver attribution framework. Given the high spatio-temporal resolution, this approach was able to map land cover dynamics that are broadly absent from Colombia’s Corine Land Cover mapping framework, which records changes on a five-year timescale. Though the current study was constrained to the Picachos foothills, the disturbances and drivers in Picachos are typical across the rest of the Andes–Amazonia transition belt region [7,64], especially pasture expansion and intensification [65]. In Picachos, as well as the rest of the Andes-Amazon region, conversion to pasture is highly related to cattle expansion [40] and has detrimental effects on landslide susceptibility in high relief regions [7,66], as well as ecosystem connectivity between the Andes and Orinoco and Amazon ecosystems [2]. Given the common challenges to disturbance detection and driver attribution across the Andes-Amazonia region, this study offers a proof-of-concept for a robust ecological monitoring system that has applicability well beyond the borders of Picachos, including the surrounding protected areas of Tinigua and Macarena.
Though Picachos is a strategic regional corridor between the Andes and Amazonia, there remains limited biophysical and sociocultural information to support effective decision-making processes. Detecting patterns of pasture intensification with remote sensing-derived spatio-temporal information is a critical tool in designing and evaluating land management policies that support desired, sustainable land use strategies [67]. The results of this study record in great detail the short- and long-term effects of human influence on regional forest cover dynamics, facilitating data-driven dialogue with the community and developing land-use planning such as better management zonification in the foothills.
We detected, characterized, and attributed disturbances with relatively high agreement, but limitations of this approach may be summarized as four aspects. First, we were unable to identify multiple, sequential disturbances at a given pixel over the time-series, nor identify causal linkages between specific drivers, such as non-stand replacing disturbance preceding agricultural conversion [16,36] or pasture expansion [62,63]. In future research, engaging local community knowledge could be relevant for establishing linkages between localized drivers and disturbances and identifying the most relevant metric for disturbance detection or driver attribution [68,69]. Second, the training data used in building the RF classifier were not randomly sampled but rather sampled based on the availability of reference datasets (e.g., CLC, RapidEye, and high-resolution images). Despite limitations in our sampling approach, the out-of-bag sample provided by the RF model gave us an unbiased estimate of the classification error, which achieved a value of 13%. Third, the discrimination and attribution of drivers was more accurate for conversion to pasture and non-stand replacing disturbance than for conversion to agriculture, and distinguishing between non-stand replacing disturbance and conversion to agriculture remains a challenge. Fourth, although overall agreement was accurate for both disturbance and drivers’ attribution detection, omission errors are significantly high, especially for disturbances (21%) and conversion to pasture (30%), which reflects how low image availability after cloud removal may affect area estimates.
Future work should thus include additional information on disturbance detection and driver attribution, for instance, the previous number of breaks before a given disturbance event [70], the incorporation of change duration [21] such as greening periods (positive trends) or browning periods (negative trends), and spatial context (pixel neighborhood), which has shown improvement in the mapping of changes [38]. Finally, ensemble methods which account for the contributions of many different change detection algorithms have enhanced forest change mapping [71].

6. Conclusions

This study integrated spectral, trend, pattern, and topographic metrics derived from BFAST-Monitor and the Landsat archive to detect disturbances and identify related drivers in the Picachos foothills. We have been able to consistently differentiate the major drivers of forest disturbance within Picachos National Park. We confirmed that conversion pasture is the main cause of change, especially from 2007–2014, followed by non-stand replacing disturbance and conversion to agriculture. The results of this study should encourage the remote sensing community to further develop robust frameworks using dense time-series metrics to analyze the drivers and consequences of human-induced changes, especially in regions with low image availability or regular cloud cover, in support of sustainable resource use strategies that support people and forests alike.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4907/9/5/269/s1, Table S1: Confusion matrix at pixel-level for disturbances and stable forest, Table S2: Full area-weight confusion matrix for disturbances and stable forest, Table S3: Confusion matrix at object-level for drivers of disturbance using 40% data training from Random Forest, Table S4: Full area-weight confusion matrix for drivers of forest disturbance.

Author Contributions

P.J.M.-S., T.H. and J.V.D.H. conceived and designed the study. P.J.M.-S. performed the data processing. P.J.M.-S., J.V.D.H., and M.A.K. analyzed and interpreted the data, and wrote the paper.

Acknowledgments

Oregon State University and Universidad del Tolima (Colombia) support P.J.M.-S. Cormacarena and Parques Nacionales de Colombia provided the Corine Land Cover dataset and Rapideye Images. We also thank the United States Geological Survey (USGS), which maintains the Landsat archive and makes it available.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McNeely, J. Conserving forest biodiversity in times of violent conflict. Oryx 2003, 37, 142–152. [Google Scholar] [CrossRef]
  2. Etter, A.; Mcalpine, C.; Phinn, S.; Pullar, D.; Possingham, H. Characterizing a tropical deforestation wave: A dynamic spatial analysis of a deforestation hotspot in the Colombian Amazon. Glob. Chang. Biol. 2006, 12, 1409–1420. [Google Scholar] [CrossRef]
  3. Ministerio de Ambiente. Estrategia integral de control a la deforestacion y gestion de los bosques (EICDGB); Ministerio de Ambiente: Bogota, Colombia, 2017. [Google Scholar]
  4. Etter, A.; McAlpine, C.; Possingham, H. Historical Patterns and Drivers of Landscape Change in Colombia Since 1500: A Regionalized Spatial Approach. Ann. Assoc. Am. Geogr. 2008, 98, 2–23. [Google Scholar] [CrossRef]
  5. Stibig, H.-J.; Achard, F.; Carboni, S.; Raši, R.; Miettinen, J. Change in tropical forest cover of Southeast Asia from 1990 to 2010. Biogeosciences 2014, 11, 247–258. [Google Scholar] [CrossRef] [Green Version]
  6. Murillo-Sandoval, P.J.; Van Den Hoek, J.; Hilker, T. Leveraging Multi-Sensor Time Series Datasets to Map Short- and Long-Term Tropical Forest Disturbances in the Colombian Andes. Remote Sens. 2017, 9, 179. [Google Scholar] [CrossRef]
  7. Armenteras, D.; Rodríguez, N.; Retana, J.; Morales, M. Understanding deforestation in montane and lowland forests of the Colombian Andes. Reg. Environ. Chang. 2011, 11, 693–705. [Google Scholar] [CrossRef]
  8. Unidad Administrativa Especial del Sistema de Parques Nacionales Naturales Convenio de Asociación Tripartita P.E. GDE.1.4.7.1.14.022 Suscrito entre Parques Nacionales Naturales, Cormacarena y Patrimonio Natural Fondo para la Diversidad y Áreas Protegidas; UAESPNN–Dirección Territorial Costa Orinoquia: Bogota, Colombia, 2015; ISBN 9789563589610.
  9. Rodríguez, N. Deforestación y cambio en la cobertura del suelo en colombia: Dinámica espacial, factores de cambio y modelamiento. Ph.D. Thesis, Universitat Autònoma de Barcelona, Barcelona, Spain, 2011. [Google Scholar]
  10. Rodríguez, N.; Armenteras, D.; Retana, J. Effectiveness of protected areas in the Colombian Andes: Deforestation, fire and land-use changes. Reg. Environ. Chang. 2013, 13, 423–435. [Google Scholar] [CrossRef]
  11. Dávalos, L.M.; Bejarano, A.C.; Hall, M.A.; Correa, H.L.; Corthals, A.; Espejo, O.J. Forests and drugs: coca-driven deforestation in global biodiversity hotspots, Supporting Information. Environ. Sci. Technol. 2011, 45, 1219–1227. [Google Scholar] [CrossRef] [PubMed]
  12. Sanchez-Cuervo, A.; Aide, T.M. Identifying hotspots of deforestation and reforestation in Colombia (2001-2010): Implications for protected areas. Ecosphere 2013, 4, 1–20. [Google Scholar] [CrossRef]
  13. Sánchez-Cuervo, A.M.; Aide, T.M.; Clark, M.L.; Etter, A. Land Cover Change in Colombia: Surprising Forest Recovery Trends between 2001 and 2010. PLoS ONE 2012, 7, e43943. [Google Scholar] [CrossRef] [PubMed]
  14. Unidad Administrativa Especial del Sistema de Parques Nacionales Naturales Plan de Manejo Coordillera de los Picachos (under review); UAESPNN–Dirección Territorial Costa Orinoquia: Neiva, Huila Colombia, 2016; p. 192.
  15. Cohen, W.B.; Healey, S.P.; Yang, Z.; Stehman, S.V.; Brewer, C.K.; Brooks, E.B.; Gorelick, N.; Huang, C.; Hughes, M.J.; Kennedy, R.E.; et al. How Similar Are Forest Disturbance Maps Derived from Different Landsat Time Series Algorithms? Forests 2017, 8, 98. [Google Scholar] [CrossRef]
  16. Lambin, E.F. Monitoring forest degradation in tropical regions by remote sensing: some methodological issues. Glob. Ecol. Biogeogr. 1999, 8, 191–198. [Google Scholar] [CrossRef]
  17. Cortes, A.; Sergio, D. Assessing the Ground Truth Drivers of Land Cover and Land Use Changes at a Local Scale in Cundinamarca and Tolima Departments in Colombia. Master’s Thesis, University of Southampton, Southampton, UK, 2012. [Google Scholar]
  18. Salazar Villegas, M.H. Effectiveness of Colombia’ s Protected Areas in Preventing Evergreen Forest Loss: A Study Using Terra-i Near Real-Time Monitoring System. Master’s Thesis, Technische Universitat Dresden, Dresden, Germany, 2013. [Google Scholar]
  19. Pasquarella, V.J.; Holden, C.E.; Kaufman, L.; Woodcock, C.E. From imagery to ecology: leveraging time series of all available Landsat observations to map and monitor ecosystem state and dynamics. Remote Sens. Ecol. Conserv. 2016, 1–19. [Google Scholar] [CrossRef]
  20. Fjeldsa, J.; Ehrlich, D.; Lambin, E.; Prins, E. Are biodiversity “hotspots” correlated with current ecoclimatic stabil ity? A pilot study using the NOAA-AVHRR remote sensing data. Biodivers. Conserv. 1997, 6, 401–422. [Google Scholar] [CrossRef]
  21. Kennedy, R.E.; Yang, Z.; Braaten, J.; Copass, C.; Antonova, N.; Jordan, C.; Nelson, P. Attribution of disturbance change agent from Landsat time-series in support of habitat monitoring in the Puget Sound region, USA. Remote Sens. Environ. 2015, 166, 271–285. [Google Scholar] [CrossRef]
  22. 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]
  23. Hermosilla, T.; Wulder, M.A.; White, J.C.; Coops, N.C.; Hobart, G.W. Regional detection, characterization, and attribution of annual forest change from 1984 to 2012 using Landsat-derived time-series metrics. Remote Sens. Environ. 2015, 170, 121–132. [Google Scholar] [CrossRef]
  24. Pickell, P.D.; Hermosilla, T.; Coops, N.C.; Masek, J.G.; Franks, S.; Huang, C. Monitoring anthropogenic disturbance trends in an industrialized boreal forest with Landsat time series. Remote Sens. Lett. 2014, 5, 783–792. [Google Scholar] [CrossRef]
  25. Frazier, R.J.; Coops, N.C.; Wulder, M.A.; Kennedy, R. Characterization of aboveground biomass in an unmanaged boreal forest using Landsat temporal segmentation metrics. ISPRS J. Photogramm. Remote Sens. 2014, 92, 137–146. [Google Scholar] [CrossRef]
  26. Forkel, M.; Carvalhais, N.; Verbesselt, J.; Mahecha, M.D.; Neigh, C.S.R.; Reichstein, M. Trend Change detection in NDVI time series: Effects of inter-annual variability and methodology. Remote Sens. 2013, 5, 2113–2144. [Google Scholar] [CrossRef]
  27. Hilker, T.; Coops, N.C.; Coggins, S.B.; Wulder, M. a.; Brownc, M.; Black, T.A.; Nesic, Z.; Lessard, D. Detection of foliage conditions and disturbance from multi-angular high spectral resolution remote sensing. Remote Sens. Environ. 2009, 113, 421–434. [Google Scholar] [CrossRef]
  28. Vogelmann, J.E.; Gallant, A.L.; Shi, H.; Zhu, Z. Perspectives on monitoring gradual change across the continuity of Landsat sensors using time-series data. Remote Sens. Environ. 2016, 185, 258–270. [Google Scholar] [CrossRef]
  29. Pflugmacher, D.; Cohen, W.B.; Kennedy, R.E. Using Landsat-derived disturbance history (1972-2010) to predict current forest structure. Remote Sens. Environ. 2012, 122, 146–165. [Google Scholar] [CrossRef]
  30. Alonzo, M.; Van Den Hoek, J.; Ahmed, N. Capturing coupled riparian and coastal disturbance from industrial mining using cloud-resilient satellite time series analysis. Sci. Rep. 2016, 1–13. [Google Scholar] [CrossRef] [PubMed]
  31. Jamali, S.; Jonsson, P.; Eklundh, L.; Ardo, J.; Seaquist, J. Detecting changes in vegetation trends using time series segmentation. Remote Sens. Environ. 2015, 156, 182–195. [Google Scholar] [CrossRef]
  32. Verbesselt, J.; Zeileis, A.; Herold, M. Near real-time disturbance detection using satellite image time series. Remote Sens. Environ. 2012, 123, 98–108. [Google Scholar] [CrossRef]
  33. Schultz, M.; Clevers, J.G.P.W.; Carter, S.; Verbesselt, J.; Avitabile, V.; Quang, H.V.; Herold, M. Performance of vegetation indices from Landsat time series in deforestation monitoring. Int. J. Appl. Earth Obs. Geoinf. 2016, 52, 318–327. [Google Scholar] [CrossRef]
  34. DeVries, B.; Verbesselt, J.; Kooistra, L.; Herold, M. Robust monitoring of small-scale forest disturbances in a tropical montane forest using Landsat time series. Remote Sens. Environ. 2015, 161, 107–121. [Google Scholar] [CrossRef]
  35. Dutrieux, L.P.; Verbesselt, J.; Kooistra, L.; Herold, M. Monitoring forest cover loss using multiple data streams, a case study of a tropical dry forest in Bolivia. ISPRS J. Photogramm. Remote Sens. 2015, 107, 112–125. [Google Scholar] [CrossRef]
  36. Devries, B.; Pratihast, A.K.; Verbesselt, J.; Kooistra, L.; Herold, M. Characterizing forest change using community-based monitoring data and landsat time series. PLoS ONE 2016, 11, 1–25. [Google Scholar] [CrossRef] [PubMed]
  37. Hamunyela, E.; Reiche, J.; Verbesselt, J.; Herold, M. Using Space-Time Features to Improve Detection of Forest Disturbances from Landsat Time Series. Remote Sens. 2017, 9, 515. [Google Scholar] [CrossRef]
  38. Hamunyela, E.; Verbesselt, J.; Herold, M. Using spatial context to improve early detection of deforestation from Landsat time series. Remote Sens. Environ. 2016, 172, 126–138. [Google Scholar] [CrossRef]
  39. Hansen, M.C.; Potapov, P.V.; Moore, R.; Hancher, M.; Turubanova, S.A.; Tyukavina, A. High-resolution global maps of forest cover change. Science 2013, 342, 850–853. [Google Scholar] [CrossRef] [PubMed]
  40. Etter, A.; McAlpine, C.; Wilson, K.; Phinn, S.; Possingham, H. Regional patterns of agricultural land use and deforestation in Colombia. Agric. Ecosyst. Environ. 2006, 114, 369–386. [Google Scholar] [CrossRef]
  41. Dudley, N. Guidelines for Applying Protected Area Management Categories; IUCN: Gland, Switzerland, 2008; p. 86. ISBN 9782831710860. [Google Scholar]
  42. Vásquez, T. El Papel Del Conflicto Armado En La Construcción Y Diferenciación Territorial De La Región De El Caguán Amazonía Occidental Colombiana. Ágora U.S.B. 2013, 14, 147–175. [Google Scholar] [CrossRef]
  43. Sánchez-Cuervo, A.M.; Aide, T.M. Consequences of the Armed Conflict, Forced Human Displacement, and Land Abandonment on Forest Cover Change in Colombia: A Multi-scaled Analysis. Ecosystems 2013, 16, 1052–1070. [Google Scholar] [CrossRef]
  44. Zhu, Z.; Woodcock, C.E. Object-based cloud and cloud shadow detection in Landsat imagery. Remote Sens. Environ. 2012, 118, 83–94. [Google Scholar] [CrossRef]
  45. Gao, B.C. NDWI—A normalized difference water index for remote sensing of vegetation liquid water from space. Remote Sens. Environ. 1996, 58, 257–266. [Google Scholar] [CrossRef]
  46. Hayes, D.J.; Cohen, W.B. Spatial, spectral and temporal patterns of tropical forest cover change as observed with multiple scales of optical satellite data. Remote Sens. Environ. 2007, 106, 1–16. [Google Scholar] [CrossRef]
  47. Devries, B.; Verbesselt, J. Tracking disturbance-regrowth dynamics in tropical forests using structural change detection and Landsat time series. Remote Sens. Environ. 2015, 169, 320–334. [Google Scholar] [CrossRef]
  48. Masek, J.G.; Hayes, D.J.; Joseph Hughes, M.; Healey, S.P.; Turner, D.P. The role of remote sensing in process-scaling studies of managed forest ecosystems. For. Ecol. Manag. 2015, 355, 109–123. [Google Scholar] [CrossRef]
  49. Jin, S.; Sader, S.A. Comparison of time series tasseled cap wetness and the normalized difference moisture index in detecting forest disturbances. Remote Sens. Environ. 2005, 94, 364–372. [Google Scholar] [CrossRef]
  50. Souza, C.M.; Roberts, D.A.; Cochrane, M.A. Combining spectral and spatial information to map canopy damage from selective logging and forest fires. Remote Sens. Environ. 2005, 98, 329–343. [Google Scholar] [CrossRef]
  51. Verbesselt, J.; Hyndman, R.; Newnham, G.; Culvenor, D. Detecting trend and seasonal changes in satellite image time series. Remote Sens. Environ. 2010, 114, 106–115. [Google Scholar] [CrossRef]
  52. Olofsson, P.; Foody, G.M.; Herold, M.; Stehman, S.V.; Woodcock, C.E.; Wulder, M.A. Good practices for estimating area and assessing accuracy of land change. Remote Sens. Environ. 2014, 148, 42–57. [Google Scholar] [CrossRef]
  53. Cohen, W.B.; Yang, Z.; Kennedy, R. Detecting trends in forest disturbance and recovery using yearly Landsat time series: 2. TimeSync—Tools for calibration and validation. Remote Sens. Environ. 2010, 114, 2911–2924. [Google Scholar] [CrossRef]
  54. Cochran, W.G. Sampling Techniques; John Wiley Sons: New York, NY, USA, 1977. [Google Scholar]
  55. Kennedy, R.E.; Andréfouët, S.; Cohen, W.B.; Gómez, C.; Griffiths, P.; Hais, M.; Healey, S.P.; Helmer, E.H.; Hostert, P.; Lyons, M.B.; et al. Bringing an ecological view of change to landsat-based remote sensing. Front. Ecol. Environ. 2014, 12, 339–346. [Google Scholar] [CrossRef] [Green Version]
  56. Healey, S.P.; Cohen, W.B.; Zhiqiang, Y.; Krankina, O.N. Comparison of Tasseled Cap-based Landsat data structures for use in forest disturbance detection. Remote Sens. Environ. 2005, 97, 301–310. [Google Scholar] [CrossRef]
  57. Bogaert, J.; Rousseau, R.; Van Hecke, P.; Impens, I. Alternative area-perimeter ratios for measurement of 2D shape compactness of habitats. Appl. Math. Comput. 2000, 111, 71–85. [Google Scholar] [CrossRef]
  58. McGarigal, K.; Tagil, S.; Cushman, S.A. Surface metrics: An alternative to patch metrics for the quantification of landscape structure. Landsc. Ecol. 2009, 24, 433–450. [Google Scholar] [CrossRef]
  59. Krummel, J.R.; Gardner, R.H.; Sugihara, G.; O’Neill, R.V. Landscape patterns in a distubed environment. OIKOS 1987, 48, 321–324. [Google Scholar] [CrossRef]
  60. Reents, C. Detection and Characterization of Forest Disturbances in California. Master’s Thesis, University of Illinois at Urbana-Champaign, Champaign, IL, USA, 2016. [Google Scholar]
  61. Mitchell, S.W.; Remmel, T.K.; Csillag, F.; Wulder, M.A. Distance to second cluster as a measure of classification confidence. Remote Sens. Environ. 2008, 112, 2615–2626. [Google Scholar] [CrossRef]
  62. Walker, R.; Browder, J.; Arima, E.; Simmons, C.; Pereira, R.; Caldas, M.; Shirota, R.; Zen, S. de Ranching and the new global range: Amazonia in the 21st century. Geoforum 2009, 40, 732–745. [Google Scholar] [CrossRef]
  63. Fernández, G.; Obermeier, W.; Gerique, A.; Sandoval, M.; Lehnert, L.; Thies, B.; Bendix, J. Land Cover Change in the Andes of Southern Ecuador—Patterns and Drivers. Remote Sens. 2015, 7, 2509–2542. [Google Scholar] [CrossRef]
  64. Achard, F.; Beuchle, R.; Mayaux, P.; Stibig, H.J.; Bodart, C.; Brink, A.; Carboni, S.; Desclée, B.; Donnay, F.; Eva, H.D.; et al. Determination of tropical deforestation rates and related carbon losses from 1990 to 2010. Glob. Chang. Biol. 2014, 20, 2540–2554. [Google Scholar] [CrossRef] [PubMed]
  65. Piquer-Rodríguez, M.; Butsic, V.; Gärtner, P.; Macchi, L.; Baumann, M.; Gavier Pizarro, G.; Volante, J.N.; Gasparri, I.N.; Kuemmerle, T. Drivers of agricultural land-use change in the Argentine Pampas and Chaco regions. J. Appl. Geogr. 2018, 91, 111–122. [Google Scholar] [CrossRef]
  66. Restrepo, J.D.; Kettner, A.J.; Syvitski, J.P.M. Recent deforestation causes rapid increase in river sediment load in the Colombian Andes. Anthropocene 2015, 10, 13–28. [Google Scholar] [CrossRef]
  67. Rufin, P.; Müller, H.; Pflugmacher, D.; Hostert, P. Land use intensity trajectories on Amazonian pastures derived from Landsat time series. Int. J. Appl. Earth Obs. Geoinf. 2015, 41, 1–10. [Google Scholar] [CrossRef]
  68. Miranda, J.J.; Corral, L.; Blackman, A.; Asner, G.; Lima, E. Effects of Protected Areas on Forest Cover Change and Local Communities: Evidence from the Peruvian Amazon. World Dev. 2016, 78, 288–307. [Google Scholar] [CrossRef]
  69. Davies, T.; Everard, M.; Horswell, M. Community-based groundwater and ecosystem restoration in semi-arid north Rajasthan (3): Evidence from remote sensing. Ecosyst. Serv. 2016, 21, 20–30. [Google Scholar] [CrossRef]
  70. Watts, L.M.; Laffan, S.W. Effectiveness of the BFAST algorithm for detecting vegetation response patterns in a semi-arid region. Remote Sens. Environ. 2014, 154, 234–245. [Google Scholar] [CrossRef]
  71. Healey, S.P.; Cohen, W.B.; Yang, Z.; Kenneth Brewer, C.; Brooks, E.B.; Gorelick, N.; Hernandez, A.J.; Huang, C.; Joseph Hughes, M.; Kennedy, R.E.; et al. Mapping forest change using stacked generalization: An ensemble approach. Remote Sens. Environ. 2018, 204, 717–728. [Google Scholar] [CrossRef]
Figure 1. Picachos National Park foothills study area, outlined in red.
Figure 1. Picachos National Park foothills study area, outlined in red.
Forests 09 00269 g001
Figure 2. Conceptual flowchart for disturbance detection and driver attribution.
Figure 2. Conceptual flowchart for disturbance detection and driver attribution.
Forests 09 00269 g002
Figure 3. Example of trajectories and disturbance detection using all available observations and associated drivers of change. This example illustrates the presence of abrupt and gradual changes: (a) Forest-to-agriculture conversion with a cyclical pattern after the breakpoint; (b) Non-stand replacing disturbance where NDMI trend exhibits a subtle decrease over time; and (c) Forest to pasture conversion where an abrupt breakpoint is followed by low, stable NDMI values. Vertical red lines indicate a breakpoint, blue lines indicate a seasonal + trend fitted model, and M denotes the magnitude of change.
Figure 3. Example of trajectories and disturbance detection using all available observations and associated drivers of change. This example illustrates the presence of abrupt and gradual changes: (a) Forest-to-agriculture conversion with a cyclical pattern after the breakpoint; (b) Non-stand replacing disturbance where NDMI trend exhibits a subtle decrease over time; and (c) Forest to pasture conversion where an abrupt breakpoint is followed by low, stable NDMI values. Vertical red lines indicate a breakpoint, blue lines indicate a seasonal + trend fitted model, and M denotes the magnitude of change.
Forests 09 00269 g003
Figure 4. Metric importance ranking for discriminating drivers based on the mean decrease accuracy of the RF classification.
Figure 4. Metric importance ranking for discriminating drivers based on the mean decrease accuracy of the RF classification.
Forests 09 00269 g004
Figure 5. Boxplots of the relationships between the eight most important metrics for classifying drivers of change showing median (black line) and interquartile range.
Figure 5. Boxplots of the relationships between the eight most important metrics for classifying drivers of change showing median (black line) and interquartile range.
Forests 09 00269 g005
Figure 6. Temporal distribution of drivers of change at patch-level. Year 2003 is not included due to the lack of Landsat imagery.
Figure 6. Temporal distribution of drivers of change at patch-level. Year 2003 is not included due to the lack of Landsat imagery.
Forests 09 00269 g006
Figure 7. Spatial distribution of drivers from 1999–2015 (upper). Disturbances are concentrated in close proximity to rivers. Temporal distribution of disturbance (below). Inset map represents disturbance locations of in Picachos’ highest elevations (A) and Platanillo settlement (B).
Figure 7. Spatial distribution of drivers from 1999–2015 (upper). Disturbances are concentrated in close proximity to rivers. Temporal distribution of disturbance (below). Inset map represents disturbance locations of in Picachos’ highest elevations (A) and Platanillo settlement (B).
Forests 09 00269 g007
Table 1. Summary of patch-level metrics used to characterize drivers of change. SRTM DEM = Shuttle Radar Topography Mission Digital Elevation Model (90 m).
Table 1. Summary of patch-level metrics used to characterize drivers of change. SRTM DEM = Shuttle Radar Topography Mission Digital Elevation Model (90 m).
Metrics EmployedSource
Spectral summary metrics
Average spectral value pre-changeNIR, SWIR1, SWIR2, NDMI
Average spectral value post-changeNIR, SWIR1, SWIR2, NDMI
Standard deviation spectral value pre-changeNIR, SWIR1, SWIR2, NDMI
Standard deviation spectral value post-changeNIR, SWIR1, SWIR2, NDMI
Trend time-series metrics
Aggregated annual trend (mean, median)NDVI
Average and median change magnitudeNDMI
Standard deviation change magnitudeNDMI
Kurtosis and skewness change magnitudeNDMI
Pattern metrics
Area
Perimeter
Shape index
Fractal dimension
Topographic indicators
ElevationSTRM DEM
SlopeSTRM DEM
Topographic position indexSTRM DEM
NIR: near infrared, SWIR1: shortwave infrared 1, SWIR2: shortwave infrared 2, NDMI: Normalized Difference Moisture Index
Table 2. Area-weighted error matrix comparing BFAST-Monitor-and TimeSync-derived classification at validation sites.
Table 2. Area-weighted error matrix comparing BFAST-Monitor-and TimeSync-derived classification at validation sites.
TimeSync Validation (Reference)
DisturbanceStable ForestProportion of Area Mapped (Wi)User’s AccuracyProducer’s AccuracyTotal Accuracy
BFAST-Monitor (Map)Disturbance0.1340.0040.140.96 ± 0.0240.79 ± 0.0800.96 ± 0.017
 
Stable Forest0.030.820.860.96 ± 0.0200.99 ± 0.004
Total0.1860.8111
Table 3. Area-weighted error matrix for disturbance driver attribution based on Random Forest output.
Table 3. Area-weighted error matrix for disturbance driver attribution based on Random Forest output.
Reference
Conversion to PastureConversion to AgricultureNon-Stand ReplacingProportion of Area Mapped (Wi)User’s AccuracyProducer’s AccuracyTotal Accuracy
MapConversion to pasture0.853 0.0040.8570.99 ± 0.0050.99 ± 0.0050.98 ± 0.008
Conversion to agriculture 0.0210.0010.0220.96 ± 0.0570.7 ± 0.121 
 
Non-stand replacing0.0060.0090.1070.1210.88 ± 0.0520.95 ± 0.040 
 
Total0.8590.030.1121

Share and Cite

MDPI and ACS Style

Murillo-Sandoval, P.J.; Hilker, T.; Krawchuk, M.A.; Van Den Hoek, J. Detecting and Attributing Drivers of Forest Disturbance in the Colombian Andes Using Landsat Time-Series. Forests 2018, 9, 269. https://doi.org/10.3390/f9050269

AMA Style

Murillo-Sandoval PJ, Hilker T, Krawchuk MA, Van Den Hoek J. Detecting and Attributing Drivers of Forest Disturbance in the Colombian Andes Using Landsat Time-Series. Forests. 2018; 9(5):269. https://doi.org/10.3390/f9050269

Chicago/Turabian Style

Murillo-Sandoval, Paulo J., Thomas Hilker, Meg A. Krawchuk, and Jamon Van Den Hoek. 2018. "Detecting and Attributing Drivers of Forest Disturbance in the Colombian Andes Using Landsat Time-Series" Forests 9, no. 5: 269. https://doi.org/10.3390/f9050269

APA Style

Murillo-Sandoval, P. J., Hilker, T., Krawchuk, M. A., & Van Den Hoek, J. (2018). Detecting and Attributing Drivers of Forest Disturbance in the Colombian Andes Using Landsat Time-Series. Forests, 9(5), 269. https://doi.org/10.3390/f9050269

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