Integration of Hydrogeological Factors for Identification of Groundwater Potential Zones Using Remot...

Lazarus G. Ndatuwong, G. S. Yadav

  Open Access OPEN ACCESS  Peer Reviewed PEER-REVIEWED

Integration of Hydrogeological Factors for Identification of Groundwater Potential Zones Using Remote Sensing and GIS Techniques

Lazarus G. Ndatuwong1, 2,, G. S. Yadav1

1Department of Geophysics, Banaras Hindu University, Varanasi, India

2Department of Physics, Adamawa State University, Mubi, Nigeria

Abstract

Remote sensing data can be used as a reconnaissance and features identification tool for identifying surface and sub-surface water potential zone. The present study has been carried out to evaluate the potential zones for groundwater targeting using an integrated remote sensing data, Survey of India (SOI) topographical sheets and field verification. Four features (geomorphologic units, slope, drainage density and lineaments density) that influence groundwater occurrence were extracted and integrated to evaluate the hydrogeomorphological characteristics of the study area and demarcate the groundwater potential zones. Thematic maps of the extracted features were prepared and integrated through geography information system (GIS) environment. The groundwater potential map was prepared by overlaying the thematic layers. Weightage percentages were assigned to the different parameters according to their relative importance to groundwater potentiality. The integrated map of the area shows different zones of groundwater prospects, viz. very high (0.77% of the area), high (35.57% of the area), moderate (54.53% of the area), while poor and very poor are made up of 9.13% of the area.

At a glance: Figures

Cite this article:

  • Ndatuwong, Lazarus G., and G. S. Yadav. "Integration of Hydrogeological Factors for Identification of Groundwater Potential Zones Using Remote Sensing and GIS Techniques." Journal of Geosciences and Geomatics 2.1 (2014): 11-16.
  • Ndatuwong, L. G. , & Yadav, G. S. (2014). Integration of Hydrogeological Factors for Identification of Groundwater Potential Zones Using Remote Sensing and GIS Techniques. Journal of Geosciences and Geomatics, 2(1), 11-16.
  • Ndatuwong, Lazarus G., and G. S. Yadav. "Integration of Hydrogeological Factors for Identification of Groundwater Potential Zones Using Remote Sensing and GIS Techniques." Journal of Geosciences and Geomatics 2, no. 1 (2014): 11-16.

Import into BibTeX Import into EndNote Import into RefMan Import into RefWorks

1. Introduction

Groundwater is a form of water held under the ground in the saturated zone that fills all the pore space of soils and geologic formations. It is formed by rainwater or snowmelt water that seeps down through the soil and into the underlying rocks [1]. It is the major resource of water supply for about half of the nations. It plays a key role in Nature. It provides more than half of humanity’s freshwater for everyday uses such as drinking, cooking, and hygiene, as well as thirty percent of irrigated agriculture and industrial development [2]. Groundwater potential zones can be said to be water bearing formation of the earth crust that acts as conduits for transmission and as reservoirs for storing water. Its identification and location is based on indirect analysis of some observable terrain features such as geologic, geomorphic, landforms and their hydrologic characteristics.

Conventional approach of groundwater exploration using geological, hydrogeological, and geophysical methods are expensive due to high cost of drilling and time consuming investigation [3, 4]. Furthermore, these methods of surveys do not always account for the different factors that control the occurrence and movement of groundwater [5]. Remote sensing technique is popular due to several advantages of spatial and spectral data, having access to large coverage and inaccessible areas with regular revisit capability [6]. It is also rapid and cost-effective tool that integrates valuable data on geology, geomorphology, lineaments, slope, morphometric etc. that helps in evaluating groundwater potential zone.

The use of remote sensing and GIS techniques is a fast emerging field in groundwater resource identification, mapping and planning. Remote sensing provides an opportunity for better observation and more systematic analysis of various geomorphic units, landforms, lineaments and drainages, due to its synoptic and multi-spectral coverage of a terrain [7]. The techniques is increasingly used in prospecting for groundwater potential zones, because of their ability to identify and outline different ground features that may serve as direct or indirect indicators of the presence of groundwater resources [8].

