Evaluation of habitat suitability for Rhino ( Rhinoceros unicornis ) in Orang National Park using Geospatial tools

The Orang National Park (Orang NP) is one of the important conservation areas in the Brahmaputra valley within North East India Biogeographic Zone covering an area of 78.8 km2. It is one of the prime habitats of one horned rhino in its distribution range in south Asia. Satellite imagery of November 2008 was used to evaluate the rhino habitat pattern in the park. A habitat suitability model for greater one horned rhino Rhinoceros unicornis was prepared using criteria that were derived from ground observations, coefficient of correlations and also from the forest department sources. Result indicates that out of total geographical area of the park 25.85% is covered by woodland mainly composed with moist mixed deciduous forest. About 26.06 % of the total geographical area of the park is covered by wet alluvial grassland and 17.97 % is covered by dry savannah grassland. Similarly degraded grassland is covering 15.23 % and eastern seasonal swamp forest is covering 1.72 % of the total geographic area of the park. About 8.22 % of the park is covered by water body and 6.83% is covered by sandy area. The habitat suitability model for rhino shows that 25.13 % of the total geographic area of the park is most suitable habitat for rhino, 13.62% of the total geographic area is moderately suitable and 61.23% is less suitable habitat for rhino in the park. This information derived out of the analysis reverberates the need of adaptive management of rhino habitat to offer conducive habitats to the rhinos to build up its population to ensure a long term future of rhinos in the Orang NP. Introduction: Wilderness areas for Indian rhinos (Rhinoceros unicornis) continue to shrink and fragment due to multiplicity of natural phenomena as well as ever increasing anthropogenic pressures. Rhinos are in critical demographic crisis; primarily by over-exploitation through poaching for rhino horn and other products and secondarily by loss of habitat due to expanding and developing human populations (Foose & Strien, 1997). Revised IUCN categories and criteria, approved by the 40th meeting of the IUCN council has rated onehoned rhino as vulnerable based upon the parameters viz. population reduction, population estimate and probability of extinction. With the decreasing size of rhino habitat and increasing fragmentation, it has become essential to develop species-specific habitat suitability maps. Habitat is a place occupied by a specific population within a community population (Smith, 1974). Habitat selection is important part of organism’s life history pattern. Roy et al. 2001 states that preservation of wildlife requires a complete knowledge of their spatial requirements commonly referred to as habitat. Habitat evaluation is the assessment of the suitability of land or water as habitat for specific wildlife species. A wildlife habitat suitability map is defined as a map displaying the suitability of land or water as a habitat for a specific wildlife species (Lui et al. 2001). To achieve this one need a model to predict the suitability of land in a given particular set of land conditions. Such model is called a habitat (environmental) suitability model (De Leeuw and Albricht , 1996). The species habitat information is required to be known by the wildlife managers in order to prepare proper habitat suitability analysis of a species in an integrated scientific manner (Parihar et. al. 1986). Conservation biologists and managers need a range of both classical analyses and specific modern tools to face the increasing threats to biodiversity (Caughley, et al. 1996). Among these tools, habitat-suitability modeling using geo-spatial tools has recently emerged as a relevant technique to assess global impacts, for example, those due to climate change, (Berry, et al. 2002) to define wide conservation priorities (Margules, et al. 1994) and to evaluate the completeness of regional nets of protected areas (Araujo, et al. 2002). During the last four decades, development of remote sensing and GIS techniques has made significant contribution in the management of natural resources (Marble et al. 1983 ; Gugan 1993) and environmental monitoring (Kushwaha 1990, 1997). Remote sensing and GIS have been widely used in wildlife habitat studies (Roy et al. 1995; Porwal et al. 1996; Kushwaha et al. 2000, 2004; Hazarika et al. 2005). Remote sensing and GIS technologies together provide vital geo-information support for relevant, reliable and timely information needed for conservation planning (Nellis et al. 1990). GIS has assumed a central role over the years in numerous species-specific applications but there are more scope for GIS in modeling species assemblages, scale-dependent habitat preferences and geographical fragmentation of population, habitat heterogeneity and ecological integrity (Duncan et al. 1995). However these techniques have not been widely used for wildlife habitat studies in Assam. Here in this research geo-spatial tool with a modeling approach was adopted to address the questions of habitat suitability for rhino in Orang National Park (Orang NP). The population of rhino in Orang NP is fluctuating from 35 rhinos in the year 1972 to 97 rhinos in the year 1991 and which is again reduce to 64 rhinos in the year 2009. This unpredictable population fluctuation in the floodplain ecosystem of Orang NP demands habitat suitability evaluation for identifying the key habitat factors and total suitable area for determining fate of rhinos in the park. Materials and Methods: Study Area: Orang NP covering an area of 78.8 km2 is a prime habitat for rhino, located in the north bank of river Brahmaputra and within the administrative boundary of Darrang and Sonitpur districts of Assam, India. This park enjoys a flood plain ecosystem and is a prime habitat for other important species of conservation importance like Royal Bengal Tiger, Asiatic Elephant and different Deer species. The park has been often regarded as the man made forest and it lies within the geographical limits of 26 ̊ 29' N to 26 ̊ 40' N latitude to 92 ̊ 16' E to 92 ̊ 27' E longitude. Figure – 1 show the geographical location of RG Orang NP. Fig. 1. Location of Study Area


