Site Suitability Analysis For Willow Tree Plantation In Uttarakhand Using GIS

The objective of a present study was to evaluate the land of Uttarakhand for suitability of Willow. The study area located between latitude 28o43’N and 31o27’N and longitude 77o34’E and 81o02’E. and covers approximately 53,485 sq km. A suitability resulting from the overlay process of the identified theme layers has unique information of tree requirement on which the suitability is based. The identified theme layers include temperature (maximum, minimum and average temperature), precipitation and soil texture. All thematic layers with their associated attribute data were encoded in GIS database. Overlay operation was performed on these layers through the model developed using Model Builder module of ARC-view 3.2. The model was a result of the crop requirements and the range of input variables in the state. Model was applied to formulate the suitability classes. The results indicate that 38%, 29.3%, 13.2% and 19.5% area of the total geographical area of the state fall under the most suitable, suitable less suitable and not suitable categories, respectively for Willow plantation. On the basis of the results, it can be recommended that Willow can be grown on the larger area of Uttarakhand state, especially the wastelands can be utilized for cultivation of Willow, which will not only be beneficial for socio-economic point of view but also sequester large amount of the carbon from the atmosphere.


INTRODUCTION
Most of the plant species in India are exploited for their commercial values with least or no consideration for their sustainable exploitation.Willow is one such species on which several local level livelihoods are dependent.As a plant, willow is used for in the manufacture of boxes, broom, cricket bats, cradle broads, chairs and other furniture, dolls, flutes, poles, sweat lodges, toys, turnery, tool handles, veneer, wands and whistles.In addition, tannin, fibre, paper, rope and string can be produced from the wood.The leaves and barks of the willow tree have been mentioned in ancient texts from Assyria, Sumer and Egypt as a remedy for aches and fever, willows contains salicin a substance that chemically resembles aspirin (Breasted, 2007).Willow is also used for biofiltration, ecological wastewater treatment systems, hedges, landscaping, phytoremediation streambank stabilization (bioengineering), slope stabilization, soil erosion control, shelterbelt and windbreak, soil reclamation tree and wildlife habitat.Willow is grown for biomass or biofuel, in energy forestry systems, as a consequence of its high energy in-energy out ratio, large carbon mitigation potential and fast growth (Matthew, 2008).
In particular, GIS has been used extensively for spatial analysis and land suitability as GIS functions could be employed for several forms of information including point, line and area.The geospatial technology can be very efficiently used in analyzing the suitability of agricultural crops, which will not only save the huge cost of experimentation but also save precious time (Harasheh, 1994).Keeping this in view present study aimed to use GIS to classify the suitability of land with integrated information for Willow plantation.The spatial information resulted from this study could be utilized the wastelands for plantation which will not only help to increased income of villages but also minimize environmental problems.

Study area
The area of present study has been considered the state of Uttarakhand which is surrounded by Himachal Pradesh in the northwest and Uttar Pradesh in the south and shares its international borders with Nepal and China.This area is located between latitude 28º43'N and 31º27'N and longitude 77º34'E and 81º02'E.

Soil
The soil information and maps of the region were acquired from NBSSLUP, Nagpur.Entire state of Uttarakhand is covered in two maps at a scale of 1:500000.The maps were scanned using roller scanner and were mosaicked (combined) using the pixel base algorithm embedded in ENVI image processing software after rotating, resizing and contrast enhancement.The mosaicked map were georeferenced with existing georeferenced district boundary map using map to map registration option available in ENVI image processing software after collection of ample Ground Control Points (GCPs).The registered raster soil map was exported in imagine image (.img) compatible format so that it could be directly viewed in GIS software (ArcView 3.2a).The map was imported in ArcView software and soil polygons were digitized on line.The soil unit boundaries of the map were digitized very carefully using the polygon feature.Then the soil attribute data were added to the project.The soil attribute table