Reviewed papers on the application of remote sensing and GIS in groundwater potential mapping can be found in [9, 10, 11, 12, 13] More recent studies have been carried out on the use of remote sensing data and GIS tools for defining the spatial distribution of different groundwater prospecting classes on the basis of hydrogeomorphological units and other associated parameters [14-20][14]. In this study the technique of remote sensing data and GIS were used to evaluate the hydrogeomorphological characteristics of the study area and demarcate the groundwater potential zones for investigation and exploration of the resource in shallow aquifer using thematic maps of geomorphologic units, slope, drainage density and lineaments density. Result obtained from this analysis will help to narrow down the targeting area for further geophysical exploration for the resource.

2. Study Area

Geologically the study area is situated in the Vindhyan basin in Central India. It falls within the survey of India toposheets 63L/9, 63L/10, 63L/13, 63L/14 located in the southern part of Uttar Pradesh state, bounded by 82° 40’ – 83° 00’ east of longitude and 24° 40’ – 24° 50’ north of latitude (Figure 1). The relief ranges from 250 m over the plains to the highest elevation of 389 m. It is underlain by the Kaimur Group of the Vindhyan Supergroup. The series comprises of sandstones and quartzites which are prominent horizons with wide distribution. All the major sandstone horizons form the scrap while the shale horizons form gentle slopes [21]. The area is covered by thin surface soil with varying thickness ranging from 1-4m and few exposures of sandstones rocks [22]. The bedrock (Vindhyan Super Group of rocks) is expected below the surface soil cover since the study area is lying in close vicinity of the Vindhyan exposures. The area is mainly flat and gently undulating terrain except in few parts. The occurrence and movement of groundwater is mainly restricted within the weathered and fractured sandstone/shale [23]. The amount of water that can be extracted from the fractured zones depends on the size and location of fractures [24].

The climatic condition comprises of a mean temperature ranging from 4.4°C in January to above 46.6°C in May. The rainfall is uncertain and erratic with about 84 percent of the annual rainfall in the monsoon period, with an annual average of about 1,530 mm. The area is covered by the Kaimur hills which constitute the greater part of the tableland. The Soils of the area are excessively drained, highly permeable and have low to medium productivity. The red soil is comprised mainly of red loam and is very easy for cultivation purposes and responds to good manure and other treatments. This soil is particularly suited for growing vegetables. Patches of laterite soil are found all over the area.

Figure 1. Location and topographic map of the study area

3. Material and Methods

Remote sensing data which comprises of 90 m interval Shuttle Radar Topographic Mission (SRTM), Digital Elevation Model (DEM) and Landsat image along with Survey of India (SOI) toposheets and field observation data were used in this study.

The methodology includes generation of thematic maps showing drainage pattern, lineament, slope and geomorphologicic units through processing and visual interpretations of remote sensed data. The identification and delineation of various units on the thematic maps was based on tone, texture, shape, colour etc. Thematic layers of geomorphologic units, slope, lineament and drainage density were generated and integrated in GIS environment to determine suitable zones for groundwater prospecting.

4. Hydrogeomorphological Analysis

The relief, slope, depth of weathering, type of weathered material, thickness of alluvium, nature of the deposited material and the overall assemblage of different landforms play an important role in defining the ground water regime, especially in the hard rocks and the unconsolidated sediments [7]. Hydrogeomorphological investigations include the delineation and mapping of various landforms, drainage characteristics and structural features that could have a direct control on the occurrence and flow of groundwater. Many of these features are favourable for the occurrence of groundwater and are classified in terms of groundwater potentiality as poor, moderate, good and very good prospecting zones [25, 26, 27].

5. Result and Discussion

The result of the study have brought out five distinct geomorphic units (Figure 2), which are discussed systematically here with respect to groundwater prospect;

5.1. Piedmont Plain

Piedmont plain are gently sloping longitudinal strip of land running parallel to foot hills traversed by innumerable rivulets with parallel to sub-parallel drainage; comprising of unconsolidated sediments of sand, silt and clay with boulders and pebbles. This unit covers a total area of 39.43 km2, which is about 8.5 % of the study area and are mainly found in the southern part of the study area at an elevation of above 313 m above msl. The groundwater prospect of this unit is generally poor, but moderate prospect can be expected where fracture/lineament are located.

5.2. Pediments

