Close Home
iForest - Biogeosciences and Forestry
vol. 7, pp. 126-131
Copyright © 2014 by the Italian Society of Silviculture and Forest Ecology
doi: 10.3832/ifor1014-007

Research Articles

Monitoring of damage from cedar shoot moth Dichelia cedricola Diakonoff (Lep.: Tortricidae) by multi-temporal Landsat imagery

H Oguz ÇobanCorresponding author, Ramazan Özçelik, Mustafa Avci


Taurus cedar (Cedrus libani A. Rich) is economically and ecologically one of the most important tree species in Turkey, occurring in both naturally regenerated stands and plantations ([6]). This species has reached its largest natural distribution in the Mediterranean region of Anatolia, covering a large area on the Taurus Mountains and providing significant environmental and social benefits as well as high-quality timber and other forest products. Taurus cedar forests cover a total area of 417 188 ha in Turkey, with a current standing volume of approximately 27.4 million m3 and an annual growth of almost 600 000 m3 ([2]). Hughes et al. ([15]) indicated that precipitation is the most important limiting factor on the growth of Taurus cedar trees in the Mediterranean region.

The most common pests of Taurus cedar forests in Turkey are the cedar leaf moth, Acleris undulana Walsingham (Lep.: Tortricidae), and the cedar bark beetle, Orthotomicus tridentatus Egger (Col.: Curculionidae). Another important pest is the cedar shoot moth, Dichelia cedricola Diakonoff (Lep.: Tortricidae), an insect spread all over the western Mediterranean region that can cause significant outbreaks - albeit rarely - in some cedar forests of this area ([3], [4], [5]). D. cedricola is a monovoltine species with adults flying in mid-June ([3]). Initially the larvae feed on needles but then sprouts are also remarkably damaged. The larvae hatch in July and reach their maximum activity in two periods: (a) in summer until November (autumn feeding period); and (b) between April and May (spring feeding period). In the latter period, insects may affect trees more seriously, and especially during the outbreaks trees may be completely defoliated by late May / June.

We monitored the population dynamics of D. cedricola since 1998, when it was first discovered in Kapidag cedar forest, which is also the area studied in this investigation. The infestation area was surveyed regularly in May and June every year, recording population density of the moth, outbreak status and conditions of the infested stands. For 4 consecutive years, between 1998 and 2001, the moth caused outbreaks of varying intensity, whereas the population reached a steady state in 2003. In 2001, when the most intense outbreak occurred, some parts of the forest showed heavily-damaged stands with totally defoliated trees. However, low-damaged and undamaged stands were also detected in the same forest. In the following years, the insects continued to affect the same areas - though to a lesser extent - also during non-epidemic periods.

During an outbreak, larvae feed on tree needles and sprouts making possible to analyse the impact of defoliation using remote sensing data. Severe damage to leaf tissues or total defoliation significantly decrease the reflectance rate in the near-infrared region of the electromagnetic spectrum ([17]). Furthermore, it is possible to identify the mass damage in the stand by means of the Landsat Thematic Mapper (TM), which has a medium spatial resolution ([20], [12], [22], [13], [9]). Besides, the multispectral satellite data such as Spot HRV (High Resolution Visible - [19], [26]) and Modis ([11]), hyperspectral and high spatial resolution remote sensing data such as Compact Airborne Spectrographic Imager (CASI-2 - [25]), and aerial photographs ([28], [1]) were also used in scientific studies for insect damage analysis.

In this study, defoliation damage caused by D. cedricola in cedar forests was examined by using multi-temporal Landsat data. Undamaged, low-damaged and heavily-damaged areas were located by assessing the variation of the Normalized Difference Vegetation Index (NDVI) derived from satellite data. The reliability of the classification based on damages was statistically tested by comparing mean annual ring widths measured on increment cores extracted from sample trees exposed to different damage levels.

Material and Methods 

Study area

