Next Article in Journal
Surface-Enhanced Raman Spectroscopy Based on a Silver-Film Semi-Coated Nanosphere Array
Next Article in Special Issue
Surface Heterogeneity-Involved Estimation of Sample Size for Accuracy Assessment of Land Cover Product from Satellite Imagery
Previous Article in Journal
Epoxy Composites with Reduced Graphene Oxide–Cellulose Nanofiber Hybrid Filler and Their Application in Concrete Strain and Crack Monitoring
Previous Article in Special Issue
A Novel Performance Assessment Approach Using Photogrammetric Techniques for Landslide Susceptibility Mapping with Logistic Regression, ANN and Random Forest
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 19
Issue 18
10.3390/s19183965
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

Climate-Triggered Insect Defoliators and Forest Fires Using Multitemporal Landsat and TerraClimate Data in NE Iran: An Application of GEOBIA TreeNet and Panel Data Analysis

by
Omid Abdi
Institute for Cartography, Department of Geosciences, Faculty of Environmental Sciences, TU Dresden, 01069 Dresden, Germany
Sensors 2019, 19(18), 3965; https://doi.org/10.3390/s19183965
Submission received: 6 August 2019 / Revised: 9 September 2019 / Accepted: 12 September 2019 / Published: 14 September 2019
(This article belongs to the Special Issue Remote Sensing of Earth Observation: Special Focus on Natural Resources and the Environment)
Browse Figures
Versions Notes

Abstract

:
Despite increasing the number of studies for mapping remote sensing insect-induced forest infestations, applying novel approaches for mapping and identifying its triggers are still developing. This study was accomplished to test the performance of Geographic Object-Based Image Analysis (GEOBIA) TreeNet for discerning insect-infested forests induced by defoliators from healthy forests using Landsat 8 OLI and ancillary data in the broadleaved mixed Hyrcanian forests. Moreover, it has studied mutual associations between the intensity of forest defoliation and the severity of forest fires under TerraClimate-derived climate hazards by analyzing panel data models within the TreeNet-derived insect-infested forest objects. The TreeNet optimal performance was obtained after building 333 trees with a sensitivity of 93.7% for detecting insect-infested objects with the contribution of the top 22 influential variables from 95 input object features. Accordingly, top image-derived features were the mean of the second principal component (PC2), the mean of the red channel derived from the gray-level co-occurrence matrix (GLCM), and the mean values of the normalized difference water index (NDWI) and the global environment monitoring index (GEMI). However, tree species type has been considered as the second rank for discriminating forest-infested objects from non-forest-infested objects. The panel data models using random effects indicated that the intensity of maximum temperatures of the current and previous years, the drought and soil-moisture deficiency of the current year, and the severity of forest fires of the previous year could significantly trigger the insect outbreaks. However, maximum temperatures were the only significant triggers of forest fires. This research proposes testing the combination of object features of Landsat 8 OLI with other data for monitoring near-real-time defoliation and pathogens in forests.

1. Introduction

Despite prosperous traditional approaches such as dendrological assessment and field observations for identifying driving forces of insect outbreaks from individual tree to stand scales [1], remotely sensed approaches are extensively progressing either for delineation insect-infested objects or for the mensuration of infestations induced by abiotic and biotic agents throughout forest biomes [2,3,4,5]. However, some novel algorithms for data mining and machine learning such as TreeNet [6] for delineation insect-infested objects from non-insect-infested objects of images, some high-resolution climate data such as TerraClimate [7] for assessing drought and climate hazards dimensions, and some associations such as interactions between insect outbreaks, forest fires, and climate hazards have received less attention in earlier studies.
Although monitoring the bark beetle infestation and coniferous defoliation are dependent on high-resolution and multi-spectral images [8,9], detecting broadleaved defoliation has been predestined by the spectral–temporal information of images, even by single near-infrared-derived vegetation indices of images with high-temporal resolutions [10]. Moreover, Landsat images have indicated high accuracy for detecting forest-infested patches using either classification algorithms in a specific date [11] or by applying multitemporal spectral-derived indices [12,13,14]. Moderate stages of tree infestation are significantly discernible through multitemporal spectral indices, while severe infestation is highly discernible through classification approaches [15].
Several number of studies exerted data mining and machine learning algorithms such as random forest [4,16,17] and decision tree [18] or maximum likelihood classification [19,20] to discern insect-defoliated from non-insect-defoliated forests. However, numerous remote sensing-derived indices are growing for quantifying the insect-induced defoliation intensity in terms of the long-term archive of Landsat products [13,14,18]. For example, Townsend et al. [13] claimed that Landsat-derived near-infrared (NIR) and short-wave infrared (SWIR) indices such as the normalized different infrared index (NDII) and the moisture stress index (MSI) were superior to the Landsat-derived red and NIR indices such as the normalized difference vegetation index (NDVI) for estimating the defoliation induced by Lymantria dispar in the oak forests at five different times. Similarly, Rullán-Silva et al. [14] examined the efficiencies of 10 Landsat-derived vegetation indices for estimating the defoliation induced by Rhynchaenus fagi in the beech forests. They concluded that the MSI showed the highest correlation with the intensity of defoliation with respect to the mixed-effects model. However, the attitude of analyzing meaningful image-objects for the classification of infested and non-infested forests [21] by the cooperation of image-derived features and GIS-based methods and databases is developing, that is, Geographic Object-Based Image Analysis (GEOBIA) [22]. In addition, earlier results of studies have demonstrated that the combination of object features derived from medium-resolution images such as Landsat and SPOT with ancillary data such as topography has increased the precision of discerning forest infestations induced by bark beetle [8]. Machine learning algorithms such as random forest improved the classification efficiency with image-derived object features for mapping forest infestations [23]. However, evaluating the efficiency of combination GEOBIA and other machine learning methods such as TreeNet for mapping infested forests induced by defoliators has not been addressed up to now.
Forests are mainly degraded by interactions between abiotic agents such as high temperatures, drought, moisture variability, anthropogenic interventions, and forest fires and biotic agents such as invasive species, tree mortality, insect infestations, and pathogens [24,25,26,27]. The abiotic agents can trigger the effects of biotic agents, and their synchronization significantly results in tree species’ hydraulic deterioration and carbon starvation [11,24,25,27]. Climate change studies demonstrated that defoliators’ population rate during a growing season is correlated to the temperature condition of its hydrological year [28]. Temperature increases can enable insects to survive during the cold season [29,30,31] and provide rich sources of nutrition from the mature leaves by changing the trees phenology cycle; along with that, they can also increase the risk of insect outbreaks during the growing season [32]. However, there is existing uncertainty about the effects of droughts on the insect outbreaks in the forest biomes [33]. Increasing the droughts’ dimensions may result in providing conditions for insect outbreaks [11], tree mortality [34], or increasing the forest fires severity [35]. However, serious effects of insect defoliators emerge during the moderate drought condition or wetter seasons following droughts occurring [11,24,33,36]. Changes in moisture capacity regarding either high moisture availability [1,37] or low moisture availability [38] are identified as the main driver forcing of some insect outbreaks. The synchronization of droughts following the conditions of above-average moisture availability may result in providing an appropriate condition for insect outbreaks [24] as well.
Mutual interactions between insect infestations and forest fires were documented in some studies [39]. Insect infestations may affect the fuel loads of the tree species and increase the severity of forest fires at landscape-level scales [34,39,40]. However, some studies reported declining the severity of forest fires by increasing the mortality induced by insect attacks [41] or neutral effects of insects and pathogens on the fire characteristics [42,43] particularly in the coniferous forests. Additionally, there is evidence of the probability of increasing the risk of insect outbreaks in those trees that were damaged or weakened following a low severity of forest fires [44] or in fire-induced larger patches of canopy cover [45].
The Caspian Hyrcanian broadleaves and mixed forests have been degrading, as a virgin ecoregion of the temperate forests’ biome, by a variety of biotic and abiotic agents such as deforestation [46], forest fires [47], drought [48], and climate hazards with consequences of massive tree mortality [49] in northeast (NE) Iran during recent decades. Droughts’ dimensions could significantly affect the water content and greenness properties of Hyrcanian forests based on the MODIS-derived normalized difference water index (NDWI) and NDVI [48,49]. The stages of moderate to extreme tree mortality events showed a significant association with the high intensity of forest water content deficit derived from the MODIS–NDWI; however, the severe defoliation only showed a significant relationship with the intensity of forest greenness loss derived from NDVI in NE Iran [49]. In addition to climate hazards, drought, and forest fires, there has been rising concerning reports about the outbreaks of some insect defoliators such as Lymantria dispar, Erannis defoliaria, and Operophtera brumata [50,51] as well as pathogens [52,53,54] throughout Hyrcanian forests during recent years. Therefore, this study used TreeNet to delineate insect-infested forests from non-insect-infested forests using numerous Landsat 8 OLI-derived object features, topographic-derived features, and tree species types in Hyrcanian forest, NE Iran. Moreover, it will explore the mutual relationships between the intensity of insect infestation and the severity of forest fires in the presence of TerraClimate-derived climate hazard variables for the period of time of insect outbreaks and forest fires within the TreeNet-derived insect-infested forest objects.

