Analysis of Temporal and Spatial Changes in the Vegetation Density of Similipal Biosphere Reserve in Odisha ( India ) Using Multitemporal Satellite Imagery

National parks and protected areas require periodic monitoring because of changing land cover types and variability of landscape contexts within and adjacent to their boundaries. In this study, remote sensing andGIS techniques were used to analyse the changes in the vegetation density particularly in the zones of higher anthropogenic pressure in the Similipal Biosphere Reserve (SBR) of Odisha (India), using Landsat imagery from 1975 to 2005. A technique for the detection of postclassification changes was followed and the change in vegetation density as expressed by normalized difference vegetation index was computed. Results indicate that high dense forest in the core zone has been conserved and the highest reforestation has also occurred in this zone of SBR.The results also reveal that anthropological interventions are more in the less dense forest areas and along the roads, whereas high dense forest areas have remained undisturbed and rejuvenated. This study provides baseline data demonstrating alteration in land cover over the past three decades and also serves as a foundation for monitoring future changes in the national parks and protected areas.


Introduction
Promoting and facilitating the use of remotely sensed images for the monitoring and management of national parks and other protected areas have been the focus of many research and governmental agencies.Monitoring landscape dynamics can answer questions of what the spatial extent of land cover types within and adjacent to the parks and the protected areas and is how they have changed over time.Vegetation mapping is a very important application of remotely sensed data.Spectral profiles of vegetation clearly show that the peak reflectance can be found in the near-infrared wavelengths mainly because of the internal structure of "green" leaves, and the greatest absorption is in the red wavelengths because of the presence of chlorophyll pigments.The normalized difference vegetation index (NDVI) is one of the most successfully used vegetation indices for land cover classification [1][2][3][4].It is also used as an environmental indicator to monitor temporal and spatial variation in vegetation density as well as the health and viability of plant cover over time [5][6][7][8] for the derivation of biophysical properties [9][10][11] and for estimation of net primary production [12][13][14].The NDVI is correlated with many biophysical parameters of the vegetation canopy such as leaf area index (LAI), fractional vegetation cover, vegetation condition, and biomass.
The NDVI values increase as green cover density increases [15]; the values are calculated from red and nearinfrared reflection values [5,6,16], which have a wide range used in vegetation indices [17].Vegetated landscapes typically have NDVI values ranging from 0.1 in the desert to 0.8 in dense tropical rain forest [18,19].It was demonstrated that when NDVI values are between 0.0 and 0.2, the pixel is considered bare soil; for NDVI values greater than 0.5, the pixel is considered fully vegetated.
The Similipal is a densely forested hill range in the heart of the erstwhile Mayurbhanj State (now Odisha), lying close to the eastern-most end of the Eastern ghats (India).Located in the Mahanadian Biogeographical Region and within the Biotic Province, Chhotanagpur Plateau, Similipal is the richest watershed in Odisha (India) giving rise to many perennial rivers like the Budhabalanga, Hadkei, Khairi, Bhandan, West Deo, Salandi, East Deo, Sanjo, and Palpalá.This biosphere reserve is unique in many respects; notable among these are its flora, fauna, landscapes, tribes, and waterfalls.
Most of the people within the Biosphere Reserve have a very low level of skill.Agriculture is not well developed and employment opportunities are very poor.Most of the local inhabitants of the reserve derive their livelihood and income from the forest wood and products.This has led to an undue pressure on wildlife too.Therefore, there is an urgent need to monitor the temporal and spatial status of the forest reserves.With this background, the present investigation was carried out with the broad objective to study the changes taking place in the vegetation density of different areas of high anthropological pressure in the Similipal Biosphere Reserve (SBR) for over a period of 30 years.