The study was carried out in the Kapidag cedar forest, a protected natural area located 70 km from Isparta province (38.08°-38.10° N, 30.67°- 30.75° E), west of Egirdir Lake (Fig. 1). The forest grows between 1300-1900 m a.s.l., mainly north-faced, with a mean slope of about 40%. The total area of the forest, excluding the masked area, was 580 hectares. Taurus cedar is the dominant tree species in the study area. Pure cedar stands have a crown closure of at least 70%.

Fig. 1 - Location of the study area.


Multi-temporal Landsat dataset downloaded from United States Geological Survey website (⇒ http:/­/­earthexplorer.­usgs.­gov) and forest inventory data were used. Such inventory data consists of graphical layers (e.g., stand boundaries, forest roads and rivers, and attributes such as tree species, period of development, and crown closure) re-presenting the forest status in 1997. Dataset, was available in a digital format fully compatible with the ArcGIS software.

The Landsat7-ETM+ (Enhanced Thematic Mapper) image (Path 178, Row 34) dated 5 June 2001 (when the outbreak occurred) was corrected geometrically with the root-mean square error (RMSE) value less than 0.3 pixels using a 1/25 000 scale topographic map. The Erdas Imagine software (version 2011) was used to process raster data such as satellite image. Based on the corrected image, Landsat4-TM (01 June 1988) and Landsat5-TM (19 June 2009) images were georeferenced using the image-to-image registration technique. In order to analyze the vegetation dynamics, NDVI data were used. Song et al. ([24]) suggested to apply the atmospheric correction algorithm to reduce the atmospheric effects on the NDVI signals. Atmospheric dispersion had a non-linear impact with high effect on visible bands and low effect on infra-red bands. In order to minimize the atmospheric effects, the improved Image-Based DOS (Dark Object Subtraction) Model was applied to the image dated 5 June 2001. Using this method the digital information of the image may be converted to surface reflectance values with high accuracy ([8], [18]). “Multiple-date Image Normalization Using Regression” method was applied for the radiometric synchronization between the atmospherically corrected image and the other image data ([16]).

The NDVI data were obtained by using the red (band 3rd) and near-infrared (band 4th) bands of the radiometrically corrected images ([21]). Healthy green plants usually reflect 40-50% of the infra-red wavelengths (0.7-1.1 μm), while plant chlorophyll absorbs approximately 80-90% of the energy in the visible (0.4-0.7 μm) part of the electromagnetic spectrum. The reflectance values fall down in the infrared region if the plants either have partially-damaged green tissues or are dead ([16]). By using NDVI such reflectance variations of the plants may be detected more easily.

A subset of NDVI data, which was large enough to cover the study area, was considered. In the subset images, open areas such as roads, streams and rocks and mixed stands were masked using the stand type and topographic maps.

Image data processing

Stands damaged by the insect outbreak were clearly detectable from the NDVI data (NDVI5) obtained from the image dated 5 June 2001 (Fig. 2). Pixels representing the damaged areas are indicated by near-black gray because the NDVI values fell down. Moreover, a transition zone was observed between the undamaged and the heavily- damaged areas. Based on field observation, insect population started to increase in 1998 and reached its peak in 2001. It was also observed that during the outbreak stands were affected to varying extents.

Fig. 2 - NDVI5 sub-areas and sample points where the increment cores were extracted.

Based on ground-truth and NDVI data, stands in the study area were classified into three categories: undamaged (UD), low-damaged (LD) and heavily-damaged (HD) areas. Damage-level thresholds were set as follows: the NDVI5 data that had pixel values ranging from 0.26 to -0.23 was divided into three groups: greater that 0; from 0 to -0.115; and from -0.115 to -0.23 (see Fig. 3). Some 40 points were randomly selected over the studied area, with a minimum of 10 points in each group. For each point, an NDVI sub-area covering a total of 9 pixels was considered, including the selected point and its 8 surrounding pixels (Fig. 2). The mean NDVI values were calculated for the NDVI sub-areas in the ArcGIS environment. At that stage, some points were displaced manually in order to keep the NDVI sub-areas within the same pixel groups. The NDVI sub-areas identified on NDVI5 were transferred onto the other NDVI data used in this study. Changes occurring in the green vegetation of trees exposed to the outbreak were also analyzed (Fig. 3).