2. Materials and Methods

2.1. Study Area

The eastern forests of the Hyrcanian ecoregion were selected for this research. This area is extended from Gorgan to Galikesh in the Golestan province, NE Iran (Figure 1). These forests comprise a mixture of broadleaved tree species such as Quercus castaneifolia, Fagus orientalis, Carpinus betulus, Acer spp.,Tilia platyphyllos, and Parrotia persica [55]. The western parts were infested by the defoliators of Erannis defoliaria and Operophtera brumata, while the eastern parts were affected by Lymantria dispar (Figure 1). Moreover, the frequency of forest fire events has been increasing in this region during recent years (Figure 1).

2.2. Data and Field Mensuration

This study identified the insect-infested regions from available reports, local media, and field observations. The attributes of defoliated spots including the type of defoliators, the position, the dominant type of host tree-species, and the evidence of current and previous egg laying were documented. The accurate spatial extent areas of the defoliation were delineated using time-series composite bands of Landsat 8 OLI [56] and the images of Google Earth (Figure 2) coincided with the advanced larval stage of insect defoliators (Figure 1). A peak of attack was recognized in 2014; the infested objects for this time were delineated from the healthy forest using GEOBIA and TreeNet. The forest fires data were achieved from field surveying, local media, and the available historical database that were provided by the Department of Natural Resources and Watershed Management of Golestan province from 2010 to 2017. The burnt area, duration, and frequency of forest fires were used during a fire season to reach the forest fire severity of a specific location.
We used monthly TerraClimate data for calculating [7] long-term dimensions of anomalies of drought, temperature, and soil moisture from 1987 to 2017. We derived the annual intensity hazard of climate variables from their dimensions as causes of triggering insect attack and forest fires for modeling by the panel data approaches (see Section 2.3.3). The ancillary data such as the Topographic Position Index (TPI), Terrain Ruggedness Index (TRI), and Topographic Wetness Index (TWI) were derived from the ALOS PALSAR data elevation [57] for the study area. Forest types were vectorized concerning the scanned maps of the forest management plans [58]. These data were used for mapping insect-infested forests along with Landsat 8 OLI data.

2.3. Methodology

2.3.1. TreeNet-Based Insect Infestation Mapping

GEOBIA was used to delineate the insect-defoliated areas from the healthy forests through image segmentation and TreeNet classification. We derived image objects from a set combination of the main spectral and panchromatic channels of Landsat 8 OLI for the peak times of defoliation (May–June) through the multiresolution segmentation algorithm. To minimize the mean heterogeneity of image objects, we assigned optimal scale parameters by trial and error, with higher weights for the red, NIR, and SWIR bands; the compactness value of 1; and the shape value of 0 in the eCognition Developer 9 [60].
• Object features
Various object features (95 features) were derived from the main channels of Landsat 8 OLI, topography data, and forest types, as shown in Table 1. The object features were classified into four main groups including spectral features, surface texture features, geometric features, and the features derived from ancillary data in the GIS. A single database was created including all the derived features of segment objects, and was utilized for classification using TreeNet.
• TreeNet classification
The stratified random sampling method was used for selecting samples of insect-infested objects (defoliation >50%) and non-insect-infested objects for assessing the object features that control insect outbreaks and discriminate them from the non-insect-affected forests. The TreeNet algorithm was applied for determining influential variables depending on the test sets and generalizing the obtained scores to all feature objects to distinguish between the insect-infested and non-insect-infested objects. Classification in TreeNet is a particular form of regression that produces a possibility of response for a variable and accurately ranks the predictor variables based on their importance from “most likely” to “least likely” to the target variable. TreeNet creates boosting regression models through sequentially fitting a very small tree in several stages. Accordingly, each stage learns from the available training data by a specific learning rate [6].
We applied TreeNet in Salford Predictive Modeler 8.3 [75] to build our model with respect to its classification precision, working with both parametric and non-parametric variables, handling big data and missing datasets, connecting to GIS, plotting the univariate and bivariate relationships between the response and predictor variables through partial dependency plots (PDPs), yielding reliable results despite existing non-stationary in data, and ranking the predictor variables in terms of their importance in the model [6,76,77]. One-third of the samples were randomly assigned for the testing set, and the remaining were considered as the learning samples. We set the number of trees and maximum nodes per tree to 500 and 6, respectively. The shrinkage method was chosen as “auto” to eliminate the complexity of overfitting and to set the learning rate of the model, which was calculated to be 0.01. Moreover, we evaluated the number of optimal trees using the criterion of the area under the receiver operating characteristic (ROC) curve and the misclassification rate for the test samples. The confusion matrix was used to assess the performance of the classification model with respect to the test samples using four measures: (1) Sensitivity: the proportion of the insect-infested objects that are correctly classified, (2) Specificity: the proportion of the non-insect-infested objects that are correctly classified, (3) Precision: the proportion of the actual classified insect-infested objects divided by total number of the insect-infested class testing objects, and (4) an F1 statistic derived from the sensitivity and precision metrics as shown in Equation (1), in which its values close to one indicate a high sensitivity of the model:
F 1   s t a t i s t i c = 2 ( S e n s i t i v i t y × P r e c i s i o n ) S e n s i t i v i t y + P r e c i s i o n

2.3.2. Intensity of Insect Infestation, Severity of Forest Fire, and Climate Hazards

The intensity of defoliation was retrieved by analyzing the long-term deficit of NDWI-derived forest-water content within the TreeNet-derived insect-infested polygons. We calculated the dimensions of defoliation including the severity, frequency, and duration of the yearly defoliations from the anomalies of the Landsat NDWI from 1987 to 2017. The yearly hazard intensity of defoliation was mapped concerning the introduced approach by Abdi et al. [49], which is a combination of standardized values of defoliation dimensions through the fuzzy gamma operator [78].
We obtained the spatial variations and characteristics of daily forest fires from diverse resources from 2010 to 2017. Several studies verified the strength of Landsat-derived burn severity indices such as the differenced Normalized Burn Ratio (dNBR) [79] and the Relativized dNBR (RdNBR) [80] for large fires; however, some fires’ characteristics such as size area, duration, and the type of fire may restrict these indices applications for post-fire monitoring in northeast Iran. Therefore, the severity of forest fires was calculated from the combination of the ground-based characteristics of fires including the frequency, size, and duration of fires within segment objects. The fire characteristics’ values were standardized between zero (low) and one (high) by exerting the membership function of fuzzy linear [81]. The severity of forest fires was obtained by overlaying the standardized layers of fire characteristics using the fuzzy gamma operator [78] during a fire season.
TerraClimate data were applied to model the hazard intensity of climate variables throughout the hydrological years (October to September) from 1987 to 2017. The features (i.e., severity, frequency, and duration) of yearly anomalies of the climate variables were combined to obtain the hazard intensity indices of drought, maximum temperatures, and soil moisture deficit using the fuzzy gamma operator [78] as well.

2.3.3. Relationships among Insect Infestation, Forest Fires, and Climate Hazards

We examined mutual relationships between insect infestation and forest fires in the presence of the hazard intensity of climate variables within the insect-affected objects using the panel data models. To include both the spatial and time effects of the data, we performed the panel data models [82] for the estimation of the intensity of defoliation affected by forest fire severity and climate hazards (Equation (2)), as well as the severity of forest fires induced by insect defoliation and the climate hazards (Equation (3)):
I n s e c t   i n f e s t a t i o n = f n ( F t − 1 ,   D t ,   D t − 1 ,   T t ,   T t − 1 ,   S t ,   S t − 1 )
F i r e   s e v e r i t y = f n ( I t ,   I t − 1 ,   D t ,   D t − 1 ,   T t ,   T t − 1 ,   S t ,   S t − 1 )
where It, Ft, Dt, Tt, and St are the averages of insect infestation intensity, fire severity, drought intensity, high temperature, and soil moisture deficits in the current year (t); and It − 1, Ft − 1, Dt − 1, Tt − 1, and St − 1 are the averages of these variables for the previous year (t − 1), respectively.