Pediments are gently sloping rock flooring areas with erosional bedrock of low relief between hills and plains developed by the process of weathering and consisting of a veneer of detritus, alluvial sand and gravel that thickens downslope. They show massive and compact rocky nature surfaces, with numerous fractures/joints permitting infiltration and storage of groundwater [26]. In the continuous process of pedimentation, when pediment gets buried by detritus and regolithic cover, it is termed as a buried pediment in which the sub-surface rock subsequently undergoes weathering. Groundwater prospecting in pediments can be considered as moderate to poor, but the presence of lineaments or fractures can provide some good scope for groundwater prospecting. This unit occupies 25.9 % of the study area with a total area of about 129.21 km2 found in the southern part with some patches towards the north that flanks from west to east of the study area.

5.3. Alluvial Plain

This is a gently sloping plain that consist of thick unconsolidated sediments comprising mostly of loamy soil formed by fluvial action. It shows smooth texture and extensive agricultural activity. It is classified as shallow and deep alluvium plain depending on the thickness of the unconsolidated material. It acts as a good recharge zone and as such represent a good to very good potential zone for groundwater exploration depending on the thickness of the alluvium and the weathered column underlying the formation.

The unit occupies an area of about 157.82 km2 which represent 34.0 % of the total study area. It mainly occupied the central portion of the map

5.4. Flood Plain

This consists of unconsolidated sediments of sand, silt, clay with kankar layers mixed with cobbles and pebbles, gravels and alluvial materials derived from surrounding uplands and deposited by fluvial action along the floor of river channels. It is highly permeable zone helping the partial bank recharge, while the flow of subsurface groundwater occurs under semi-confine to perched water table condition [28]. The unit occupies the central part of the survey area just along the course of the main river (Bellan Nadi) that cuts across the study area and some occupying an area up north. The unit covers a total area of about 118.83 km2 which is about 25.6% of the study area and is considered as very good groundwater potential zone.

5.5. Pediplain

This is a generally flat and smooth surface of weathered pediplain of granite gneiss weathered material usually covered with red soil. In hard rocks, this landform forms very good recharge and storage zones depending upon the thickness of the weathering/accumulated material and its composition [29]. Moderate to good yields are expected in this geomorphic unit. Faults/ fracture zones passing through this unit act as conduits for movement and occurrence of ground water. Good yields are expected along these fractures/lineaments. The unit is made of 5.9 % of the total study area covering an area of about 27.44 km2 at the north-western end of the map.

6. Lineaments Density Map

A lineament is usually defined as a straight or somewhat-curved feature in an image. In a satellite image, lineaments can be the result of man-made structures such as transportation networks (roads, canals, etc.), or natural structures such as geological structures (faults/fractures, lithological boundaries, unconformities) or drainage networks (rivers). Generally lineaments are underlain by zones of localized weathering and increased permeability and porosity. Previous studies have revealed a close relationship between lineaments and groundwater flow and yield [30, 31]. Meanwhile, some researchers studied relationships between groundwater productivity and the number of lineament within specifically designated areas or lineament density rather than the lineament itself [32]. Therefore mapping of lineaments closely related to groundwater occurrence and yield is essential to groundwater surveys, development, and management.

Lineaments are usually extracted and interpreted from satellite imagery manually or automatically. PCI-Geomatica Version 10 software was used in this work to automatically extra lineament from satellite image of the survey area. The extracted linear features were further screen for non-geological features such as roads, fences, field boundaries, by comparing the lineament map with the corresponding toposheet map of the survey area and field verification, thereby deleting non-geological features and leaving only possible geological lineament (Figure 3a). Lineament density map was prepared by dividing the study area into 1 km/1 km grids. The total length of lineament in each grid was measured and plotted in the respective grid centres. These values are joined by isolines to prepare a lineament density map as shown in Figure 3b. The extension of large lineaments representing a shear zone or a major fault and can extend from hilly terrain to alluvial terrain [27]. It may form a productive groundwater reserve. Similarly intersection of lineaments can also be probable sites of groundwater accumulation. Therefore, areas with high lineament density may have important groundwater prospects even in hilly regions which otherwise have nil groundwater prospects than area with low density.

7. Slope

Slope analysis is an important parameter in geomorphic studies. The slope elements, in turn are controlled by the climatomorphogenic processes in the area having the rock of varying resistance [33]. Slope plays a key role in groundwater occurrence as infiltration is inversely related to slope [34]. Slope map of the survey area was generated from DEM data using spatial analyst tool in ArcGIS. The slope varies between zero and 24 degree (Figure 4). Area of high slope value will cause more runoff and less infiltration thus, have poor groundwater prospect compared to low slope region.