Introduction:
Wilderness areas for Indian rhinos (Rhinoceros unicornis) continue to shrink and fragment due to multiplicity of natural phenomena as well as ever increasing anthropogenic pressures.Rhinos are in critical demographic crisis; primarily by over-exploitation through poaching for rhino horn and other products and secondarily by loss of habitat due to expanding and developing human populations (Foose & Strien, 1997).Revised IUCN categories and criteria, approved by the 40th meeting of the IUCN council has rated onehoned rhino as vulnerable based upon the parameters viz.population reduction, population estimate and probability of extinction.With the decreasing size of rhino habitat and increasing fragmentation, it has become essential to develop species-specific habitat suitability maps.Habitat is a place occupied by a specific population within a community population (Smith, 1974).Habitat selection is important part of organism's life history pattern.Roy et al. 2001 states that preservation of wildlife requires a complete knowledge of their spatial requirements commonly referred to as habitat.
Habitat evaluation is the assessment of the suitability of land or water as habitat for specific wildlife species.A wildlife habitat suitability map is defined as a map displaying the suitability of land or water as a habitat for a specific wildlife species (Lui et al. 2001).To achieve this one need a model to predict the suitability of land in a given particular set of land conditions.Such model is called a habitat (environmental) suitability model (De Leeuw and Albricht , 1996).The species habitat information is required to be known by the wildlife managers in order to prepare proper habitat suitability analysis of a species in an integrated scientific manner (Parihar et. al. 1986).Conservation biologists and managers need a range of both classical analyses and specific modern tools to face the increasing threats to biodiversity (Caughley, et al. 1996).
Among these tools, habitat-suitability modeling using geo-spatial tools has recently emerged as a relevant technique to assess global impacts, for example, those due to climate change, (Berry, et al. 2002) to define wide conservation priorities (Margules, et al. 1994) and to evaluate the completeness of regional nets of protected areas (Araujo, et al. 2002).
During the last four decades, development of remote sensing and GIS techniques has made significant contribution in the management of natural resources (Marble et al. 1983 ;Gugan 1993) and environmental monitoring (Kushwaha 1990(Kushwaha , 1997)).Remote sensing and GIS have been widely used in wildlife habitat studies (Roy et al. 1995;Porwal et al. 1996;Kushwaha et al. 2000Kushwaha et al. , 2004;;Hazarika et al. 2005).Remote sensing and GIS technologies together provide vital geo-information support for relevant, reliable and timely information needed for conservation planning (Nellis et al. 1990).GIS has assumed a central role over the years in numerous species-specific applications but there are more scope for GIS in modeling species assemblages, scale-dependent habitat preferences and geographical fragmentation of population, habitat heterogeneity and ecological integrity (Duncan et al. 1995).However these techniques have not been widely used for wildlife habitat studies in Assam.Here in this research geo-spatial tool with a modeling approach was adopted to address the questions of habitat  The methodologies for assessment of habitat suitability condition for rhino in Orang NP was consisted of extensive survey and direct monitoring of rhino, habitat assessment of the park through IRS P6 LISS III satellite imagery of 2008, relationship between rhino presence and habitat parameters and preparation of suitability model using GIS tool based upon habitat parameters.