Panel data models

Analyzing the panel data was established based on the combination of multiple observations for the same cross-sections and time series, which is written in a standard model as shown in Equation (4) [83]:
y i , t = β ′ X i t + Z i α + ε i t
For i = 1 ,   2 , 3 , … , N and t = 1 ,   2 , 3 , … , T ; where N stands for the number of individuals (cross-sections), and T is the number of times. The vector X i t refers to the k regressors (Equations (2) and (3)). The vector β ′ refers to unobserved coefficients, which need to be estimated. The term Z i α expresses the cross-section specific effects. The error terms of the model were indicated by ε i t . Three panel data models are defined regarding the nature of included variables in the vector Z i , including the common effects, fixed effects, and also random effects models (Equations (5)–(7)).
The common effects model does not consider time and cross-sectional effects. The vector Z i contains only one constant coefficient: α . The coefficients of this model are estimated by applying the ordinary least squares (OLS) approach:
y i , t = β ′ X i t + α + ε i t
The fixed effects model assumes specific cross-sectional effects from different intercepts to obtain unobserved heterogeneity. The parameter of fixed effect α i is constantly estimated over time with the estimators’ normality assumption. This model assumes that cross-sectional effects are correlated with the included regressors X i t :
y i , t = β ′ X i t + α i + ε i t
The random effects model assumes that unobserved cross-sectional heterogeneity is not correlated to the included regressors X i t . The coefficients of this model are estimated using the generalized least squares (GLS) estimator:
y i , t = β ′ X i t + α + μ i t
μ i t = α i + ε i t
The random effects model is preferable to the fixed effect model if samples were taken from a larger population and the estimations’ generalization to the population was aimed. The higher accuracy of estimates and greater flexibility of this model in comparison with the fixed effects model were reported [84,85] as well.

Testing for fixed effects and random effects

Tests were accomplished to make a distinction between panel data models. The Chow test was used for testing fixed effects against common effects [83], which determines whether a single regression model (common effects) is superior to the two separate regression models (fixed effects). The Hausman’s specification test [86] was used for selecting whether the fixed effects model or the random effects model is appropriate [87]. The Hausman test makes a distinction between a model that assumes its cross-sectional effects are not correlated with its regressors (random effects), and a model that assumes that these relationships are established (fixed effects) [83].

3. Results

3.1. Insect Defoliation Mapping

The highest performance of TreeNet was obtained after the formation of 333 trees with the optimal ROC value of 93.4% for discriminating insect-infested and non-insect-infested forest objects. The model results using 83 out of 95 predictor variables yielded an average correctness of 87% for testing samples to predict the insect-infested and non-insect-infested forests (Table 2).
Relative variable importance values describe that the top predictors are the mean of PC2, tree species, and the mean of the red channel derived from the gray-level co-occurrence matrix (GLCM), NDWI, and global environment monitoring index (GEMI) (Figure 3). The positive log-odds values of the insect-infested class were significantly increased when the mean of the PC2 values had passed the point of 0.70. Tree species with the domination of Carpinus betulus-Quercus castaneifolia and Carpinus betulus-Acer spp.-Tilia platyphyllos exhibited the highest partial relationships with the insects’ presence. The average log-odds values of the insect-infested class were increased from –0.287 to 0.314, where the mean of the red band derived from GLCM ranged from 122 to 130. The average log-odds dropped at values of 0.264 and 0.635 of the NDWI and GEMI, respectively, where the probability of infestation also reduced (Figure 4).
The TreeNet model was rebuilt with respect to the top influential object features that have gained a minimum score of 15% importance relative to the most important variable. The insect-infested and non-insect-infested map was created with the contribution of the top 22 predictor variables over the study area (Figure 5).
The hazard intensity of defoliation derived from the time-series anomaly of Landsat–NDWI was mapped as an insect infestation indication within the delineated insect-infested forest objects from 2010 to 2017 (Figure 6). The maps demonstrated that the infestation intensity was considerable in 2011, 2014, and 2015 (Figure 6b,e,f).

3.2. Insect Infestation, Forest Fires, and Climate Hazards Modelling

The results of the Hausman test indicate that the random effects model was superior to the fixed effects model for the variables expressing that control either the intensity of insect infestation (X2 = 11.87, df = 7, p= 0.105) or severity of forest fires (X2 = 6.72, df = 8, p = 0.567). Furthermore, the results of the Chow test showed that the fixed effects model was superior to the common effects in both models (X2 = 37.39, df = 18, p = 0.004; X2 = 43.31, df = 18, p = 0.000).
The summary of the random effects model reveals that the drought condition (SPI), maximum temperature (Tmax), and the deficit of soil moisture of the current year—along with the maximum temperature, high soil moisture availability, and the severity of forest fires of the previous year—were the significant variables, which are controlling the insect infestation intensity during 2010–2017 in the study area (Table 3). However, the intensity of drought of the previous year was not significant, but its coefficient is positive. Although the majority of climate hazards demonstrated positive coefficient with the severity of forest fires, Tmax of the current year (β = 0.330, p < 0.05) and Tmax of the previous year (β = 0.196, p < 0.01) were the only significant variables. The insect infestation (II) of the previous year showed a positive relationship with the severity of forest fire as well; however, it is not statistically significant (β = 0.106, p > 0.05) in NE Iran (Table 4).

4. Discussion