8. Drainage Density Map

The drainage is expressed as density map for easy characterization of runoff and infiltration of rainwater in the survey area. The drainage density map was prepared in similar way as the lineament density map by dividing the study area into 1 km/1 km grids. The total length of stream in each grid was measured and plotted in the respective grid centres. These values are joined by isolines and the resultant density map was plotted in arcgis (Figure 5).

Drainage density is an inverse function of permeability, and therefore it is an important parameter in evaluating the groundwater zone [34]. Area of high drainage density indicates less infiltration which favours runoff and hence acts as poor groundwater prospect, because major part of the rainwater over the area is lost as surface runoff with little infiltration for recharging the groundwater reservoir. On the other hand low drainage density areas permit more infiltration and recharge to the groundwater reservoir, hence can be described as a good potential zone for groundwater prospect.

9. Mapping Groundwater Potential Map

In this study the weighted index overlay analysis (WIOA) using a multi-criteria approach was used to investigate the potential zone for groundwater prospecting through integration of hydrogeomorphological map, slope, lineament and drainage density map. All the thematic maps generated for the study area were integrated and projected to a common co-ordinate system through spatial analysis in ArcGis and reclassified for weighted index overlay analysis (WIOA). The weighted index overlay analysis is a technique that applies a common measurement scale of values to diverse and dissimilar inputs to create an integrated analysis using a multi-criteria approach. This analysis according to [35] has no standard scale, but incorporates human judgment for its efficiency. For this work, each of the thematic maps was assigned a weight that represents its importance in respect to groundwater reserve based on criteria of previous work [35, 36, 37].

Theme weight and class rank assigned to the different parameters considered for groundwater prospects evaluation in this study are given in Table 1.

Table 1. Weightage and class rank of different parameters used for groundwater prospects

The groundwater potential map of the study area is shown in Figure 6. The groundwater potential zones were grouped into five different potential zones viz; very high, high, moderate, poor and very poor. Analysis of the potential zones shows that the very high groundwater potential zones constitute just less than 1 % of the study area. Few patches of this zone were observed at areas associated flood plain geomorphologic unit in the central and north-eastern end of the map. High groundwater potential zones are seen at the central towards the northern part and constitute 35.51 % of the study area. This zone is associated with geomorphologic flood and alluvial units and part of the pediments with high lineament density in the southern part f the map. Moderate groundwater potential zones occupy the largest area of about 54.53 % of the total study area, scattered all over the map. The poor and very poor potential zone constitute about 9.13 % of the study area and are mostly seen along denudational hills with steep slope and part of the alluvial and pediments units that are associated with high drainage density.

Figure 6. Map of groundwater potential zone.

10. Conclusion

The study area was able to reveal the use of remotely sensed data in identifying favourable zones for groundwater exploration in hard rock terrain. Groundwater prospects zonation based on this study clearly indicated that the combination of flood plain, alluvial plains, lineament and gentle slope areas are the favourable terrain conditions having good groundwater potential. Pediplain zones have much better prospects as compared to pediments and piedmonts regions due to their colluvial material composition, moderate drainage density, moderate slope and the existence of lineaments/fractures.

References

[1]  Banks, D., Robins, N., “An introduction to groundwater in crystalline bedrock”. Norges geologiske undersøkelse, Trondheim, p 64. 2002
In article      
 
[2]  Zuppi, G.M., “The groundwater challenge. In: Clinic C, Musu I and Gullino M.L (eds). Sustainable development and environmental management. Experience and case studies”, Dordrecht, Springer, 49-64. 2007.
In article      
 
[3]  Singh, A.K. and Prakash, S.R., “An integrated approach of remote sensing, geophysics and GIS to evaluation of groundwater potentiality of Ojhala subwatershed, Mirzapur district, UP, India”, Map India conference. 2003.
In article      
 
[4]  Sander, P., Chesley, M.M., Minor, T.B.,“Groundwater assessment using remote sensing and GIS in a rural groundwater project in Ghana: lessons learned”, Hydrogeol J., 4(3). 40-49. 1996.
In article      CrossRef
 
[5]  Oh, H.J., Kim, Y.S., Choi, J.K., Park, E. and Lee, S., GIS mapping of regional probabilistic groundwater potential in the area of Pohang City, Korea” J. Hydrol. 399, 158-172. 2011.
In article      CrossRef
 
