GIS for the assessment of the groundwater recharge

Environ Geol (2009) 58:185–195
DOI 10.1007/s00254-008-1504-9
ORIGINAL ARTICLE
GIS for the assessment of the groundwater recharge
potential zone
Hsin-Fu Yeh Æ Cheng-Haw Lee Æ
Kuo-Chin Hsu Æ Po-Hsun Chang
Received: 3 May 2008 / Accepted: 27 July 2008 / Published online: 12 August 2008
Ó Springer-Verlag 2008
Abstract Water resources in Taiwan are unevenly distributed in spatial and temporal domains. Effectively
utilizing the water resources is an imperative task due to
climate change. At present, groundwater contributes 34%
of the total annual water supply and is an important fresh
water resource. However, over-exploitation has decreased
groundwater availability and has led to land subsidence.
Assessing the potential zone of groundwater recharge is
extremely important for the protection of water quality and
the management of groundwater systems. The Chih-Pen
Creek basin in eastern Taiwan is examined in this study to
assess its groundwater resources potential. Remote sensing
and the geographical information system (GIS) are used to
integrate five contributing factors: lithology, land cover/
land use, lineaments, drainage, and slope. The weights
of factors contributing to the groundwater recharge are
derived using aerial photos, geology maps, a land use
database, and field verification. The resultant map of the
groundwater potential zone demonstrates that the highest
recharge potential area is located towards the downstream
regions in the basin because of the high infiltration rates
caused by gravelly sand and agricultural land use in these
regions. In contrast, the least effective recharge potential
area is in upstream regions due to the low infiltration of
limestone.
H.-F. Yeh C.-H. Lee K.-C. Hsu P.-H. Chang
Department of Resources Engineering,
National Cheng Kung University,
Tainan, Taiwan
C.-H. Lee (&) K.-C. Hsu
Sustainable Environment Research Center,
National Cheng Kung University,
Tainan, Taiwan
e-mail: [email protected]
Keywords Geographical information system Groundwater recharge Chih-Pen Creek basin
Introduction
Groundwater recharge refers to the entry of water from the
unsaturated zone into the saturated zone below the water
table surface, together with the associated flow away from
the water table within the saturated zone (Freeze and
Cherry 1979). Recharge occurs when water flows past the
groundwater level and infiltrates into the saturated zone. It
is an extremely important water component of the circulation cycle in nature. Many factors affect the occurrence
and movement of groundwater in a region including
topography, lithology, geological structures, depth of
weathering, extent of fractures, primary porosity, secondary porosity, slope, drainage patterns, landform, land use/
land cover, and climate (Mukherjee 1996; Jaiswal et al.
2003). On-site hydrogeology experiments and geophysics
surveys help to explain the process of groundwater
recharge and evaluate the spatial–temporal difference in
the study region. However, these surveys often focus on a
single affecting factor or an indirect site-specific experiment for groundwater recharge, reducing the reliability
of the explanation. Recently, remote sensing has been
increasingly employed to replace on-site exploration or
experiments. Remote sensing not only provides a widerange scale of the space–time distribution of observations,
but also saves time and money (Murthy 2000; Leblanc
et al. 2003; Tweed et al. 2007). Sener et al. (2005) pointed
out that remote sensing can effectively identify the characteristics of the surface of the earth (such as lineaments
and geology) and can also be used to examine groundwater recharge. Bierwirth and Welsh (2000) applied remote
123
186
sensing to determine the preferential path of groundwater
recharge in an area.
The National Remote Sensing Agency (NRSA 1987) in
India was the first to integrate information from remote
sensing and the technology of the geographical information
system (GIS) for delineating the groundwater recharge
potential zone. GIS is used to manage, utilize, and classify
the results of remote sensing, to explore sites, to combine
the factors of groundwater recharge potential, and to provide appropriate weight relationships (Krishnamurthy et al.
1996; Saraf and Choudhury 1998; Sener et al. 2005).
Salama et al. (1994) used aerial photos and information
from a satellite to derive the lithology, topography, and
geological characteristics of the Salt River in Western
Australia. These properties can be used to determine the
mechanism of groundwater flow and the groundwater
recharge zone. Their analytical results demonstrate that the
sandy plain is the major recharge zone. Edet et al. (1998)
classified groundwater recharge potential zones in southeast Nigeria as high, medium, or low. They found
that linear features, drainage, lithology, temperature of
groundwater, vertical hydraulic conductivity, yield, and
transmissivity closely control the recharge potential zones.
Singh and Prakash (2002) plotted a groundwater recharge
potential map of a sub-watershed from the geology, lineament maps, drainage, slope, and the thickness of the soil
covered. Their results show that the well-yield data in India
is closely related to the groundwater recharge potential
zone. Jaiswal et al. (2003) concluded that there is a need to
adjust the information of satellites and GIS to agree with
the on-site geology, particularly in typical hard rock terrain, where groundwater occurrence is complex and
restricted. Shaban et al. (2006) explored the recharge
potential map of the Occidental Lebanon, and found that
the regions of hard, fractured, and karstified limestone were
excellent potential areas for groundwater recharge, while
the least effective recharge potential was in high-populated
areas and in relatively flat areas covered by soft materials.
The diffusion of pollutants in groundwater is fastest in the
most efficient recharge zones. Tweed et al. (2007) verified
that the integration of remote sensing and GIS reduces the
uncertainty of hydro-geological data in terms of both
macroscopic (climate, change of land utilization) and
microscopic (preferential flow) factors. Their data can be
used to analyze groundwater numerical models or water
balance.
Groundwater is an important source of water for
industrial, agricultural, and domestic uses in Taiwan. Taiwan is narrow, with a small area and a high-population
density. The elevation of the terrain of basins is great and
steep, meaning that most precipitation becomes runoff and
drains directly to the ocean very quickly. Water resources
and water demand are unevenly distributed spatially and
123
Environ Geol (2009) 58:185–195
temporally. Water shortage has recently become an
important issue due to climate change. The government
needs to regulate the usage of water resources in order to
solve the problem of water shortages. The continuous
development of the economy has led to an increase in water
consumption, and has consequently resulted in shortages of
surface water. Therefore, the reliance on groundwater
resources has increased, leading to the over-consumption
of groundwater, and causing ecological problems such as
decreased groundwater levels, water exhaustion, water
pollution, deterioration of water quality, and seawater
intrusion. These cause serious problems and threaten both
people’s livelihoods and overall national development.
Therefore, it is important to thoroughly understand the
groundwater resources of Taiwan in order to enhance the
efficiency and performance of their planning, utilization,
administration, and management. The government of Taiwan has invested significant labor and financial resources
to survey five main groundwater areas (Cho-Shui River
alluvial fan, Pingtung Plain, Chia-Nan Plain, Lan-Yang
Plain, and Hsin-Chu and Miaoli Region) to construct a
database of hydrogeology and groundwater. Additionally,
groundwater monitoring stations have been established, but
only on the western plain of Taiwan. Very limited information, such as precipitation, river flux, hydro-geological
properties, groundwater consumption, and groundwater
recharge, is available for the eastern mountain area of
Taiwan. The assessment and planning of groundwater
resources is particularly difficult for the mountain area of
Taiwan.
Groundwater resources in the eastern mountainous
region of Taiwan are becoming increasingly insufficient
due to economic development. The declining groundwater
level in the mountain area is an indication of the decrease
in groundwater resources (Yeh et al. 2007). Therefore, this
study uses the Chih-Pen Creek basin in eastern Taiwan as
the study domain. The Chih-Pen Creek basin was selected
because it is compact but important for eastern Taiwan.
The groundwater recharge potential zones are assessed for
the Chih-Pen Creek basin.
Methodology
Study approach
Remote sensing technology, such as aerial photos, was
used in the present study to identify the geological features,
topography, and distribution of the rivers in the region.
Additionally, the Land Utilization Survey Database, geologic maps, and on-site investigation were adopted to
quantitatively and qualitatively describe the hydro-geological conditions of the area. The different polygons in the
Environ Geol (2009) 58:185–195
thematic maps were labeled separately. The influence of
the factors of groundwater recharge and the interaction
between the factors were examined. Weighting values were
assigned according to the on-site situation. The distribution
of the groundwater recharge potential zone was determined
by coordinating it with the space integrating function of the
geographical information system (GIS). Figure 1 illustrates
the flowchart of this investigation.
187
Table 1 Factors influencing groundwater recharge classified criteria
Factor
Basis of categorization
Lithology
Rock type, weathering character, joints, fractures
Land cover/land
use
Type, areal extent, associated vegetation
Lineaments
Lineament-density value
Drainage
Drainage-density value
Slope
Slope gradient
The Infrastructure of the Groundwater Recharge
Potential Model
The groundwater recharge potential zone has been assessed
in many countries (Krishnamurthy et al. 1996; Saraf and
Choudhury 1998; Shahid et al. 2000; Jaiswal et al. 2003;
Sener et al. 2005; Shaban et al. 2006). However, the
groundwater recharge potential zone has not yet been
assessed in Taiwan, especially in mountain watersheds. In
this study, the weights of different factors for groundwater
recharge potential and the score under various characteristics were assessed based on the characteristics of the
Chih-Pen Creek basin. The factors influencing groundwater
recharge, and their relative importance, are compiled from
previous literature. Duplicate factors were combined and
only representative factors were extracted. This study uses
lithology, land use/cover, lineaments, drainage, and slope
as the five significant factors affecting groundwater
recharge potential. The factors influencing groundwater
recharge potential, which are listed in Table 1. GIS technology was used to digitize the hydrologic and geographic
information, and a fundamental database was constructed.
Appropriate scores were set for different factors. Finally,
the spatial analysis function was used to demonstrate the
groundwater recharge potential zone of the research area.
Establishment of groundwater recharge
potential-related factors
Lithology
Shaban et al. (2006) pointed out that the type of rock
exposed to the surface significantly affects groundwater
recharge. Lithology affects the groundwater recharge by
controlling the percolation of water flow (El-Baz and
Himida 1995). Although some investigations have ignored
this factor by regarding the lineaments and drainage characters as a function of primary and secondary porosity, this
study includes lithology to reduce uncertainty in determining lineaments and drainage.
Land use/cover
Fig. 1 Methodology flowchart for the groundwater potential zone
Land use/cover is an important factor in groundwater
recharge. It includes the type of soil deposits, the distribution of residential areas, and vegetation cover. Shaban
et al. (2006) concluded that vegetation cover benefits
groundwater recharge in the following ways. (1) Biological
decomposition of the roots helps loosen the rock and soil,
so that water can percolate to the surface of the earth easily.
(2) Vegetation prevents direct evaporation of water from
soil. (3) The roots of a plant can absorb water, thus preventing water loss. Leduc et al. (2001) estimated the
difference in the amount of groundwater recharge due to
changes of land utilization and vegetation from changes in
the groundwater level. Land use/cover was included in this
study as an important factor affecting the groundwater
recharge process.
123
188
Environ Geol (2009) 58:185–195
Lineaments
The analysis of lineaments has been applied extensively to
explain geological status since geological images were
first utilized in the 1930s. Lineaments are generally
referred to in the analysis of remote sensing of fractures or
structures. Lineament photos from satellites and aerial
photos have similar characteristics but the results of the
explanation in on-site may be different. Lineaments are
currently not fully defined. O’Leary et al. (1976) has
defined lineaments as the simple and complex linear
properties of geological structures such as faults, cleavages, fractures, and various surfaces of discontinuity, that
are arranged in a straight line or a slight curve, as detected
by remote sensing. Many non-geological structures, such
as roads and channels, cause errors in the analysis of
lineaments. Therefore, geologic maps and on-site investigations must be used to eliminate possible errors.
Lineaments may be used to infer groundwater movement
and storage. Lattman and Parizek (1964) were the first to
adopt a lineaments map to exploit groundwater. Thereafter, many scholars have applied this approach in
complicated geological regions (Solomon and Quiel
2006). The present study used lineament-length density
(Ld, L-1) (Greenbaum 1985), which represents the total
length of lineaments in a unit area, as:
iP
¼n
Li
Ld ¼ i¼1
A
iP
¼n
Li denotes the total length of lineaments (L) and
where
i¼1
A denotes the unit area (L2). A high lineament-length
density infers high secondary porosity, thus indicating a
zone with high groundwater potential.
of groundwater recharge. Many studies have integrated
lineaments and drainage maps to infer the groundwater
recharge potential zone (Edet et al. 1998; Shaban et al.
2006).
Slope
Rainfall is the main source of groundwater recharge in
tropic and subtropic regions. The slope gradient directly
influences the infiltration of rainfall. Larger slopes produce
a smaller recharge because water runs rapidly off the surface of a steep slope during rainfall, not having sufficient
time to infiltrate the surface and recharge the saturated
zone.
Interrelationships between the factors
of the groundwater recharge potential
There might be interactions between the factors of
groundwater recharge. This study used five factors of
groundwater recharge potential, namely lithology, land
use/cover, lineaments, drainage, and slope. A plot of the
interrelationship between these factors is shown in Fig. 2.
Figure 2 illustrates the primary and secondary interrelationships among the factors. Each relationship is weighted
according to its strength. The representative weight of a
factor of the recharge potential is the sum of all weights
from each factor. A factor with a higher weight value
shows a larger impact on groundwater recharge. Spatial
integration and analysis was performed using GIS technology to demonstrate the groundwater recharge potential
zone as depicted in Fig. 3.
Drainage
The structural analysis of a drainage network helps assess
the characteristics of the groundwater recharge zone. The
quality of a drainage network depends on lithology, which
provides an important index of the percolation rate. The
drainage-length density (Dd, L-1), as defined by Greenbaum (1985), indicates the total drainage-length in a unit
area, and is determined by:
iP
¼n
Dd ¼
where
Si
i¼1
A
iP
¼n
Si denotes the total length of drainage (L) and A
i¼1
denotes the unit area (L2). The drainage-length density is
significantly correlated with the groundwater recharge; a
zone with a high drainage-length density has a high level
123
Fig. 2 The interactive influence of factors concerning the recharge
property (modified from Shaban et al. 2006)
Environ Geol (2009) 58:185–195
189
and it lies between longitudes 121°050 –121°500 E and latitudes 22°350 –22°450 N. Figure 4 shows the geographical
location of the Chih-Pen Creek basin. The research region
belongs to the tropical marine climate, with a mean annual
temperature of 24.5°C and an average annual precipitation
for 1971–2006 of 1,800 mm year-1. During the summer,
southwest monsoons occur and typhoons bring heavy
rainfall. The northeast monsoon brings vapor from the
Pacific Ocean during the winter. Because water vapor is
blocked by the Central Mountains of Taiwan, there is little
rainfall in the winter. Therefore, the wet and dry seasons
are very distinct in this region. The wet season is from May
to October, and the dry season is from November to April.
Evapotranspiration is approximately 750 mm year-1. The
maximum streamflow on annual hydrographs occurs during
August and September, and the minimum flow occurs
during January and February (Yeh et al. 2007).
Fig. 3 GIS technology used in spatial integration and analysis to
demarcate basin groundwater recharge potential zone
Study area
Geographical position and meteorological hydrology
The study area, the Chih-Pen Creek basin, is in the
southeast of Taiwan. The basin encompasses an area of
about 198.4 km2. The length of the river is about 39.3 km
Topography and geology
The Chih-Pen Creek basin can be divided into two topographic units, namely the mountainous terrain and the
alluvial plain. The Chih-Pen Fault splits the basin into east
and west parts. The terrain in the west of the fault is precipitous and consists of metamorphic rock, covering part of
the Central Mountain Range. The east area of the fault is
the Taitung alluvial fan-delta. Because the end of the
alluvial fan has already entered the coastline, the delta
Fig. 4 Location of the study
area
123
190
Environ Geol (2009) 58:185–195
Fig. 5 Three-dimensional
topographic map of the
Chih-Pen Creek basin
plain gradually changes to an alluvial fan delta. The relief
declines by 70 m km-1 from east to west, which is faster
than the north–south decline of about 30 m km-1. Figure 5
shows the three-dimensional topographic map of the ChihPen Creek basin.
Taiwan can be divided into three major geological
terrains, namely the Tananao Schist terrain, the slate formation terrain, and the Neogene-Quaternary sedimentary
strata terrain (Yen 1970). The only rock stratum, the ChihPen Formation, in the Chih-Pen Creek basin belongs to the
slate formation terrain. In contrast to the other terrains,
which are composed of a variety of rock types, the slate
formation terrain has monotonous lithology. The rock types
in the Chih-Pen Formation are mainly slate with subordinate meta-sandstone. However, lithofacies was determined
from the interpretation of the deposits as a submarine-slope
sequence. Additionally, the area contains folds and cleavage, indicating that the whole region is suffering very
strong compression (Lin and Lin 1998).
Groundwater recharge potential factor
establishment and spatial analysis
Factor establishment
This study analyzed the hydrologic and geographic attributes of the Chih-Pen Creek basin, and identified five
major factors influencing groundwater recharge potential,
namely lithology, land use/cover, lineaments, drainage, and
slope. Each factor was examined and was assigned an
appropriate weight.
Each recharge potential factor may influence the
groundwater recharge process to a different degree.
Moreover, the factors are interdependent. Figure 2 illustrates the conceptual graph of the interrelationships
among the groundwater recharge potential factors. A major
123
Table 2 Relative rates for each factor
Factor
Calculation process
Proposed relative
rates
Lithology
3 9 1.0 = 3.0
3.0
Land cover/land
use
1 9 1.0 + 3 9 0.5 = 2.5
2.5
Lineaments
2 9 1.0 = 2.0
2.0
Drainage
1 9 1.0 + 1 9 0.5 = 1.5
1.5
Slope
1 9 1.0 + 1 9 0.5 = 1.5
1.5
R10.5
interrelationship between two factors is assigned a weight
of 1.0. A minor interrelationship between two factors is
assigned a weight of 0.5. Finally, the total weight of each
factor is the representing weight of the recharge potential.
For instance, major interrelationships exist for lithology on
lineaments, drainage and, land use/cover. Therefore, its
evaluated weight is 3.0. This high weight value means
that the factor significantly influences the groundwater
recharge. Table 2 shows the process for determining the
relative rate of each factor.
The extent of the influence of every factor on groundwater recharge was assessed from the interrelationships
among the factors (major and minor). Analytical results
demonstrate that the factors influencing the groundwater
recharge potential of the Chih-Pen Creek basin, in
descending order, are lithology, land use/cover, lineaments,
drainage, and slope. Lithology and land use/cover are the
major factors influencing the basin groundwater recharge
potential. The score of each recharge potential factor was
calculated as 100 multiplied by the weight of the recharge
potential divided by the total weight of each recharge
potential factor. Table 3 shows the calculation approach.
A 1 km 9 1 km grid was used to quantify the score of
each recharge potential factor based on the characteristics
of lithology, land use/cover, lineament-length density,
Environ Geol (2009) 58:185–195
191
Table 3 Score of each recharge potential factor
Factor
Calculation process
Proposed score
of each influencing
factor
Lithology
100 9 (3/10.5) = 29
29
Land cover/land use
100 9 (2.5/10.5) = 24
24
Lineaments
Drainage
100 9 (2/10.5) = 19
100 9 (1.5/10.5) = 14
19
14
Slope
100 9 (1.5/10.5) = 14
14
Analysis of land use/cover
R100
drainage-length density, and slope of the Chih-Pen Creek.
The upper threshold of the score of each recharge potential
factor was set to be the score of the corresponding recharge
potential factor. For example, the highest value of lithology
was 29. Table 4 lists the scores of the recharge potential
factors of the study region.
Spatial analysis
Analysis of the types of lithology
The 1/250,000 geological map of Taiwan reveals that the
upstream region of the Chih-Pen Creek basin is mainly
Table 4 Categorization of factors influencing recharge potential in
the Chih-Pen Creek basin
Factor
Slope gradient
Domain of effect
55–90°
Lineament
density
4
7
15–35°
11
0–15°
14
4
1.5–3.0
7
3.0–4.5
11
[4.5
14
0.0–0.4(lineament per 1 km2)
6
0.4–0.8
13
0.8–1.2
19
Land cover/land Building
use
Forest
Agricultural land
Surface water body or river
channel
Lithology
Proposed
weight
of effect
35–55°
Drainage density 0.0–1.5(segment per 1 km2)
Shale, slate, Phyllite black schist
composed of metamorphic limestone. In the midstream
region, the basin mainly consists of the Tananao Schist, the
slate and phyllite of the Pilushan Formation, and the slate
of the Lushan Formation. The downstream region is mainly
composed of gravelly sand. The gravel sand in the downstream and the metamorphic limestone in the upstream are
excellent regions for percolation. Figure 6 shows the distribution of lithology in the study area.
6
12
18
24
7
Phyllite intermixed with quartz
sandstone
15
Marble
22
Gravelly sand
29
The land use/cover of the Chih-Pen Creek basin was
assessed in accordance with the Territory Utilization Status
Survey Database established by the Land Administration
Bureau of the Taiwan Provincial Government. The information was obtained from the original map of the territory
utilization survey drawn by the Land Administration Office
of the counties and cities. In the urban planning districts,
the 1/500, 1/600, 1/1,000 or 1/1,200 cadastral map was
used as the base map. In the non-urban planning districts,
the 1/5,000 farming map from the Agriculture and Food
Agency was used as the base map. The total area of
the Chih-Pen Creek basin is approximately 200 km2. The
mountain area occupies around 185 km2, i.e. 92.5% of the
basin. The area of level ground is around 15 km2, which is
7.5% of the basin. Most of the mountain area, which is
forestland covered by vegetation is located in the midstream and upstream. This reduces runoff and increases the
recharge. The level ground is mainly distributed in the
downstream alluvial plain and fragmentary small terrace
along the riverbank. The left bank is mostly paddy fields,
with some orchards. The right bank is mainly dry farming
and orchards. Figure 7 illustrates the land use/cover distribution diagram of the study.
Analysis of lineament-length density
For the assessment of the lineaments of Chih-Pen Creek
basin, a stereoscope was used to interpret the aerial photos
from an agricultural aerial survey from 1996 to 2000,
which was verified on-site. The lineaments of the Chih-Pen
Creek basin are mainly distributed in the mid- and downstream. The lineament-length density of unit grids through
which the lineaments passed through was approximately
0.8–1.2 km (km2)-1. Some unit grids had densities of 0.4–
0.8 km (km2)-1. Figure 8 depicts the lineament-length
density diagram of the study area.
Analysis of drainage-length density
The distribution of drainage in the Chih-Pen Creek basin
was determined using the aerial photos of the agricultural
aerial survey undertaken in 1996–2000. The Chih-Pen
123
192
Environ Geol (2009) 58:185–195
Fig. 6 Lithology map of the
study area
Fig. 7 Land use/cover map of
the study area
Creek basin can be roughly divided into dendritic and
grid drainage patterns. The dendritic drainage pattern is
distributed mostly in the northwest region of the basin,
while the southeast is mainly grid drainage. The dendritic
drainage pattern is influenced more by lithology than by
its structure. For example, the mudstone area in southern Taiwan always results in a high density dendritic
123
drainage pattern. The drainage of the Chih-Pen Creek
basin is distributed mainly downstream. The drainagelength density is often larger than 3.0 km (km2)-1,
making the region an excellent percolation recharge zone.
The up- and midstream are mostly 1.5–3.0 km (km2)-1.
Figure 9 shows the drainage length density of the study
area.
Environ Geol (2009) 58:185–195
193
Fig. 8 Lineament density map
of the study area
Fig. 9 Drainage density map of
the study area
Analysis of slope
The slope analysis function in GIS was used to assess the
variation of slope in the Chih-Pen Creek basin using data
from the Digital Terrain Model (DTM) database in Taiwan.
The Council of Agriculture authorized the Aerial Survey
Office of Forestry Bureau to measure and produce the
digital terrain information. Three-dimensional aerial photos
with regular 40 m sampling intervals were applied. The
Chih-Pen Creek basin belongs to a mountain area. The
gradient of its slope is greater than 35°. The precipitous
terrain causes rapid runoff and does not store water easily.
A small part of the downstream is fairly gentle, so the time
for percolation is increased. Since the terrain varies
123
194
Environ Geol (2009) 58:185–195
Fig. 10 Slope gradient map of
the study area
Fig. 11 Groundwater potential
zones map prepared using GIS
significantly, the factor slope varies inversely with the
groundwater recharge. Figure 10 illustrates the distribution
of slopes in the study area.
Demarcation of the groundwater recharge zone
Our results demonstrate that the groundwater recharge
potential zone of this basin can be divided into five grades,
123
namely very good, good, moderate, low, and poor, based
on the analysis of the five factors of groundwater recharge
potential. Analytical results demonstrate that the excellent
groundwater recharge potential zone is concentrated in the
downstream region due to the distribution of gravelly
stratum and agricultural land with a high infiltration ability.
Additionally, the concentration of drainage also helps
the streamflow to recharge the groundwater system. The
Environ Geol (2009) 58:185–195
upstream region is less important and is influenced by
metamorphic limestone. Figure 11 shows the groundwater
recharge potential diagram of the study area.
Conclusions and recommendations
This study produced a groundwater recharge potential map
of the Chih-Pen Creek basin in Taiwan. The results indicate
that the most effective groundwater recharge potential zone
is located downstream. In this region, the gravelly stratum
and agricultural land have a high infiltration ability. Additionally, the concentration of drainage also indicates the
ability of streamflow to recharge the groundwater system.
The upstream region is least effective for groundwater
recharge, mainly due to its metamorphic limestone. This
study has established the interrelationships between the
groundwater recharge potential factors and the groundwater
recharge potential scores from the general hydrology
characteristics of Taiwan. Since the groundwater recharge
potential is directly correlated with percolation, the established scores may be more accurate and objective if the rate
of percolation and hydraulic conductivity of each recharge
potential factor can be measured in a laboratory or on-site.
The groundwater recharge potential zones are demonstrated using the grid model, which can be partially
modified to study groundwater recharge potential factors (such as changes in terrain and river courses caused by
an earthquake, or changes in land utilization) in a small area
in the future, thus avoiding extensive re-estimation, which
requires a lot of time and labor. The isotopic tracer technique will be applied to verify the model results in future
research. Additionally, the quantity of pumped groundwater
should be considered in the groundwater recharge.
Acknowledgments The authors would like to thank the Water
Resources Agency of the Ministry of Economic Affairs of the
Republic of China, Taiwan, for financially supporting this research
under Contract No. MOEAWRA0950416.
References
Bierwirth PN, Welsh WD (2000) Delineation of recharge beds in the
Great Artesian Basin using airborne gamma-radiometrics and
satellite remote sensing. Report for the National Landcare
Program, Bureau of Rural Sciences, Canberra, Australia, 33 p
Edet AE, Okereke CS, Teme SC, Esu EO (1998) Application of remote
sensing data to groundwater exploration: a case study of the Cross
River State, southeastern Nigeria. Hydrogeol J 6:394–404
El-Baz F, Himida I (1995) Groundwater potential of the Sinai
Peninsula. Egypt, Project Summery, AID, Cairo
Freeze RA, Cherry JA (1979) Groundwater. Prentice-Hall, Inc.,
Englewood Cliffs
Greenbaum D (1985) Review of remote sensing applications to
groundwater exploration in basement and regolith. Br Geol Surv
Rep OD 85:36
195
Jaiswal RK, Mukherjee S, Krishnamurthy J, Saxena R (2003) Role of
remote sensing and GIS techniques for generation of groundwater prospect zones towards rural development—an approach.
Int J Remote Sens 24:993–1008
Krishnamurthy J, Venkatesa Kumar N, Jayaraman V, Manivel M
(1996) An approach to demarcate groundwater potential zones
through remote sensing and geographic information system. Int
J Remote Sens 17:1867–1884
Lattman LH, Parizek RR (1964) Relationship between fracture traces
and the occurrence of groundwater in carbonate rocks. J Hydrol
2:73–91
Leblanc M, Leduc C, Razack M, Lemoalle J, Dagorne D, Mofor L
(2003) Application of remote sensing and GIS for groundwater
modeling of large semiarid areas: example of the Lake Chad
Basin, Africa. In: Hydrology of Mediterranean and semiarid
regions conference, Montpieller, France. Red Books Series, vol
278. IAHS, Wallingford, pp 186–192
Leduc C, Favreau G, Schroeter P (2001) Long-term rise in a Sahelian
water table: the Continental Terminal in southwest Niger.
J Hydrol 243:43–54
Lin MS, Lin WH (1998) Explanatory text of the geologic map of
Taiwan. Scale 1:50,000 Sheet 64. Central Geological Survey,
Taiwan
Mukherjee S (1996) Targeting saline aquifer by remote sensing and
geophysical methods in a part of Hamirpur–Kanpur, India.
Hydrogeol J 19:53–64
Murthy KSR (2000) Groundwater potential in a semi-arid region of
Andhra Pradesh-a geographical information system approach. Int
J Remote Sens 21:1867–1884
NRSA (1987) Land and water resources survey of drought affected
district-Kolar, Karnataka. Technical Report, National Remote
Sensing Agency, Hyderabad
O’Leary DW, Friedman JD, Poh HA (1976) Lineaments, linear,
lineations: some standards for old terms. Geol Soc Am Bull
87:1463–1469
Salama RB, Tapley I, Ishii T, Hawkes G (1994) Identification of areas of
recharge and discharge using Landsat-TM satellite imagery and
aerial photography mapping techniques. Hydrogeol J 162:119–141
Saraf AK, Choudhury PR (1998) Integrated remote sensing and GIS
for groundwater exploration and identification of artificial
recharge sites. Int J Remote Sens 19:1825–1841
Sener E, Davraz A, Ozcelik M (2005) An integration of GIS and
remote sensing in groundwater investigations: a case study in
Burdur, Turkey. Hydrogeol J 13:826–834
Shaban A, Khawlie M, Abdallah C (2006) Use of remote sensing and
GIS to determine recharge potential zone: the case of Occidental
Lebanon. Hydrogeol J 14:433–443
Shahid S, Nath SK, Roy J (2000) Groundwater potential modeling in
a soft rock area using a GIS. Int J Remote Sens 21:1919–1924
Singh AK, Prakash SR (2002) An integrated approach of remote
sensing, geophysics and GIS to evaluation of groundwater
potentiality of Ojhala sub-watershed, Mirjapur district, U. P.,
India, http://www.GISdevelopment.net
Solomon S, Quiel F (2006) Groundwater study using remote sensing
and geographic information system (GIS) in the central
highlands of Eritrea. Hydrogeol J 14:1029–1041
Tweed SO, Leblanc M, Webb JA, Lubczynski MW (2007) Remote
sensing and GIS for mapping groundwater recharge and
discharge areas in salinity prone catchments, southeastern
Australia. Hydrogeol J 15:75–96
Yen TP (1970) Structural analysis of the slate formation of Taiwan.
Bull Geol Surv Taiwan 21:1–51
Yeh HF, Chang PH, Hsu KC, Lee CH (2007) Assessment of
groundwater recharge in Chih-Pen and Jin-Lin Creek basins. In:
Proceedings of 3rd conference on resources engineering in
Taiwan, pp 1–15
123