Fig. 1: Flow diagram of the method for plantation suitability assessment
comprised of different columns of fields like surface form, parent material, soil depth, mineralogy, particle size, calcareous, soil temperature, soil reaction, drainage, surface texture, slope class, erosion class, salinity, surface stone, flooding etc.

Climatic data
The data of all weather stations falling within the geographical boundaries of Uttarakhand were collected.Additionally, the data of meteorological observatories of adjoining states of Uttarakhand were also used in the present study.Data collected from different weather stations which located in Uttarakhand like VPKAS Almora, ARIES Nainital, DEBER Haldwani, CRC Pantnagar, College of Forestry and Hill Agriculture, Ranichauri, FRI Dehradun and data of some stations covered in UPROBE project of IIT, Roorkee were also considered, while some other stations, data were taken from published IMD periodicals.The complete list of stations of Uttarakhand with number of years, mean and CV% of available parameters has been appended in the Table .1,while the climate normal's of other stations computed by IMD on the basis of thirty years weather data and published in IMD periodicals have been presented in Table 2.The geocordinates of all the weather stations were collected recorded with the help of GPS or were collected from literature / published maps.Weather data with geocordiantes were stored in table and saved as a txt file.Text file was added in GIS environment and thereafter was added to GIS view using "add event theme"

Soil and climatic requirement of Willow
The optimum climatic requirement (i.e.maximum, minimum and optimum temperature, rainfall) and soil type for fast growing tree species are essential in order to check the spatial suitability of plant in a given region.This information was collected from the published literatures and is presented in Table 3.