Fig. 3 - Mean NDVI values of the NDVI sub-areas and threshold boundaries for different damage levels.

The NDVI data varied depending on the density of the green-leaf vegetation in the stand. Holben ([14]) stated that high and low-density changes affect the NDVI data similarly. Based on a similar approach, we classified the pixels in the range 0.13-0.26 as high green-leaf vegetation density, and pixels in the range 0.06-0.13 as low green-leaf vegetation density. As some of the values lower than 0.06 still belonged to the green-leaf vegetation, various thresholds were tested with an interval of 0.01. At that stage, the categories classified according to each threshold were overlapped with the inventory data, the NDVI data from other years without an outbreak (1988 and 2009) and the high spatial resolution images on Google Earth®, and then they were checked. Furthermore, information obtained from ground observations was also assessed. In summary, it was found that the range from 0.06 to 0.02 consisted of the pixels that had very low green-leaf vegetation density. The pixels were adjacent to the masked roadside and streamside as well as the open areas such as the bare soil and rocks. The rest of the NDVI5 values were mostly negative. As the study was carried out only on pixels representing the totally green-leaf vegetation, it was assumed that pixels with negative values were exposed to the insect outbreak. For that reason, the NDVI values’ range from 0.02 to -0.23 was divided into two parts: the range from 0.02 to -0.12 represented the low-damaged area (LD), the range from -0.12 to -0.23 represented the heavily-damaged area (HD), while the range from 0.26 to 0.02 represented the undamaged areas (UD). According to that classification, the mean NDVI values of the non-epidemic years 1988 and 2009 fell in the undamaged level (Fig. 3). NDVI5 values for year 2001 were re-coded according to the above thresholds and all the pixels were consequently classified in the above categories.

Four different sample plots were randomly distributed in each damage class. The coordinates of sample plots were installed in the Garmin GPSMAP® 62s hand-held GPS (Global Positioning System) receiver.

Field and laboratory studies

In autumn 2012, an increment core was extracted at breast height from each of ten cedar trees sampled in each sample plot. Cores were labeled and sealed in plastic bags to preserve moisture and prevent shrinkage.

Habitat characteristics of the chosen sample plots were assumed to be similar. As suggested by Carus & Avci ([7]), if non-host and host areas responded in a similar manner to climatic variations, then the differences between standardized ring chronologies of non-host and host trees would primarily reflect non-climatic environmental variations, such as the effects of the insect damage.

Plots were arbitrarily selected over the whole studied area in order to span the whole range of variation of cedar in terms of age, stand density and site characteristics. Sample trees were selected within each plot from both dominant and co-dominant crown class. Trees having multiple stems, broken tops, cankers or crooked boles were excluded from sampling. Total height of each sampled tree was measured to the nearest 0.05 m. Diameter outside bark (dob) at breast height (1.3 m) was measured and recorded to the nearest 0.3 cm. Descriptive statistics for sampled trees in sample plots are reported in Tab. 1.

Tab. 1 - Descriptive statistics for sampled trees. (SD): standard deviation; (dob): diameter outside bark.

Annual ring widths of 40 increment cores extracted from 4 sample plots for each damage level identified on the NDVI5 image were measured to the nearest 0.01 mm using a stereomicroscope (Fig. 4). Overall, 120 increment cores were measured for the three different damage classes. All cores consisted of approximately 40 annual growth rings, the most reliable chronologies spanning over the period 1970-2012. In this study, only the mean annual ring widths in the year 2001 were taken into consideration.

Fig. 4 - Annual rings on increment core from heavily-damaged sample plot. Ring width was extremely narrow for the year 2001.

Differences among defoliation classes in annual ring widths in 2001 were tested by one-way analysis of variance. Mean values were compared using the Tukey’s range test (α = 0.05) to detect possible differences between different defoliation levels and non-defoliation group. All statistical analyses were performed using the software package SPSS® version 11.0.