[6]  Teeuw, R.M., “Groundwater exploration using remote sensing and a low cost GIS”, Hydrogeol J;3(3):21-30. 1999.
In article      CrossRef
 
[7]  Obi Reddy, G.P., Chandra Mouli, K., Srivastava, S.K., Srinivas, C.V. and Maji, A.K., “Evaluation of ground water potential zones using remote sensing data—a case study of Gaimukh watershed, Bhandra district, Maharashtra”, J Indian Soc. Remote Sens., 28. 19-32. 2000.
In article      
 
[8]  Gopinath, G. and E Seralathan, E., “Identification of groundwater prospective zones using IRS-ID LISS III and pump test methods”, J. Indian Soc. Remote Sensing, 32(4), 329-342. 2004.
In article      CrossRef
 
[9]  Meijerink, A.M.J., “Groundwater. In: Schultz, G.A, Engman ET (eds) Remote sensing in hydrology and water management” Springer, Berlin, pp 305-325. 2000.
In article      CrossRef
 
[10]  Hoffman, J. and Sander, P., “Remote sensing and GIS in hydrogeology”, Hydrogeol J., 15(1). 2007.
In article      CrossRef
 
[11]  Jha, M.K., Chowdhury, A. Chowdary, V.M. and Peiffer, S., “Groundwater management and development by integrated remote sensing and geographic information systems: prospects and constraints”, Int Ser Prog Wat Res., 21. 427-467. 2007.
In article      
 
[12]  Jasmin, I. and Mallikarjuna, P., “Review: satellite-based remote sensing and geographic information systems and their application in the assessment of groundwater potential, with particular reference to India” Hydrogeol J., 19(4). 729-740. 2011.
In article      CrossRef
 
[13]  Teixeira, J., Chamin, H.I., Carvalho, J.M., Prez-Alberti, A. and Rocha, F., “Hydrogeomorphological mapping as a tool in groundwater exploration” Journal of Maps, 9(2). 263-273. 2013.
In article      CrossRef
 
[14]  Elmahdy, S.I. and Mohamed, M.M., “Groundwater potential modelling using remote sensing and GIS: a case study of the Al Dhaid area, United Arab Emirates”, Geocarto International. Geocarto International, 2013
In article      
 
[15]  Nag, S.K. and Ghosh, P., “Delineation of groundwater potential zone in Chhatna Block, Bankura District, West Bengal, India using remote sensing and GIS techniques”, Environ Earth Sci., 2012.
In article      
 
[16]  Magesh, N.S., Chandrasekar, N., Soundranayagam, J.P., “Morphometric evaluation of Papanasam and Manimuthar watersheds, parts of Western Ghats, Tirunelveli district, Tamil Nadu, India: a GIS approach” Environ Earth Sci., 64, 373-381.2011.
In article      
 
[17]  Rashid, M., Ahmad Lone, M. and Ahmed, S., “Integrating geospatial and ground geophysical information as guidelines for groundwater potential zones in hard rock terrains of south India”, Environ Monit Assess., 2011.
In article      
 
[18]  Javeed, A. and Wani, M.H., “Delineation of Groundwater Potential Zones in Kakund Watershed, Eastern Rajasthan, Using Remote Sensing and GIS Techniques”, Journal Geological Society of India, 73. 229-236. 2009.
In article      CrossRef
 
[19]  Rao, P.J., Harikrishna, P., Srivastav, S.K., Satyanarayana. P.V.V. and Vasu Deva Rao, B., “Selection of groundwater potential zones in and around Madhurawada Dome, Visakhapatnam District - A GIS approach”, J. Ind. Geophys. Union, 13(4). 191-200. 2009.
In article      
 
[20]  Bose P K, Sarkar S, Chakrabarty S, and Banerjee S., (2001), Overview of Meso- to Neoproterozoic evolution of the Vindhyan basin, Central India. J. Sediment. Geol., 142, pp 395-419.
In article      CrossRef
 
[21]  Krishnan M S.,” Geology of India and Burma”. 6th Ed. CBS publisher, New Delhi, 1982.
In article      
 
[22]  Amaresh K. S & Prakash S. R.., “An integreted approach of Remote Sensing, Geophysics and GIS to evaluation of Groundwater potentiality of Ojhala subwatershed, Mirzapur district, U.P., India”, Map India Conference, 2003, Available at http:// www.GISdevelopment.net.
In article      
 
