Development Team: Paper No: 6
Development Team: Paper No: 6
Development Team: Paper No: 6
Module: 31 Applications of Remote Sensing and GIS in Groundwater and Water Pollution
Studies
Development Team
Principal Investigator Prof. R.K. Kohli
& Prof. V.K. Garg & Prof. Ashok Dhawan
Co- Principal Investigator Central University of Punjab, Bathinda
Dr. Puneeta Pandey
Assistant Professor
Paper Coordinator Centre for Environmental Sciences and Technology
Central University of Punjab, Bathinda
Dr. Jitendra Kumar Pattanaik
Content Writer Central University of Punjab, Bathinda
Content Reviewer Dr. Puneeta Pandey
Central University of Punjab, Bathinda
Module Name/Title Applications of Remote Sensing and GIS in Groundwater and Water Pollution Studies
Module Id EVS/RSGIS-EVS/31
Pre-requisites Basic knowledge of school- level physics and fundamentals of remote sensing and GIS
CONTENTS
1. Aim of the Module
2. Introduction
3. Remote sensing and GIS techniques for Groundwater studies
4. Sensing of Ground Water Fluxes
5. Ground water and Land Surface
6. Image interpretation for ground water study
7. Remote sensing and GIS in water pollution studies
8. Conclusions
9. References
2. Introduction
This module is divided into two parts. The first part deals with the application of remote sensing and
geographic information system (GIS) techniques for the groundwater studies and the second part
(section 6) explains the application on water pollution studies.
Groundwater is a vital natural resource for mankind. This resource is extensively being used for
drinking and household utilization; irrigation and industrial purpose. It plays an important role for the
economic development and food security of the country. Only 2.5% of the earth’s water is available as
freshwater. Groundwater comprises nearly 30% of freshwater resource of the world whereas glaciers
and ice cap consist of 68.7% (http://water.usgs.gov/edu/earthwherewater.html) which is difficult to
3
Based on the topographic driving force Tόth (1963) conceptualized a model (Fig. 1) of
groundwater flow system for local and regional scale. This model shows that the ground water
recharged at higher elevation in the regional scale tend to move deeper compare to local scale
recharge. So based on the topographical information from remote sensing data predicting local or
regional scale groundwater flow will be more effective (Becker, 2005). The rate and behavior of flow
depend up on the geology and it can be expressed by the Darcy’s Law. Darcy’s law defines the flow of
fluid in a porous medium and also states that there is a linear relationship between flow velocity and
hydraulic gradient (I) for any given saturated soil or medium under steady laminar flow conditions. It
can be stated as
q = K. I
Where ‘q’ is the specific discharge vector representing flow per unit area (flux of ground
water), ‘K’ is the hydraulic conductivity which is a function of geology, ‘I’ is the hydraulic gradient
which is a function of surface forcing (Becker, 2005). Geological maps prepared in combination with
remotely sensing data and ground verification provides useful information about the hydraulic
conductivity, water bearing formations, lineaments such as faults, fractures in the hard-rock terrain.
This information is used for groundwater prospecting as evident from literature.
For preparation of groundwater model of an area, surface water treated as the boundary
conditions for the subsurface flow equation which is based Darcy’s law. Remotely sensed imageries
are used to define boundary conditions such as streams, lakes, wetlands, seepage areas, recharge
zones, or evapotranspiration zones for prediction of ground water flow. The important mathematical
boundary conditions are hydraulic head (Fig. 2), flux or discharge (Fig. 3), mixed (both head and
discharge). So the remote sensing applications in ground water studies can be structured into the
sensing of hydraulic potential (heads) and hydraulic flux (or discharge).
3.2 Surface Water Elevations: Generally the elevation of surface water depicts the possible
groundwater head of that region. Therefore the spring or first-order stream originates at an elevation
where water table intersects the slope. In the catchment scale it provides an opportunity for dynamic
monitoring of water table. Sometime it may be difficult if hydraulic conductivity changes
considerably. Use of satellite based altimetry and interferometry for obtaining surface water elevations
provides higher accuracy (Becker, 2005) compare to digital elevation models (DEMs) and topographic
digital line graphs (DLGs).
Figure 1: Conceptual model (Tόth) of topographically driven ground water flow systems (after
Becker, 2005; Fetter 2001).
3.3 Water Column Mass: Water storage in the subsurface and hydraulic head in an aquifer can be
estimated using satellite or aerial gravity surveys. This method is useful for studying very large aquifer
system due to very coarse spatial resolution. Additionally this method does not have the vertical
resolution; hence the influences from water present in the atmosphere and vegetation, and unsaturated
water content i.e. soil moisture has to be removed to enhance the accuracy in estimates of saturated
water mass (Becker, 2005). Determination of ground water potential gradient is difficult using this
method because of coarse spatial resolution which requires aquifer to be continuous over hundreds of
kilometer. Data obtained from the NASA Gravity Recovery and Climate Experiment (GRACE)
satellite proved to be an asset for estimating ground water storage.
3.4 Heat Capacity: Heat capacity of saturated soil is higher than the dry soil. Using this property of
soil, depth of water table can be estimated from remotely sensed thermal image. This technique is very
useful to locate shallow water tables and this was proposed by Cartwright (1968) and Chase (1969) in
the early days of remote sensing application. Some researchers have found that night time thermal
images are more useful to predict depth of shallow water table compare to the day time thermal image.
Annual variation in soil temperature like heat sink during summer and heat source in the winter should
7
3.5 Land Subsidence: Ground water generally occupies the pore space of the sediments. In case of
unconsolidated sediment, addition (recharge) or withdrawal (depletion) of water from the pore space
will change the net volume. During recharge effective pressure will be high in the pore space hence it
will increase the volume as well as water level (Becker, 2005). Withdrawal of water will cause
reduction in the pore pressure which leads to decrease in volume and compaction of unconsolidated
sediment. This will result in land subsidence. This volume change will be reflected as variation in the
surface elevation. Although this variation is small, it can be measured by interferometric synthetic
aperture radar (InSAR). For this analysis image of a location is taken from different angle and time.
Accuracy of elevation change estimated by In SAR analysis is control by topography and
concentration of water vapor in the atmosphere. This accuracy is different for humid (10cm) and dry
(1mm) region (Galloway et al 1998). Using surface elevation change storativity of a porous medium
can also be estimated with the help of other data such as geodetic controls from GPS, water level,
hydrological flux and strain measurement of the study area. Integration of numeric model with InSAR
analysis widens the application of this method and provides better resolution and spatial extent of land
subsidence than ground base measurement (Becker, 2005).
3.6 Soil Moisture: Presence of shallow water table can be predicated based on the soil moisture
content. Different remote sensing method has been applied extensively to delineate shallow water
table from soil moisture content and vegetation stress or proliferation (Becker, 2005). Visible and
near-infrared sensor is also used to monitor the change of vegetation cover/agricultural performance
which can be linked with the water logging or change in the soil moisture content. Passive and active
microwave sensor can be used to monitor flood and ground water recharge. Predicting water table
depth from soil moisture content is conditional because it requires surface soil should be continuous as
drying of surface soil may decouple from the subsurface soil moisture.
The nature of connection between ground water and surface water differ in the arid and humid climate
or during wet and dry season. For example: water table is well separated from the surface water by a
large vadose zone in the arid climate (Becker, 2005). Water deprived vegetation from dry season or
arid climate regime indicates the position and flow of the ground water. During wet season or humid
climate region the ground water perched below the surface water demarcating the interface and hence
it can be used for inferring ground water condition. In this region shallow water tables are more
common and generally it is easy to interpret ground water condition compare to arid region.
Surplus water after evaporation infiltrates the surface and recharges the ground water table. In this
case movement of water is reverse while comparing with the discharge of ground water discussed in
the previous section. Vegetation cover is also control the infiltration of surface water hence it regulate
the ground water recharge (Becker, 2005). Growth, speciation and abundance of vegetation of a region
is depend up on the availability of water and nutrients, atmospheric moisture content, salinity and
acidity/alkalinity. Different plant species provides clue to the occurrence of ground water but the link
between them varies considerably in different climatic regime. In the arid environment ground water
discharge or shallow water table is the only source of water for vegetation whereas in the humid
climate region it is more complex. Soil chemistry also plays a major role to support selective plant
species, hence while considering the vegetation species assemblage to study the ground water
condition, data on soil chemistry should be taken into consideration (Becker, 2005). Distribution,
growth and type of plant species are used as indicator in the remote sensing and GIS method to
determine the ground water conditions. However, it is difficult to estimate the change of ground water
flow from vegetation cover due to its late response to the flow change. Some researchers (Batelaan et
al., 1998) classified vegetation cover using principal component analysis to study the ground water
discharge in a wetland in Belgium. They also estimated the flow rates and travel time using combined
data of remote sensing and hydrochemistry in a GIS GRASS environment. Although based on type of
vegetation indicator classifying interaction of ground water with the surface i.e. recharge/ discharge
area, appears simple but in practice it is tricky.
11
12
13
Figure 4: A section of alluvial fan schematically representing change of facies from the upper to the
lower fan and the groundwater flow system. Arrows indicate groundwater flow. (Modified
after Meijerink et al. 2007)
14
Figure 5: Alluvial fans in an arid region of Death Valley, USA with distinct difference between upper
(U), middle (M) and lower (L) parts, but with different provenance. Fans on the left: minor seepage on
lower parts, Fans on the right: middle part of the fan slightly dissected and has pheareophyte
vegetation. F.P.: flood plains. (Source: Meijerink et al. 2007)
Another case study from Tibet (Fig. 6) shows large fan complex developed in the tectonic regime
where outflows are perennial and mainly recharged by the snowmelt and occasional flash flow. Fans
are highly permeable hence surface runoff infiltrates and groundwater leaves the fan through few
springs. In the head of the fan few channels with dense vegetation are found and differential upliftment
along the southern flank of the area is evident of dissected older fans.
15
16
Using remotely sensed image interpretations of the various terrains with their specific flow model,
recharge and storage properties can be distinguished and lineaments can be traced out. However,
hydrogeological assessment relies mostly on geological and geophysical survey and on analysis using
experiential data from the existing wells.
17
For sustainable management and development of water resource the monitoring water quality and
quantity is very essential (Sharma et al. 2015). Remote sensing and GIS techniques are used directly or
indirectly for studying water quality and quantity with temporal changes mainly for the river, lake,
snow/glacier or ground water resources. Concentrations of specific parameters and specific properties
of water are being monitored using suitable sensors for assessment of water quality.
Emitted energy (reflectance) of surface water will change due to turbidity, presence of
phytoplankton/algae, specific chemical constituents or dissolve organic matter, oil spill etc. (Sharma et
al., 2015). Based on the emitted energy from the water surface, which is recorded by different sensor,
water quality is monitored and the change in the energy is studied by the remote sensing tools. In India
water quality of river, lake and pond etc. has been analyzed by many researchers using remote sensing
tool. Suspended materials are common pollutant in the surface water. Different sensors carried by
18
R = X + YZ or R = XYZ
Where R= reflectance, Z= water quality parameter, X & Y = empirically derived factors (Sharma et
al., 2015; Ritchie et al., 1974). This equation can be rewritten (Sharma et al., 2015; Schiebe et al.,
1992) based on the physical relationship model between spectral and physical properties of surface
water.
Here Ri = reflectance of surface water for specific wave band i, Cs = concentration of suspended
sediments, Si = reflectance saturation level at high suspended sediment concentration for wave band i,
Pi = concentration parameter, which is equal to the concentration for reflectance of 63% of saturation
level in wave band i (Sharma et al. 2015). Similar to turbidity presence of chlorophyll will change the
reflectance of surface water hence various algorithm and wavelength are used to monitor surface water
bodies and eutrophication of in lake. Seasonal change of chlorophyll concentration can be estimated
using following equation (applied for Chesapeake Bay by Harding et al. 1995):
Here, x & y = empirical value derived from in situ measurement, R1= radiance at 460 nm, R2 =
radiance at 490 nm, R3 = radiance at 520 nm (Sharma et al. 2015). Now various satellite sensors like
IKONOS, OCTS (ocean color and temperature scanner), MOS (Modular optical scanners) are used for
measuring chlorophyll in the surface water (Sharma et al. 2015).
Remote sensing and GIS techniques are also used for preparing map showing spatial variation of
groundwater quality parameter such as arsenic, fluoride, chloride, TDS (total dissolve solid), TH (total
hardness) nitrate, iron, SAR (Sodium absorption ratio) and bacterial contamination to identify affected
area and for risk assessment. In addition to this salt water intrusion is another issue in the coastal areas.
For this study groundwater samples are collected from predetermined area followed by chemical
19
Generation of study
Lab analysis data tabulation map
Interpretation of Thematic
Spatial distribution of groundwater
map
quality for drinking and irrigation
purpose as per the national /
international standard
Recommendation
Figure 8: Flow chart for application of remote sensing and GIS technique for water pollution study
20
The first part of this module deals with the application of remote sensing and geographic information
system (GIS) techniques for the groundwater studies and the second part (section 6) explains the
application on water pollution studies. In the first part (Section 1 -5) is written for those who are
interested in applying these techniques for groundwater studies, be it exploration, evaluation of
resources, management or required data processing. Inadequate experience in image interpretation and
lack of knowledge will hinder the appropriate uses of satellite images and aerial photographs. Here
emphasis is given on the interpretation of aerial photograph and satellite images of diverse geological
terrain having different climatic settings for extraction of information about groundwater. Map of the
principle aquifer system of India is given in the last part of the section 5. At last, this module
discusses about the determination of water pollution using remote sensing and GIS techniques.
9. References
MurryY. B., 2013, Geospatial Modeling for Assessing Ground Water Resources: A Study in Dimapur
Area, NE India, M.Tech. Thesis, IIRS, ISRO, Dehradun.
Aquifer Systems of India, Central Ground Water Board (CGWB), Ministry of Water Resources,
Government of India, 2012
Report of Central Ground Water Board (CGWB), 2009
NRSA (National Remote Sensing Agency) (2008) Ground Water prospects mapping for Rajiv Gandhi
National Drinking Water Mission, manual. p.3-5.
Becker M. W. (2006): Potential for Satellite Remote Sensing of Ground Water. GROUND WATER,
Vol. 44, No. 2. pp: 306 – 311.
Meijerink A.M.J., BannertD., Batelaan O., Lubczynski M. W., Pointeter T. , (2007). Remote Sensing
Application to Groundwater, IHP-VI, Series on Groundwater No.16, Published by the United
Nations Educational, Scientific and Cultural Organization, France.
Sharma B., Tyagi S., Singh P., Dobhal R., Jaiswal V., (2015) Application of Remote Sensing and GIS
in Hydrological Studies in India: An Overview. Natl. Acad. Sci. Lett., Vol- 38 (1), PP:1–8.
21
22
Web resources:
http://water.usgs.gov/edu/earthwherewater.html
http://www.kcse-online.info/geog/5.html
23