This study exerted novel remote sensing-based data collections (e.g., Landsat and TerraClimate) and approaches (e.g., TreeNet and panel data models) in order to discern insect-infested forests and quantify the multitemporal intensity of forest defoliation, climate hazards, drought, and the severity of forest fires along with modeling their temporal and spatial interactions in NE Iran.
We obtained reliable discrimination between insect-infested forests and non-insect-infested forests (Table 2) using the important collections of Landsat 8 OLI-derived and ancillary object features using the TreeNet algorithm. The summary of variable importance represented that the object features extracted from the Landsat channels show higher performance in comparison with the main spectral channels of Landsat (Figure 3). For example, the highest value recorded for the mean of PC2, and after that the mean of red channel derived from GLCM and the mean values of vegetation indices, including the NDWI and GEMI, respectively. Analyzing one partial dependency revealed that the probability of the presence of defoliation increases along with the increase of average values of the PC2 and the mean of the red channel derived from GLCM, while it decreases by the increase of values of the NDWI and GEMI (Figure 3). The PC2 was positively loaded on the visible bands of Landsat 8 OLI with a higher coefficient for the red and green bands, which high-defoliated forest objects have obtained high values on PC2 (Figure 4b). The increase in the reflectance of some ranges of visible wavelengths due to the vegetation stress was reported in earlier studies [88,89], which are near to the green and red bands of Landsat 8 OLI. In contrast, PC2 was loaded negatively on the NIR and SWIR bands, as the objects with high biomass showed lower values on the PC2. While our analysis indicated that the PC2–Landsat 8 OLI is the top variable predictor, the third principal component (PC3) was positively loaded on the red band, and has showed considerable importance value in discriminating between insect-infested forest objects and healthy forest objects (Figure 3) as well. In addition, the ability of PC3 was reported for identifying year-to-year forest defoliation by Lymantria dispar L. using the time series of SPOT [90]. Moreover, the importance of Landsat-derived indices was presented for detecting insect-affected forests for different types of insects in other biomes [13,14,18]. This study results confirm that vegetation indices with respect to the NIR–SWIR (i.e., NDWI) are superior to those vegetation indices that are depending on the visible NIR channels for the delineation of insect-infested forest objects [13,14] in the Hyrcanian ecoregion as well.
Insect defoliators can affect the structure of vegetation [91]; therefore, image-derived textural attributes are appropriate for detecting insect-defoliating forest objects. However, the mean of the red channel derived from GLCM was among the top five predictors; Figure 3 indicates that about one-third of effective predictor variables are categorized in the first-order or second-order textural attributes for detecting the defoliation induced by insects in NE Iran. The superiority of image textures derived from GLCM and gray-level occurrence matrices (GLOM) [91,92] for detecting vegetation degradation induced by insect defoliation has been particularly demonstrated for high-spatial resolution images.
The results of variable importance describe that the type of tree species has a significant effect on the insect outbreaks as well. However, the partial dependence plot shows that mixed forest stands such as Carpinus betulus and Quercus castaneifolia with the other tree species types were highly correlated with the positive values of probability infestation in the Hyrcanian forests (Figure 4a). However, some studies reported insect defoliation in the monospecific stands in the other forest biomes [9,14,93].
The results of interactions between the intensity of insect defoliation with the intensity of climate hazards demonstrated that Tmax, drought, and soil moisture significantly increased the intensity of insect infestation. Both current and previous temperatures presented significant coefficients (Table 3). The high temperatures of the previous year can enable insects to survive during the winter [29,30,31] and change the cycle of forest phonology with the insect outbreaks consequences during the growing season [32].
Although some studies reported uncertainty about the effect of drought on the insect eruption [33], this study’s analyses confirmed that drought condition is a key trigger of increasing insect outbreaks in Hyrcanian broadleaved forests (Table 3). This area has experienced severe droughts from 2010 to 2011 [48]; however, the wetter seasons during 2012 and 2013 provided the condition for eruption insects during moderate drought seasons. Also, these fluctuations in the wet and dry seasons have been reported as a trigger of insect outbreaks [24]. The panel data models using random effects verified that the soil moisture availability [1,37] in the previous year significantly increased the intensity of insect infestation (Table 3). However, the effect of the deficit of soil moisture in the current year was not significant; its positive coefficient indicates that it was more likely that it intensified the insect eruption [38]. Therefore, the conditions of the previous year regarding its high temperatures and availability of moisture have supplied the sources of nutrients for insects and caused massive outbreaks in the current year with existing high temperatures and drought conditions. Earlier studies emphasized that high temperature is the main cause of increasing the intensity of tree mortality in the Hyrcanian forests [49].
Moreover, the random effects model indicated that the severity of forest fires of the previous year was a significant driver in increasing the intensity of insect infestation (Table 3). The trends of forest fires were dramatically increased following 2010 in the study area (Figure 1). Since the majority of the type of fires are classified as surface fires and single-tree burning in the study area, the likelihood of damaging trees [44] and opening fire-induced spaces [45] have increased the risk of insect infestation in these spots. The presence insect defoliators was higher in the locations with a high density of fires (Figure 1). In contrast, the intensity of insect defoliation either in the current year or in previous year was not a significant trigger for increasing the severity of forest fires in the Hyrcanian ecoregion (Table 4). The neutral effects of insect infestation on fire severity were demonstrated in the studies accomplished by Meigs et al. [42] and Kane et al. [43] as well. However, some researchers concluded that insect-induced tree mortality has decreased fuels and associated fire proneness, resulting in declining fire severity [41]. We can justify these results for two reasons in Hyrcanian forests: (1) insect infestations have not yet led to such a massive tree mortality that could load extensive fuels for fires [34,39,40] in this area, and (2) human activities are identified as the main causes of fires in this area, and mostly occurred in the condition of high temperatures [47]. This study confirmed that the maximum temperature of the current year and previous year are the only significant variables controlling forest fires in the Hyrcanian ecoregion (Table 4).
Landsat 8 OLI-derived features have shown reasonable efficiency for mapping insect-defoliated Hyrcanian mixed broadleaved forests for the investigated periods. Alternatively, this study proposes testing the capabilities of new data such as Sentinel-2 with higher spectral, spatial, and swat width that can produce dense time series, which are more appropriate for monitoring near-real-time insect infestation. There have been increasing reports about the outbreaks of other biotic agents such as bark beetles and pathogens in the Hyrcanian forests [52,53,54]; therefore, further investigations are required for identifying the spatial extent of these agents and exploring their triggers. Furthermore, separating mortality induced by insect defoliators from bark beetles or pathogens should be scrutinized.

5. Conclusions

This study has designed an object-based TreeNet framework to discern insect-defoliated forests using Landsat 8 OLI-derived and ancillary object features. Furthermore, we performed panel data models for quantifying the effects of TerraClimate-derived climate hazards on insect outbreaks and forest fires along with mutual associations of insect infestation and forest fires in the Hyrcanian broadleaved mixed forests, NE Iran. According to the analyses, we drew the following conclusions:
  • GEOBIA TreeNet indicated excellent performance with the contribution of Landsat 8 OLI-derived and ancillary object features for discriminating insect-defoliated forests from healthy forests.
  • Although the object features of Landsat 8 OLI recorded a higher importance for discriminating insect-defoliated objects, tree species has obtained the second rank of importance following the mean of PC2. In addition, other top image object features were the mean of the red channels derived from GLCM, the mean of NDWI, and the mean of GEMI, respectively.
  • The random effects model demonstrated higher performance in comparison with the fixed effects and common effects models to model mutual interaction of the intensity of insect defoliation and the severity of forest fire and their associations with the TerraClimate-derived climate hazards.
  • Maximum temperatures significantly triggered both insect outbreaks and forest fires. Although the drought conditions of the current year and the availability of soil moisture of the previous year were significant regarding the intensity of insect infestation, they have indicated neutral effects on the severity of forest fires.
  • The severity of forest fires of the previous year has triggered the intensity of insect infestation; however, the insect infestation was not effective for the forest fires.
  • Future studies will be required to explore the application of novel satellite images such as Sentinel-2 or the combination of Landsat 8 and Sentinel-2 for monitoring near-real-time insect-induced defoliation, identifying infestations resulting from bark beetles and pathogens, and discriminating between them.

Author Contributions

O.A. implemented all parts of this paper, including design of the work, the mapping of insect infestation using eCognition® and Salford Predictive Modeler®, analyzing panel data in R, and writing the manuscript.

Funding

This research received no external funding.

Acknowledgments