[23]  L. G. Ndatuwong and G. S. Yadav., “ Quantitative estimation of natural recharge due to monsoon rainfall using the principle of information theory in the area of Ghorawal block of Sonebhadra district, U.P., India” International Journal of Environmental Sciences, 3, (3), 976-985, 2012.
In article      
 
[24]  Ganapuram, S., Vijaya Kumar, G.T., Murali Krishna, I.V., Ercan K. and Demirel, C.M., “Mapping of groundwater potential zones in the Musi basin using remote sensing data and GIS”, Advances in Engineering Software, 40. 506-518. 2009.
In article      CrossRef
 
[25]  NRSA., National rural drinking water mission methodology manual for preparation of groundwater prospects maps. National Remote Sensing Agency, Govt. of India, Hyderabad, 1999
In article      
 
[26]  Deepika, B., Avinash, K. and Jayappa, K.S., “Integration of hydrological factors and demarcation of groundwater prospect zones: insights from remote sensing and GIS techniques”, Environ Earth Sci., 2013.
In article      
 
[27]  Mondal, Md.S., Pandey, A.C. and Garg, R.D., “Groundwater Prospects Evaluation Based on Hydrogeomorphological Mapping using High Resolution Satellite Images: A Case study in Uttarakhand”, J. Indian Soc. Remote Sens., 36. 69-76. 2007.
In article      CrossRef
 
[28]  Sankar, K.,“Evaluation of groundwater potential zones using remote sensing data in Upper Vaigai river basin, Tamil Nadu, India”, J Indian Soc Remote Sensing, 30(30). 119-129. 2002.
In article      CrossRef
 
[29]  Raghu, V. and Reddy, K.M., “Hydrogeomorphological Mapping at Village Level Using High Resolution Satellite Data and Impact Analysis of Check Dams in Part of Akuledu Vanka Watershed, Anantapur District, Andhra Pradesh”, J. Ind. Geophys. Union, 15(1). 1-8. 2011.
In article      
 
[30]  Magowe, M., Carr, J.R., “Relationship between lineaments and ground water occurrence in western Botswana”, Ground Water, 37 (2). 282-286. 1999.
In article      CrossRef
 
[31]  Fernandes, A. and Rudolph, D., “The influence of Cenozoic tectonics on the groundwater production capacity of fractured zones: a case study in Sao Paulo, Brazil” Hydrogeol J., 9. 151-167. 2001.
In article      CrossRef
 
[32]  Hardcastle, K., “Photolineament factor: a new computer-aided method for remotely sensing the degree to which bedrock is fractured. Photogramm” Eng Remote Sensing, 61(6). 739-747. 1995.
In article      
 
[33]  Sreedevi, P.D., Owais, S., Khan, H.H. and Ahmed, S., “Morphometric Analysis of a Watershed of South India Using SRTM Data and GIS”, J. Geological Society of India 73, 543-552. 2009.
In article      CrossRef
 
[34]  Agarwal, E., Agarwal, R., Garg, R.D. and Garg, P.K., 2013. Delineation of groundwater potential zone: An AHP/ANP approach. J. Earth Syst. Sci. 122(3), 87-898.
In article      CrossRef
 
[35]  Nag, S.K., “Application of lineament density and hydrogeomorphology to delineate groundwater potential zones of Baghmundi block in Prulia district, West Bengal”, Jour. Indian Soc. Remote Sensing, 33, 521-529. 2005.
In article      CrossRef
 
[36]  Narendra, K., Nageswara Rao, K.N. and Swarna latha, P. “Integrating Remote Sensing and GIS for Identification of Groundwater Prospective Zones in the Narava Basin, Visakhapatnam Region, Andhra Pradesh”, Journal Geological Society of India, 81. 248-260. 2013.
In article      CrossRef
 
[37]  IMSD., 1995 Integrated mission for sustainable development technical guidelines. National Remote Sensing Agency, Department of Space, Govt. of India.
In article      
 
  • CiteULikeCiteULike
  • Digg ThisDigg
  • MendeleyMendeley
  • RedditReddit
  • Google+Google+
  • StumbleUponStumbleUpon
  • Add to DeliciousDelicious
  • FacebookFacebook
  • TwitterTwitter
  • LinkedInLinkedIn