Direct monitoring of rhino:
A year-long field survey was conducted in Orang NP from September 2008 to September 2009 to understand the habitat utilization pattern of rhino in different seasons.GPS locations of the direct evidences like sighting and indirect evidences like dung piles, foot print and wallowing sign of rhino was taken and plotted over the boundary layer of the park, which was digitized from the map available with the state forest department of Assam using Arc GIS 9.3 software.Throughout the different seasons of the year, rhino was sighted in 183 times in the park.

Habitat assessment using satellite imagery:
Satellite imagery of IRS P6 LISS III of 8 th November 2008 having path and row no.110 and 52 was used to assess the land cover or habitat types in Orang NP.The satellite imagery of 2008 was procured from National Remote Sensing Centre (NRSC), Hyderabad.The imagery was then projected to UTM -WGS 84 projection system using Landsat ETM imagery as reference which was downloaded freely from National Aeronautics and Space Administration's (NASA) Global Land Cover Facilitator's (GLCF) website (www.glcfapp.umiacs.umd.edu).Sub-pixel image to image registration accuracy was achieved through repeated attempts.A radiometric correction of the image was done using dark pixel subtraction technique (Lillesand, et al. 2004).
Subset operation of satellite imagery was carried out by creating an area of interest (AOI) layer of the vector layer of forest boundary of Orang NP, which was digitized from the published maps of department of forest and environment, Govt. of Assam at 1:50,000 scale.After sub setting, the image of the study area was processed through spectral enhancement technique using ERDAS Imagine 9.2 software.Principal component analysis (PCA) was carried out to the image.The imagery was then converted into three principal components.PCA is often used as a method of data compression.It allows redundant data to be compacted into fewer bands-that is, the dimensionality of the data is reduced.The bands of PCA data are non-correlated and independent, and are often more interpretable than the source data (Jensen, 1996).After generating the hybrid PCA image for the year, a supervised classification technique was used using maximum likelihood algorithm to assess the habitat pattern of Orang NP for the year 2008.Since supervised classification is a process where the image analyst supervised the pixel categorization process by specifying to the computer algorithm, numerical descriptors of the various land cover types present in a scene (Lillesand, et al. 2004).Many researchers have been using supervised classification technique to extract the features from the remotely sensed imagery, as it demonstrates the classification that can incorporate both the spectral and spatial features of the pixels in the image resulting in better defined categories in terms of its homogeneity (Fortain, et al., 1999, Dubeni, et al. 2008)

Relationship between rhino presence and habitat parameters:
Wildlife habitat suitability analysis is considered as most important criteria for the conservation and management of wildlife and its habitat (Kushwaha et al. 2000).Such suitability analysis includes a wide variety of factors like habitat pattern, habitat quality, distance from road, availability of water, topography, land cover characteristics including human interferences.It is very essential to understand the relationship between these controlling factors and the species distribution, to make an assessment of the species habitat suitability in a landscape.A variety of analytical techniques have been used to investigate species-environment relationships.These include logistic regression (Pereira et al. 1991;Buckland et al. 1993;Osborne et al. 1992;Walker 1990), discriminant analysis (Haworth, et al. 1990), classification and regression trees (Walker et al. 1988;Skidmore et al. 1996), correlation analysis (Andries et al. 1994) and artificial neuron network (Skidmore et al. 1997).
Here in this research a correlation analysis method was used to understand the relationship between rhino and their habitat.Based upon the direct monitoring of rhino and its relations with habitat types, some habitat parameters were identified.These are habitat types, availability of water, location of human settlement, elevation, distance from roads and their impacts on distribution of rhinoceros in Orang NP.

Habitat suitability modeling:
Based upon the co-efficient of correlation analysis habitat suitability for rhino in the park was categorized in to three categories namely most suitable, moderately suitable and less suitable.The table 1 shows the habitat suitability classes and their respective parameters.