Study Area.
The Similipal is a densely forested hill range in the heart of District Mayurbhanj (Odisha, India) and located between 20 ∘ 17  and 22 ∘ 34  north latitude and 85 ∘ 40  and 87 ∘ 10  east longitude.Similipal is unique in many respects, notable among which are its flora, fauna, landscapes, waterfalls, and native tribal population.The landscape is beautifully scattered with numerous small and high hills densely covered with vegetation.The highest mountain peak is Khairiburu which is approximately 1168 m, above the sea level.Similipal is the richest watershed in the state of Odisha giving rise to many perennial rivers.Natural colour composite of the SBR along with core, buffer, and transition boundaries is presented in Figure 1.People residing in the Reserve area are largely tribal.Due to low level of skills, lower educational levels, and sociocultural traits, they are mainly dependent on local resources.There are four villages inside the core, namely, Kabatghai, Jenabil, Jamunagarh, and Bakua.
There are 61 villages in the buffer zone.The climate of the area is monsoonal and divisible into three seasons: summer (March-June), rainy (July-October), and winter (November-February).Premonsoon showers are received during May and June.Postmonsoon showers are received during November and December.The natural vegetation is moist deciduous type [20] and is dominated by Shorea robusta, Anogeissus latifolia, Buchanania lanzan, Dillenia pentagyna, Syzygium cumini, and Terminalia alata.The area and status of SBR since 1979 are presented in Table 1.Similipal was declared as a biosphere reserve by Government of India on June 22 1994, due to its rich biodiversity and natural heritage.There are 1076 plants recorded from the area including 60 species of ferns, 92 species of orchids, and two gymnosperms [21,22].
Our study sites were located in places of high anthropological pressure zones in order to evaluate the impact of human induced disturbances on the vegetation density.Four core villages along with buffer region surrounding roads and higher elevation points were selected for this study.Topographic maps of Survey of India with scale of 1 : 50,000 were used mainly for geometric correction of the satellite images and for some ground truth information.

Geometric Correction.
After atmospheric correction and elimination of offset value of satellite data, in order to prepare two or more satellite images for an accurate detection of change and comparison, it is imperative to geometrically rectify the imagery.For accurate registration of the different time periods, the images were georeferenced to the Universal Transverse Mercator (UTM) projection and WGS84 datum using 40-45 ground control points (GCPs) with a root mean square error not exceeding 0.5 pixels [23].Accordingly, the 2005 image was used as the base image to more precisely georeference the other scenes.The correction was made using first order polynomial transformation model and nearest neighbor method for resampling [23].The geo-referenced images were then clipped with vector masks for the transition, buffer and core boundary to generate the areas of interest.leaf biomass [24] and calculated by the following formula: NDVI = (IR − R)/(IR + R), where IR refers to infrared light and R refers to red light [25].For its measurement, scientists use satellite sensors that observe the distinct wavelengths of visible and near-infrared sunlight which is absorbed and reflected by the plants.Then the ratio of visible and nearinfrared light reflected back up to the sensor is calculated.The ratio gives a number from minus one (−1) to plus one (+1).An NDVI value of zero means no green vegetation and close to +1 (0.8-0.9) indicates the highest possible density of green leaves.The group of pixels having NDVI values from 0 to 0.2 were categorized under canopy density class of <10% and 0.2-0.4 as canopy density class of 10-40% and, the group  of pixels having NDVI value >0.4 were kept under the canopy density class of 40%.Normalized difference vegetation index (NDVI) was used to prepare a forest density map that was categorized into three canopy density classes: <10% (nonforest), 10-40% (open), and >40% (medium/high/very high) depending on the NDVI value.Image elements like tone, texture, shape, size, shadow, location, and association were also evaluated to aid in the class delineations.

Results and Discussion
There are many techniques for the detection of change(s) ranging from visual comparison to detailed quantitative approaches [26].The fact that sums and differences of bands are used in the NDVI rather than absolute values makes the NDVI more appropriate for use in studies where comparisons over time for a single area are involved.This is because the NDVI might be expected to be influenced to a lesser extent by variation in atmospheric conditions.An assessment of change of vegetation cover in terms of NDVI between 1975, 1990, and 2005 was conducted using an analysis involving detection of postclassification change(s), and the results are presented in Tables 3 and 4 and Figures 4 and 5.