Image pixels of the studied area were classified as the undamaged, low-damaged and heavily-damaged classes according to the prefixed thresholds based on NDVI data (Fig. 5). In summary, 73% of the forest was found to be undamaged, 14% to be low-damaged, while 13% was found to be heavily-damaged. Misclassification was observed for some pixels belonging to the low-damaged area, previously classified as undamaged. Combining the visual interpretation of images and ground-level information data, misclassified pixels were located near high rocky slopes without vegetation, where the stand crown closure was poor and the soil exposed to erosion because of steep inclination. The above pixels had a total surface area of 18 hectares, representing about 23% of the pixels classified in this category. The error rate of the heavily-damaged area was approx. 0.2%.

Fig. 5 - NDVI5 data classified according to the density of the insect damage. (M): undamaged; (S): low-damaged; (K): heavily-damaged sample plots.

The effect of the insect defoliation on wood increments before, during and after the outbreak is shown in Fig. 6. Mean annual ring widths were 0.54 mm for heavily-damaged areas, 1.04 mm for low-damaged areas and 1.64 mm for undamaged areas (Tab. 2). The above mean values are fairly in accordance to the degree of tree defoliation inferred from satellite images.

Fig. 6 - Mean values of annual tree ring width calculated from 40 trees at each damage level.
Tab. 2 - Summary statistics and comparison of severity of defoliation for radial growth. (SD): Standard deviation; (df): degree of freedom.

Significant differences in mean annual ring width among the three defoliation classes (UD, LD, HD) were found after ANOVA (p<0.001 - Tab. 2). As expected, defoliated cedar trees showed the greatest decreases in growth.

When compared with the year 2001, the mean annual ring widths of the consecutive 11 years increased by 14% in the low-damaged trees and by 112% in the heavily-damaged trees and decreased by 3% in the undamaged trees.


Damage caused by the D. cedricola outbreak in Taurus cedar stands in Turkey was identified by means of the Landsat TM bands. According to the NDVI data obtained from the satellite image of 2001, defoliation of the stands was analyzed by classifying the pixels into three different categories. The slight decrease in the mean annual ring widths of the undamaged trees in 11 years after defoliation may be due to an age-related decrease in the increment, as well as to changes in the climatic conditions. On the other hand, increased growth as the insect population diminished was observed especially for trees that were more heavily affected by the outbreak. However, considering that the general ecological and climate conditions affected the entire study area in the same way, such changes should also affect the increment rates of the damaged stands. Between 1970 and 2012 the wood increment in the damaged stands was always lower than those of the undamaged stands; this may be due to the presence of insect populations in these areas.

Extensive damages to forest trees caused by insects have been previously assessed using Landsat TM data ([20], [22]). In this study, the accuracy of classification based on NDVI5 data diminished in isolated locations and in mixed forests. High spatial resolution and hyperspectral data can also help in locating isolated damaged areas, as well as increase the accuracy of damage classification ([25]).

Previous studies used vegetation indexes or other measures to examine canopy defoliation caused by a variety of insects ([13], [1]). Landsat multi-temporal images have been proven to help in the detection of actual insect defoliation, although their use in mixed stands still present some uncertainty. Heterogeneity of forest tree species also limits the ability of medium resolution remote sensors in determining insect damage ([27], [23], [10]). Nevertheless, in this study Landsat images provided detailed and accurate information of the insect epidemic occurring in the monitored homogenous forest. In general, high spatial resolution data, such as QuickBird, WorldviewII etc, are preferred for the identification of small clusters of defoliated trees. However, the relative costs, availability and processing requirements of the various data sources are all important considerations for the end-user ([27]). Indeed, the use of Landsat data is more cost-effective and less time-consuming than other conventional methods.

In conclusion, the use of Landsat bands may be considered a useful and reliable tool to map and monitor insect damage occurring in cedar stands. With respect to specific thresholds to map insect damage, we found that pixels representing damaged areas had NDVI values lower than 0. Based on the analysis of the increment cores, significant differences between the annual ring widths of trees exposed to different damage levels were observed. These findings indicate that damage mapping can be performed remotely by using NDVI values.