iii) Dry Savannah Grassland:
This type of grassland is dominated by grasses such as Narenga porphyrocoma, Imperata cylindrica, Phragmites karka, Arundo donax, Saccharum spontaneum, Themeda arundinacea, etc.This land cover type is mainly concentrated in the transitional zones of high and low lying areas of the park.This land cover type covers an area of 14.17 km² of the Orang NP.

iv)
Wet Alluvial Grassland: The area under this category of habitat type is 20.54 km².This habitat type is scattered all over the park area.Pure patches of grassland and presence of water characterize it during the rainy season.This wet alluvial grassland plays a critical role in rhinoceros habitat utilization pattern, as rhinoceros prefer to use this habitat throughout the year.This grassland type is mainly composed of Alpinia allughas, Mikania scandens, Saccharum procerum, Pharagmites karka etc.

v) Eastern Seasonal Swamp Forest:
This land cover type occupies an area of 1.36 km² in Orang NP.This type is mainly found along the river Brahmaputra, Dhansiri and Panchnoi and also in and around the wetlands of the park.This is mainly composed of Barringtonia type of vegetation.

vi) Degraded Grassland:
This type of grassland covers an area of 12 km² of Orang NP.These areas are mainly found near the easternmost and westernmost boundary of the park.The leading factor behind the formation of degraded grassland in Orang NP is the over grazing by the domestic cattle that comes from the fringe villages of the park.The impact of invasive species like Mimosa invesa is also a major factor in the formation degraded grassland in Orang NP.

vii) Water Body:
The area under water bodies in Orang NP is 6.48 km².The mighty river Brahmaputra flows along the southern boundary of Orang NP which covers an area of 5.78 km² in the park.Besides this, there are several wetlands in the park, which are also recognized as good habitat for rhinoceros in the park.

viii) & ix) Sandy Area (Dry and Moist):
The area covered by sand in Orang NP is 3.

Rhino-Habitat relationship:
Out of total 183 rhinos sighted in all the seasons of the park, 109 (59.56%) rhinos were found in wet alluvial grassland habitat, 45 (24.59%) rhinos were found in dry savannah grassland habitat and rest 29 (15.84%)rhinos were found in woodlands and wetland habitats.The co-efficient of correlation shows a positive correlation between rhino sighting and wet alluvial grassland (r = 0.582).This indicates that rhino prefers wet alluvial grassland more than dry savannah grassland, woodland and wetland habitats in the park.Thus in selection of habitat suitability parameters, areas with wet alluvial grassland were considered as most suitable for rhino.The fig. 4 shows the relationship

Rhino-Road network relationship:
The correlation coefficient between rhino distribution and road network shows a positive relationship between distance from road and rhino sighting (r = 0.9).
It indicates that with the increase of distance from the road the number of rhino sighting has also increases.Fig. 5 shows the trend of relationship.decreases with the increase of distance from the wetlands.Fig. 6 shows the relationship between rhino sighting and distance from wetlands.During the field survey period it was observed that out of total rhino sighted, 129 (70.49%) rhinos were present within 50 meter elevation in Orang NP.
Total 41 (22.40%) rhinos were found in more than 50 meter height but less than 60 meter height and rest 7 (3.82%)rhinos were found in more than 60 meter elevation zone.The fig. 7 shows the graphical representation of height range and rhino sighting.In the selection of habitat parameters in suitability modeling for rhino, areas less than 50 meter elevation were considered as most suitable for rhino.Wildlife suitability maps and their underlying suitability models have been criticized because of their assumed poor accuracy (Norton et al. 1992).The maps produced by these models have rarely been (Stoms et al. 1992;Williams 1988), although this was clearly advised in the habitat evaluation procedures (USFWS 1981).An accuracy assessment of habitat suitability map was done based on rhino presence and absence and relationship with habitat suitability and rhino sighting.
The result shows that out of 183 rhino sighted, 100 (54.64%) rhinos were found in most suitable habitat of RG Orang NP, 68 (37.16%) rhinos were found in moderately suitable area and rest 15 (8.20%) were found in less suitable habitat of the park.It