The author acknowledges support by the German Research Foundation and the Open Access Funding by the Publication Fund of the TU Dresden. The author expresses his acknowledgment to the Department of Natural Resources and Watershed Management of Golestan for providing some forest information and fire attributes during the research.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Swetnam, T.W.; Lynch, A.M. Multicentury, Regional-Scale Patterns of Western Spruce Budworm Outbreaks. Ecol. Monogr. 1993, 63, 399–424. [Google Scholar] [CrossRef]
  2. Rullan-Silva, C.D.; Olthoff, A.E.; La Delgado de Mata, J.A.; Pajares-Alonso, J.A. Remote Monitoring of Forest Insect Defoliation. A Review. For. Syst. 2013, 22, 377–391. [Google Scholar] [CrossRef]
  3. Hall, R.J.; Castilla, G.; White, J.C.; Cooke, B.J.; Skakun, R.S. Remote sensing of forest pest damage: A review and lessons learned from a Canadian perspective. Can. Entomol. 2016, 148, S296–S356. [Google Scholar] [CrossRef]
  4. Dash, J.P.; Watt, M.S.; Pearse, G.D.; Heaphy, M.; Dungey, H.S. Assessing very high resolution UAV imagery for monitoring forest health during a simulated disease outbreak. ISPRS J. Photogramm. Remote Sens. 2017, 131, 1–14. [Google Scholar] [CrossRef]
  5. Senf, C.; Seidl, R.; Hostert, P. Remote sensing of forest insect disturbances: Current state and future directions. Int. J. Appl. Earth Obs. Geoinf. 2017, 60, 49–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Friedman, J.H. Stochastic gradient boosting. Comput. Stat. Data Anal. 2002, 38, 367–378. [Google Scholar] [CrossRef]
  7. Abatzoglou, J.T.; Dobrowski, S.Z.; Parks, S.A.; Hegewisch, K.C. TerraClimate, a high-resolution global dataset of monthly climate and climatic water balance from 1958–2015. Sci. Data 2018, 5, 170191. [Google Scholar] [CrossRef] [PubMed]
  8. Latifi, H.; Fassnacht, F.E.; Schumann, B.; Dech, S. Object-based extraction of bark beetle (Ips typographus L.) infestations using multi-date LANDSAT and SPOT satellite imagery. Prog. Phys. Geogr. Earth Environ. 2014, 38, 755–785. [Google Scholar] [CrossRef]
  9. Meddens, A.J.H.; Hicke, J.A. Spatial and temporal patterns of Landsat-based detection of tree mortality caused by a mountain pine beetle outbreak in Colorado, USA. For. Ecol. Manag. 2014, 322, 78–88. [Google Scholar] [CrossRef]
  10. De Beurs, K.M.; Townsend, P.A. Estimating the effect of gypsy moth defoliation using MODIS. Remote Sens. Environ. 2008, 112, 3983–3990. [Google Scholar] [CrossRef]
  11. Senf, C.; Wulder, M.A.; Campbell, E.M.; Hostert, P. Using Landsat to Assess the Relationship between Spatiotemporal Patterns of Western Spruce Budworm Outbreaks and Regional-Scale Weather Variability. Can. J. Remote Sens. 2016, 42, 706–718. [Google Scholar] [CrossRef]
  12. Babst, F.; Esper, J.; Parlow, E. Landsat TM/ETM+ and tree-ring based assessment of spatiotemporal patterns of the autumnal moth (Epirrita autumnata) in northernmost Fennoscandia. Remote Sens. Environ. 2010, 114, 637–646. [Google Scholar] [CrossRef]
  13. Townsend, P.A.; Singh, A.; Foster, J.R.; Rehberg, N.J.; Kingdon, C.C.; Eshleman, K.N.; Seagle, S.W. A general Landsat model to predict canopy defoliation in broadleaf deciduous forests. Remote Sens. Environ. 2012, 119, 255–265. [Google Scholar] [CrossRef]
  14. Rullán-Silva, C.; Olthoff, A.E.; Pando, V.; Pajares, J.A.; Delgado, J.A. Remote monitoring of defoliation by the beech leaf-mining weevil Rhynchaenus fagi in northern Spain. For. Ecol. Manag. 2015, 347, 200–208. [Google Scholar] [CrossRef]
  15. Meddens, A.J.H.; Hicke, J.A.; Vierling, L.A.; Hudak, A.T. Evaluating methods to detect bark beetle-caused tree mortality using single-date and multi-date Landsat imagery. Remote Sens. Environ. 2013, 132, 49–58. [Google Scholar] [CrossRef]
  16. Vastaranta, M.; Kantola, T.; Lyytikäinen-Saarenmaa, P.; Holopainen, M.; Kankare, V.; Wulder, M.; Hyyppä, J.; Hyyppä, H. Area-Based Mapping of Defoliation of Scots Pine Stands Using Airborne Scanning LiDAR. Remote Sens. 2013, 5, 1220–1234. [Google Scholar] [CrossRef] [Green Version]
  17. Adelabu, S.; Mutanga, O.; Adam, E. Evaluating the impact of red-edge band from Rapideye image for classifying insect defoliation levels. ISPRS J. Photogramm. Remote Sens. 2014, 95, 34–41. [Google Scholar] [CrossRef]
  18. Goodwin, N.R.; Coops, N.C.; Wulder, M.A.; Gillanders, S.; Schroeder, T.A.; Nelson, T. Estimation of insect infestation dynamics using a temporal sequence of Landsat data. Remote Sens. Environ. 2008, 112, 3680–3689. [Google Scholar] [CrossRef]
  19. Paritsis, J.; Veblen, T.T.; Smith, J.M.; Holz, A. Spatial prediction of caterpillar (Ormiscodes) defoliation in Patagonian Nothofagus forests. Landsc. Ecol. 2011, 26, 791–803. [Google Scholar] [CrossRef]
  20. Thayn, J.B. Using a remotely sensed optimized Disturbance Index to detect insect defoliation in the Apostle Islands, Wisconsin, USA. Remote Sens. Environ. 2013, 136, 210–217. [Google Scholar] [CrossRef]
  21. Heurich, M.; Ochs, T.; Andresen, T.; Schneider, T. Object-orientated image analysis for the semi-automatic detection of dead trees following a spruce bark beetle (Ips typographus) outbreak. Eur. J. For. Res. 2010, 129, 313–324. [Google Scholar] [CrossRef]
  22. Hay, G.J.; Castilla, G. Geographic Object-Based Image Analysis (GEOBIA): A new name for a new discipline. In Object-Based Image Analysis: Spatial Concepts for Knowledge-Driven Remote Sensing Applications; Blaschke, T., Lang, S., Hay, G.J., Eds.; Springer: Berlin, Germany; London, UK, 2008; pp. 75–89. ISBN 978-3-540-77057-2. [Google Scholar]
  23. Latifi, H.; Schumann, B.; Kautz, M.; Dech, S. Spatial characterization of bark beetle infestations by a multidate synergy of SPOT and Landsat imagery. Environ. Monit. Assess. 2014, 186, 441–456. [Google Scholar] [CrossRef] [PubMed]
  24. Flower, A.; Gavin, D.G.; Heyerdahl, E.K.; Parsons, R.A.; Cohn, G.M. Drought-triggered western spruce budworm outbreaks in the interior Pacific Northwest: A multi-century dendrochronological record. For. Ecol. Manag. 2014, 324, 16–27. [Google Scholar] [CrossRef]
  25. Millar, C.I.; Stephenson, N.L. Temperate forest health in an era of emerging megadisturbance. Science 2015, 349, 823–826. [Google Scholar] [CrossRef] [PubMed]
  26. Kautz, M.; Meddens, A.J.H.; Hall, R.J.; Arneth, A. Biotic disturbances in Northern Hemisphere forests—A synthesis of recent data, uncertainties and implications for forest monitoring and modelling. Glob. Ecol. Biogeogr. 2017, 26, 533–552. [Google Scholar] [CrossRef]
  27. Ward, S.F.; Aukema, B.H. Anomalous outbreaks of an invasive defoliator and native bark beetle facilitated by warm temperatures, changes in precipitation and interspecific interactions. Ecography 2019, 42, 1068–1078. [Google Scholar] [CrossRef] [Green Version]
  28. Ryan, M.G.; Vose, J.M. Effects of Climatic Variability and Change on Forest Ecosystems: A Comprehensive Science Synthesis for the U.S. Forest Sector. PNW-GTR-870; 2012. Available online: https://www.fs.usda.gov/treesearch/pubs/42651 (accessed on 2 August 2019).
  29. Paradis, A.; Elkinton, J.; Hayhoe, K.; Buonaccorsi, J. Role of winter temperature and climate change on the survival and future range expansion of the hemlock woolly adelgid (Adelges tsugae) in eastern North America. Mitig. Adapt. Strateg. Glob. Chang. 2008, 13, 541–554. [Google Scholar] [CrossRef]
  30. Trần, J.K.; Ylioja, T.; Billings, R.F.; Régnière, J.; Ayres, M.P. Impact of minimum winter temperatures on the population dynamics of Dendroctonus frontalis. Ecol. Appl. 2007, 17, 882–899. [Google Scholar] [CrossRef]
  31. Pureswaran, D.S.; Roques, A.; Battisti, A. Forest Insects and Climate Change. Curr. For. Rep. 2018, 4, 35–50. [Google Scholar] [CrossRef] [Green Version]
  32. Jepsen, J.U.; Kapari, L.; Hagen, S.B.; Schott, T.; Vindstad, O.P.L.; Nilssen, A.C.; Ims, R.A. Rapid northwards expansion of a forest insect pest attributed to spring phenology matching with sub-Arctic birch. Glob. Chang. Biol. 2011, 17, 2071–2083. [Google Scholar] [CrossRef]
  33. Kolb, T.E.; Fettig, C.J.; Ayres, M.P.; Bentz, B.J.; Hicke, J.A.; Mathiasen, R.; Stewart, J.E.; Weed, A.S. Observed and anticipated impacts of drought on forest insects and diseases in the United States. For. Ecol. Manag. 2016, 380, 321–334. [Google Scholar] [CrossRef]
  34. Fettig, C.J.; Mortenson, L.A.; Bulaon, B.M.; Foulk, P.B. Tree mortality following drought in the central and southern Sierra Nevada, California, U.S. For. Ecol. Manag. 2019, 432, 164–178. [Google Scholar] [CrossRef]
  35. Littell, J.S.; Peterson, D.L.; Riley, K.L.; Liu, Y.; Luce, C.H. Fire and Drought. In Effects of Drought on Forests and Rangelands in the United States: A Comprehensive Science Synthesis; U.S. Department of Agriculture, Forest Service: Washington, DC, USA, 2016; pp. 135–154. [Google Scholar]
  36. Rouault, G.; Candau, J.-N.; Lieutier, F.; Nageleisen, L.-M.; Martin, J.-C.; Warzée, N. Effects of drought and heat on forest insect populations in relation to the 2003 drought in Western Europe. Ann. For. Sci. 2006, 63, 613–624. [Google Scholar] [CrossRef]
  37. Ryerson, D.E.; Swetnam, T.W.; Lynch, A.M. A tree-ring reconstruction of western spruce budworm outbreaks in the San Juan Mountains, Colorado, U.S.A. Can. J. For. Res. 2003, 33, 1010–1028. [Google Scholar] [CrossRef]
  38. Thomson, A.J.; Shepherd, R.F.; Harris, J.W.E.; Silversides, R.H. Relating weather to outbreaks of western spruce budworm, Choristoneura occidentalis (Lepidoptera: Tortricidae), in British Columbia. Can. Entomol. 1984, 116, 375–381. [Google Scholar] [CrossRef]
  39. Chen, G.; He, Y.; de Santis, A.; Li, G.; Cobb, R.; Meentemeyer, R.K. Assessing the impact of emerging forest disease on wildfire using Landsat and KOMPSAT-2 data. Remote Sens. Environ. 2017, 195, 218–229. [Google Scholar] [CrossRef] [Green Version]
  40. Hicke, J.A.; Johnson, M.C.; Hayes, J.L.; Preisler, H.K. Effects of bark beetle-caused tree mortality on wildfire. For. Ecol. Manag. 2012, 271, 81–90. [Google Scholar] [CrossRef]
  41. Agne, M.C.; Woolley, T.; Fitzgerald, S. Fire severity and cumulative disturbance effects in the post-mountain pine beetle lodgepole pine forests of the Pole Creek Fire. For. Ecol. Manag. 2016, 366, 73–86. [Google Scholar] [CrossRef] [Green Version]
  42. Meigs, G.W.; Campbell, J.L.; Zald, H.S.J.; Bailey, J.D.; Shaw, D.C.; Kennedy, R.E. Does wildfire likelihood increase following insect outbreaks in conifer forests? Ecosphere 2015, 6, 1–24. [Google Scholar] [CrossRef]
  43. Kane, J.M.; Varner, J.M.; Metz, M.R.; van Mantgem, P.J. Characterizing interactions between fire and other disturbances and their impacts on tree mortality in western U.S. Forests. For. Ecol. Manag. 2017, 405, 188–199. [Google Scholar] [CrossRef]
  44. Kulakowski, D.; Jarvis, D. Low-severity fires increase susceptibility of lodgepole pine to mountain pine beetle outbreaks in Colorado. For. Ecol. Manag. 2013, 289, 544–550. [Google Scholar] [CrossRef]
  45. McCarley, T.R.; Kolden, C.A.; Vaillant, N.M.; Hudak, A.T.; Smith, A.M.S.; Kreitler, J. Landscape-scale quantification of fire-induced change in canopy cover following mountain pine beetle outbreak and timber harvest. For. Ecol. Manag. 2017, 391, 164–175. [Google Scholar] [CrossRef] [Green Version]
  46. Shirvani, Z.; Abdi, O.; Buchroithner, M.F.; Pradhan, B. Analysing Spatial and Statistical Dependencies of Deforestation Affected by Residential Growth: Gorganrood Basin, Northeast Iran. Land Degrad. Dev. 2017, 28, 2176–2190. [Google Scholar] [CrossRef]
  47. Abdi, O.; Kamkar, B.; Shirvani, Z.; Teixeira da Silva, J.A.; Buchroithner, M.F. Spatial-statistical analysis of factors determining forest fires: A case study from Golestan, Northeast Iran. Geomat. Nat. Hazards Risk 2018, 9, 267–280. [Google Scholar] [CrossRef]
  48. Abdi, O.; Shirvani, Z.; Buchroithner, M.F. Spatiotemporal drought evaluation of Hyrcanian deciduous forests and semi-steppe rangelands using moderate resolution imaging spectroradiometer time series in Northeast Iran. Land Degrad. Dev. 2018, 29, 2525–2541. [Google Scholar] [CrossRef]
  49. Abdi, O.; Shirvani, Z.; Buchroithner, M.F. Forest drought-induced diversity of Hyrcanian individual-tree mortality affected by meteorological and hydrological droughts by analyzing moderate resolution imaging spectroradiometer products and spatial autoregressive models over northeast Iran. Agric. For. Meteorol. 2019, 275, 265–276. [Google Scholar] [CrossRef]
  50. Setayeshgar, F. The Attack of Lymantria Dispar to 420 Hecatres of Hyrcanian Forests: Outbreak of Defoliators in the Forests. Available online: http://khorasannews.com/?nid=18724&type=0 (accessed on 1 July 2014).
  51. Kia-Daliri, H.; Kazemi-Najafi, S.; Ahangaran, Y. The effect of leaf feeder moth (Erannis defoliaria & Operophtera brumata) on radial growth of 3 tree species in north of Iran (case study: Mashalak, Nowshahr). Iran. J. For. Poplar Res. 2007, 15, 301–309. [Google Scholar]
  52. Ahmadi, A.; Kavosi, M.R.; Soltanloo, H. Zelkova carpinifolia reservoir from Hyrcanian Forests, Northern Iran, a new sacrifice of Ophiostoma novo-ulmi. Biodiversitas 2014, 15, 48–52. [Google Scholar] [CrossRef]
  53. F.A.O. Food Chain Crisis Early Warning Bulletin. Alerts on Threats to the Food Chain Affecting Food Security in Countries and Regions No. 18. 2016. Available online: www.fao.org/food-chain-crisis (accessed on 20 July 2019).
  54. Mirabolfathy, M.; Groenewald, J.Z.; Crous, P.W. The Occurrence of Charcoal Disease Caused by Biscogniauxia mediterranea on Chestnut-Leaved Oak (Quercus castaneifolia) in the Golestan Forests of Iran. Plant Dis. 2011, 95, 876. [Google Scholar] [CrossRef]
  55. Sagheb Talebi, K.; Sajedi, T.; Pourhashemi, M. Forests of Iran. A Treasure from the Past, a Hope for the Future; Springer: Dordrecht, The Netherlands, 2014; ISBN 978-94-007-7371-4. [Google Scholar]
  56. U.S. Geological Survey, Earth Resources Observation Science Center. Collection-1 Landsat OLI Level-2 Surface Reflectance (SR) Science Product; U.S. Geological Survey, Earth Resources Observation Science Center: Reston, VA, USA, 2019.
  57. ASF DAAC. PALSAR Radiometric Terrain Corrected High Res; Includes Material © JAXA/METI 2007; ASF DAAC: Fairbanks, AK, USA, 2019. [Google Scholar]
  58. Department of Natural Resources and Watershed Management of Golestan. The Booklets of Forest Management Plans: Jafarabad, Golriz, Tavir, Kohmian, Zarringol, Sorkhdari, Rezaeian, Qalemorian, Daland, Vatan, Nodeh-Bon-Qale, Minudasht, Takht, Loveh, Sijan-Kiaram, Korankaftar, Terjenli; Department of Natural Resources and Watershed Management of Golestan: Gorgan, Iran, 1981–2011. [Google Scholar]
  59. Google Earth Pro V 7.3.2.5491. Sousara, Golestan Province. 37° 00’ 47.58”N, 55° 16’ 30.88”E, Eye alt 10.17 km; DigitalGlobe 2019; Google LLC: Mountain View, CA, USA, 2019. [Google Scholar]
  60. eCognition Developer. T. 9.0 User Guide; Trimble Germany GmbH: Munich, Germany, 2014. [Google Scholar]
  61. Crippen, R. Calculating the vegetation index faster. Remote Sens. Environ. 1990, 34, 71–73. [Google Scholar] [CrossRef]
  62. Pinty, B.; Verstraete, M.M. GEMI: A non-linear index to monitor global vegetation from satellites. Vegetatio 1992, 101, 15–20. [Google Scholar] [CrossRef]
  63. Kaufman, Y.J.; Tanre, D. Atmospherically resistant vegetation index (ARVI) for EOS-MODIS. IEEE Trans. Geosci. Remote Sens. 1992, 30, 261–270. [Google Scholar] [CrossRef]
  64. Gitelson, A.A.; Kaufman, Y.J.; Merzlyak, M.N. Use of a green channel in remote sensing of global vegetation from EOS-MODIS. Remote Sens. Environ. 1996, 58, 289–298. [Google Scholar] [CrossRef]
  65. Tucker, C.J. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens. Environ. 1979, 8, 127–150. [Google Scholar] [CrossRef] [Green Version]
  66. Jiang, Z.; Huete, A.R.; Didan, K.; Miura, T. Development of a two-band enhanced vegetation index without a blue band. Remote Sens. Environ. 2008, 112, 3833–3845. [Google Scholar] [CrossRef]
  67. 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]
  68. Singh, A.; Harrison, A. Standardized principal components. Int. J. Remote Sens. 1985, 6, 883–896. [Google Scholar] [CrossRef]
  69. Kauth, R.J.; Thomas, G.S. The Tasselled Cap—A Graphic Description of the Spectral-Temporal Development of Agricultural Crops as Seen by Landsat; Purdue University: West Lafayette, USA, 1976. [Google Scholar]
  70. Haralick, R.M.; Shanmugam, K.; Dinstein, I.H. Textural Features for Image Classification. IEEE Trans. Syst. Man Cybern. 1973, 3, 610–621. [Google Scholar] [CrossRef]
  71. Navulur, K. Multispectral Image Analysis Using the Object-Oriented Paradigm; CRC Press: Boca Raton, FL, USA, 2006; ISBN 9780429146305. [Google Scholar]
  72. Moore, I.D.; Grayson, R.B.; Ladson, A.R. Digital terrain modelling: A review of hydrological, geomorphological, and biological applications. Hydrol. Process. 1991, 5, 3–30. [Google Scholar] [CrossRef]
  73. Wilson, J.P.; Gallant, J.C. Terrain Analysis. Principles and Applications; John Wiley & Sons: New York, NY, USA, 2000; ISBN 978-0-471-32188-0. [Google Scholar]
  74. Riley, S.J.; DeGloria, S.D.; Elliot, R. Index that quantifies topographic heterogeneity. Intermt. J. Sci. 1999, 5, 23–27. [Google Scholar]
  75. Salford Systems Ltd. Salford Predictive Modeler 8.3; Salford Systems Ltd.: San Diego, USA, 2018. [Google Scholar]
  76. Hastie, T.; Tibshirani, R.; Friedman, J. The elements of statistical learning: Data mining, inference and prediction. Math. Intell. 2005, 27, 83–85. [Google Scholar] [CrossRef]
  77. Kandel, K.; Huettmann, F.; Suwal, M.K.; Ram Regmi, G.; Nijman, V.; Nekaris, K.A.I.; Lama, S.T.; Thapa, A.; Sharma, H.P.; Subedi, T.R. Rapid multi-nation distribution assessment of a charismatic conservation species using open access ensemble model GIS predictions: Red panda (Ailurus fulgens) in the Hindu-Kush Himalaya region. Biol. Conserv. 2015, 181, 150–161. [Google Scholar] [CrossRef]
  78. ESRI. How Fuzzy Overlay Works. Available online: http://desktop.arcgis.com/en/arcmap/10.5/tools/spatial-analyst-toolbox/how-fuzzy-overlay-works.htm (accessed on 1 May 2019).
  79. Key, C.H.; Benson, N.C. Landscape Assessment: Remote Sensing of Severity, the Normalized Burn Ratio. FIREMON: Fire Effects Monitoring and Inventory System RMRS-GTR-164-CD: LA 1-51; 2006. Available online: http://pubs.er.usgs.gov/publication/2002085 (accessed on 2 August 2019).
  80. Miller, J.D.; Thode, A.E. Quantifying burn severity in a heterogeneous landscape with a relative version of the delta Normalized Burn Ratio (dNBR). Remote Sens. Environ. 2007, 109, 66–80. [Google Scholar] [CrossRef]
  81. ESRI. How Fuzzy Membership Works. Available online: http://desktop.arcgis.com/en/arcmap/10.5/tools/spatial-analyst-toolbox/how-fuzzy-membership-works.htm (accessed on 4 April 2019).
  82. Elhorst, J.P. Spatial Panel Data Models. In Spatial Econometrics: From Cross-Sectional Data to Spatial Panels; Elhorst, J.P., Ed.; Springer: Berlin/Heidelberg, Germany, 2014; pp. 37–93. ISBN 978-3-642-40340-8. [Google Scholar]
  83. Greene, W.H. Econometric Analysis, 5th ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2002; ISBN 0-13-066189-9. [Google Scholar]
  84. Elhorst, J.P. Spatial Panel Data Models. In Handbook of Applied Spatial Analysis; Fischer, M.M., Getis, A., Eds.; Springer: Berlin/Heidelberg, Germany, 2010; pp. 377–407. ISBN 978-3-642-03646-0. [Google Scholar]
  85. Bell, A.; Jones, K. Explaining Fixed Effects: Random Effects Modeling of Time-Series Cross-Sectional and Panel Data. PSRM 2015, 3, 133–153. [Google Scholar] [CrossRef]
  86. Hausman, J.A. Specification Tests in Econometrics. Econometrica 1978, 46, 1251. [Google Scholar] [CrossRef] [Green Version]
  87. Baltagi, B.H. Econometric Analysis of Panel Data, 4th ed.; John Wiley & Sons: Chichester, UK, 2008; ISBN 0-470-01456-3. [Google Scholar]
  88. Carter, G.A. Responses of Leaf Spectral Reflectance to Plant Stress. Am. J. Bot. 1993, 80, 239. [Google Scholar] [CrossRef]
  89. Carter, G.A.; Cibula, W.G.; Miller, R.L. Narrow-band Reflectance Imagery Compared with ThermalImagery for Early Detection of Plant Stress. J. Plant Physiol. 1996, 148, 515–522. [Google Scholar] [CrossRef]
  90. Muchoney, D.M.; Haack, B.N. Change Detection for Monitoring Forest Defoliation. Photogramm. Eng. Remote Sens. 1994, 60, 1243–1252. [Google Scholar]
  91. Moskal, L.M.; Franklin, S.E. Relationship between airborne multispectral image texture and aspen defoliation. Int. J. Remote Sens. 2004, 25, 2701–2711. [Google Scholar] [CrossRef]
  92. Lottering, R.; Mutanga, O.; Peerbhay, K.; Ismail, R. Detecting and mapping Gonipterus scutellatus induced vegetation defoliation using WorldView-2 pan-sharpened image texture combinations and an artificial neural network. J. Appl. Remote Sens. 2019, 13, 014513. [Google Scholar] [CrossRef]
  93. Olsson, P.-O.; Jönsson, A.M.; Eklundh, L. A new invasive insect in Sweden—Physokermes inopinatus: Tracing forest damage with satellite based remote sensing. For. Ecol. Manag. 2012, 285, 29–37. [Google Scholar] [CrossRef]