Computation of PDI
We used the Precipitation Distribution Index developed by Nain et al. (2010).The precipitation distribution index (PDI) was calculated with the help of following formula: Where, PDI = Precipitation Distribution Index 1 = the value of a particular month when the ration between P/PET is >0.3, the value of month is 0 when the ratio between P/PET is < 0 P = Precipitation (mm) PET= Potential Evapotranspiration (mm) PET was calculated using Thornthwaite method.The formula was found suitable for calculation of PET by different researchers (  The unadjusted values of "e" are corrected for actual day light hours and days in a month.For daily computation, the formula is modified as: mm/day Where, k = adjustment factor for which table values are given by Michael (1978).
The PDI was computed on point data and later the spatial surface covering entire state of Uttarakhand was generated by developing a relationship between point precipitation and point PDI and later applying the model on the precipitation surface (Fig. 4).PDI surface was divided in 3 classes, which are more suitable, moderate suitable and less suitable (Table . 4).If PDI value is high then it suitable for planting and if it is low then it is less suitable.

Model for suitability analysis
The methodology of Suitability analysis included the digitization of the map.The soil unit boundaries of the map were digitized very carefully using the polygon feature.Then the soil attribute data were added to the project.The layers of different meteorological parameters were generated using the spatial analyst Extension of ARC-view.Ssuitability model was developed using model builder module available in Arc-View 3.2a GIS.The interpolated data related to climate and soil were= converted in raster format.The reclass function was added and reclasses of data was done.Then all these reclass data were combined to the weighted overlay function and a weighted overlay model was developed.After model was run and suitability map was generated.
Weightes given to the different parameters are shown in Table 5 on the basis of their effect on growth of tree.A value has been assigned to the different ranges of parameters on the basis of optimum, minimum and maximum range of parameter as mentioned in the Table 3.The area with value of parameter less than minimum value required by plant and higher value of the parameter than the maximum required by the plant were restricted as no growth of willow is possible in those regions.The optimum range has been assigned value of three (3), while suboptimum range has been assign value of two (2) and subsuboptimum range has been assigned value of one (1).A suitability map with 5 classes: most suitable, suitable, moderately suitable, less suitable and not suitable was generated and most suitable class was found when all parameters weightage values are highest and other classes, classified on basis of descending values.

Thematic map of soil
Loamy soils are found in most of the regions of Uttarakhand.However, sandy soils are limited to small pockets of Pithoragarh, Chamoli, Uttarkashi, Rudraprayag, Pauri, Tehri, Haridwar, Dehradun, Udham Singh Nagar, and Champawat and Bageshwar districts.A considerable region of Upper hills of Uttarakhand is snow bound and primarily has frigid soils.The spatial distribution of soil in Uttarakhand has been given in Fig. 2.

Thematic map of climatic parameter
Maximum temperature occurs in the month of June.However, a large area of the state has average maximum temperature of 16.4 oC.The low altitudinal areas of the state including Udham Singh Nagar, Haridwar and parts of Pauri and Nainital record maximum temperature in the range of 33.6 to 39.3ºC during June.Minimum temperature over the state prevails during January.During this month, upper hilly regions of the state experience minimum temperature of -31.5 to -6.7 oC.The low lying areas of the state including, Udham Singh Nagar, Haridwar and parts of Pauri, Dehradun and Nainital recorded a minimum temperature of 8.2 to 13.2oC during January.The upper hilly region records average temperature of -12.6 to 3.8 ºC.The upper part of Uttarkashi, Chamoli, Pithoragarh and a small portion of Tehri-Garhwal, Bageshwar and Rudraprayag constitute this region.However, a large area of the state has an average temperature range of 19.9 to 25.3 ºC.In districts like Udham Singh Nagar, Haridwar and some parts of Nainital, Pauri-Garhwal and Tehri-Garhwal, situated in low altitudinal areas, average temperature ranges between 25.3 to 30.7 ºC (Fig. 3).The maximum amount of rainfall occurs during June to September over Uttarakhand due to southwest monsoon.The precipitation reached its peak in the month of August.In some part of Dehradun precipitation is received in the range of 922 to1035 mm during August.In this month the upper hilly region receives rainfall of 312 to 612 mm.Soil Texture 10 In the state, the lowest precipitation occurs in the month of November.

Thematic map of precipitation distribution index (PDI)
The values of PDI for each station were calculated and three distinct regions were identified from the map (Fig. 4 According to the PDI surface almost entire region of Uttarakhand has the capability to support the growth of fast growing tree species because most part of the state comes under the high and moderately high PDI values except upper hilly regions which has low PDI values and Skelton and frozen soil surface.High and medium water demanding tree species can be grown in high and moderately high PDI values region, while low water demanding tree species in low PDI region.

Suitability classes of Willow
By using the all parametric information (precipitation, PDI, soil types, average, minimum and maximum temperature) crop suitability model was developed.The model was used to produce the suitability map of willow (Fig. 5) for Uttarakhand.

Most suitable and suitable zones
Suitable area covers 15654 km 2 , which has 29.3% of total geographical area of the state.The zone includes parts of U. S. Nagar, Haridwar, Tehri Garhwal, Almora, Champawat, Chamoli, Pithoragarh, Uttarkashi and some part of Dehradun, Bageshwar, Nainital and Pauri Garhwal.Most suitable areas for willow have estimated at 20343 km 2 (38%), which includes most portions of Dehradun, Bageshwar, Almora, Pauri Garhwal, Tehri Garhwal, Champawat, Udham Singh Nagar, Pithoragarh and Some part of Uttarkashi and very small portion of Haridwar.

Less suitable zones
There was a very weak relation between climatic conditions and requirements of Willow in these areas compared to the first zone.In all climatic requirements was not favourable, so these zones are less suitable for Willow.The area under the zone has 7035 km 2 or 13.2%, which includes central and some upper parts of Chamoli, Pithoragarh, Tehri Garhwal, Uttarkashi, Bageshwar, Rudraprayag and some parts of Dehradun districts.

Not suitable zones
This zone has considered as areas that are not suitable for Willow because of noncompatibility of climatic and soil conditions in respect to requirement of Willow.The most important factors that limits planting of Willow in these areas have minimum temperature (5-10 0 C) and soil type (Skelton surface).The zone includes upper most hilly region of Uttarkashi, Rudraprayag, Chamoli, Pithoragarh, Bageshwar and Tehri Garhwal.This zone has 10453 km 2 (19.5%) area.
Willow suitability map shows that most of the central parts are most suitable, lower part suitable and less suitable area falls in middle and lower Himalayan hilly region of the state.More weight to temperature has been assign on the basis of literature.This map indicate that some region of temperate areas fall under the suitable and moderate suitable zone.It has possible that some area shows suitability where willow currently not found but it can florish because the climatic conditions have favorable for its growth.Willow has 400 species and climatic requirements (like precipitation, elevation maximum, minimum and optimum temperature) are difference, hence suitable species can be selected for a particular place.However due to adverse climatic conditions in less suitable areas, willow plantations can be avoided.
Fischer et al. (2005) used GIS database of climate, soil, terrain and vegetation cover to assess the potential productivity for miscanthus, willow and poplars in countries of Eastern Europe and Northern and Central Asia.The result shows that, the potential for bio-energy production (miscanthus, willow and poplar), both absolute and percentage areas suitable for poplar and willow alone and by poplar, willow and miscanthus combined.In a few countries such as Lithuania, Latvia, Romania, Georgia, Belarus and Azerbaijan, the potential for producing energy from these species amounts to more than one-third of current commercial energy use even when current forests and land potentially highly suitable for cereals are excluded from the assessment.

CONCLUSION
The present study demonstrates use the available geo-spatial information for the identification of suitable areas for Willow plantations in Uttarakhand.According to the results of the study, it can be concluded that the area classified as most suitability, suitability, less suitability and not suitability occupied 38%, 26.4%, 13.2% and 19.5% of the total geographical area respectively.Uttarakhand has wastelands on considerable geographical area, which is 26% of the total geographical area of state.Thus, it can be recommended that on the basis of the suitability map; the wastelands can be used for Willow plantation.It will not only be beneficial for socio-economic point of view but also enhance the carbon sequestration potential of the land.

Fig. 2 : 3 :Fig. 4 : 5 :
Fig. 2: Soil Surface texture of Uttarakhand Fig. 3: Thematic map of average temperature over Uttarakhand Temperature corrected in °C T int = Temperature interpolated in °C.Alt act = Actual altitude in m Alt int = Interpolated altitude in m 2008) in the various parts of the world.Thornthwaite proposed the following formula for monthly potential evapotranspiration: unadjusted potential evapo-transpiration (cm/ month) (month of 30 days each and 12 hours day time).t= mean air temperature (ºC) I= annual or seasonal heat index, the summation of 12 values of monthly heat indices (i) when, an empirical exponent computed by the equation, a = 0 .0 0 0 0 0 0 6 7 5 1 * ( I 3 ) -0 .0 0 0 0 7 7 1 * ( I 2 ) +0.01792*I+0.49239 option.The "spatial extension" was loaded in order to use the interpolation functions available in GIS environment.Inverse Distance Weightage (IDW) interpolation technique was used to construct the thematic layers of the different weather parameters.Interpolated weather data provides values of weather parameters on flat surface, however the topography of Uttarakhand is highly variables.Therefore weather surface was corrected using Digital Elevation Model (DEM) of Uttarakhand.DEM provides a digital representation of a portion of the earth's surface terrain over a two dimensional surface.The DEM was used for constructing temperature surfaces keeping in view the environmental lapse rate.The environmental lapse rate which describes that with increasing elevation by 1 km the temperature will decrease by 6.5ºC was used to correct the interpolated temperature surface using following equation.

Table 1 : Climatic normals with no. of years and CV% of different stations of Uttarakhand Station Max. Temp. (C) Min. Temp. (C) Aver. Temp. (C) Rainfall (mm) No. of
Michael,