Conclusion:
It is evident from the present study that the distribution of rhino in Orang NP completely depends upon the habitat parameters like availability of food plant, distance from water body, distance from road, distance from human settlement, elevation, etc.The habitat parameters have tremendous impact over the habitat utilization and suitability pattern of rhino in Orang NP.From this study it is also clear that most suitable habitat for rhino in Orang NP is only 19.81 km² which is 25.13% of the total geographical area of the park.This indicates that most suitable habitat for rhino in the park is not sufficient for the rhino population that have in the park.
Immediate attention should be taken to conserve the existing suitable habitat for rhino in the park and measures should also be taken to expand the most suitable habitat of the park from 25% to at least 40% to 45% of the total geographical area of the park.The park managers should also take the initiative to increase the wet alluvial grassland habitat in the park, which rhino prefer most in different seasons throughout the year.Finally from this study it is evident that geo-spatial technology has the capability to evaluate the habitat suitability condition of wild animals.
Through spatial modeling in GIS environment it is quite possible to understand the wildlife habitat suitability condition of any wildlife species.
suitability for rhino in Orang National Park (Orang NP).The population of rhino in Orang NP is fluctuating from 35 rhinos in the year 1972 to 97 rhinos in the year 1991 and which is again reduce to 64 rhinos in the year 2009.This unpredictable population fluctuation in the floodplain ecosystem of Orang NP demands habitat suitability evaluation for identifying the key habitat factors and total suitable area for determining fate of rhinos in the park.an area of 78.8 km² is a prime habitat for rhino, located in the north bank of river Brahmaputra and within the administrative boundary of Darrang and Sonitpur districts of Assam, India.This park enjoys a flood plain ecosystem and is a prime habitat for other important species of conservation importance like Royal Bengal Tiger, Asiatic Elephant and different Deer species.The park has been often regarded as the man made forest and it lies within the geographical limits of 26˚ 29' N to 26˚ 40' N latitude to 92˚ 16' E to 92˚ 27' E longitude.Figure -1 show the geographical location of RG Orang NP.

Fig
Fig. 1.Location of Study Area

Fig
Fig.2.Rhino Habitat Suitability Model 87 km².River sand banks devoid of any vegetation are mainly concentrated around the dried river bed of Brahmaputra.The change in course of river Brahmaputra along with excessive siltation during the rainy season has resulted the expansion of such areas.The table 2 shows the area covered by each habitat types and their percentage of area cover.

Fig. 5 .
Fig.5.Correlation between distance from road and rhino sighting

Fig. 6 .
Fig.6.Correlation between distance from wetlands and number of rhino sighting

Fig. 7 .
Fig. 7. Distribution of rhino in Orang NP in different elevation zones

Fig. 8 .
Fig. 8. Correlation between distance from forest camps and number of rhino sighting

Fig. 9 .
Fig. 9. Rhino habitat suitability status of RG Orang NP indicates that the model prepared for the assessment of rhino habitat suitability in RG Orang NP has its validation with the reality.Similarly the relationship between habitat suitability and rhino sighting shows that there is a positive correlation between most suitable habitat and number of rhino sighting (r = 0.682).This indicates that the number of rhino sighting increases with the increase of most suitable habitat in Orang NP.Similarly correlation coefficient was done between less suitable habitat and number of rhino sighting.The result shows a negative correlation between these two variables (r = -0.525).It indicates that the number of rhino sighting decreases with the increase of less suitable habitat in Orang NP.It shows that the output of the model is well correspond to reality, hence it can be assumed that the model is accurate and well validated and this model can be implemented in other rhino bearing areas of Assam like Pabitora Wildlife Sanctuary and Kaziranga National Park.
. Ground truth verification was made during the period from September 2008 to September 2009 and based on the ground verification data, classes were assigned in the PCA based images.Nine land cover types were identified from the field observation and training sets of the land cover classes were gathered using handheld GPS receiver.The nine classes are as follows:

Table : 1. The parameters used for rhino habitat suitability model
Based upon the above mentioned parameters a habitat suitability model for rhino was generated in Arc GIS 9.3 environment.A new tool box in arc tool box was generated under which a model for rhino habitat suitability was designed.The habitat parameters for rhino were placed in the model and spatial analysis tools namely select, buffer, erase, union and intersect were used to get the final habitat suitability map of Orang NP showing the most suitable, moderately suitable and less suitable habitats for rhino in the park.