We thank Selçuk Kütük and Umut Akbulut, students of Süleyman Demirel University Faculty of Forestry, for providing assistance during the field studies; and master’s students Serdar Göktepe and Mustafa Önder Ersin. We have been indebted to Prof. M. Dogan Kantarci (Faculty of Forestry, Istanbul University) for his valuable recommendations.


Adelabu S, Mutanga O, Cho MA (2012). A review of remote sensing of insect defoliation and its implications for the detection and mapping of Imbrasia belina defoliation of Mopane woodland. The African Journal of Plant Science and Biotechnology 6: 1-13.
::Online::Google Scholar::
Anonymous (2006). Forest Resources. Turkish Forest General Directorate Press, Ankara, Turkey.
::Google Scholar::
Avci M (2000). Biology, damage and natural enemies of a new Cedar pest in Turkey Dichelia cedricola (Diakonoff) 1974 (Lep.: Tortricidae). In: Proceedings of the “4th Turkish National Congress of Entomology”. Aydin (Turkey), 12-15 September 2000, pp. 447-454.
::Google Scholar::
Avci M (2004). Two important harmful insect species for Turkish cedar forests. Forest Engineering Periodical, Chamber of Forest Engineers publishing, Ankara, Turkey, vol. 41, pp. 7-9.
::Google Scholar::
Avci M, Sarikaya O (2009). Orthotomicus tridentatus Eggers: distribution and biology in cedar forests of Turkey. Turkish Journal of Agriculture and Forestry 33: 277-283.
::Online::Google Scholar::
Boydak M (2003). Regeneration of Lebanon cedar (Cedrus libani A. Rich.) on karstic lands in Turkey. Forest Ecology and Management 178: 231-243.
::CrossRef::Google Scholar::
Carus S, Avci M (2005). Growth loss of Lebanon cedar (Cedrus libani) stands as related to periodic outbreaks of the cedar shoot moth (Dichelia cedricola). Phytoparasitica 31: 118-123.
::CrossRef::Google Scholar::
Chavez PS (1996). Image-based atmospheric corrections-revisited and improved. Photogrammetric Engineering and Remote Sensing 62: 1025-1036.
::Online::Google Scholar::
Çoban HO, Eker M, Avci M (2011). Identification and monitoring the damage of Dichelia cedricola (Diakonoff) by using multi-temporal Landsat imagery. In: Proceedings of the “1st Forest Entomology and Pathology Symposium”. Antalya (Turkey), 23-25 October 2011, pp. 178.
::Google Scholar::
DeRose RJ, Long JN, Ramsey RD (2011). Combining dendrochronological data and the disturbance index to assess Engelmann spruce mortality caused by a spruce beetle outbreak in southern Utah, USA. Remote Sensing of Environment 115: 2342-2349.
::CrossRef::Google Scholar::
Eklundh L, Johansson T, Solberg S (2009). Mapping insect defoliation in Scots pine with MODIS time-series data. Remote Sensing of Environment 113: 1566-1573.
::CrossRef::Google Scholar::
Hall RJ, Fernandes RA, Hogg EH, Brandt JP, Butson CR, Case BS, Leblanc SG (2003). Relating aspen defoliation to changes in leaf area derived from field and satellite remote sensing data. Canadian Journal of Remote Sensing 29: 299-313.
::CrossRef::Google Scholar::
Hall RJ, Skakun RS, Arsenault E (2006). Remotely sensed data in the mapping of insect defoliation. In: “Understanding Forest Disturbance and Spatial Pattern: Remote Sensing and GIS Approaches” (Wulder MA, Franklin SE eds). Taylor & Francis, CRC Press, Boca Raton, FL, USA.
::Google Scholar::
Holben BN (1986). Characteristics of maximum-value composite images from temporal AVHRR data. International Journal of Remote Sensing 7: 1417-1434.
::CrossRef::Google Scholar::
Hughes MK, Kuniholm PI, Eischeid JK, Garfin G, Griggs CB, Latini C (2001). Aegean tree-ring signature years explained. Tree-ring Research 57: 67-73.
::Online::Google Scholar::
Jensen JR (1996). Introductory digital image processing (2nd edition). Prentice-Hall Publisher, Series in Geographic Information Science, USA.
::Google Scholar::
Lillesand TM, Kiefer RW (1999). Remote sensing and image interpretation (4th edn). John Wiley & Sons Inc., Hoboken, NJ, USA, pp. 804.
::Google Scholar::
Lu D, Mausel P, Brondizio E, Moran E (2002). Assessment of atmospheric correction methods for Landsat TM data applicable to amazon basin LBA research. International Journal of Remote Sensing 23: 2651-2671.
::CrossRef::Google Scholar::
Muchoney DM, Haack BN (1994). Change detection for monitoring forest defoliation. Photogrammetric Engineering and Remote Sensing 60: 1243-1251.
::Online::Google Scholar::
Price KP, Jakubauskas ME (1998). Spectral retrogression and insect damage in lodgepole pine successional forests. International Journal of Remote Sensing 19: 1627-1632.
::CrossRef::Google Scholar::
Rouse JW, Haas RH, Shell JA, Deering DW (1974). Monitoring vegetation systems in the great plains with ERTS. In: Proceedings of the “Third Earth Resources Technology Satellite-1 Symposium” (Freden SC, Mercanti EP, Becker MA eds). Washington (DC, USA) 10-14 December 1973. NASA, Goddard Space Flight Center, MD, USA, pp. 309-317.
::Online::Google Scholar::
Skakun RS, Wulder MA, Franklin AE (2003). Sensitivity of the thematic mapper enhanced wetness difference index to detect mountain pine beetle red-attack damage. Remote Sensing of Environment 86: 433-443.
::CrossRef::Google Scholar::
Somers B, Verbesselt J, Ampe EM, Sims N, Verstraeten WW, Coppin P (2010). Spectral mixture analysis to monitor defoliation in mixed-aged Eucalyptus globulus Labill plantations in southern Australia using Landsat 5-TM and EO-1 Hyperion data. International Journal of Applied Earth Observation and Geoinformation 12: 270-277.
::CrossRef::Google Scholar::
Song C, Woodcock CE, Seto KC, Lenney MP, Macomber SA (2001). Classification and change detection using Landsat TM data: When and how to correct atmospheric effects? Remote Sensing of Environment 75: 230-244.
::CrossRef::Google Scholar::
Stone C, Coops NC (2004). Assessment and monitoring of damage from insects in Australian eucalypt forests and commercial plantations. Australian Journal of Entomology 43: 283-292.
::CrossRef::Google Scholar::
White JC, Wulder MA, Grills D (2006). Detecting and mapping mountain pine beetle red-attack damage with SPOT-5 10-m multispectral imagery. BC Journal of Ecosystems and Management 7: 105-118.
::Online::Google Scholar::
Wulder MA, White JC, Bentz B, Alvarez MF, Coops NC (2006). Estimating the probability of mountain pine beetle red-attack damage. Remote Sensing of Environment 101: 150-166.
::CrossRef::Google Scholar::
Wulder MA, White JC, Grills D, Nelson T, Coops NC, Ebata T (2009). Aerial overview survey of the mountain pine beetle epidemic in British Columbia: communication of impacts. BC Journal of Ecosystems and Management 10: 45-58.
::Online::Google Scholar::


Çoban HO, Özçelik R, Avci M (2014).
Monitoring of damage from cedar shoot moth Dichelia cedricola Diakonoff (Lep.: Tortricidae) by multi-temporal Landsat imagery
iForest - Biogeosciences and Forestry 7: 126-131. - doi: 10.3832/ifor1014-007
First Previous Next Last
© iForest

Download Reference

Paper ID# ifor1014-007
Title Monitoring of damage from cedar shoot moth Dichelia cedricola Diakonoff (Lep.: Tortricidae) by multi-temporal Landsat imagery
Authors Çoban HO, Özçelik R, Avci M
Close Download