Sensors 19 03965 g001 550
Figure 1. Study area in the Hyrcanian forest ecoregion in the southeast Caspian Sea. Spatial scatter of the insect defoliators and forest fires samples of 2012–2017 in northeast (NE) Iran. The eastern part of the study area was affected by Lymantria dispar, while the western area was affected by Erannis defoliaria and Operophtera brumata.
Figure 1. Study area in the Hyrcanian forest ecoregion in the southeast Caspian Sea. Spatial scatter of the insect defoliators and forest fires samples of 2012–2017 in northeast (NE) Iran. The eastern part of the study area was affected by Lymantria dispar, while the western area was affected by Erannis defoliaria and Operophtera brumata.
Sensors 19 03965 g001
Sensors 19 03965 g002 550
Figure 2. Site C (Figure 1): The images of Google Earth (GE) [59] and Landsat 5 TM (SWR, NIR, R) show that the forest was in a healthy condition before the attack of Lymantria dispar in 2010 (a, b); the symptoms of the infestation emerged on the Landsat 8 OLI (SWR1, NIR, R) as “olivenite green” in 2013 (c); the insect significantly infested the region based on the GE (d) and Landsat 8 OLI (e) images in 2014, while it declined in 2015 (f) and ended in 2016 (g).
Figure 2. Site C (Figure 1): The images of Google Earth (GE) [59] and Landsat 5 TM (SWR, NIR, R) show that the forest was in a healthy condition before the attack of Lymantria dispar in 2010 (a, b); the symptoms of the infestation emerged on the Landsat 8 OLI (SWR1, NIR, R) as “olivenite green” in 2013 (c); the insect significantly infested the region based on the GE (d) and Landsat 8 OLI (e) images in 2014, while it declined in 2015 (f) and ended in 2016 (g).
Sensors 19 03965 g002
Sensors 19 03965 g003 550
Figure 3. The importance of top predictor variables of Landsat 8 OLI and ancillary features for discerning insect-infested from the non-insect-infested objects using TreeNet. The most important predictor variable (i.e., the mean of the second principal component, or PC2) has gained a value of 100%, and other features were rescaled based on their importance relative to the PC2.
Figure 3. The importance of top predictor variables of Landsat 8 OLI and ancillary features for discerning insect-infested from the non-insect-infested objects using TreeNet. The most important predictor variable (i.e., the mean of the second principal component, or PC2) has gained a value of 100%, and other features were rescaled based on their importance relative to the PC2.
Sensors 19 03965 g003
Sensors 19 03965 g004 550
Figure 4. Univariate partial dependency graphs for the top five-predictor variables for classification of insect-infested forests in NE Iran. Positive values of the outputs indicate a direct association of a class of the categorical variables or values of the continuous variables with the focus class. Eight tree species showing positive relationship with the infested class (a), the mean of PC2 and the mean of red band derived from GLCM show a positive association at values of greater than 0.711 (b) and 126.30 (c), and the mean values of the NDWI and GEMI show a positive relationship with the presence of infestation until the values of 0.267 (d) and 0.697 (e).
Figure 4. Univariate partial dependency graphs for the top five-predictor variables for classification of insect-infested forests in NE Iran. Positive values of the outputs indicate a direct association of a class of the categorical variables or values of the continuous variables with the focus class. Eight tree species showing positive relationship with the infested class (a), the mean of PC2 and the mean of red band derived from GLCM show a positive association at values of greater than 0.711 (b) and 126.30 (c), and the mean values of the NDWI and GEMI show a positive relationship with the presence of infestation until the values of 0.267 (d) and 0.697 (e).
Sensors 19 03965 g004
Sensors 19 03965 g005 550
Figure 5. The insect-affected and non-insect-affected forest-objects derived from the influential object features of Landsat 8 OLI and ancillary data using Geographic Object-Based Image Analysis (GEOBIA) and TreeNet in the Hyrcanian region, NE Iran.
Figure 5. The insect-affected and non-insect-affected forest-objects derived from the influential object features of Landsat 8 OLI and ancillary data using Geographic Object-Based Image Analysis (GEOBIA) and TreeNet in the Hyrcanian region, NE Iran.
Sensors 19 03965 g005
Sensors 19 03965 g006 550
Figure 6. The yearly hazard intensity of defoliation for the time period of 2010 to 2017 (a–h) within the insect-infested forest objects. The monthly anomalies of NDWI were obtained from the time series of Landsat 5, 7, and 8 from 1987 to 2017. The dimensions of forest water content deficit including severity, frequency, and duration were derived from Landsat–NDWI anomalies for the time period of defoliation (2010–2017). The values of these dimensions were standardized by the membership functions of fuzzy linear and fuzzy large between zero and one. The hazard intensity of defoliation was obtained by overlaying the standardized layers of dimensions of the NDWI anomalies using the fuzzy gamma operator within the insect-infested segment objects during a growing season in Hyrcanian forests, NE Iran.
Figure 6. The yearly hazard intensity of defoliation for the time period of 2010 to 2017 (a–h) within the insect-infested forest objects. The monthly anomalies of NDWI were obtained from the time series of Landsat 5, 7, and 8 from 1987 to 2017. The dimensions of forest water content deficit including severity, frequency, and duration were derived from Landsat–NDWI anomalies for the time period of defoliation (2010–2017). The values of these dimensions were standardized by the membership functions of fuzzy linear and fuzzy large between zero and one. The hazard intensity of defoliation was obtained by overlaying the standardized layers of dimensions of the NDWI anomalies using the fuzzy gamma operator within the insect-infested segment objects during a growing season in Hyrcanian forests, NE Iran.
Sensors 19 03965 g006
Table
Table 1. Object features derived from Landsat 8 OLI channels and ancillary data (topographic and forest data) for discriminating defoliated forests from healthy forests in NE Iran.
Table 1. Object features derived from Landsat 8 OLI channels and ancillary data (topographic and forest data) for discriminating defoliated forests from healthy forests in NE Iran.
Object featuresInput dataFeatures 1 (No.)
Spectral features (32)Blue, Green, Red, NIR, SWIR1, SWIR2Mean (6), StdDev (6) of the spectral bands
Mean and StdDev of spectral indices (14) (IPVI [61] GEMI [62], ARVI [63], GVI [64], NDVI [65], EVI2 [66], NDWI [67]
Principal components (6) [68] Greenness (2), Wetness (2) [69], Brightness (1), Max. diff. (1)
Surface texture-features (56)Single bands and all bands in all directionsGLCMall dir. (Homogeneity (7), Contrast (7), Dissimilarity (7), Entropy (7), Mean (7), Angle second moment (7), StDev (7), Correlation (7)) [70,71]
Geometric features (3)ObjectsArea (1), Compactness (1), Asymmetry(1) [71]
Ancillary data (4)ALOS PALSAR
Forest data		Topographic Wetness Index [72], Topographic Position Index [73], Terrain Ruggedness Index [74], Forest types
1 StdDev: Standard deviation; IPVI: Infrared percentage vegetation index; GEMI: Global environment monitoring index; ARVI: Atmospherically resistant vegetation index; GVI: Green vegetation index; NDVI: Normalized difference vegetation index; EVI2: Enhanced vegetation index 2; NDWI: Normalized difference water index; GLCM: gray-level co-occurrence matrix.
Table
Table 2. Classification correctness of test samples for the TreeNet analysis to discriminate the insect-infested from the non-insect-infested forests in NE Iran.
Table 2. Classification correctness of test samples for the TreeNet analysis to discriminate the insect-infested from the non-insect-infested forests in NE Iran.
MeasureAverageOverall accuracySpecificitySensitivityPrecisionF1 statistic
Percent87.1586.7680.5693.7581.0886.96
Table
Table 3. Results of panel data models for relationships between the intensity of insect infestation with the intensity of climate hazards of the current year and previous year as well as forest fires of the previous year in NE Iran.
Table 3. Results of panel data models for relationships between the intensity of insect infestation with the intensity of climate hazards of the current year and previous year as well as forest fires of the previous year in NE Iran.
ModelConstantSPItSPIt − 1TmaxtTmaxt − 1SoilMtSoilMt − 1Firet − 1R2
Common effects0.2040.184 *0.070 ns0.463 *0.165 *0.144 *−0.134 *0.171 *0.680
Fixed effects0.290 **−0.039 ns−0.080 ns0.762 **0.367 *0.152 *−0.072 ns0.154 *0.798
Random effects0.210 **0.153 *0.048 ns0.718 **0.321 *0.146 *−0.126 *0.194 *0.706
** p value < 0.01, * p value < 0.05, and ns: not significant. SPI: standardized precipitation index; Tmax: maximum temperature; SoilM: soil moisture deficit; Fire: forest fire; t: current time; and t − 1: previous time.
Table
Table 4. Results of panel data models for relationships between the severity of forest fires with the intensity of climate hazards and insect infestation of the current year and the previous year in NE Iran.
Table 4. Results of panel data models for relationships between the severity of forest fires with the intensity of climate hazards and insect infestation of the current year and the previous year in NE Iran.
ModelConstantSPItSPIt − 1TmaxtTmaxt − 1SoilMtSoilMt − 1IItIIt − 1R2
Common effects−0.0169 ns−0.085 ns0.0300.385 *0.177 *0.032 ns−0.163 ns−0.059 ns0.105 ns0.210
Fixed effects−0.032 ns−0.113 ns0.094 ns0.254 ns0.213 **0.049 ns−0.126 ns−0.031 ns0.106 ns0.550
Random effects−0.017 ns−0.101 ns0.117 ns0.330 *0.196 **0.041 ns−0.148 ns−0.020 ns0.106 ns0.236
** p value < 0.01, * p value < 0.05, and ns: not significant. SPI: standardized precipitation index; Tmax: maximum temperature; SoilM: soil moisture deficit; II: insect infestation intensity; t: current time; and t − 1: previous time.

Share and Cite

MDPI and ACS Style

Abdi, O. Climate-Triggered Insect Defoliators and Forest Fires Using Multitemporal Landsat and TerraClimate Data in NE Iran: An Application of GEOBIA TreeNet and Panel Data Analysis. Sensors 2019, 19, 3965. https://doi.org/10.3390/s19183965

AMA Style

Abdi O. Climate-Triggered Insect Defoliators and Forest Fires Using Multitemporal Landsat and TerraClimate Data in NE Iran: An Application of GEOBIA TreeNet and Panel Data Analysis. Sensors. 2019; 19(18):3965. https://doi.org/10.3390/s19183965

Chicago/Turabian Style

Abdi, Omid. 2019. "Climate-Triggered Insect Defoliators and Forest Fires Using Multitemporal Landsat and TerraClimate Data in NE Iran: An Application of GEOBIA TreeNet and Panel Data Analysis" Sensors 19, no. 18: 3965. https://doi.org/10.3390/s19183965

APA Style

Abdi, O. (2019). Climate-Triggered Insect Defoliators and Forest Fires Using Multitemporal Landsat and TerraClimate Data in NE Iran: An Application of GEOBIA TreeNet and Panel Data Analysis. Sensors, 19(18), 3965. https://doi.org/10.3390/s19183965

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