High Elevation
Points.The higher elevation points/places are supposed to be sacred places of worship and other activities for tribal people living inside the forest.For this purpose, 4 km buffer area surrounding the higher elevation points was selected for analysis.It is observed that the area under <0.2-0.4 class has decreased greatly over the years (Table 3, Figures 2 and 3).During the same period, the area under >0.4 NDVI class has increased though by approximately 1,000 ha, yet only marginally if examined in terms of percent increase.This shows that less dense forest area has undergone deforestation to a greater extent, whereas  4 and Figures 4 and 5).This shows that less dense forest area has undergone deforestation whereas medium dense forest areas have remained almost constant.However, high dense forest areas surrounding the major roads inside the biosphere reserve have decreased marginally over a period of 30 years.In the area surrounding the roads, anthropological intervention seemed to have increased threefold in less dense forest areas since 1975 through 2005.However, it remained constant in medium dense forest areas and marginally increased in highly dense forest areas.

Buffer Area
Surrounding Four Core Villages.Four core villages, namely, Bakua, Jenabil, Jamunagarh, and Kabatghai, are located inside the core zone of the biosphere reserve.Five km buffer area surrounding each of these villages has been analysed for studying the change in vegetation density as expressed in terms of NDVI over a period of 30 years.The results are presented in Table 5 and Figures 6 and 7.
The five km buffer areas surrounding all the four core villages present a very good picture of human-nature coexistence and irrespective of human settlement; there has been an increase in the vegetation density.In all the villages, there is tremendous decrease in the area under <0.2 NDVI class and in fact zero in three villages.These areas are mostly covered by dense vegetation (>0.4) throughout the decadal period.The buffer area surrounding all the four villages shows a sign of greening up as the area under >0.4 NDVI class increased to a maximal extent which constitutes the maximal proportion of area of analysis.

Conclusion
Based on the analysis of postclassification changes observed from the time series landsat images it was inferred that the high dense forest in the core zone has been conserved and the highest reforestation had also occurred in this zone of SBR.Over the years, forest density in terms of NDVI, as evidenced from the time series Landsat images, has increased as the area under >0.4 NDVI class has shown a positive trend throughout the zones of analysis.This also indicates that anthropological interventions are more in the less dense forest areas and along the roads, whereas high dense forest areas have remained undisturbed and rejuvenated.Further study involving extensive field work and survey is suggested for the analysis of species level vegetation change(s) based on high resolution satellite images.

Figure 1 :
Figure 1: Natural colour composition of the SBR along with core, buffer, and transition boundaries.

Figure 2 :
Figure 2: Vegetation density map in terms of NDVI classes in 5 km buffer area surrounding the highest elevation points for 1975.

3. 1 .
Generation of NDVI Image.The normalized difference vegetation index (NDVI) is strongly correlated to the green

Figure 3 :Figure 4 :
Figure 3: Vegetation density map in terms of NDVI classes in 5 km buffer area surrounding the highest elevation points for 2005.

Figure 5 :
Figure 5: Vegetation density map in terms of NDVI classes in 3 km buffer area surrounding dense road points for 2005.

Figure 6 :
Figure 6: Vegetation density map in terms of NDVI classes in 5 km buffer area surrounding Bakua, Jenabil, Jamunagarh, and Kabatghai villages in SBR for 2005.

Figure 7 :
Figure 7: Vegetation density map in terms of NDVI classes of 5 km buffer area surrounding Bakua, Jenasbil, Jamunagarh, and Kabatghai for 2005.

Table 1 :
Area and status of SBR since 1979.

Table 2 :
Details of the Landsat images used in the study.

Table 3 :
Area occupied by different NDVI classes in different years in the 4 km buffer zones surrounding high elevation points in SBR.

Table 4 :
Area occupied by different NDVI classes in different years in the 3 km buffer area surrounding dense road points in SBR.

Table 5 :
Area occupied by different NDVI classes in different years in the 5 km buffer area surrounding four villages in SBR.