Jodhpur Final Report PDF
Jodhpur Final Report PDF
FINAL REPORT
Sponsored by:
Ground Water Department
Government of Rajasthan
Jodhpur
Studied By:
National Institute of Hydrology,
Roorkee – 247 667
May, 2011
Project Data
ii
Study Team
iii
Foreword
The problem of rising groundwater level in some parts of the Jodhpur city
located in the Thar Desert of the western Rajasthan that is known for a water scarce
region is not only a matter of concern but also an eye opening issue from the hydrology
point of view. As such, groundwater is a hidden resource; the fact findings from the
study of such resource for a problem that has been originated from the human
interventions into the associated hydrological and hydro-geological components are, in
fact, a challenging task.
In order to find the source and causes of groundwater table rise and to develop
an appropriate management plan to revert back the rising trend and to contain the
groundwater at a safe level, the project entitled “Study on rising groundwater table in
Jodhpur City, and to evolve a management plan for containing the rising trend” has
been entrusted to the National Institute of Hydrology, Roorkee. I put on records my
thanks to Government of Rajasthan, particularly to the Ground Water Department,
Jodhpur for considering NIH capable for this challenging task.
The study has been carried out by the group of scientists and scientific staff of
NIH under the coordination of Dr. N. C. Ghosh, Scientist-F and Principal Investigator of
the study. The study team deserves appreciation and thanks for their sincere and
commendable efforts. I put on record my appreciation to Prof. (Retd.) G. C. Mishra,
consultant of this project, for his contribution.
iv
Acknowledgement
I thankfully acknowledge the help and data support provided by the officials
of the Ground Water Department and Public Health Engineering Department,
Jodhpur in general, and particularly to Shri K. N. Mathur then Chief Engineer,
PHED; Shri U. K. Mathur then Chief Engineer, GWD; Shri C. S. Sekhawat, Chief
Engineer, GWD; Shri G. S. Marwaha, Superintending Hydrogeologist; Shri G. L.
Bora, Senior Hydrogeologist; Shri Niranjan Mathur, TA–Hydrogeologist; Dr. P. S.
Rathore, TA-Hydrogeologist for their keen interest to the study and constant
support as and when requested during the course of the study. Thanks are due to
the members of the ‘Task Force’ for their keen interest and critical review of the
activities from time to time. The data support provided by the SRRSC, RRSC, and
CGWB are also duly acknowledged.
Last but not the least, the assistance provided by Shri Avdesh Sharma, PA
and Shri Haridas, Messanger-Sr. Gr. is also duly acknowledged.
v
Executive Summary
The rise of groundwater table near to the ground surface in some parts of the
Jodhpur city has resulted in a hazard to the people living in the affected areas. Controlling
the rising trend of groundwater level in the urban area of Jodhpur city has emersed as a
challenging task to find the cause and an immediate problem resolving scientific solution
for the Public Health Department and Ground Water Department of the Government of
Rajasthan. In order to find the exact source and causes of groundwater table rise and to
develop an appropriate management plan to revert back the rising trend to contain the
groundwater at a safe level, the present study had been entrusted to the National Institute
of Hydrology, Roorkee, with the following objectives:
For addressing the issues and to reach to logical conclusions and solution, a
systematic in depth analysis of the data/information related to topography, demography,
geological formations, hydrometeorology-hydrology and hydrogeology, groundwater
quantity and quality, sewage flows, inflows and outflows of waters to/from the Jodhpur
city and from the Kailana-Takhatsagar Reservoir have been carried out. The data have
mostly been supplied by the Ground Water Department (GWD) and the Public Health
Engineering Department (PHED) of Jodhpur; and some data have been collected by the
executing agency from the field investigations/surveys and from the other sources (IMD,
NRSC, local agencies) during course of the study.
vi
To analyze the data, all spatially varying databases have been geo-referenced with
reference to the geographic coordinate system (WGS-1984) by their latitudes and
longitudes. Therefore, the analyzed data and results presented in this report can be
verified with the in situ field truth.
The data analyses and results have been reported in 17 Sections. The Section-1
presents an introduction to the problem; Section-2 elaborates geography and topographic
of the study area; Section-3 describes the geological formations of the area; Section-4
illustrates hydro-meteorological data analyses; Section-5 explains demography and water
usages and requirements; Section-6 brings out the Stage-Area-Capacity analyses and
curves for the Kailana-Takhatsagar and Umaid Sagar Reservoirs; Section-7 provides
analyses of Inflows-outflows data of the Kailana-Takhatsagar Reservoir and their water
balances; Section-8 contains analyses of the groundwater data; Section-9 illustrates
sewages flow measurements and data analyses; Section-10 describes the surface water
balance of the study area; Section-11 elaborates analyses of aquifer parameters
estimation; Section-12 focusses on groundwater quality data and scenario analyses;
Section-13 demonstrates the inputs data preparation for groundwater modelling of the
study area; Section-14 provides description of the modelling aspects and analyses of
results obtained from the modelling of water table evolution in the study area; Section-15
illustrates the possible management strategies to remediate the problem; Section-16
brings out the summary and conclusions of the study; and finally Section-17 states the
recommendations as to what actions would be necessary to resolve the problem.
The study area comprising of 76 sq. km. encompasses the old and the sprawled
Jodhpur city area including the waterlogged area. The Kailana and the Takhatsagar
Reservoirs, which are two naturally formed cascading type geological faults and are the
source of water supply to the city area, are located outside the boundary of the study area
as they are located in a different geological entity. The topography of the Jodhpur city
area, which has been analyzed making use of surveyed data supplied by the GWD,
Jodhpur in conjunction with the ASTER data, showed a general slope of the Jodhpur city
towards south-west, south, and south-east directions, except for the hilly terrain in the
north-western side. The southern terrain has flatter slope than that of the south-west, and
vii
south-east terrain. The topography level of the city area varies between 202 m and 360 m
above MSL (Mean Sea Level) with a small stretch near to the fort area having higher
elevation of 310 m, whereas most of the area in the city is largely flat terrain having
elevation below 250 m above MSL.
The geological formations of the Jodhpur city area have been analyzed making use of
the bore logs data of 93 locations supplied by the GWD, Jodhpur in ‘ROCKWORKS’
software. Based on the analysis one could infer that the Kailana and Takhatsagar area is
located on the Rhyolite formation, whereas Jodhpur city is mainly located on the
sandstone and shale formation. These two formations have different hydro-geological
properties and cannot be considered as a single system. The geological formations of the
Jodhpur city area are primarily composed of Shales, Sandstones and Rhyolites with the
Quaternary alluvium formations at the top particularly in the plain areas. The Quaternary
alluvium formations vary in thickness from few centimeters to about 75 m and form an
unconfined aquifer.
Analysis of thirty six years (1971-2006) daily rainfall data of the Jodhpur city
indicated that about 68% and 86% of the annual rainfall occurred respectively during
July-August and during the monsoon months (June through September) with average
annual rainfall of 378 mm. The maximum water surface evaporation is observed to be in
the month of May (14.127 mm/day) followed by June (13 mm/day), and minimum is in
the month of December (4.059 mm/day) followed by January (4.255 mm/day) every year.
The demographic data analysis showed a population of 11, 08,950 in the Jodhpur
city area by the year 2010. The population censuses of previous four decades indicated a
growth rate of 3.21% per year. The projected population by the year 2015 is 12, 74,830.
The water supply to the population in the city area has been made on the basis
of a thumb rule in accordance to the supply-demand norm. The quantities of water
supplied per capita per day to the population in the city area from the Kailana-
Takhatsagar Reservoir in different years were 17% to 60% higher than the quantity of
140 lcpd prescribed by the Ministry of Welfare and Housing (MoWH), Govt. of India .
viii
The Stage-Area-Capacity relationships for the Kailana, Takhatsagar and Umaid
Sagar Reservoirs have been formulated which provide functional relationships of the
water spread areas and the reservoir capacities with depths of water in the respective
reservoirs.
The groundwater level contours prepared for the pre and post monsoon of the
years 1996 to 2008 indicated that groundwater flow from the Kailana-Takhatsagar side is
not causing water logging in the waterlogged area as the flow direction is not towards the
water logged area. The source of water causing water logging is generated locally.
To assess the daily sewages outflow from the city area, field investigations and
measurements have been carried out in the three sewerages drains; one near to the
Jodhpur Airport, other one near to the Jodhpur Polytechnic Institute, and the third one at
the Nandri sewage treatment site. The analyses of the measured sewages data show that
the total discharge of sewages from the city areas through these three sewerage systems is
37% of the daily water supplied. The wastewater generated is approximately 65% of the
water supplied during a year. Thus about 28% of the water supplied is joining the
unconfined aquifer below the city.
The groundwater quality data supplied by the PHED, Jodhpur have been analyzed
The spatial variation of the parameters pH, Cl2, TDS, NO3 and SAR in the study area
indicated that the source of water logging and rise in groundwater level in the
ix
problematic area are due to the return flow of water from water supply system and from
the source other than the sewage waters originating from domestic supply. In some
pockets, the seepage from sewage system cannot be ruled out. The quality of
groundwater showed that the groundwater can safely be used for irrigation purposes.
The groundwater simulation modeling has been carried out using visual
MODFLOW software. The responses of the aquifer for different stress periods have been
simulated setting transient state model. Based on the analyses of data and modeling,
different management options have been analyzed, and the best one of those has been
recommended accordingly.
Based on the analysis of satellite imageries, it is found that Jodhpur city is mainly
located on the sandstone and shale formations, which are relatively pervious. Kailana and
Takhatsagar area is located on the Rhyolite formation. As Rhyolite formation exhibits
low permeability (0.058 m/day), the seepage losses from the lake would be very small.
Considering two extreme situations ( i.e. in the first situation, the reservoirs are
hydraulically connected with the under lying aquifer implying that a rise in water table
position in the aquifer in the vicinity of the reservoir boundary would reduce the seepage
from the reservoir, and in such situation the water level in the reservoir forms as a known
water level contour line; in the second situation, in which any change in water level in
the aquifer near the reservoirs do not affect the seepage from the reservoir, and water
level in the reservoir is not linked to the water level in the aquifer), contours are drawn
and flow directions ascertained. As seen from the water level contours maps, the flow
directions are not converging towards the water logged area. Therefore, the seepage from
the reservoir irrespective of its quantity is not entering to the water logged area. The flow
directions are towards the south-east region only; and in the south east region, the water
table contours are having less value than those in the waterlogged area. From the
consideration of hydraulic principle, and from the consideration of the direction of flow,
it is evident that the seepage from the reservoirs is not entering to the waterlogged areas.
x
The contours of pH of groundwater in the study area exhibit three concentration
peaks, two in the water logged area, and another in the south west region. The contours
around each peak show evidence of mechanical dispersion. Mechanical dispersion is
caused by groundwater flow. The two peaks in water logged area indicate that the high
concentration of pH is caused locally as the contours nearer to each peak are mostly
circular. The gradient directions of the contours and decreasing trend around the peak
exhibited in the south west region indicate that groundwater is not entering to the water
logged area from the West or South-West.
Besides three peaks, there are two troughs one towards North and the other in
South East. The contours around each trough show evidence of dispersion. From the
gradient direction of the contours and decreasing trend of pH towards North in the
northern trough, it could be argued that no ground water flow is occurring from north to
the waterlogged area.
There is only one lineament, which is oriented towards the waterlogged area. The
Lineament analysis survey and the geological and geophysical study conducted
independently by the National Geophysical Research Institute (NGRI), Hederabad (2010)
surrounding the Kailana-Takhatsagar Reservoir indicated that the lineaments are oriented
in NNE-SSW to NE-SW directions with no connectivity to the city areas. A few
lineaments with ESE-WNW directions are present but these are small and do not have
connectivity to the city areas. The findings of the NGRI thus corroborate the present
finding based on hydraulic principle. The chances of seepage from the Kailana-
Takhatsagar Reservoir to the waterlogged areas through lineaments are, therefore, very
less.
As water is not entering to the part of the aquifer below the water logged area
from any side, it is evident that water logging is caused because of vertical infiltration of
sewage and return flow of domestic supply.
xi
(ii) Remedial Measures
1. As the first and foremost remedial measure, it is suggested to regulate the quantity
of water being supplied to the city area at the source itself, i.e., regulation of water
from the Kailana-Takhatsagar Reservoir. The regulation needs to be based on per
capita per day water requirement basis. The Jodhpur city being located in the arid
and water scarce region, about 110 liters per capita per day could be taken as the
guideline. The break up of 110 liters is as follows: 70 liter (drinking & toilet
flushing) + 20 liter (commercial uses) + 20% conveyance losses. Industrial water
requirements are to be included separately. For 110 lpcd supply, the quantity of
water requirement for the estimated population of 11,08,950 in the Jodhpur city
for the year 2010 is worked out to be 268.69 lac gallon per day. To meet water
requirement for domestic animals, and kitchen gardens, swimming pools, and
public parks 1/3 of the requirement for the population i.e., about 100 lac gallon
per day extra water be supplied, which will reduce the quantity by about 30%
over the quantity supplied (521.7 lac gpd) in the year 2009.
2. In the affected area, the water supply lines need to be thoroughly checked to find
the locations of leakages, and suitable remedial measures to stop the leakages
need to be taken up. The sewages/drainage lines in the affected area need to be
properly sealed to stop seepage, if any.
3. The topography of the city area in the northern and middle part is of undulating
type. The thickness of the alluvium formation varies from about 2.0 m to 13.0 m
in that part. Therefore, a single generalized safe ground water level is incorrect to
suggest. In general, considering the possibility of capillary rise in alluvium,
groundwater table should be at least 0.5m below the foundation level. In areas
where building basement floors are located, assuming that basement level is 3.5m,
and capillary rise is 0.5m, the groundwater table should be at 4.0 m below the
ground surface. The groundwater level in the affected area thus has to be lowered
down below 4.0 m from the respective ground surface elevation.
xii
4. The terrain being undulating , the area being an urban area, the requirement of
lowering the water table by 4m, minimum depth of alluvium being 2m in some
places, all these aspects do not promote provision of a usual horizontal sub
surface drainage system. However random sub surface drainage trench of 4m
depth, filled with coarse sand and gravel, where possible, can be constructed to
control the rising water. The water entering into such trenches can be led to a
collector caisson, from where water can be pumped out. Construction of such
trenches to control water table rise would depend upon the local terrain and
building locations.
5. Provision of vertical drainage system i.e. by pumping the water from the aquifer
in the problematic area looks feasible, as drainage wells can be constructed with
least interference with the urbanized area. The pumping rate and schedule can be
controlled, the number drainage wells can be increased in a locality as required,
and already such practice has been initiated in the area, all these factors favor
provision of vertical drainage. In region of low transmissivity area i.e.
transmissivity< 30m2 / day large diameter wells of 0.5m can be constructed.
However vertical drainage system would require electrical energy, and would
cause noise pollution.
6. There are three large ponds, namely; Baiji Ka Talab, Fateh Sagar, and Gulab
Sagar, located near to the problematic area. The pond beds are more or less
impervious, or if necessary these can be lined. The pumped water can be
discharged to these ponds through conveyance pipe. From these ponds surface
channel can be constructed to convey the water stored in the pond to the existing
surface drainage system through gravity.
7. A Bentonite clay grout curtain across the lineament which is terminating before
but directing towards the water logged area can be constructed proximity to the
reservoir site to check the groundwater flow, if any, from the reservoirs to the
water logged area.
xiii
Contents
xiv
Section – 9.0: Sewerage and Drainage Data Analysis for Jodhpur City……. 76 - 82
Section – 13.0: Discretization of The Study Area, And Input Data For Groundwater
Modeling………………………………………………………. 102 - 109
13.1 Discretization of the Study Area………………………………. 102
13.2 Boundaries of the Study Area……………………………………104
13.3 Water Bodies in the Study Area………………………………….105
13.4 Initial and Boundary Conditions………………………………….106
13.5 Aquifer Properties and Parameters………………………………. 106
13.6 Inflow and Outflow Stresses …………………………………….108
13.7 Time Step Size and Simulation Period……………………………109
xv
LIST OF FIGURES
2.1 A geo-referenced map of the Jodhpur city showing its sprawled area 8
and water bodies in and around the city
2.3 2-D Digital Elevation Map (DEM) for Jodhpur area showing locations 10
of Kailana & Takhatsagar, small water bodies and two contour lines
along which water table lies at 3 m and 5 m depth below ground
surface
2.4 3-D Digital Elevation Map (DEM) of the Jodhpur area showing 10
location of Kailana and Takhatsagar and area in which water table lies
up to 3 m and 5 m below ground surface
3.2 Map showing locations of 93 bore logs and their sectional lines along 15
which geological formations are shown
4.1 Annual variation of rainfall in Jodhpur area during the last 36 years for 20
the period 1971-2006
4.2 Annual variation of rainfall in Jodhpur area during the last 25 years for 21
the period 1982-2006
4.3 Distribution of daily rainfall pattern in Jodhpur area during the last 28 21
xvi
years for the period 1978-2006
5.2 Municipal wards encompassed in the study area in the Jodhpur city 24
5.3 Water requirements as per MOWH norms and actually supplied per 29
day to the Jodhpur city in different years
7.2 Variation of average inflow per day in different years to the Kailana 43
Lake from the Rajiv Gandhi Lift Canal
xvii
7.7 Variation of supply of water from Kailana Lake to Takhatsagar 48
Reservoir during the year 2009 staring from January,1 –December, 31
(red color line corresponds to without seepage component; blue color
line corresponds to with seepage component)
7.8 Variation of water balance components of the Kailana Lake during the 50
year 2003 staring from January,1 –December, 31
7.9 Variation of water balance components of the Kailana Lake during the 50
year 2005 staring from January,1 –December, 31
7.10 Variation of water balance components of the Kailana Lake during the 51
year 2008 staring from January,1 –December, 31
7.11 Variation of water balance components of the Kailana Lake during the 51
year 2009 staring from January,1 –December, 31
xviii
during the year 2005 staring from January,1 –December, 31
8.1 Pre (red color) and post (black color) monsoon groundwater table 64
contour map for 1996 with the assumption that the Kailana-
Takhatsagar is hydraulically connected to the aquifer below the
Jodhpur city ( Æ Indicates groundwater flow direction)
8.2 Pre (red color) and post (black color) monsoon groundwater table 64
contour map for 2000 with the assumption that the Kailana-
Takhatsagar is hydraulically connected to the aquifer below the
Jodhpur city ( Æ Indicates groundwater flow direction)
8.3 Pre (red color) and post (black color) monsoon groundwater table 65
contour map for 2004 with the assumption that the Kailana-
Takhatsagar is hydraulically connected to the aquifer below the
Jodhpur city ( Æ Indicates groundwater flow direction)
8.4 Pre (red color) and post (black color) monsoon groundwater table 65
contour map for 2008 with the assumption that the Kailana-
Takhatsagar is hydraulically connected to the aquifer below the
Jodhpur city( Æ Indicates groundwater flow direction)
8.5 Pre (red color) and post (black color) monsoon groundwater table 66
contour map for 1996 with the assumption that the Kailana-
Takhatsagar is not hydraulically connected to the aquifer below the
Jodhpur city( Æ Indicates groundwater flow direction)
8.6 Pre (red color) and post (black color) monsoon groundwater table 67
contour map for 2000 with the assumption that the Kailana-
Takhatsagar is not hydraulically connected to the aquifer below the
Jodhpur city( Æ Indicates groundwater flow direction)
8.7 Pre (red color) and post (black color) monsoon groundwater table 67
contour map for 2004 with the assumption that the Kailana-
Takhatsagar is not hydraulically connected to the aquifer below the
Jodhpur city( Æ Indicates groundwater flow direction)
8.8 Pre (red color) and post (black color) monsoon groundwater table 68
contour map for 2008 with the assumption that the Kailana-
Takhatsagar is not hydraulically connected to the aquifer below the
Jodhpur city( Æ Indicates groundwater flow direction)
xix
8.9 Pre (red color) and Post (black color) monsoon groundwater table 69
contour map for the year 2008 with assumption of the Kailana-
Takhatsagar reservoir to be hydraulically connected to the aquifer
below the Jodhpur city
8.10 Pre (red color) and Post (black color) monsoon groundwater table 70
contour map for the year 2008 with assumption of the Kailana-
Takhatsagar Reservoir not hydraulically connected to the aquifer
below the Jodhpur city
9.2 Sewage measurement location near the Jodhpur Airport site; (a) 78
section chosen for measurement of cross-section and flow velocity, and
(b) temporary structural arrangement made for measurement of cross-
section and flow velocity, and measurement of depth of flow
9.5 Generalized rating curve of the sewage flows in the sewerage drain 80
near to the Jodhpur Airport area; developed based on the continuous 7
days field measurements from 6:00 A.M to 11:00 P.M.(April, 2010)
9.6 Generalized rating curve of the sewage flows in the sewerage drain 81
near to the Jodhpur Polytechnic Institute; developed based on the
continuous 7 days field measurements from 6:00 A.M to 11:00
P.M.(June, 2010)
9.7 Generalized rating curve of the sewage flows in the sewerage drain at 82
Nandri Sewage Treatment Plant site; developed based on the
xx
continuous 7 days field measurements from 6:00 A.M to 11:00
P.M.(July, 2010)
xxi
12.6 Spatial variation of NO3 concentration during pre-monsoon season in 98
the study area (• indicates location of sampling points)
13.1 Discretized view of the study area including the water bodies and the 102
affected area within the Jodhpur city ( Number of grid along X
direction = 107; Number of grid along Y direction = 113; size of each
grid = 100 m x 100m)
13.2 Vertical section of the study area indicating discretized view of the 103
geological formations (Number of layers = 6 each of different
thickness)
13.3 Locations of water bodies in and around the Jodhpur city area 104
13.7 A scheme of the input stresses assigned zone wise to the modelling 108
area
14.1 Comparison of the steady-state simulated and observed groundwater 111
table for the pre-monsoon period of the year 2004
15.1 Location of ponds in the problematic area, and the alignment of 117
sewerage system in the Jodhpur city area.
xxii
Section - 1.0: INTRODUCTION
The Jodhpur city (also known as Sun City), that had spread in an area of about
14.5 sq. km. in the year 1972 (Survey of India map, 1972), is presently sprawled over 76
sq. km. between latitudes 26015’ N to 260 20’ N and longitudes 730 0’ E to 730 4’ E in the
Thar Dessert in Western Rajasthan. Being a city of historical importance and the second
largest in the state having progressive infra-structural urban facilities, the city is
expanding very rapidly, mainly towards south-east, south and south-west directions
despite the area located in the arid region having sparse rainfall, hot and dry climatic
conditions. In order to facilitate human activities and their aspirations, administratively
the Jodhpur area has been categorized into 60 wards based on the households and
population. Domestic and other needs of water for the people living in those wards of the
city are mainly supplied by the Jodhpur Municipal Corporation; while the availability of
water from the source to the City is facilitated by the Public Health Engineering
Department (PHED), Jodhpur. The city has a network of organized water supply
distribution lines (not connected to every household) of which, some distribution lines,
mainly in the old city areas, had been developed long back, and some others have been
constructed with the passage of the progression of city’s extension. As the city and its
human habitations are expanding, the demands of water supply and networking of
distribution lines are also increasing, connecting more people and human activities to the
water supply system of the city. These are resulting towards more and more usages of
supplied water to the city. Waters supplied to the city areas are used for domestic uses for
the purposes of drinking, bathing, laundering, car washing, lawn watering, gardening,
etc., and in some areas for supporting minor to medium commercial activities. As such,
waters supplied to the city areas are not used for any major agricultural irrigation
activities. Whatever quantity of water is supplied for the domestic and other activities to
the city area, some fraction of the applied water in excess of the consumptive uses of
water, will be returned back in the form of wastewater from each household and from
water users after their usages. The wastewaters of different qualities, which depend on
To support the water supply requirement to the city area for domestic, commercial
and industrial purposes, since the year 1997 the PHED had switched over from the earlier
mixed supply pattern of surface water and groundwater to the present completely canal
water based supply system through storages in the Kailna-Takhatsagar Reservoir.
According to the data of PHED, Jodhpur (GWD et al., 1999), the supplies of water to the
city were about 183.5 lacs gpd (i.e.,117 lpcd) in the year 1994, 278.6 lacs gpd (i.e., 164.8
lpcd) in the year 1998, 341.7 lacs gpd( i.e., 185.1 lpcd) in the year 2000, 462.4 lacs gpd
(i.e, 211.9 lpcd) in the year 2006, and 521.7 lacs gpd (i.e., 220 lpcd) in the year 2009. The
waters supplied to the city areas in the respective years were liable to generate
wastewaters. To maintain a perfect balance between water supply-disposal of
wastewaters of such quantities, and resulting shallow water table position, a fitting water
supply and drainage systems would be necessary. The conveyance losses during
transportation due to leakage/seepage or otherwise would cause accumulation of water in
the subsurface system.
The city is devoid of any perennial natural drainage system. The Jojri River,
which routes through the outer periphery about 10 km away from the city, is the only
natural ephemeral stream that flows intermittently during the monsoon. There are number
of small water bodies inside the city, some of those are manmade, constructed long back
to store monsoon water; while others are formed by the natural topographic depressions.
Prior to the current arrangement of water supply from the transit storage in the
Kailana lake, Takhatsagar Reservoir and Umaid sagar Reservoir , which are fed from the
IGNP linked Rajiv Gandhi Lift canal, the domestic and municipal water supply including
drinking water was met from hand-pumps, tube-wells, step-wells, baories, and from
surface water storage in the city. After the existing arrangement of the water supply
scheme to the city in the year 1996-‘97, i.e., feeding the Kailana Lake and Takhatsagar
Reservoir from the IGNP linked Rajiv Gandhi Lift canal and transferring the water from
the Kailana and Takhatsagar Reservoir largely by pumping and partially by gravity flow
for treatment of water and then supplying to the city, almost all the previous provisions of
water supply from the groundwater storages have been put into hold. The enhanced-
urban-water-supply-return-flows and the seepage from the water bodies would cause rise
in ground water level. The existing hydrogeological condition and water usages
transformation have consequently given rise to the problem of groundwater level
increase. The resulting impact of rise in groundwater level over the years was such that a
number of pockets and stretches in the city area have been experiencing waterlogged
conditions. The rising trend of groundwater level was noticed since the year 1997, and
over the years the situation had continued to be so aggravated that considerable area,
mainly in the old city area had come under the grim of water logging conditions. Many
depressed land surfaces along the built up groundwater flow direction had also
experienced the surface water logging conditions. To cope up with the situation of rising
groundwater levels and water logging conditions, the State Groundwater Department has
made several pumping stations to lower the groundwater levels in a number of affected
areas. Even people having basement floors in the affected areas have also regularly
drained out the accumulated groundwater in the basement floors. Pumping out of
To find the sources and causes of groundwater rise in the Jodhpur city, a number
of organizations investigated the problem; these organizations are: GWD, Jodhpur
(1998), joint study by GWD, PHED, MBM Engineering College, and RRSSC, Jodhpur
(1999); BARC, Mumbai (2000); GSI, Western Region (2001); CGWB, Western Region
(2001); GWD, Jodhpur (2006); RRSSC, Jodhpur (2007), etc. Different and mixed
opinions have been stated regarding the causes of groundwater rise in the above
investigations. For example, GWD et al., (1999) indicated the cause of groundwater rise
is due to seepage from the used water and the Kailana-Takhatsagar Reservoir ; BARC
(2000) indicated the reason to be direct seepage from the Kailana-Takhatsagar Reservoir
through rock fractured aligned to the city areas or seepage from pipelines carrying lake
water; CGWB (2001) explained the reasons as the seepage from the open drains and
sewerage lines, RRSSC (2007) strongly expressed the source as the Kailana-Takhatsagar
Reservoir and the cause of rise is due to seepage through lineaments/joints present in the
rocks. Different investigators had suggested several remedial measures to contain the rise
in groundwater level. However, the problem resolving issues are yet to be focused.
In order to find the exact source and causes of groundwater level rise in the
Jodhpur city and to develop an appropriate management plan to revert back the rising
trend or containing the groundwater at a safe level, the task has been entrusted to the
National Institute of Hydrology, Roorkee with the following objectives:
Based on the preliminary analyses of the pertinent data supplied by GWD, PHED,
Jodhpur and collected till December, 2009, an interim report with the following prima-
facie observations had been provided in the month of April, 2010:
¾ The Kailana-Takhatsagar Reservoir and Jodhpur city are located in two different
geologic formations. The geological formation on which Kailana-Takhatsagar
Reservoir is located is of Malani group having low hydraulic conductivity.
Therefore, the seepage loss from the lake will not be very significant.
¾ The source of water causing water logging is getting originated locally in the
problematic area.
¾ The seepage from the Kailana-Takhatsagar Reservoir does not enter to the
waterlogged area.
¾ A lineament is noted orienting towards the waterlogged area. Near to the Kailana
Lake the connection of this lineament with the lake is not clear. If at all the flow
from the Kaikana Lake to the waterlogged area is occurring through this
lineament, it can be prevented by intercepting wells.
• Two to three lines of drainage wells each line consisting of few staggered drainage
wells along equipotential lines will be required to restrict water logging.
• The pumping rates and number of wells can only be ascertained using groundwater
modelling.
This report is the final outcome of the detailed analysis of different hydrologic and
hydro-geologic components and groundwater flow modeling. In order to analyze the data,
all spatially varying databases, such as locations, have been first geo-referenced with
reference to the geographic coordinate system (WGS-1984) by their latitudes and
longitudes. Therefore, the analyzed data and results can easily be verified with the field
conditions.
The report brings out the results of different remedial options, and suggests the
suitable one, fitting to the field conditions. The analyses in the report primarily deal with
the following:
(i) Detailed description of the problematic and the study area.
(ii) Geological formations and aquifer characterization.
(iii) Hydro-meteorological data and analysis.
(iv) Demography and water requirement,
(v) Stage-area-capacity data of the Kailana-Takhatsagar Reservoir.
(vi) Inflow-outflow data analysis of the Kailana-Takhatsagar Reservoir.
(vii) Groundwater data analysis.
(viii) Sewerage and city’s drainage data analysis.
(ix) Water balance of the study area.
(x) Aquifer parameters estimation.
(xi) Groundwater quality data and analysis.
(xii) Discritization of the study area, and input data for groundwater modeling.
(xiii) Modeling scenario for different remedial options.
(xiv) Remedial options and groundwater management plan.
***************
2.1 Geography
The Jodhpur city located in the Thar Desert is the second largest and one of the
fastest growing urban areas in the state of Rajasthan even though it has the characteristics
of arid, hot and dry climatic conditions. Presently, the city is on an area of about 76 sq.
km located between latitudes 26015’ N to 260 20’ N and longitudes 730 0’ E to 730 4’ E,
and it is sprawled mainly towards south-east, south and south-west directions (Figure
2.1). The old walled city area, which is thickly populated, is located on the hill slope area
and on the base of the fort hill ridge, whose surface topography is of undulating type. The
sloping area gradually turns to plain alluvial terrain towards south, east and south-west
sides. The north-western side of the city is largely comprised of hillocks covered with
shale and sandstone. The city area has been sprawled partially on the Malani group of
rock formations, mainly the old city area and the major part has been expanded on the
Quaternary Alluvium formations. As such, there is no river/stream passing through the
city except the ephemeral Jojri River, that routes through the outer periphery about 10 km
away from the city. There are number of small water bodies inside the city, some of those
are man-made, constructed long back to store monsoon water; while others are formed by
the natural topographic depressions. The Kailana and Takhatsagar Reservoir , which are
two naturally formed cascading type geological faults providing large storage volume
aligned towards north-south west direction, are located about 5 km away in the western
side of the city. The Kailana and the Takhatsagar Reservoir are located in a single
geological unit and form one continuous depression storage. For operational benefits,
these have been divided into two separate storage units by constructing an embankment
in-between them. The upper part is known as the Kailana Lake and the lower part is
named as Takhatsagar Reservoir. The Kailana Lake has larger water spread area than the
Takhatsagar Reservoir. The Takhatsagar Reservoir has been connected to the Kailana
Lake through gate operated pipes concealed inside the embankment at different operating
depths of water in the Kailana Lake. The flow from the Kailana Lake is also diverted
through open channel constructed at a particular height on the right bank of the Kailana
Lake. The Takhatsagar Reservoir storage has been created by constructing a dam across
the natural depression. The Kailana water level is maintained at higher level than the
Takhatsagar Reservoir level so as to enable the Kailana water to flow to Takhatsagar
Reservoir by gravity flow. The storage capacity of the Kailana Lake is 4.814MCM at its
maximum elevation level of 273.7 m, and that of the Takhatsagar Reservoir is 6.524
Jodhpur City
Figure 2.1: A geo-referenced map of the Jodhpur city showing its sprawled area and
water bodies in and around the city.
Figure 2.2: Topographic map of Jodhpur city showing topographical contour lines,
Kailana & Takhatsagar Reservoirs including locations of some of the
important places, and contours of water table up to a depth of 3m and 5 m
below the ground surface.
Figure 2.4: 3-D Digital Elevation Map (DEM) of the Jodhpur area showing location
of Kailana and Takhatsagar and area in which water table lies up to 3 m
and 5 m below ground surface.
Study on Rising Groundwater Table in Jodhpur City - NIH 10
The level of topography within the city area varies between 202 m and 360 m
above msl with a small stretch near to the fort area having higher elevation with peak
value of 310 m, whereas most of the area in the city is largely flat terrain having elevation
below 250 m above msl. The Kailana Lake and the Takhatsagar Reservoir, which have
the minimum bottom elevation of 254 m and 251 m above msl, have the maximum water
elevation level of 273.7 m and 269.75 m respectively, above msl. Grossly, the bottom
elevations of the Kailana Lake and the Takhatsagar Reservoir are above the elevation of
the most part of the Jodhpur city. The storage conditions in the the Kailana Lake and the
Takhatsagar Reservoir, apparently indicate that the Kailana Lake and the Takhatsagar
Reservoir have the potential of groundwater flow movement to the city area provided the
areas are hydraulically connected. However, the rising groundwater level and the
waterlogged conditions in the city were observed along the extreme northern side in the
old city area down below the Jodhpur fort area where the area has ground elevation
between 230m and 310m, and the problematic area has an alignment almost at the same
horizontal line if drawn from the Kailana and Takhatsagar Reservoir, which are separated
by two distinct watersheds divider by a distance of about 5 km from the city area
(Figures 2.2, 2.3, and 2.4). In Figures 2.2, 2.3 and 2.4, two contour lines of groundwater
table; one indicating depth at 3 m below ground level (bgl) and other one indicating depth
at 5m bgl have been shown. The orientation of the groundwater contour lines in the
problematic area, which have been developed from the groundwater level data, indicates
sprawling tendency towards south-easterly and south-westerly direction. The positional
setups of the Kailana and Takhatsagar Reservoir, and the problematic area do not show
any signature of homogeneous hydraulic connectivity rather has disjointed formations. If
at all, both the setups have any linkage between them; that may be possible by sub-
surface conduit created by rock fractures and lineaments.
(i) Whether the Kailana Lake and the Takhatsagar Reservoir could be the source for
cause of rise in groundwater level in the area facing the water logging
conditions?
(ii) If yes, how? And if no, then what are other possible causes of such rise in
groundwater level?
(iii) What are the remedies, and how to contain the rising groundwater level in a
sustainable manner?
***************
The Malani suite of Igneous rocks is comprised of grey, buff and brown colored
volcanic flows of Rhyolite. These rocks constitute the basement rocks in the area. The
volcanic rocks are exposed as bold ridges with respect to regional ground level and
separate Jodhpur group of rocks into two distinct outcrops. These rocks are exposed
mainly in the western and northern part of the area.
From the analysis of satellite imageries, it is noted that Kailana and Takhatsagar
area is located on the Rhyolite formation, whereas Jodhpur city is mainly located on the
sandstone and shale formation (Figure 3.1). These two formations have different hydro-
and cannot be considered as a single system. Rhyolite is relatively impervious, whereas
sandstone is reasonably pervious.
The litholog data of 93 bore wells (details are given in Table-A3.1 in the
Annexure), supplied by the Groundwater Department, Jodhpur, have been analyzed using
ROCKWORKS 2006 software. The geo-referenced locations of the borelogs are shown
in Figure 3.2. From the analysis of the lithologs data, it is observed that the Jodhpur city
area is characterized by a number of flat-topped hills trending in N-S or NE –NW
direction and are primarily composed of Shales, Sandstones and Rhyolites (Figure 3.3),
while the plain terrain represents the Quaternary alluvium formations at the top, which is
formed by weathering of Rhyolites and Sandstones. The Quaternary alluvium formation
varies in thickness from few centimeters to about 75 m in depth. The alluvium formation
is underlain mainly by the Sandstone and Shale of Jodhpur group (Figure 3.4), which is
further underlain by the Rhyolites.
Figure 3.3 : Stratigraphic model of the geological formations of the Jodhpur city.
Figure 3.5: Geological formations of Jodhpur area along line B-B’ of Figure-3.2.
Figure 3.6: Geological formations of Jodhpur area along lines C-C’ of Figure-3.2.
It can be seen from Figure 3.1 that, there are number of faults and lineaments
mostly oriented towards N-E to S-W direction and few from North to south direction.
There is only one fault and lineament that is aligned in west to east direction
perpendicular to the Kailana-Takhatsagar Reservoir and extended to the Jodhpur city
close to the problematic area, apparently showing not exactly connected near to the
proximity of the Kailana-Takhatsagar Reservoir, but are separated by two distinct
geological formations; one is comprised of Malani Group of rocks on which Kailana-
Takhatsagar Reservoir is situated, and the other one is Quaternary alluvium formation
encompassing the problematic area. However, whether these two formations are
connected to each other through fractures can be studied by hydraulic analysis of flow
in addition to the geophysical survey.
The geological formations underneath the Jodhpur city can be seen in details from
the cross-sectional views along different directions. The cross-sectional view along A-
A’ direction of Figure 3.2, that indicates formations along N-S direction, is shown in
Figure 3.4. It can be seen from Figure 3.4 that the top formation is of Quaternary
alluvium having thickness varying from about 7 m along northern side to a maximum
of about 75 m as one move along south direction. The Quaternary alluvium formation
is underlain by a thin layer of Sandstone formation, which is further underlain
respectively by; Shale, Sandstone, and finally by Rhyolite of varying thickness. As
such, the borelogs data representing formations up to a depth of 75 m below ground
surface do not represent existence of any other aquifer at larger depth except the one at
The cross-sectional view along B-B’ direction of Figure 3.2, that indicates
formations along NW towards SE direction, is shown in Figure 3.5. It can be seen from
Figure 3.5 that the top formation is of Quaternary alluvium having thickness varying
from zero near the NW side to a maximum of about 21 m as one move to the SE
direction. The Quaternary alluvium formation is underlain by a thin layer of Sandstone
formation, which is further underlain respectively by; Shale, Sandstone, and finally by
Rhyolite of varying thickness. The general slope of the alluvium and Shale formations
are towards south-east direction; whereas the thickness of the Sandstone formation
gradually reduces in the south-east direction. The formations along this direction also
do not indicate presence of any deeper aquifer underneath the Rhyolite formation.
The cross-sectional view along C-C’ direction of Figure 3.2, that indicates
formations along NE towards SW direction, is shown in Figure 3.6. It can be seen from
Figure 3.6 that the top formation is of Quaternary alluvium having by and large
relatively less thickness varying from zero along NE direction to a maximum of about 6
m towards the SW direction. The Quaternary alluvium formation is underlain by a thin
layer of Shale formation, which is further underlain respectively by; Sandstone, and
finally by Rhyolite of varying thickness. The general slope of the alluvium and
Rhyolite formations are towards south-east; whose thickness varies from location to
location. The formations along this direction also do not indicate presence of any
deeper aquifer underneath the Rhyolite formation.
The cross-sectional view along D-D’ direction of Figure 3.2, that indicates
formations along W to E direction, is shown in Figure 3.7. It can be seen from Figure
3.7 that the top formation is of Quaternary alluvium formation having zero thickness
near to the western side and gradually increases to a large thickness of about 22 m as
one move towards east direction. The Quaternary alluvium formation is underlain by a
very thin layer of Sandstone formation that disappears in many locations, which is
further underlain respectively by; Share, Sandstone, and finally by Rhyolite of varying
****************
4.1 Rainfall
Thirty six years (1971-2006) daily rainfall data of the Jodhpur city obtained from
the IMD showed that average annual rainfall of the city is 378 mm with minimum of 91
mm in the year 2002 and maximum of 821 mm in the year 1990 (Figure 4.1), while the
last 25 years (1982-2006) annual rainfall data indicate an average annual rainfall of 357
mm (Figure 4.2). The main rainfall months are July and August. July is the most wetted
month, followed by August. About 68% and 86% of the annual rainfall occurs
respectively during July and August during the monsoon months (June through
September). The number of intensified rainy days in a year is reported to be about 15
days. Distribution of daily rainfall pattern in a year, analyzed based on the data of last 36
years (1971-2006) is shown in Figure 4.3. The variation of annual rainfalls in the
Jodhpur city during the year 1978-2006 is shown in Figure 4.3.
900
821
800
691
673
650
700
577
572
Rainfall (mm)
554
538
600
527
526
521
451
500
396
389
386
374
359
357
400
322
304
283
279
276
273
269
250
300
233
232
221
210
194
197
182
175
165
200
91
100
0
1971
1972
1973
1974
1975
1976
1977
1978
1979
1980
1981
1982
1983
1984
1985
1986
1987
1988
1989
1990
1991
1992
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
Year
Figure 4.1: Annual variation of rainfall in Jodhpur area during the last 36 years for the
period 1971-2006.
The average maximum temperature and the average minimum temperature during
summer are reported to be 42.2°C, and 27.3°C respectively, and during winter the
average maximum temperature is 27.5°C and the average minimum temperature is 9.5°C.
Figure 4.2: Annual variation of rainfall in Jodhpur area during the last 25 years for the
period 1982-2006
Figure 4.3: Distribution of daily rainfall pattern in Jodhpur area during the last 28 years
for the period 1978-2006
4.2 Evaporation
Twenty five years (1978-2002) monthly Pan-A evaporation data of the Jodhpur
station (station no. 42339) obtained from IMD shown in Table-4.1 and Figure 4.4
16
Monthly average daily evaporation
14 Yearly average daily evaporation
Evaporation (mm/day)
12
10
8
6
4
2
0
Jan Feb. March April May June July Aug. Sept. oct. Nov. Dec.
Month
Figure 4.4: Average monthly variation of water surface evaporation in the Jodhpur area.
Table-4.1: Average monthly evaporation from Class 'A' Pan considering 25 years (1978-
2002) IMD data for Jodhpur city
5.1 Demography
Because of its historical and commercial importance, strategic location and
infrastructural facilities, Jodhpur city has grown very fast in the past and is also growing
continuously. It being the Headquarters of the Western Air Command, it has a strategic
importance for the country. The Aerodrome is located in the city in the south-eastern
side. The city is known to be the second largest city in Rajasthan after Jaipur.
The city that had a population of about 3,65,000 in the year 1971, was expected to
have a population of 11,00,000 by the year 2010. According to the censuses of the years
1981, 1991, and 2001, the population of the Jodhpur city had been reported to be
5,05,000; 6,66,280 and 8,51,051, respectively. The population censuses of previous four
decades indicated a growth rate of 3.21% per year. The growth rate of 3.21% per year
gives the projected population of about 11,08,950 at the end of the year 2010, and
12,74,830 by the year 2015. The variation of population growth over different years since
the year 1971 is shown in Figure 5.1. The fitted population data give an exponential rise
satisfying the following equation:
P = 5.1 x 10 −19 Exp (0.02788 X ) ……………. (5.1)
in which, P is the population in a particular calendar year; X is the calendar year in its full
term, e.g.,1990, for which population is to be calculated. Equation (5.1) fits to the data
with a correlation coefficient, R2 = 0.9976.
1200000
0.02788 x
1000000
y = 5.1 E-19 e
2
R = 0.9976
Population
800000
600000
400000
Population growth rate = 3.21% per year
200000
0
1965 1970 1975 1980 1985 1990 1995 2000 2005 2010 2015
Year
Figure 5.2 : Municipal wards encompassed in the study area in the Jodhpur city
Estimated Population
Position Approx. area
Wards Population Total area in the area
in study within study
No. (2001) (sq.km) encompassed by
area area (sq. km)
study area (2001)
1 Partially 11,857 39.79 3.75 1117
2 Partially 21,231 20.90 4.81 4883
3 Partially 22,339 3.43 0.79 5138
4 Outside 40,590 17.23 0.00
5 Partially 17,153 1.38 0.32 3945
6 Fully 12,905 1.60 1.60 12905
7 Partially 22,410 15.98 8.75 12268
8 Fully 14,313 0.30 0.30 14313
9 Fully 10,791 0.60 0.60 10791
10 Fully 15,842 1.10 1.10 15801
11 Fully 16,941 2.55 2.55 16910
12 Partially 9,556 0.85 0.20 2198
13 Partially 8,669 0.39 0.20 4437
14 Fully 11,756 0.12 0.12 11935
15 Fully 13,488 1.14 1.14 13519
16 Fully 8,841 0.39 0.39 8783
17 Fully 21,734 0.86 0.86 21732
18 Fully 11,012 0.50 0.50 10941
19 Fully 10,513 2.94 2.94 10498
20 Fully 12,586 1.00 1.00 12552
21 Fully 8,021 0.73 0.73 8015
22 Fully 11,555 0.80 0.80 11564
23 Fully 10,579 0.19 0.19 10447
24 Fully 9,063 0.16 0.16 9019
25 Fully 11,969 0.26 0.26 11816
26 Fully 13,117 0.23 0.23 13178
27 Fully 10,623 0.20 0.20 10497
28 Partially 11,702 0.91 0.25 3230
29 Partially 12,061 0.38 0.22 6935
30 Partially 7,115 0.40 0.35 6240
31 Fully 9,289 0.26 0.26 9115
32 Fully 8,104 0.36 0.36 8076
33 Fully 12,297 0.40 0.40 12326
34 Fully 8,161 0.49 0.49 8233
35 Fully 10,995 0.39 0.39 10863
36 Fully 12,492 0.66 0.66 12411
37 Fully 12,970 0.07 0.07 12251
38 Fully 10,282 0.40 0.40 10306
39 Fully 7,065 0.66 0.66 7014
40 Partially 10,499 1.74 1.55 9362
Table 5.2: Projected population of each ward encompassed by the study area; calculated
on pro-rata basis with growth rate of 3.21% per year.
Out of the total area of 303.44 sq. km of 60 wards covered by the Jodhpur city,
the area encompassed by 51 wards (34 fully, and 17 partially) within the study area is
nearly 55.55 sq. km. The total population in the 60 wards as per the year 2001 was about
8,41,505, out of which nearly 4,80,731 ( estimated on pro-rata basis in accordance with
the proportion of area of each ward encompassed by the study area) is estimated to be the
population in the study area. The ward-wise details of area and population can be seen in
The water supply requirements to the population in the city area for their various
designated uses are primarily met by the Jodhpur Municipality through water supply &
distribution systems managed by the Public Health Engineering Department (PHED),
Jodhpur. As such, no specific guidelines towards quantity of water being supplied are
followed; it is based on thumb rule in accordance to the supply-demand norm. As the city
has expanded over the years, more areas and more number of people have been brought
under the coverage of the water supply distribution; and the water supply from the intake
has been increased accordingly on pro-rata basis. The city, by and large, has a large
network of water supply distribution lines not connecting to all households. In some
areas, the distribution lines are very old, constructed long back, and in some parts
people/houses are not connected to organized water supply lines, use water from the
hydrants constructed along the roadsides or from the nearby water supply points. Owing
to the unregulated and uncontrolled supply of water, people in the city area practice water
uses round the clock without understanding to its far reaching consequences.
Water requirement in lac gallon/day Water actually supplied to the city area in lac gallon/day
600
Water Supplied to the city
500
400
300
200
100
0
1994
1996
1998
2000
2002
2004
2006
2008
2010
Year
Figure 5.3 : Water requirements as per MOWH norms and actually supplied per day to
the Jodhpur city in different years
As such, no specific data were available on the area-wise water supplied to the
different locations in the city. Based on the ward-wise population and their projections in
different years, a conservative estimate of water requirement (MOWH norm) and water
supplied in different years on pro-rata basis has been made, and the estimated values are
given in Table- 5.4. For projection of ward-wise population, the criterion of uniform
growth rate of 3.21% per year on the basic data of population (year, 2001) in the
respective ward has been considered. The quantities of water supplied in different years
to the respective wards are also considered for further analysis. As the sum total of water
distribution in the different wards satisfies the total quantity of water supplied in the city
area, hence, it assures the overall water balance.
*******************
0.8
Area, A (Sq.KM)
0.6
0.4
A = 0.0064 V3 - 0.0496 V2 + 0.2756 V + 0.0224
0.2 R2 = 0.9991
0
0 1 2 3 4 5
Capacity, V (MCM)
4
Capacity, V (MCM)
1 2
V = 2.585 A + 3.791 A - 0.102
R² = 0.998
0
0 0.2 0.4 0.6 0.8
Area, A (Sq.KM)
The functional relationships that hold good to estimate the capacities for varying
stages and capacities for different wetted areas of the Kailana Lake are as follow:
The Takhatsagar Reservoir , which is fed by the Kailana Lake through regulated
gates, is situated at an altitude of 241 m above msl between the latitude ranges from 26°
16´ 51˝ N to 26° 17´ 58˝ N and longitude ranges from 72° 57´ 51˝ E to 72° 58´ 25˝E . The
Takhatsagar Reservoir has been created by constructing a masonry dam across the
depressed valley formed by natural geological faults. In fact, the Kailana-Takhatsagar
Reservoir is situated on the same geological faults, but separated by a man-made earthen
dam in between. The bed levels as well as the water surface levels of the Kailana and
Takhatsagar Reservoirs are different, and the bed level and water surface level of the
Kailana Lake are at considerably higher elevation than those of the Takhatsagar
Reservoir. The full level of the Takhatsagar Reservoir is 269.75 m above msl. The area
and capacity of the Takhatsagar Reservoir at its full level are estimated to be 0.601 sq.
km and 6.523 mcm, respectively. The stage-capacity and the area-capacity curve of the
including their functional relationships are given in Figures 6.4, 6.5 and 6.6, respectively
and the stage-area-capacity table is given in Table-A6.2 in the annexure. The water is
The functional relationships that hold good to estimate the capacities for varying
stages and capacities for different wetted areas of the Takhatsagar Reservoir are as
follow:
(a) Stage-Capacity relationship:
V = 0.019 h1.89 ………………………….. (6.4)
(c) Area-Capacity relationship
V = 14.75 A 2 + 1.673 A − 0.12 .....………….. (6.5)
or
0.70
0.60
Area, A (Sq.KM)
0.50
0.40
0.30
2
0.20 A = ‐0.006 V + 0.129 V + 0.071
R² = 0.972
0.10
0.00
0.000 1.000 2.000 3.000 4.000 5.000 6.000 7.000 8.000 9.000
Capacity, V (MCM)
12
9
Capacity, V (MCM)
3 2
V = 14.75 A + 1.673 A - 0.120
R² = 0.996
0
0.0 0.2 0.4 0.6 0.8
-3
Area, A (Sq.KM)
All three equations fitted to the field data have correlation coefficients, R2 >
0.972, and hence can satisfactorily be used for estimation of storage capacity of the
Takhatsagar Reservoir from known stage or water spread area. Alternately, when the
storage capacity of the Reservoir is known from the measurement of stages, the
corresponding water spread area can be computed using equation (6.6).
The functional relationships that hold good to estimate the capacities for varying
stages of the Umaidsagar are as follow:
The Kailana-Takhatsagar Reservoir receives its inflow from the Indira Gandhi
Canal through Rajiv Gandhi Lift Canal by pumping at PS-8 pumping station. There are six
pumps to lift water from the Indira Gandhi canal; out of which normally, 3-4 pumps
operate almost round the clock. There is no gate that regulates the flow between the
pumped out water till it reaches to the Kailana Lake. The regulation of flow is based on
some pre-decided thumb rules, which are based on the supply of water to the city area and
to the downstream reservoirs, namely Takhatsagar and Umaidsagar; i.e., a constant rate of
inflow is maintained, and as and when the water level in the Kailana Lake goes down
below a specific level, the inflow rate is increased by operating more number of pump. A
line diagram of Canal-Lake-Reservoir-City water supply management practice is shown in
Figure 7.1.
Indira Gandhi Canal
R
A
J
I
V Inflow through Pumping
G
A
N
D
H
I
C
A
N
A
L
TAKHATSAGAR
RESERVOIR
Gravity Flow
140.52
0.0997 x
y = 58.802 e
119.14
160
106.17
Inflow to Kailana Lake
2
140 R = 0.9318
90.21
120
81.10
80.88
(MCM/day)
72.04
68.31
100
80
60
40
20
0
2002 2003 2004 2005 2006 2007 2008 2009
Year
Figure 7.2 : Variation of average inflow per day in different years to the Kailana Lake
from the Rajiv Gandhi Lift Canal
The water balance equation connecting the hydrological components involved in the
Kailana Lake shown in Figure 7.3 can be written as:
in which, Qik is the quantity of inflow from the Rajiv Gandhi lift canal to the Kailana
Lake in L3T-1; Qrk is the quantity of rainfall over the Kailana Lake in L3T-1; Qsrk is the
The components, Qik is known from the daily inflow data of the Kailana Lake
obtained from the PHED; Qrk is known from the daily rainfall data obtained from the
IMD as shown in Figures 4.1, 4.2 & 4.3 in section 4; Qsrk can be estimated from the
rainfall data using suitable rainfall-runoff relationship; Qok1 and Qok2 although are
regulated flow, however, are not known in exact quantities as to how these are regulated
and hence are to be ascertained; Qevk can be estimated externally from the evaporation rate
data given in Figure 4.4 and Table 4.1 in section 4.0 for the corresponding water spread
area at a particular head of water above the lake bottom; Qgwk can be estimated for a
particular head of water above the lake bottom using seepage theory; ∆t can be chosen
suitably for which the water balance is to be performed; and ∆H is known from the
measurement of heads of water in the lake = H2 – H1, in which H1 is the head of water at
time T1, and H2 is the head of water at time T2; and ∆t = T2 – T1.
Without GW component
With GW component
0.20
0.18
0.16
0.14
Flow (MCM/day)
0.12
0.10
0.08
0.06
0.04
0.02
0.00
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
16
31
46
61
76
91
1
Day
Figure 7.4 : Variation of supply of water from Kailana Lake to Takhatsagar Reservoir
during the year 2003 staring from January,1 –December, 31 (red color line
corresponds to without seepage component; blue color line corresponds to
with seepage component).
With GW component
0.35
0.30
0.25
Flow (MCM/day)
0.20
0.15
0.10
0.05
0.00
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
16
31
46
61
76
91
1
Day
Figure 7.5 : Variation of supply of water from Kailana Lake to Takhatsagar Reservoir
during the year 2005 staring from January,1 –December, 31 (red color line
corresponds to without seepage component; blue color line corresponds to
with seepage component).
Without GW component
With GW component
0.40
0.35
0.30
Flow (MCM/day)
0.25
0.20
0.15
0.10
0.05
0.00
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
16
31
46
61
76
91
1
Day
Figure 7.6 : Variation of supply of water from Kailana Lake to Takhatsagar Reservoir
during the year 2008 staring from January,1 –December, 31 (red color line
corresponds to without seepage component; blue color line corresponds to
with seepage component).
With GW component
0.50
0.45
0.40
0.35
Flow (MCM/day)
0.30
0.25
0.20
0.15
0.10
0.05
0.00
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
16
31
46
61
76
91
1
Day
Figure 7.7 : Variation of supply of water from Kailana Lake to Takhatsagar Reservoir
during the year 2009 staring from January,1 –December, 31 (red color line
corresponds to without seepage component; blue color line corresponds to
with seepage component).
It can be seen from Figs. 7.4-7.7 that, the supplies of water from the Kailana to the
Takhatsagar Reservoir varied during a year, and even year to year. In some of the days in
a year, virtually there was no supply of water from the Kailana to the Takhatsagar
Reservoir. The variation of seepage below the lakebed has been estimated to be varied
very marginally following the same trend as that of the outflow pattern from the Kailana
Lake, i.e., higher the rate of outflow from the Kailana Lake to the Takhatsagar Reservoir,
larger is the rate of seepage below the lakebed. It can, in other way, be explained in the
following way: large inflow to the Takhatsagar Reservoir means increased outflows from
the accumulated storages in the Kailana Lake; larger storages designate higher potential of
head, and thereby more seepage through the lakebed. Moreover, larger the inflow rates
from the Rajiv Gandhi lift canal to the Kailana Lake, higher is the outflows from the
Kailana to the Takhatsagar Reservoir, which can be seen from Figs.7.8-7.11. Figs. 7.8-
7.11 show year-wise variation of inflows to the Kailana Lake from the Rajiv Gandhi lift
canal, supply of water from the Kailana Lake to the city, and the outflows from the
Kailana Lake to the Takhatsagar Reservoir. It can be seen from Figs.7.8-7.11 that, the
inflows to and outflows from the Kailana followed almost a similar trend except during
certain period, while the supply to the city area remained largely uniform in a year. The
The year-wise sum of each water balance component which is involved in the
water balance equation of the Kailana Lake is given in Table-7.1. It can be seen from
Table-7.1 that the differences between the change in storage estimated from the inflow-
outflow balance and the change in storage calculated independently are very minor. These
minor differences may be due to the errors in the fitting of the rating curve.
Table 7.1 : Water balance components of the Kailana Lake in different years.
Inflow
Change in Change in
Components Outflow Components (MCM) storage = storage
(MCM) ∑ Inflow - Calculated
Diversion Diversion
Year
Pumping from Kailana from Kailana
∑ Outflow ( MCM)
Ground ( MCM)
Inflow to Rainfall from to to
water Evaporation
Kailana Volume Kailana Takhatsagar Takhatsagar
Component
to City without GW with GW
Component Component
2003 69.31 0.38 25.03 45.86 34.27 11.56 1.48 - 2.68 -1.57
(- 2.65)*
2005 85.46 0.24 32.49 54.16 40.91 12.20 1.40 - 2.35 -1.70
(- 2.30)*
2008 125.94 0.10 49.13 78.26 64.78 13.50 1.47 - 2.82 -1.07
(- 2.84)*
2009 144.15 0.00 46.23 96.02 81.75 14.25 1.66 0.24 0.00
(0.26)*
(*) indicates change in storage with ground water component.
Diversion to Takhatsagar
0.25
0.20
Flow (MCM/day)
0.15
0.10
0.05
0.00
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
16
31
46
61
76
91
1
Figure 7.8 : Variation of water balance components of the Kailana Lake during the year
2003 staring from January,1 –December, 31.
Water balance of the Kailana Lake for the year 2005 without groundwater Inflow to Kailana
component Supply to City
Diversion to Takhatsagar
0.45
0.40
0.35
Flow (MCM/day)
0.30
0.25
0.20
0.15
0.10
0.05
0.00
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
16
31
46
61
76
91
1
Figure 7.9 : Variation of water balance components of the Kailana Lake during the year
2005 staring from January,1 –December, 31.
Diversion to Takhatsagar
0.50
0.45
0.40
Flow (MCM/day)
0.35
0.30
0.25
0.20
0.15
0.10
0.05
0.00
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
16
31
46
61
76
91
1
Figure 7.10 : Variation of water balance components of the Kailana Lake during the year
2008 staring from January,1 –December, 31.
Water balance of the Kailana Lake for the year 2009 without groundwater Inflow to Kailana
component Supply to City
Diversion to Takhatsagar
0.50
0.45
0.40
Flow (MCM/day)
0.35
0.30
0.25
0.20
0.15
0.10
0.05
0.00
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
16
31
46
61
76
91
1
Figure 7.11 : Variation of water balance components of the Kailana Lake during the
year 2009 staring from January,1 –December, 31.
The diverted outflow from the Kailana Lake to the Takhatsagar Reservoir through
gate operated pipes and regulated diversion canal, which flow under gravity, form the
inflow to the Takhatsagar Reservoir. The stored water in the Takhatsagar reservoir is
further diverted to the downstream Chopasni Filter House by gravity flow for supply to
some of the areas in the city after necessary treatment. The regulation of the flow from the
Takhatsagar Reservoir depends on the downstream requirement, which is decided based
on the gauge height in the reservoir. The Takhatsagar Reservoir has a small catchment
encompassed by its surrounding sloping topography. Whatever rainfall occurs on its
catchment area, some fraction of it in the form of runoff, in addition to the direct rainfall
on the water surface area also accumulates in the reservoir as its storage water. During the
course of accumulation of water in the Takhatsagar Reservoir, a component of the
accumulated water is lost as water surface evaporation, and another component may
seep/infiltrate-percolate below the Takhatsagar Reservoir bed to the underneath
aquifer/groundwater whose rate depends on the potential head, reservoir bed material;
geological formations & their hydraulic properties. A schematic diagram showing the
components involved in the water balance of the reservoir is given in Figure 7.12.
The water balance equation for the Takhatsagar reservoir can be written as:
[( Q it + Qrt + Qsrt ) (
− Qot + Qevt + Q gwt ) ] Δt = A(H ) ΔH …….. (7.4)
in which, Qit is the quantity of inflow from the Kailana Lake to the Takhatsagar Reservoir
in L3T-1; Qrt is the quantity of rainfall over the Takhatsagar Reservoir in L3T-1; Qsrt is the
quantity of surface runoff from the Takhatsagar catchment in L3T-1; Qot is the quantity of
outflow from the Takhatsagar Reservoir to the Chopasni Filter House through regulated
pipe flow in L3T-1; Qevt is the quantity of water surface evaporation from the Takhatsagar
Reservoir in L3T-1; Qgwt is the quantity of groundwater seeping (recharge) through the
Takhatsagar Reservoir bed in L3T-1; ∆t is the time step in unit of time(T); and ∆H is the
change in water level in the Takhatsagar Reservoir in unit of length (L), A(H ) is water
spread area which depends on reservoir stage.
The components, Qit is known from the daily diverted flow from the Kailana Lake;
Qrt is known from the daily rainfall data obtained from the IMD as shown in Figures 4.1,
4.2 and 4.3 in Section 4; Qsrt can be estimated from the rainfall data using suitable
rainfall-runoff relationship; Qot although is regulated flow, however, is not known in exact
quantities as to how this is regulated and hence is to be determined; Qevt can be estimated
externally from the evaporation rate data given in Figure 4.4 and Table 4.1 in Section 4.0
for the corresponding water spread area at a particular head of water above the reservoir
bottom; Qgwt can be estimated for a particular head of water above the reservoir bottom
using seepage theory; ∆t can be chosen suitably for the water balance; and ∆H is known
from the measurement of heads of water in the reservoir = H2 – H1, in which H1 is the
head of water at time T1, and H2 is the head of water at time T2; and ∆t = T2 – T1, A(H ) =
water spread corresponding to height, H .
The components, Qevt and Qgwt have been calculated following the similar
approach as explained in the case of Kailana water balance. As the daily outflows from
the Takhatsagar Reservoir to the Chopasni Filter House, Qct is not exactly known though
this is regulated flow; however, this quantity is computed by developing suitable rating
The flow of water from the Takhatsagar Reservoir to the Chopasni Filter House is
primarily regulated through the imbedded pipe of diameter 1m whose inlet point is placed
near to the bed of the reservoir at an elevation of 255.42 m above m.s.l (Fig.7.13). The
pipe is of about 80 m long aligned obliquely with a slope of 0.91 m vertical for a length of
80 m horizontal. After this slanting, the pipe line has moderate alignment up to the
Chopasni Filter House. In addition to this gravity pipe flow, there is an additional
diversion arrangement located at an elevation of 266.56 m (above msl) i.e., about 11 m
above the bottom pipe, that operates when water level in the Takhatsagar Reservoir
exceeds the level of 266.56 m.
Using the following formula of flow through pipe, the discharges of flow for
varying heads of water in the reservoir have been calculated.
Q1 = C d A 2 g (H 1 + h f ) ……………………………………… (7.5)
in which, Q1 is the discharge of flow through the pipe ( L3T-1 ); Cd is the coefficient of
discharge that ranges 0.60 to 0.75; A is the cross-sectional area of the pipe (L2); g is the
2
f L Vm
hf = ……………………………… (7.6)
2gd
in which, L is the length of the pipe = 80 m; Vm is the mean flow velocity in the pipe =
2 g H 1 ; d is the diameter of the pipe = 1.0 m.
As the diversion of flow through the upper pipe takes place when the water level in
the reservoir exceeds the level of 266.56 m, the combined affects of discharge of flows are
calculated for varying heads of water in the reservoir by taking into account both the flows
separately, and then added to have complete rating curve. The rating curve so develop is
shown in Figure 7.14.
3 2
8.0 Q = 0.0075 h - 0.18 h + 1.6615 h - 0.9045
2
R = 0.9658
6.0
4.0
2.0
0.0
10.36
17.37
0.00
0.50
1.00
1.52
2.44
3.05
3.66
4.57
6.10
7.62
8.84
9.75
Head (m)
Diversion of flow from the Takhatsagar to Chopasni Filter House during year 2003
0.80
0.50
0.40
0.30
0.20
0.10
0.00
16
31
46
61
76
91
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
1
Day
Figure 7.15 : Variation of supply of water from the Takhatsagar Reservoir to the
Chopasni Filter House during the year 2003 staring from January,1 –
December, 31 (red color line corresponds to without seepage component;
blue color line corresponds to with seepage component).
0.35
0.20
0.15
0.10
0.05
0.00
1
16
31
46
61
76
91
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
Day
Figure 7.16 : Variation of supply of water from the Takhatsagar Reservoir to the
Chopasni Filter House during the year 2005 staring from January,1 –
December, 31 (red color line corresponds to without seepage component;
blue color line corresponds to with seepage component).
Diversion of Flow from the Takhatsagar to Chopasni Filter House during year 2008
1.80
1.60
1.40 Without GW component
Flow (MCM/day)
1.00
0.80
0.60
0.40
0.20
0.00
1
16
31
46
61
76
91
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
Day
Figure 7.17: Variation of supply of water from the Takhatsagar Reservoir to the Chopasni
Filter House during the year 2008 staring from January,1 –December, 31
(red color line corresponds to without seepage component; blue color line
corresponds to with seepage component).
3.50
3.00
Without GW component
Flow (MCM/day)
2.50
With GW component
2.00
1.50
1.00
0.50
0.00
1
16
31
46
61
76
91
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
Day
Figure 7.18: Variation of supply of water from the Takhatsagar Reservoir to the Chopasni
Filter House during the year 2009 staring from January,1 –December, 31
(red color line corresponds to without seepage component; blue color line
corresponds to with seepage component).
It can be seen from Figs. 7.15-7.18 that, the supplies of water from the
Takhatsagar to the Chopasni Filter House varied during a year, and even year to year. In
some of the days in a year, virtually there was no supply of water from the Takhatsagar
Reservoir, and the variation of flows followed the similar pattern as that of the supplies
from the Kailana Lake to the Takhatsagar Reservoir. The variation of seepage below the
reservoir bed has been estimated to be varied very marginally following the same trend as
that of the outflow pattern from the Takhatsagar Reservoir, i.e., higher the rate of outflow
from the Takhatsagar Reservoir to the Chopasni Filter House, larger is the rate of seepage
below the reservoir bed.
0.80
Inflow to Takhatsagar
0.70
Takhatsagar to Chopasni Filter House without GW component
0.60 Takhatsagar to copasni with GW component
Flow (MCM/day)
0.50
0.40
0.30
0.20
0.10
0.00
1
16
31
46
61
76
91
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
Figure 7.19 : Variation of water balance components of the Takhatsagar Reservoir during
the year 2003 staring from January,1 –December, 31.
0.25
Flow (MCM/day)
0.20
0.15
0.10
0.05
0.00
1
16
31
46
61
76
91
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
Day ( Ye ar : 2005)
Figure 7.20 : Variation of water balance components of the Takhatsagar Reservoir during
the year 2005 staring from January,1 –December, 31.
1.80
1.60
Inflow to Takhatsagar
1.40
Takhatsagar to Chopasni Filter House without GW component
Flow (MCM/day)
1.20
Takhatsagar to copasni with GW component
1.00
0.80
0.60
0.40
0.20
0.00
1
16
31
46
61
76
91
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
Figure 7.21 : Variation of water balance components of the Takhatsagar Reservoir during
the year 2008 staring from January,1 –December, 31.
2.50
2.00
1.50
1.00
0.50
0.00
1
16
31
46
61
76
91
106
121
136
151
166
181
196
211
226
241
256
271
286
301
316
331
346
361
Day ( Year : 2009)
Figure 7.22 : Variation of water balance components of the Takhatsagar Reservoir during
the year 2009 staring from January,1 –December, 31.
******************
The groundwater levels data, supplied by the GWD, had been measured from the
ground surface at the respective locations, and all those levels thus represent the depth to
groundwater table. Identifying the latitude and longitude of each observation location by
GPS and making use of the ground level contour map, the groundwater level at each
observation point with respect to a common datum i.e., the mean sea level (MSL) is
obtained. After this conversion, considering the observed data to be consistent, the pre
monsoon and post monsoon water table contour maps for the years; 1996, 2000, 2004,
and 2008 have been prepared.
Reservoirs are prone to seepage depending upon the hydraulic conductivity of the
underlying geologic formation. The Kailana and the Takhatsagar Reservoirs, which
approximately conform to a rectangular strip of long length, may be prone to seepage.The
Presently, there are no observation wells in the close vicinity of the Kailana and
Takhatsagar Reservoirs. Therefore, contours of equal water table height above MSL have
been drawn for two limiting situations. In situation 1, it is envisaged that, the Kailana and
Takhatsagar Reservoirs are hydraulically connected to the aquifer under the Jodhpur city.
With this assumption, the reservoirs perform as a line of observation wells. Accordingly
contour lines incorporating the reservoir water level have been drawn. The parallel line
contours near the Kailana and the Takhatsagar Reservoir conform to this assumption.
This is the most favorable condition for seepage from the Kailana-Takhatsagar to flow to
the aquifer below the Jodhpur city to cause water logging.
Figure 8.2: Pre (red color) and post (black color) monsoon groundwater table contour
map for 2000 with the assumption that the Kailana-Takhatsagar is
hydraulically connected to the aquifer below the Jodhpur city ( Æ Indicates
groundwater flow direction).
Figure 8.4: Pre (red color) and post (black color) monsoon groundwater table contour
map for 2008 with the assumption that the Kailana-Takhatsagar is
hydraulically connected to the aquifer below the Jodhpur city( Æ Indicates
groundwater flow direction).
Figure 8.5 Pre (red color) and post (black color) monsoon groundwater table contour
map for 1996 with the assumption that the Kailana-Takhatsagar is not
hydraulically connected to the aquifer below the Jodhpur city( Æ Indicates
groundwater flow direction).
Figure 8.7 : Pre (red color) and post (black color) monsoon groundwater table contour
map for 2004 with the assumption that the Kailana-Takhatsagar is not
hydraulically connected to the aquifer below the Jodhpur city( Æ Indicates
groundwater flow direction).
In the year 2008, water level observations have been made at more number of
observation points. Therefore, the observed data in the year 2008 have been considered for
comparison of the pre and post monsoon groundwater table contours. The comparison of
the pre and post monsoon groundwater table contours maps for the cases of Kailana-
Takhatsagar Reservoir hydraulically connected and hydraulically not connected with the
aquifer underneath of the Jodhpur city are shown in Figs.8.9 and 8.10, respectively. It can
be seen from the Figs. 8.9 and 8.10 that the groundwater table contour lines in both the
cases follow the similar trend with variation in the level differences between the pre and
post monsoon contours, viz., at a particular location the elevation of the post monsoon
groundwater table has higher elevation than that of the pre-monsoon groundwater table. It
indicates that the observed groundwater data and the generated groundwater contour maps
are consistent, and the aquifer below the Jodhpur city area has a good response to the
monsoon rainfall recharge.
A 3-dimensional view of the pre and post monsoon groundwater table elevation
map for the year 2008 has been shown in Figs. 8.11 and 8.12, respectively. These figures
clearly depict the groundwater table mounds in the uppermost region of the study area.
The difference between highest (in the northern part) and the lowest (in the southernmost
part) groundwater table elevetation in the study area is of the order of about 90-95 m.
7
Area between level 0 - 3m
6 5.8101
Area between level 3 - 5m
4 3.735
3.074 3.028
3 2.7821
2.208
2.0586 1.95
2
1.527
0.9816 1.077
1
0
1996 2000 2004 2008
Year
The following inference is drawn analyzing the water table contour maps:
(a) The directions of groundwater flow during all the years from 1996 to 2008 (Figs.
8.1-8.8), the pre as well as post monsoon seasons, are observed to be broadly
from the north to the south and southeasterly direction with some minor
deviations at local pockets,
(c) The possibility of seepage from the Kailana-Takhatsagar Reservoir entering the
waterlogged area is investigated making use of the contour maps. The water table
contours shown in Figs 8.5-8.8 are based on the assumption that Kailana-
Takhatsagar Reservoir is not hydraulically connected with the aquifer underneath
the Jodhpur city. As seen from these figures, the water table contours in the
waterlogged area are having higher values than those of the contours which are in
the vicinity of the region where the Kailana-Takhatsagar is located. As the
waterlogged area is located on a region having higher water table contours than
those in the region near Kailana-Takhatsagar, from hydraulic principle, the
waterlogged area will not receive any flow from the southwest side i.e., from the
reservoir side, where hydraulic head is less than that in the waterlogged area. So it
is inferred that groundwater flow from the Kailana-Takhatsagar side is not
causing water logging in the waterlogged area,
(e) From the configurations of the flow nets (Figs.8.1-8.8), it is concluded that the
Kailana-Takhatsagar Reservoir being located in the Malani group of rocks, and
the Malani group of rocks have low permeability; the seepage losses from the lake
will be very small, which can be quantified from a water balance study. This
seepage is entering towards the south-east region in which the water table
contours are having less value than those in the waterlogged area. From the
consideration of hydraulic principle (flow takes place from region with higher
hydraulic head to region with lower hydraulic head), and from the consideration
of the direction of flow, it can be postulated that the seepage from the reservoir is
not entering to the waterlogged areas.
(f) Analysis of evolution of depth to water table below ground surface in the area
affected progressively due to water logging indicates that area getting
waterlogged has increased every year, and dynamic equilibrium has not been
reached yet.
********************
Before the year 1996-’97, the water supply of the Jodhpur city for domestic and
municipal water including drinking water used to meet partly from about 1962 numbers
of hand pumps, 109 numbers of tube wells/open wells, 4 numbers of step wells &
baories, surface water bodies in the city and from the monsoon storage available in the
Kailana-Takhatsagar Reservoir. After the year 1996-‘97, all the groundwater based
supplies have been put into hold, and the water supply to the city has been fully switched
over to meet from the Kailana Lake and Takhatsagar Reservoir through continuous
feeding from the IGNP linked Rajiv Gandhi Lift canal by partly pumping and partly by
gravity flow, after treatment. Water from the Kailana-Takhatsagar Reservoir to the city
area is transported through large diameter pipes. The quantity of water supplied from the
Kailana-Takhatsagar Reservoir in different years during the period 1991-2009 is shown
in Figure 9.1. It can be seen from Figure 9.1 that, the supplies of water from the Kailana-
Takhatsagar Reservoir to the city have increased gradually over the years from 228 lacs
gallon per day in the year 1996 to 522 lacs gallon per day by the year 2009 to meet the
rising demands for water by the growing population and their allied activities.
600
522
516
506
Water Supply (lac g/day)
476
500
420
392
378
374
373
365
366
400
337
309
259
300
228
218
204
188
183
180
200
100
0
1991
1992
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
Year
Figure 9.1: Water supplied from the Kailana-Takhatsagar Reservoir to the Jodhpur city
in different years during 1991-2009.
Table 9.1: Analyzed information on the random door-to-door survey carried out in the
problematic area for ascertaining usages of waters and wastewaters by the
householders.
Nos. of family Source of water Method of disposal of Method of disposal of
responded to supply & no. of domestic wastewater & toilet water & no. of
questionnaire family no. of family family
Nagar Nigam : 21 Sewerage system : 13 Sewerage system : 22
23 PHED : 2 Nallah, etc. :5 Nallah, etc. :1
Soak Pit :5
Total : 23 Total : 23 Total : 23
In most part of the old city area, the sewerage/drainage systems are found to be
chocked by debris, and the flow velocity is sluggish. Therefore, leakage/seeping of
wastewaters and return flow to the aquifer from the households, which are not connected
to the sewerage system, can not be ruled out, in such cases. Currently, the city’s
wastewaters are drained out through three main sewerage systems, one towards the
Nandri area, and the other one to the Jodhpur Airport side, and the third one is near to the
Polytechnic Institute, which have the disposal outlets to the Jojri River. Of these, Nandri
site has an organized sewage treatment plant having capacity of 20 MLD (Million liter
per day). In the other two cases, they are directly discharged to the Jojri River without
treatment. Whether or not the wastewaters generated from the water supplied to the city
(a) (b)
Figure 9.2 : Sewage measurement location near the Jodhpur Airport site; (a) section
chosen for measurement of cross-section and flow velocity, and (b)
temporary structural arrangement made for measurement of cross-section
and flow velocity, and measurement of depth of flow.
Figure 9.3 : Sewage measurement location near the Jodhpur Polytechnic Institute; (a)
pre-calibrated V-notch installed at the site by constructing a guided
channel for measurement of flow, and (b) depiction of flow over the V-
notch.
(b)
(a)
Figure 9.4 : Sewage measurement location at the Nandri sewage treatment plant site; (a)
calibration & fixation of the data logger connected to the current meter, and
(b) measurement taken from the current meter for estimation of flow
velocity.
From the rating curves (Figs. 9.5-9.7) it can be seen that sewage flows in all the
three drains vary during 24 hours in a day. Largely, in a day the flows are found to be
minimum in morning 6-7 A.M., thereafter, the flows gradually increase and attain to a
peak value between 10:00 – 11:30 A.M; after that flows gradually reduces. A second
peak is also found common in all the three cases between 9:00 P.M. to 10:30 P. M.
Thereafter, the flows decrease to attain the minimum in the morning. The pattern of
sewage flows depend on the uses of the supplied water by the households in the city
areas. In general, the maximum uses of water are seemed to be during the morning 8:00
A.M to 12 noon, followed by evening 6:30 P.M. to 10:00 P.M.
Out of the three sewerage drains, the Jodhpur Airport drain disposed off
maximum sewages of about 71279 m3/day (157 lacs gallon per day) followed by Nandri
Treatment Plant drain of about 13235 m3/day (29.15 lacs gallon per day); while the
Polytechnic Institute drain disposed off an average of about 3325 m5/day (7.12 lacs
gallon per day). The total discharge of sewages from the city areas through these three
sewerage systems in terms of the water supplied in the year 2009 is about 37%. The
sewage flows may vary from season to season, even from month to month that depend on
the supplies and uses of water. In addition to the three major drains, there are number of
small uncounted drains, which are connected to the nearby Talabs located in the city area.
Assuming that the wastewaters generated and their outflows from the city area are
respectively, 65% and 37% of the water supplied in the respective years. The probable
quantities of the wastewaters, which have been generated and disposed off from the city
areas in different years, are given in Table 9.2.
Table 9.2: Probable quantities of wastewaters generated and disposed off from the city
areas in different years.
Year Water supplied Probable quantity of Probable quantity of
(Lac gallon per day) wastewaters generated wastewaters disposed off
(Lac gallon per day) (Lac gallon per day)
2000 378 245.7 139.86
2002 392 254.8 145.04
2004 366 237.9 135.42
2006 476 309.4 176.12
2008 518 335.4 190.92
****************
An assessment of the unknown fraction of surface water supply that joins the
aquifer system causing water level rise in the city area has been made from surface water
balance study for the city area. In water balance, one is required to confirm that the
quantitative difference between sum of all inputs and outputs of the hydrological
components, those are involved in the governing physical processes in a system for a
specific time period, is equal to the change in water storage in the system during that time
period. In the case of the present study area, the hydrological components of inputs are:
(i) water supplied to the city area, and (ii) rainfalls and the corresponding runoffs; while
the outputs are: (i) sewage outflows from the city area, (ii) evapotraspiration, and (iii)
groundwater recharge. The surface water balance considering these components can be
written as:
[(Qws ]
+ QRi + QRo ) − ( f Qws + Qsw + Qevp + Q gr ) Δt = ΔV …………………. (10.1)
in which Qws is the rate of water supply from the Kailana-Takhatsagar Reservoir to the
city area (L3T-1); QRi is the rainfall rate over the city area (L3T-1) = Ri As, Ri is the rainfall
quantity (L), and As is the surface area on which rainfall occurs ( L2); QRo is the runoff
rate to the area corresponding to the rainfall (L3T-1); f is the fraction of the water, Qws that
is lost by the consumptive use (dimensionless), Qsw is the sewage & storm water outflows
from the city area (L3T-1); Qevp is the evapotranspiration rate of water from the study area
(L3T-1); Qgr is the rate of ground recharge from the study area (L3T-1); ΔV is the change in
storage of water in the city area (L3); and Δt is the time period (T).
The daily rate of water supply from the Kailana-Takhatsagar Reservoir to the city
area, Qws has been given in section 7. The daily and annual variation of rainfall, Ri have
been discussed in section 4. The water surface evaporation rate has also been explained in
section 4. The sewage & storm water outflows rate, Qsw and the fraction of consumptive
use (assumed to be, f = 0.35) have been elaborated in section 9. The study area, As is
measured to be 76 sq. km. Making use of the respective data, and assuming the time
period to be annual, i.e., Δt = 1 year, the water balance of the study area is carried out.
The information/data used for calculation are given in Table 10.1:
Study on Rising Groundwater Table in Jodhpur City - NIH 83
Table 10.1 : Data used for calculation of components of water balance equation of the
Jodhpur study area.
The water balance components can vary during monsoon and non-monsoon days,
even during different hours in a day, which depend on water supplied to the city area,
water usages pattern, rainfall and its intensity and duration, land-uses and land cover,
storage retention, temperature, etc. However, making use of the data given in Table 10.1,
a lumped water balance of the study area has been made for the years 2000, 2002, 2004,
2006 and 2008 and is given in Table 10.2. The lumped storage volumes shown in Table
10.2, which are the differences of Inflow components and the outflow components, are
the quantities of water those might have joined to the underneath aquifer through return
flow over the respective year in addition to the normal rainfall recharge.
Table 10.2: Water balance components of the study area in different years.
Year Inflow components Outflow components Storage (ΔV)
Water supply, Runoff Consumptive Sewage + Groundwater (MCM) =
Vws (MCM) volume use compo- runoff recharge
(VRi +VRo) nent; f Vws outflow from component, ∑ Inflow -
(MCM) (MCM) the city area, Vgr (MCM) ∑ Outflow
Vsw (MCM)
2000 56.62 6.087 19.82 27.036 8.492 7.356
2002 57.47 2.075 20.11 23.339 8.621 7.476
2004 56.92 6.430 19.92 27.492 8.539 7.404
2006 76.62 3.443 26.82 31.791 11.492 9.956
2008 85.57 8.618 29.95 40.281 12.836 11.124
(MCM : Million Cubic Meter)
***************
In a hard rock region in Jodhpur city, aquifer tests have been conducted in small
radius tube wells and in large diameter wells locally known as Bowaris. In this study, for
the aquifer test conducted in the wells with small radii, the inverse problem has been
solved using the Theis’ basic solution treating the aquifer to be confined. For the tests
conducted in large diameter wells, Hantush’ basic solution for well with finite radius has
been used considering well storage effect on drawdown data. Unit pulse kernel
coefficients are generated and used in a Marquardt Algorithm as described in
Appendices. The convolution technique used to compute drawdown is quite versatile.
The parameters have also been estimated using a simpler Paul algorithm search
technique.
Table 11.2 :Transmissivity and Storage coefficient as obtained through successive iteration.
Iteration T
*
φ* ΔT Δφ Error: C(1) Error: C(2)
no
1 0.017997 0.000373 8.00E-03 7.33E-05 0.00E+00 0.00E+00
2 0.028653 0.000338 1.07E-02 -3.55E-05 0.00E+00 0.00E+00
3 0.036955 0.000276 8.30E-03 -6.21E-05 0.00E+00 0.00E+00
4 0.039092 0.00027 2.14E-03 -5.46E-06 0.00E+00 0.00E+00
5 0.039129 0.000272 3.65E-05 1.32E-06 0.00E+00 0.00E+00
6 0.039133 0.000271 4.43E-06 -9.17E-08 2.78E-17 8.88E-16
7 0.039132 0.000271 -5.05E-07 1.11E-08 0.00E+00 0.00E+00
8 0.039133 0.000271 6.01E-08 -1.31E-09 -4.34E-19 0.00E+00
9 0.039133 0.000271 -7.14E-09 1.56E-10 5.42E-20 0.00E+00
2.5
2
Drawdown (m)
1.5
1 Observed drawdown
Simulated drawdown
0.5
0
0 100 200 300 400
Time (Min)
Table 11. 5: Transmissivity and storage coefficient as obtained through successive iteration
Iteration
no T* φ* ΔT Δφ Error: C(1) Error: C(2)
1 0.024 0.01 -0.0518 0.117059 5.70E-06 3.44E-06
2 0.024 0.127059 -0.03277 0.418365 -7.69E-06 -6.00E-07
3 0.024 0.127059 -0.03485 0.423725 -7.69E-06 -6.00E-07
4 0.024 0.127059 -0.03485 0.423725 -7.69E-06 -6.00E-07
1.5
Drawdown (m)
Observed drawdown
1
Simulated drawdown
0
0 20 40 60 80 100 120 140 160 180 200
-0.5
Time (min)
2.5
1.5
Drawdown (m)
Observed drawdown
1 Simulated drawdown
0
0 20 40 60 80 100 120 140 160 180 200
-0.5
Time (min)
Figure 11.3: Observed and simulated drawdowns for T=0.0396 m 2 / min = 57 m 2 / day ,
Storage coefficient, φ = 0.0298; estimated by Paul Algorithm (Subhash
Chowk).
0.8
Observed drawdown
Drawdown (m)
0.6
Simulated drawdown
0.4
Difference between observed and
simulated
0.2
0
0 20 40 60 80 100 120 140 160 180 200
-0.2
Time (min)
Figure 11.4: Observed and simulated drawdowns for T = 0 .1156 m 2 / min =166.46
m 2 / day and Storage Coefficient, φ = 0.2337 Estimated Using Marquardt
Algorithm ( Suraj kund ).
1.2
0.8
Observed drawdown
Drawdown (m)
0.6
Simulated drawdown
0.4
Difference between observed and
simulated
0.2
0
0 20 40 60 80 100 120 140 160 180 200
-0.2
Time (min)
Figure 11.5: Observed and simulated drawdowns for T= 0.1 m 2 / min =144.00 m 2 / day
and Storage coefficient, φ = 0.121 estimated using Paul Algorithm (Suraj
kund ).
0.6
0.5
0.4
Drawdown (m)
Observed drawdown
0.3
Simulated drawdown
0.2
Difference between observed
and simulated
0.1
0
0 50 100 150 200
-0.1
Time (min)
Figure 11.6: Observed and simulated drawdowns for T = 0 .17 m 2 / min =244.8 m 2 / day ,
Storage coefficient, φ = 0.181, estimated by Marquardt Algorithm
(Sanicherji ka than).
0.5
Difference between observed and
simulated
0.4
D rawdown (m )
0.3
0.2
0.1
0
0 50 100 150 200
-0.1
-0.2
Time (min)
Figure 11.7 : Observed and simulated drawdowns for T= 0.2 m 2 / min = 288.00 m 2 / day
and Storage coefficient, φ = 0.25; estimated using Paul Algorithm
( Sanicherji ka than).
*******************
In order to find the possible causes of waterlogging in the city area and their
sources, the groundwater quality data are also considered for analysis. The PHED, Govt.
of Rajasthan has monitored groundwater quality at 53 locations in the city area during the
year 2009, and analyzed the monitored samples for determination of concentration of 13
chemical constituents, namely; pH, Cl2, TDS, NO3, EC (Electrical Conductivity),
Sodium, Potassium, Calcium, Magnesium, Carbonate and Bicarbonate ions, Thorium,
Fluorine. Out of these 13 chemical constituents monitored at 53 locations, data of pH,
TDS, Cl2, NO3, Sodium, Potassium, Calcium, Magnesium, Carbonate and Bicarbonate
The reason of selecting pH, TDS, Cl2, NO3, Sodium (Na+), Potassium (K+), Calcium
(Ca++), Magnesium (Mg++), Carbonate (CO3-) and Bicarbonate (HCO3-) ions is; they
represent dominant characteristics of water quality fate, viz., pH is one of the basic
properties of water that indicates acid and alkaline characteristics of water, TDS gives
concentration of total dissolved solids those include cations and anions, and is an indicator
of taste of water; Cl2 indicates salty taste, that is mostly associated with fresh/brackish
water; NO3 concentration is associated with the occurrence of sanitary/landfill/agricultural
wastes; Na+, K+, Ca++, Mg++, CO3- and HCO3- ions are indicators to detect fitness of water
for irrigation uses. Making use of SURFER software and employing kriging method, the
thematic maps of the spatial variation of pH, Cl2, TDS and NO3 for pre and post monsoon
within the study area have been prepared and are shown in Figures 12.1-12.7. The Sodium
Adsorption Ratio (SAR) that gives the measure of the suitability of water for use in
agricultural irrigation, has been estimated using the following formula (Richards, 1954):
Na +
SAR = …………………………….. (12.1)
(Ca 2 + + Mg 2 + ) / 2
in which all ionic concentrations are expressed in mill equivalent per litre.
A low SAR (2 to 10) indicates no danger from sodium; SAR between 7 and 18
indicates low to medium hazards, between 11 and 26 high hazards, and a SAR > 26
indicates very high hazards. The lower the Ca 2+ and Mg 2+ ionic strength of the solution,
the greater is the sodium hazards for a given SAR (Richards, 1954). The water used for the
irrigation purpose can have a SAR value maximum up to 26. The recommended
classification with respect to electrical conductivity, sodium content, Sodium Absorption
Ratio (SAR) and Residual Sodium Carbonate (RSC) is given in Table 12.1.
Figure 12.1: Spatial variation of pH of groundwater in the study area (• indicates location
of sampling points).
It can be seen from Figure 12.1 that the pH values within the study area vary from
7.4 to 8.20, with slightly higher values in the waterlogged area. The pH values are within
the permissible limit showing no sign of interference by other chemical constituents. There
are three concentration peaks, two in the waterlogged area, and another in south-west. The
contours around each peak show evidence of dispersion. Mechanical dispersions are
caused by groundwater flow. The two peaks in water logged area indicate that the high
There are two troughs one towards North and the other in South East. The contours
around each trough show evidence of dispersion. From the gradient direction of the
contours and decreasing trend to wards north corresponding to the trough in the North, it
could be seen that no ground water flow is occurring from north to the waterlogged area.
Figure 12.2: Spatial variation of Chloride (Cl2) concentration during pre-monsoon season
in the study area (• indicates location of sampling points).
Figures 12.2 and 12.3, which show distribution of Cl2 concentration in the study
area during pre and post monsoon period, respectively, indicate the variation between 50
and 3550 mg/l, with higher values at a concentrated pocket down below the waterlogged
area. Within the waterlogged area the concentrations of Cl2 range between 100 mg/l and
1450 mg/l. The permissible limits range between 100 mg/l and 600 mg/l. The peak of
highest chlorine concentration is out side water logged area. The reason of higher Cl2
concentration at this particular location needs further verification. The Cl2 concentration is
low in the north-western part of the waterlogged area and increases further as one move
towards the south-east direction. In general, the concentrations Cl2 in the waterlogged area
and within the study area do not indicate the source of water logging to be from the
leakage of sewerage lines.
Figure 12.5: Spatial variation of TDS concentration during post-monsoon season in the
study area (• indicates location of sampling points).
Figure 12.6: Spatial variation of NO3 concentration during pre-monsoon season in the
study area (• indicates location of sampling points).
Figures 12.6 and 12.7 show the spatial variation of NO3 concentration in the study
area for the pre and post monsoon period, respectively. The variation of NO3 concentration
is observed to be between 10mg/l and 850mg/l; while the permissible limit ranges between
20 mg/l and 50 mg/l. The higher values are mainly observed in the southern part of the
study area. In the waterlogged area, the values range from 50 mg/l to 350 mg/l. It indicates
that the waterlogged area has some effects of NO3 contamination. The possibilities of
sewage leakage at different pockets can not be over ruled. The southernmost part of the
study area has larger effect of NO3 contamination. As stated earlier, the source of NO3 is
mainly from the septic tanks/sewage outflows or agricultural refusal. In the southern part
of the study area where NO3 concentration is more are mostly agricultural areas which are
generally irrigated by city’s sewage waters, mostly used for producing winter crops.
The spatial variation of SAR in the study area calculated using eq.(12.1) based on
the concentration of Na+, Ca++, and Mg++ is shown in Figure 12.8. It can be seen from the
Fig. 12.8 that the SAR values vary between 0 and 14 within the study area. The value is
higher in the eastern part of the waterlogged area and gradually decreases as one move
towards the western direction. The water which can safely be used for the irrigation
purpose can have the SAR value below 26. That means, the groundwater in the study area
can safely be used for irrigation purposes without environmental hazards.
(i) The pH values of the groundwater in the Jodhpur city area are within the
prescribed permissible limit having slight leaning towards alkalinity, which
may be because of excessive Chloride concentration. The Chlorine is used as a
disinfectant chemical during treatment of water.
(iv) The higher TDS values, which are mostly observed in the southern part far way
from the problematic area, may be due to the local intrinsic mineral
composition in the aquifer.
(v) The Nitrate concentrations in the study area have been found much higher than
the permissible limit, mostly in the southern part, which are influenced by the
agricultural activities using the city’s sewages water. In the problematic area, in
some pockets concentration of Nitrate is also observed to be beyond
permissible limit. Therefore, leakage of sewerage lines or connectivity of the
septic tanks outlets to the groundwater system in those pockets can not be ruled
out.
(vi) The source of water logging and rise in groundwater level in some parts of the
city area appears to be due to return flow of water other than the sewages.
(vii) The quality of groundwater in the city area indicates that the groundwater can
safely be used for irrigation purposes.
********************
The study area of about 76 sq. km. bounded between latitude of 26015’ N to 260 20’
N and longitudes 730 0’ E to 730 4’ E has been discretized into 107 x 113 grids (along x
direction = 107 nos; along y direction = 113 nos) of each grid size 100 m x 100 m.
Vertically along z-direction, the maximum aquifer thickness of about 78.6 m which
varies from location to location has been divided into 6 layers in accordance with the
variation of geological formations and their hydraulic properties. Each layer has different
thickness. A schematic view of the discretized study area and its sectional view
indicating the vertical discretization prepared using the Visual MODFLOW software
(version 9.1) is shown in Figs. 13.1 and 13.2.
Figure 13.1: Discretized view of the study area including the water bodies and the affected
area within the Jodhpur city (Number of grid along X direction = 107; Number of
grid along Y direction = 113; size of each grid = 100 m x 100m).
Within the discretized rectangular area (Fig. 13.1) that encompasses an area of
120.91 sq. km, the study area of 76 sq. km constitutes 7600 grids out of 12091 grids. The
remaining 4491 grids constitute cells outside the ambit of the study area. For calculation
purposes, the grids which fall outside the study area are considered to be inactive cells,
and those come within the study area are considered to be active cells. The active cells
mean they respond to perturbation, while the inactive cells do not respond to
perturbation. The conceptualized domain of the study area is shown in Fig. 13.3.
General
head
boundar
y
Active cell
General
head
Active cell boundar
y
Figure 13.3 : Conceptualized setting of the modeling area showing position of the active
and inactive cells.
The boundaries of the study area are shown in Fig. 13.3. Most of the northern part
of the boundary constitutes hilly terrain and no flows are expected through this boundary,
and hence the cells encompassing the hilly terrain are considered to be no flow boundary.
The southern part of the study area has the general groundwater flow direction, and hence
considered as the General Head Boundary (GHB). The western part, which is directed
towards Kailala-Takhatsagar Reservoir, is partly considered as GHB and partly no flow
boundary depending upon the geological formations. The eastern part of the study is
considered as the GHB.
Vertically along Z-direction, the thickness of the alluvium formation varies from
location to location from few centimeters to a maximum of 78.6 m below ground surface.
The elevation difference of the surface topography is about 86.7 m between the
maximum (269.3m) and the minimum (182.6 m) elevation in the study area. To account
Figure 13.4 : Locations of water bodies in and around the Jodhpur city area.
The status of groundwater condition in the year 2004 during the pre-monsoon
period is considered as the initial condition.
Two types of boundary conditions are considered; one is General Head Boundary
(GHB) and the other one is No Flow Boundary (NFB). The GHB is the one through
which cells just adjacent to the boundary can exchange flow in either side of the
boundary depending upon the gradient of heads in the cells. The NFB is the one through
which no flows are allowed get in or out of the boundary.
The aquifer is considered to be unconfined, i.e., under water table condition. The
hydraulic properties of the geological formations are initially considered as given in
Table 13.1:
Table 13.1: Hydraulic properties of the geological formations, considered as the initial
guess values.
The spatial variation of the Transmissivity and Storativity values estimated for the
alluvium formations from the pumping tests data (given in Section 11), which are
considered as the initial guess values for each grid are shown in Figures 13.5 and 13.6,
respectively.
The inflow and outflow stresses are those based on which the responses of the
aquifer are to be ascertained. The inflows to the study are: return flow from the
accumulated wastewaters, rainfall recharge to groundwater, seepage from water bodies,
and the inflows through boundaries; while the outflows from the study area are:
evapotranspiration, withdrawal from the groundwater, and outflows through the
boundaries. The magnitude of spatial variation of inflow and outflow components for
different years, and in different period in a year have been calculated and assigned to
each grid accordingly. A scheme of the input stresses which have been assigned zone
wise is shown in Figure 13.7.
In active cell
Figure 13.7: A scheme of the input stresses assigned zone wise to the modelling area.
The time step size for simulation of the groundwater levels is considered 1 day.
This will constitute 365 time intervals in a year. The time step size of 1 day indicates that
all input and output stresses are to be assigned on per day basis with time step size, Δt = 1
day. The simulation is carried out setting a transient-state model considering uniform
stresses between pre-monsoon to post monsoon and post monsoon to pre-monsoon, and
so on. That means within a period between pre to post monsoon, the stresses are
considered constant.
With the above inputs data, values of the parameters and variables; the
mathematical model has been set to ascertain the responses of the aquifer using the visual
MODFLOW software (version 9.0 pre.) developed by the USGS .
**********************
The inputs data required for developing the simulation model have been explained
in the Section 13.0. To develop the simulation model, one has to first calibrate the model
parameters from the comparison of the computed and the corresponding observed
groundwater table contour profiles. The calibrated parameters have thereafter to be
validated with another set of computed and observed groundwater table contour profiles.
If the responses of the parameters are found to be comparable with the profiles which are
considered for the validation, the model parameters and the setting of the model are said
to be calibrated. The calibration a model is thus a trial and error approach.
Making use of the inputs stresses and boundary conditions as explained in the
Section 13.0, the model has been calibrated for steady state scenario employing the pre-
monsoon groundwater profile of the year 2004. For developing the model, the software
Visual Modflow 9.1 has been utilized. The model parameters were calibrated by
comparing the computed and the observed water table. In absence of any physical
boundary like a river or large surface water body, use has been made of the general head
boundary to simulate the flux into/out of the system. Figure 14.1 shows the comparison
between the pre-monsoon computed water table and the observed water table for the year
2004.
Taking the steady state scenario as the initial condition for the transient
simulation, the water table elevations are computed for post-monsoon 2004. The total
simulation period is taken as 180 days. During this period, the return flow from the
wastewater and rainfall recharge are also taken into account. The simulated profile is
shown in Fig. 14.2 along with the observed water table for the post monsoon of the year
2004.
Figure 14.2 : Comparison of the steady-state simulated and observed groundwater table
for the post-monsoon period of the year 2004.
Figure 14.3 : Calibrated hydraulic conductivity (m/day) zones of the Aquifer below
Jodhpur city area
For evolving the management plan to contain the rise in groundwater table, the
following options are attempted:
The topography of the city area in the problematic and adjoining part is of
undulating type. The thickness of the alluvium formation varies from about 2.0 m to 13.0
m in this region. In general, considering the possibility of capillary rise in alluvium
formation, groundwater table should be at least 0.5m below the foundation level. In areas
where building basement floors are located, assuming that basement depth is 3.5m, and
capillary rise is 0.5m, the groundwater table should be at 4.0 m below the ground surface.
The groundwater level in the affected area thus has to be lowered down to a depth of 4.0
m from the ground surface elevations. Nevertheless, the problematic area is densely
populated with concentrated residential and commercial buildings, and practically has
very thin scope to construct horizontal sub-surface drainage system at a depth of 4.0m
below the ground surface. Moreover, laying of horizontal subsurface drains will be very
difficult as the required excavation will interfere with the already existing building
foundation structures. All these aspects do not promote provision of a usual horizontal
sub surface drainage system in the problematic area. Therefore, this option is not
appeared to be a logical proposition.
The problematic area being located in the densely populated and congested area,
provision of network of vertical pumping wells looks to be a logical proposition and can
be constructed with least interference to the area. Already, a number of vertical pumping
wells are in operation in the problematic area.
It is intended here to assess the number of pumping wells required to lower the
water table to the required depth and to maintain its position at the desired level. It is also
intended to examine how to dispose the pumped water more effectively.
The aim is to lower the groundwater table up to 4.0 m below the ground surface.
The volume of water existing up to a depth of 4.0 m below ground surface in the
saturated zone is assessed using the following equation:
Vs = Ai ∗ H i * S y …………………………………………….. (15.1)
The location of the ponds in the problematic area, groundwater table contour map,
and the alignment of the major sewerage systems in the Jodhpur city area are shown in
Figure 15.1.
Figure 15.1 : Location of ponds in the problematic area, and the alignment of sewerage
system in the Jodhpur city area.
The similar approach of discharging the pumped water to the ponds located in the
problematic area through conveyance pipe can suitably be implemented. From these
ponds, surface channel can be constructed to convey the water stored in the pond to the
surface drainage system through gravity.
*******************
16.1 Summary
In order to identify the source of water, that has caused water logging, and to
suggest appropriate remedial measures to control the rising trend of ground water levels
in Jodhpur city, detailed analyses of the data related to topography, demography,
geological formations, hydrometeorology-hydrology and hydrogeology, groundwater
quantity and quality, sewage flows, inflows and outflows of waters to/from the Jodhpur
city including the Kailana-Takhatsagar Reservoir have been carried out. These data have
been supplied mostly by the Ground Water Department (GWD) and the Public Health
Engineering Department (PHED) of Jodhpur, Government of Rajasthan; some data,
which were not made available but were required for the study, have been collected from
the field investigations/surveys and from other sources during course of the study.
In order to analyze the data, all spatially varying databases have been geo-referenced
with reference to the geographic coordinate system (WGS-1984) by their latitudes and
longitudes. Therefore, the analyzed data and results presented in this report can be
verified with the field ground truth.
The study area is comprised of 76 sq. km., and encompasses the old and the sprawled
Jodhpur city area including the waterlogged area. The Kailana-Takhatsagar Reservoir is
located outside the boundary of the study area as the study area and the reservoir are
located in different geological entities. The topography of the Jodhpur city area has been
analyzed making use of surveyed data supplied by the GWD, Jodhpur in conjunction with
the ASTER data. The geological formations of the Jodhpur city area have been analyzed
making use of the borelogs data at 93 locations supplied by the GWD, Jodhpur using
‘ROCKWORKS’ software.
The demographic data of the city area have been analyzed based on the
population census of the decades from the year 1971 to 2001. The water supply
The existing arrangement of the water supply to the city area, which is based on
feeding the Kailana-Takhatsagar Reservoir from the IGNP linked Rajiv Gandhi Lift canal
and transferring the water from the Kailana-Takhatsagar Reservoir largely by pumping
and partially by gravity flow for treatment of water and then supplying to the city, has
been analyzed in details.
The Stage-Area-Capacity curves for the Kailana Lake, Takhatsagar and Umaid Sagar
Reservoirs have been developed. The water balance of the Kailana Lake and Takhatsagar
Reservoir involving the associated hydrological components has been carried out
exclusively for two situations: (i) in the first situation, it has been assumed that the lake
has a permeable bed through which water is seeping to the underneath formation in which
water table lies at a large depth, (ii) in the other situation, the lake bed is impermeable,
therefore, the seepage is negligible.
The groundwater levels data supplied by the GWD, Jodhpur for the years 1996-2009
have been analyzed considering two situations; (i) the Kailana-Takhatsagar Reservoir is
assumed to be hydraulically connected to the aquifer below the Jodhpur city, and (ii) the
Kailana-Takhatsagar is not hydraulically connected to the aquifer below the Jodhpur city.
The increase in waterlogged areas in different years has been ascertained.
The sewages generated from the refusal of water usages, and the drainage data of the
city area have been analyzed in details. In order to assess the quantity of sewages flown
out daily from the city area through the existing sewerage channels, field investigations
and measurements have been carried out in the three major sewerages drains.
Measurements of sewages’ flows from morning 6:00 A.M to evening 11:00 P.M. for a
continuous 7 days in all the three drains have been made during the months of April-July,
2010.
The groundwater quality scenarios in the study area for the year 2009 have been
analyzed based on the data supplied by the PHED, Jodhpur. The spatial variations of the
parameters pH, Cl2, TDS, NO3 and SAR have been analyzed and their possible reasons of
occurrence and fate in the study area have been analyzed in details.
16.2 Conclusions
Based on the analyses of data and detailed results reported in different chapters,
the following conclusions are drawn:
(i) The waterlogged area and it’s extent are mainly located in the Quaternary
alluvium formation below the Jodhpur city area.
(iv) The Kailana-Takhatsagar and the Jodhpur city areas are located on two distinct
geological units. The Kailana-Takhatsagar area is located on Rhyolite formation,
whereas, Jodhpur city area is mainly located on the Quaternary alluvium,
Sandstone and Shale formation. These formations have different hydro-geological
properties and hence, cannot be considered as a single system.
(v) The geological formations of the Jodhpur city area indicated presence of an
unconfined aquifer comprised of Quaternary alluvium formation having thickness
varying from few centimeters around NW and NE direction, to a large depth
ranging from 6 to 75 m in the SW to SE direction. The unconfined aquifer is
underlain by Shale, Sandstone and Rhyolite formation of different thickness
varying from location to location. As such, there is no existence of any deeper
aquifer below the Jodhpur city area up to a depth of 75 m, below the ground
surface.
(vi) Thirty six years (1971-2006) daily rainfall data of the Jodhpur city showed the
average annual rainfall of 378 mm.
(vii) The Jodhpur city that had a population of about 3,65,000 in the year 1971, was
expected to have a population of the order of 11,00,000 by the year 2010. The
population censuses of previous four decades (1971-2001) indicated a growth rate
of 3.21% per year.
(viii) The water supply requirements to the population in the city area are primarily met
by the Jodhpur Municipality through water supply & distribution systems
managed by the Public Health Engineering Department (PHED), Jodhpur. As
such, no prescribed guidelines towards quantity of water per capita per day to be
supplied are followed; it is based on thumb rule in accordance to the supply-
demand norm.
(x) For both pre and post monsoon seasons, the directions of groundwater flow are
broadly observed from the north to the south and southeasterly direction with
some deviations at local pockets.
(xi) In the waterlogged area, the water table contours are closely spaced. In the
southern region, the spacing between two consecutive contours is comparatively
wider. The region in which the contours are closely spaced is a region of low
transmissivity. The region in which the contours are widely spaced is a region of
comparatively higher transmissivity indicating larger thickness of the Quaternary
alluvium formation. The waterlogged area is located in the geological unit having
low transmissivity where the contours are closely spaced. In the waterlogged
area, the water table contours have higher values than those in the southern side.
Since water has flown from the waterlogged region and the North side is a hilly
region and no water has come to the waterlogged region from the north side, the
source of water causing water logging is thus appeared to be generated locally
which is likely to be the return flow from the water supply including sewages.
(xii) The possibility of seepage from the Kailana-Takhatsagar Reservoir entering the
waterlogged area investigated making use of the contour maps assuming that the
Kailana-Takhatsagar Reservoir is not hydraulically connected with the aquifer
underneath the Jodhpur city, and the Kailana-Takhatsagar Reservoir is
hydraulically connected with the aquifer underneath the Jodhpur city. For the
former case, it is noted that the waterlogged area does not receive any flow from
the southwest side i.e., from the reservoir side. So it is inferred that groundwater
flow from the Kailana-Takhatsagar side is not causing water logging in the
waterlogged area. For the later case, it is noted that the direction of flow lines is
not towards the waterlogged area. The direction of flow in all the cases is broadly
towards south-east direction. The seepage from the Kailana-Takhatsagar
Reservoir is not entering the waterlogged area. Thus, whether the Kailana-
Takhatsagar Reservoir is hydraulically connected with the aquifer underneath the
Jodhpur city or not, the seepage from the lake is not flowing towards the
waterlogged area in either situation.
(xiv) Evolution of the water table has caused water logging progressively. Area getting
waterlogged has increased every year, and a dynamic equilibrium has not been
reached.
(xv) As such, there is no sign of bulk groundwater seepage from the Kailana-
Takhatsagar Reservoir to the Jodhpur city area through the underneath geological
formations of the respective areas. If at all seepage had taken place, it might have
occurred through lineaments. There is only one lineament, which is oriented
towards the waterlogged area. The Lineament analysis survey and the geological
and geophysical study conducted independently by the National Geophysical
Research Institute (NGRI), Hederabad (2010) surrounding the Kailana-
Takhatsagar Reservoir indicated that the lineaments are oriented in NNE-SSW to
NE-SW directions with no connectivity to the city areas. A few lineaments with
ESE-WNW directions are present but these are small and do not have
connectivity to the city areas. The findings of the NGRI thus corroborate the
present finding based on hydraulic principle. The chances of seepage from the
Kailana-Takhatsagar Reservoir to the waterlogged areas through lineaments are,
therefore, ruled out.
(xvi) The measurements and analyses of the sewages data of the year 2009 indicated
sewages disposal of about 37% of the supplied water through three major
sewerage systems, i.e., Airport drain, Polytechnic Institute drain and Nandri
sewage treatment plant drain.
(xvii) The pH values of the groundwater in the Jodhpur city area are within the
prescribed permissible limit having slight leaning towards alkalinity, which may
be because of excessive Chloride concentration. The Chlorine is used as a
disinfectant chemical during treatment of water. The contours of pH concentration
indicate that the ground water flow direction is not towards the water logged area.
(xix) The higher TDS values, which are mostly observed in the southern part far way
from the problematic area, may be due to the local intrinsic mineral composition
in the aquifer.
(xx) The nitrate concentrations in the study area have been found much higher than the
permissible limit, mostly in the southern part, which are influenced by the
agricultural activities using the city’s sewages water. In the problematic area, in
some pockets concentration of nitrate is also observed to be beyond permissible
limit. Therefore, leakage from sewerage lines or connectivity of the septic tanks’
outlets to the groundwater system in those pockets can not be ruled out.
(xxi) The source of water logging and rise in groundwater level in some parts of the
city area appears to be due to return flow of water from water supply system and
from the source other than the sewage waters. In some pockets, the seepage from
sewage system cannot be ruled out.
(xxii) The quality of groundwater in the city area indicates that the groundwater can
safely be used for irrigation purposes.
(xxiii) Vertical drainage system is feasible and easy to construct in the urbanized area.
Constructing subsurface horizontal drain is not feasible because of the existing
basements.
(xxiv) Pumped water can also be used for irrigating gardens and parks. Therefore, the
pumped water should not be discharged into the existing surface drains which are
meant for disposal of sewages.
****************
1. As the first and foremost remedial measure, it is suggested to regulate the quantity
of water being supplied to the city area at the source itself, i.e., regulation of water
from the Kailana-Takhatsagar Reservoir. The regulation needs to be based on per
capita per day water requirement basis. The Jodhpur city being located in the arid
and water scarce region, about 110 liters per capita per day could be taken as the
guideline. The break up of 110 liters is as follows: 70 liter (drinking, bathing &
toilet flushing) + 20 liter (commercial uses) + 20% conveyance losses. Industrial
water requirements are to be included separately. For 110 lpcd supply, the
quantity of water requirement for the estimated population of 11,08,950 in the
Jodhpur city for the year 2010 is worked out to be 268.69 lac gallon per day. To
meet water requirement for domestic animals, and kitchen gardens and public
parks 1/3 of the requirement for the population i.e., about 73 lac gallon per day
extra water be supplied, which will reduce the quantity by about 35% over the
quantity supplied (521.7 lac gpd) in the year 2009. The per capita per day water
supply requirement is estimated to be 140 litre.
2. In the affected area, the water supply lines need to be thoroughly checked to find
the locations of leakages, and suitable remedial measures to stop the leakages
need to be taken up. The sewages/drainage lines in the affected area need to be
properly sealed to stop seepage, if any.
3. The topography of the city area in the northern and middle part is of undulating
type. The thickness of the alluvium formation varies from about 2.0 m to 13.0 m
in that part. Therefore, a single generalized safe ground water level is incorrect to
suggest. In general, considering the possibility of capillary rise in alluvium
formation, groundwater table should be at least 0.5m below the foundation level.
In areas where building basement floors are located, assuming that basement
height is 3.5m, and capillary rise is 0.5m, the groundwater table should be at 4.0
m below the ground surface. The groundwater level in the affected area thus has
to be lowered down below 4.0 m from the respective ground surface elevation.
The highest groundwater table elevation in the affected area is about 285 m.
5. Provision of vertical drainage system i.e., by pumping the water from the aquifer
in the problematic area looks feasible; as drainage wells can be constructed with
least interference in the urbanized area. The pumping rate and schedule can be
controlled, the number drainage wells can be increased in a locality as required,
and already such practice has been initiated in the area, all these factors favor
provision of vertical drainage. In region of low transmissivity area i.e.,
transmissivity < 30m 2 / day , large diameter wells of 0.5m can be constructed.
However, vertical drainage system would require electrical energy, and would
cause noise pollution.
6. There are three large ponds, namely, Baiji Ka Talab, Fateh Sagar, and Gulab
Sagar (details given in table below), located near to the problematic area. The
pond beds are more or less impervious, or if necessary, these can be made
impervious by lining. The pumped water can be discharged to these ponds
through conveyance pipe. From these ponds surface channel can be constructed to
convey the water stored in the pond to the existing surface drainage system
through gravity. These waters can be used for agricultural irrigation purposes. It
should not be mixed with sewage waters.
7. A Bentonite clay grout curtain across the lineament which is terminating before
but directing towards the water logged area can be constructed proximity to the
reservoir site to check the groundwater flow, if any, from the reservoirs to the
water logged area.
Study on Rising Groundwater Table in Jodhpur City – NIH 127
REFERENCES
Table A 6.1 : Area and Capacity data of the Kailana for different stages of water level in
the Lake ( Reference : Section 6.0).
Table A6.2: Area and Capacity data of the Takhatsagar Reservoir for different stages of
water level in the reservoir (Reference: Section -6.0).
Table A6.3 : Stage - Capacity data of the Umaidsagar Reservoir (Reference : Section-6.0).
R.L.of ground
S.no Well Location Type of Well Latitude Longitude
surface
1 Raj Basera Open Well 26° 15' 59.04" 73° 3' 43.80" 206.875
2 Polo ground Tube Well 26° 15' 47.99" 73° 2' 10.98" 219.8
3 New Campus Univ. Piezometer 26° 14' 49" 73° 00' 59" 224.6
4 Thorion Ki Dhani 26° 14' 18" 72° 58' 16" 228.375
5 Subhash Chowk Open Well 26° 16' 41.36" 73° 1' 54.42" 229.375
6 R. T. O. Office Hand Pump 26° 18' 23.18" 73° 3' 41.52" 230.58
7 Jetta Bera,Opp Chand Baori Open Well 26° 17' 50.99" 73° 0' 56.0" 265.5
8 Ground Water Dept. Piezometer 26° 15' 36" 73° 0' 42.0 231.8
9 Gauesala/Umed Club Open Well 26° 17' 15.79" 73° 2' 2.04" 233.8
10 Ajeet Bhawan Open Well 26° 17' 7.19" 73° 2' 21.40" 233.9
11 Laxmi Nagar Open Well 26° 18' 3.32" 73° 2' 58.56" 234.21
26° 16'
12 Nehru Park Open Well 73° 0' 50.31" 234.5
32.1456"
13 Rai Ka Bag(Near RAE) 26° 17' 21.06" 73° 2' 43.02" 234.65
14 High Court Hand Pump 26° 17' 25.90" 73° 2' 8.09" 235.28
Tar Ghar Mahavir Complex
15 26° 16' 40.80" 73° 0' 57.70" 235.28
Sardar pura
16 Jalam Singh ka Hatta Hand Pump 26° 17' 51" 73° 2' 42" 236.53
Paota"B"Road(Near Ghanshyam
17 Open Well 26° 17' 46" 73° 2' 30" 236.8
Singh)
18 Amritiya Bera(Mahamandir) Open Well 26° 18' 7" 73° 2' 44" 237
19 Khema Ka Kua Open Well 26° 15' 16.82" 72° 58' 56.69" 237.45
20 CAZRI Piezometer 26° 15' 17" 72° 59' 42" 239.52
21 Sardarpura 3 road Open Well 26° 16' 46.28" 73° 0' 29.12" 241.37
22 Mahamandir Ka Jhalra Step Well 26° 18' 22" 73° 2' 37" 241.85
23 Nai Sarak(mehro Ka Chowk) 26° 17' 27" 73° 1' 32" 243.15
24 Gorinda Baori Step Well 26° 17' 31" 73° 1' 20.39" 243.38
25 Gantaghar/Do Kotho Ke Beech Open Well 26° 17' 42" 73° 1' 31" 244.39
26 Ramanuj Kot Open Well 26° 17' 47" 73° 1' 44" 244.85
27 Deo Nagar Open Well 26° 16' 6" 72° 59' 20" 245.38
28 Siwanchi Gate Open Well 26° 17' 4" 73° 0' 44.3" 246.85
29 Moti Chowk 26° 17' 33.9" 73° 1' 11.20" 247.35
30 Katla Bazar Hand Pump 26° 17' 37" 73° 1' 13" 247.38
31 Sindhi Bhtoo Ki Maszid 26° 18' 4" 73° 1' 53" 248.39
32 Tunwar Ji Ka Jhalra Step Well 26° 17' 49" 73° 1' 26" 248.68
33 Bhram Bag(Galori Gate) Open Wel 26° 16' 59" 73° 2' 44" 248.85
34 Khanda Falsa Open Wel 26° 17' 22.24" 73° 0' 50.2194" 249.78
35 Dau Ki Dhani Open Well 26° 17' 22" 72° 59' 14" 252.41
36 Akhe Raj ji ka Talab Open Well 26° 17' 49.30" 72° 59'5" 254.63
37 Kaylana I Piezometer 26° 17' 41.89" 72° 59' 32.60" 258.85
38 Nav Chokiya 26° 17' 43" 73° 0' 52" 259.68
39 Ganwa(Pratap Nagar) Open Well 26° 17' 46" 72° 59' 38" 260.12
40 Kailana II Piezometer 26° 17' 41.87" 72° 59' 38.80" 262.85
Using Theis’ basic solution i.e. by matching the time drawdown curve with Theis’
type curves, we can determine the parameters of an aquifer which is confined,
homogeneous, isotropic, and is of infinite area and the pumping well has small radius.
The aquifer is to be initially at rest condition, and the aquifer test is conducted under
constant pumping rate. These are the assumptions on which Theis’ solution is based.
Mishra Chachadi’s type curves (Mishra, Chachadi, 1986) could be used for
determining parameters of a confined aquifer if the test is conducted in a large diameter
well. These type curves include both the pumping phase and recovery the recovery phase
To avoid human error while curve matching, the inverse problem could be solved
conveniently using Marquardt Algorithm (Marquardt, 1963). Berg (1971), and Chander
et.al.(1981) have used the algorithm to predict parameters of aquifers in different hydro
geological settings.
⎧ Q ⎫⎡ ∞ e−u ⎤
s (Q, T , φ , r , t ) = ⎨ ⎬ ⎢∫ 2 du ⎥ = F1 ( Q, r , t , T , φ ) ………. (A11.1)
⎩ 4π T ⎭ ⎣ r φ / ( 4Tt ) u ⎦
∂s ( Q, T , φ , r , t ) ⎧ −Q ⎫ ⎡ ∞ e −u ⎤ ⎧ Q ⎫ ⎡ ⎛ e − r φ /(4Tt ) ⎞ ⎛ r 2φ ⎞ ⎤
2
2 ⎬ ⎢∫ 2
=⎨ du ⎥ + ⎨ ⎬ ⎢− ⎜ 2 ⎟ ⎜ − 2 ⎟⎥
∂T ⎩ 4π T ⎭ ⎣ r φ /(4Tt ) u ⎦ ⎩ 4π T ⎭ ⎢⎣ ⎜⎝ r φ / (4Tt ) ⎟⎠ ⎝ 4T t ⎠ ⎥⎦
⎧ Q ⎫ ⎡ − r 2φ /(4Tt ) ∞ e−u ⎤
=⎨ 2 ⎬⎢
e −∫2 du ⎥ = F2 ( Q, r , t , T , φ ) ….(A11.2)
⎩ 4π T ⎭ ⎣ r φ /(4Tt ) u
⎦
Partial derivative of s (Q, T , φ , r , t ) with respect to φ is
∂s ( Q, T , φ , r , t ) ⎧ Q ⎫⎡ ⎛ e
− r 2φ / ( 4Tt ) ⎞
⎛ r 2 ⎞⎤ Q − r 2φ / ( 4Tt )
=⎨ ⎬ ⎢ − ⎜⎜ 2 ⎟⎜ ⎟ ⎥=− e = F3 ( Q, r , t , T , φ ) …. (A11.3)
∂φ ⎩ 4π T ⎭ ⎢⎣ ⎝ r φ / ( 4Tt ) ⎟⎠ ⎝ 4Tt ⎠ ⎥⎦ 4π T φ
Let T * and φ * be the approximate values of transmissivity and storage coefficient near to
the true values of transmissivity and storage coefficient. Initially, T * and φ * are to be
guessed. Let ΔT and Δφ be incremental values in transmissivity and storage coefficient
so that T * + ΔT and φ * + Δφ are nearer to the true values. ΔT and Δφ are unknown and
are to be predicted by Marquardt algorithm. Performing Taylor series expansion of
drawdown s (Q, T , φ , r , t ) about T * and φ * , and neglecting the higher order terms, the
following equation is obtained:
∂s ( Q, T , φ , r , t ) ∂s ( Q, T , φ , r , t )
s (Q, T , φ , r , t ) = s (Q, T , φ , r , t ) + ΔT + Δφ
T * +ΔT ,φ * +Δφ T * ,φ *
∂T T * ,φ *
∂φ T *φ *
( ) ( ) (
= F1 Q, r , t , T * , φ * + F2 Q, r , t , T * , φ * ΔT + F3 Q, r , t , T * , φ * Δφ ……..(A11.4) )
( ) ( ) ( ) (
F1 t , T * , φ * , F2 Q, r , t , T * , φ * by F2 t , T * , φ * and F3 Q, r , t , T * , φ * by F3 t , T * , φ * . ) ( )
(
at time ti . The predicted drawdown , sc ti , T * + ΔT , φ * + Δφ , at observation time ti from )
equation (A11.4) is
( ) ( ) (
sc ti , T * + ΔT , φ * + Δφ = F1 ti , T * , φ * + F2 ti , T * , φ * ΔT + F3 ti , T * , φ * Δφ ) ( ) …….(A11.5)
(
E ( i ) = so ( i ) − sc ti , T * + ΔT , φ * + Δφ )
{ ( ) (
= so ( i ) − F1 ti , T * , φ * + F2 ti , T * , φ * ΔT + F3 ti , T * , φ * Δφ ) ( ) } …….. (A11.6)
The Marquardt algorithm minimizes sum of the squares of error for a set of N
observations and the minimization problem is,
⎧N
{ ( ) ( ) ⎫
( ) }
2
Min ⎨∑ ⎡ s0 ( i ) − F1 ti , T * , φ * + F2 ti , T * , φ * ΔT + F3 ti , T * , φ * Δφ ⎤ ⎬ …….. (A11.7)
ΔT , Δφ
⎩ i =1 ⎣ ⎦ ⎭
Differentiating sum of the squares of the error with respect to ΔT and equating it to zero
∑ {−2 ⎡⎣ s ( i ) − {F ( t , T ) ( ) ( ) } ( )} = 0
N
*
, φ * + F2 ti , T * , φ * ΔT + F3 ti , T * , φ * Δφ ⎤ F2 ti , T * , φ * .. (A11.8)
i =1
0 1 i
⎦
⎧N ⎡ ⎤ ⎫⎬ ΔT + ⎧⎨∑ ⎡ F3 ti , T * , φ * F2 ti , T * , φ *
( ) ( ) ( ) ( )⎤⎦ ⎫⎬ Δφ
N
⎨∑ ⎣ F2 ti , T , φ F2 ti , T , φ
* * * *
⎩ i =1 ⎦⎭ ⎩ i =1 ⎣ ⎭
( ) ( ) ( )
N N
= ∑ ⎡⎣ s0 ( i ) F2 ti , T * , φ * ⎤⎦ − ∑ ⎡⎣ F1 ti , T * , φ * F2 ti , T * , φ * ⎤⎦ ……. … (A11.9)
i =1 i =1
∑ {−2 ⎡⎣ s ( i ) − {F ( t , T ) ( ) ( ) } ( )} = 0 … (A11.10)
N
*
, φ * + F2 ti , T * , φ * ΔT + F3 ti , T * , φ * Δφ ⎤ F3 ti , T * , φ *
i =1
0 1 i
⎦
⎨∑ ⎣ F2 ti , T , φ F3 ti , T , φ , φ * ⎤⎦ ⎬ Δφ
* * * * * * *
3 i 3 i
⎩ i =1 ⎭ ⎩ i =1 ⎭
( ) ( ) ( )
N N
= ∑ ⎡⎣ s0 ( i ) F3 ti , T * , φ * ⎤⎦ − ∑ ⎡⎣ F1 ti , T * , φ * F3 ti , T * , φ * ⎤⎦ ………. (A11.11)
i =1 i =1
( ) ( ) ( ) ( )
N N
a (1,1) = ∑ ⎡⎣ F2 ti , T * , φ * F2 ti , T * , φ * ⎤⎦; a (1, 2) = ∑ ⎡⎣ F3 ti , T * , φ * F2 ti , T * , φ * ⎤⎦
i =1 i =1
( ) ( ) ( )
N N
c(1) = ∑ ⎡⎣ s0 ( i ) F2 ti , T * , φ * ⎤⎦ − ∑ ⎡⎣ F1 ti , T * , φ * F2 ti , T * , φ * ⎤⎦
i =1 i =1
and
( ) ( ) ( ) ( )
N N
a (2,1) = ∑ ⎡⎣ F2 ti , T * , φ * F3 ti , T * , φ * ⎤⎦; a (2, 2) = ∑ ⎡⎣ F3 ti , T * , φ * F3 ti , T * , φ * ⎤⎦
i =1 i =1
( ) ( ) ( )
N N
c(2) = ∑ ⎡⎣ s0 ( i ) F3 ti , T * , φ * ⎤⎦ − ∑ ⎡⎣ F1 ti , T * , φ * F3 ti , T * , φ * ⎤⎦
i =1 i =1
This iteration procedure is to be repeated till ΔT and Δφ tend to very small values.
( ) (
The functions F1 Q, r , ti , T * , φ * and F2 Q, r , ti , T * , φ * are evaluated for different )
observation time ti using a convolution technique. Let the observation time span be
discretized by uniform time steps. In case of aquifer test, a convenient time step size is
one minute as observations are generally made at different intervals of minutes.
( ) (
Accordingly, in the functions F1 Q, r , ti , T * , φ * , F2 Q, r , ti , T * , φ * and F3 Q, r , ti , T * , φ * , ) ( )
the unit of pumping rate is 1 m3 per minute, and unit of transmissivity is m 2 per minute.
Let Q ( i ) be pumping rate during i th time period. During recovery period, the pumping
r 2φ * m −1 r 2φ *
m − * − *
1 1
δ (m) = ∫ e 4T ( m −τ )
dτ − ∫ e 4T ( m −1−τ )
dτ
0
4π T *
( m − τ ) 0
4π T * ( m − 1 − τ )
1 ⎛ r 2φ * ⎞ 1 ⎛ r 2φ * ⎞
= W ⎜ * ⎟ − W ⎜ ⎟⎟ …………………….. (A11.19)
⎜
4π T ⎝ 4T m ⎠ 4π T ⎝ 4T ( m − 1) ⎠
* * *
⎛ r 2φ * ⎞ r 2φ * ⎛ r 2φ * ⎞
W ⎜ * ⎟ is Theis’ Well function with argument . W ⎜ * ⎟ is an exponential
⎝ 4T m ⎠ 4T *m ⎝ 4T m ⎠
integral. For argument X, the exponential integral, W ( X ) , is computed using the following
polynomial approximation.
For X ≤ 1
W ( X ) = − ln ( X ) − 0.57721566 + 0.99999193X − 0.24991055X 2
For X > 1
X 4 +8.5733287X 3 + 18.059017X 2 + 8.6347608X +0.26777373
Xe X W ( X ) = ….(A11.21)
X 4 + 9.5733223X 3 + 25.632956X 2 + 21.099653X +3.9584969
An Example
A set of synthetic observation data generated using T = 0.1m 2 / min and storage
coefficient φ = 0.001 are as given in Table 1. Predict the T , φ making an initial guess
T * = 0.01m 2 / min and φ * = 0.003 . The piezometer is located at a distance of 20m from
Iteration T
*
φ* ΔT Δφ C(1) C(2)
no
1 0.020507 0.004138666 1.05E-02 1.14E-03 0.00E+00 3.64E-12
2 0.040832 0.003697811 2.03E-02 -4.41E-04 -1.14E-13 0.00E+00
3 0.07088 0.001491592 3.00E-02 -2.21E-03 1.42E-14 0.00E+00
4 0.092492 0.001075698 2.16E-02 -4.16E-04 0.00E+00 -1.14E-13
5 0.099488 0.001004064 7.00E-03 -7.16E-05 0.00E+00 1.42E-14
6 0.099998 0.001000018 5.10E-04 -4.05E-06 0.00E+00 0.00E+00
7 0.1 0.001 2.42E-06 -1.81E-08 0.00E+00 0.00E+00
8 0.1 0.001 5.43E-11 -4.07E-13 0.00E+00 2.12E-22
9 0.1 0.001 2.98E-16 -9.44E-18 0.00E+00 -8.08E-28
Thus using synthetic drawdown data, we have checked that Marquardt Algorithm
successfully predicts the true transmissivity and storage coefficient when the initial guess
was different from the true value.
In hard rock region, the shallow aquifers have low transmissivity ranging from 25
to 100 m 2 / day Therefore, in hard rock region, wells with diameter ranging from 1 to 2
m are constructed to have reasonable yield. An aquifer test can be conducted in a large
Solution to unsteady flow to a dug-cum-bore well, that takes well storage into
account, has been derived by Papadopulos and Cooper (1967). According to them, the
(
well storage dominates the time-drawdown curve up to a time, t, given by t = 25rc2 / T)
where rc is radius of the well casing, and T is transmissivity of the aquifer. If a short
duration aquifer test is conducted in a large diameter well, the transmissivity can be
estimated reasonably well, but the storage coefficient may differ by an order of
magnitude. This is because, the type curves presented by Papadopulos and Cooper
contain straight line portions, which are parallel, and a short duration time-drawdown
curve if matched with one of the straight lines, it could be matched as well with either of
the adjacent straight lines. Discretising the time domain by uniform time steps, and
generating unit response function coefficients from Thies’ basic solution for unsteady
flow to a well with small radius, Patel and Mishra (1983), Mishra and Chachadi (1985)
have derived simple analytical solutions to unsteady flow during pumping, and during
recovery respectively. These solutions, as well as that by Papadopulos and Cooper are
applicable for a bore well with small radius having large casing.
Hantush has derived an analytical solution to unsteady flow to a well with finite
radius assuming that all the water pumped is from aquifer storage. The effect of well
storage on time drawdown curve has not been taken into account in the solution.
Discretising the time domain by uniform time steps, and generating unit response
function coefficients from Hantush’s basic solution, we derive a simple analytical
solutions to unsteady flow to a large diameter well during pumping and recovery. The
well storage contribution during pumping, and well storage effect on drawdown have
been accounted. After solving the direct problem, the inverse problem has been solved
using the Marquardt Algorithm as described below. Pumping as well as recovery data
could be used for estimating the aquifer parameters.
for n > m, Q p ( n ) = 0 .
Drawdown, sw ( n ) ,at the well face at the end of time step n is given by (Patel and
Mishra, 1983)
n
1
sw ( nΔt ) =
π rc2
∑
γ
Q ( γ ) Δt
=1
w ………………….………. (A11.24)
in which, Qw ( γ ) is average withdrawal rate from well storage during time step γ . Qw ( γ )
storage from aquifer storage during time of recovery. rc is radius of well casing. For
Following Duhamels’ principle and method of convolution, drawdown at the well face at
the end of time step n ( t = nΔt ) due to abstraction from aquifer storage is given by
(Morel-Seytoux, 1975)
where δ ( rw , Δt , N ) is a unit pulse response function coefficient derived from unit step
response function using a time step size Δt (Morel-Seytoux, 1975) , and N is an integer.
The kernel coefficient δ ( rw , Δt , N ) is given by
δ ( rw , Δt , N ) =
1
Δt
{ }
U ( rw , N Δt ) − U ( rw , ( N − 1) Δt ) ; N > 1 …………. (A11.26)
1
For N = 1, δ ( rw , Δt ,1) = U ( rw , Δt ) .
Δt
pumping rate which could be computed using Theis’ solution. The time step size Δt can
be chosen conveniently. In an aquifer test, draw down observations are made at different
intervals of minutes. Therefore, for solving an inverse problem, it is convenient to choose
Δt = 1 minute. Accordingly, pumping rate is to be chosen in m3 per minute and
transmissivity in m 2 per minute. For a well with finite radius, the unit step response
function has been derived by Hantush as given in Annexure A11.1
Assuming that there is no surface of seepage at the well face, the drawdown in the
well is equated to the drawdown in the aquifer at the well face i.e. sw ( n ) = S a ( rw , n ) .
π rc2
∑ Qw (γ ) = ∑ Qa (γ ) δ ( rw , Δt , n − γ + 1) ……………………..
γ =1 γ =1
(A11.27)
Qw ( n ) = Q p ( n ) − Qa ( n ) ……………………… (A11.28)
Δt
{Q (γ ) − Q (γ )} = ∑ Q (γ ) δ ( r , Δt, n − γ + 1) ……………
n n
π rc2
∑
γ=1
p
γ
a
=1
a w (A11.29)
the second part the nth step, incorporating Δt = 1 minute and solving for Qa (n)
n n −1 n −1
∑
γ
Q ( γ ) − ∑ Q ( γ ) − π r ∑ Q ( γ ) δ ( r , Δt , n − γ + 1)
=1
p
γ =1 γ
a c
2
=1
a w
Q p (1)
Qa (1) = …………………. (A11.31)
1 + π rc2δ ( rw , Δt ,1)
2Q p (1)
Qa (1) = ………………………. (A11.32)
1 + π rc2δ ( rw , Δt ,1)
Qa (n) , n = 1,..., N are solved in succession. After solving Qa (n) , for n=1, 2,…,N, the
drawdown in the aquifer at any distance r is found generating the corresponding kernel
⎡ 1 ⎤
coefficients δ ( r , Δt , N ) ⎢ =
⎣ Δt
{
U ( r , N Δt ) − U ( r , ( N − 1) Δt ) ⎥
⎦
} and applying the
convolution technique.
Having solved the direct problem, the inverse problem is solved next.
⎧N
{ (
Min ⎨∑ s0 ( i ) − sc ti , T * + Δt , φ * + Δφ )} ⎫⎬
2
………………. (A11.33)
ΔT , Δφ
⎩ i =1 ⎭
∂sc ( ti , T , φ )
=
( ) (
sc ti , T * , φ * + ε 2 − sc ti , T * , φ * ) …………... (A11.34c)
∂φ T * ,φ *
ε2
where ε1 , ε 2 are small increments in transmissivity and storage coefficient.
( )
Incorporating sc ti , T * + ΔT , φ * + Δφ i.e. equation (A11.34a), in equation (A11.33), the
minimization problem reduces to
⎡N ⎧ ⎛ ⎞ ⎪⎫ ⎤
2
⎪ ∂s (
Min ⎢ ∑ ⎨ s0 ( i ) − ⎜ sc ( ti , T , φ ) T * ,φ * + c i
t , T , φ ) ∂s (
ΔT + c i
t , T , φ ) Δφ ⎟ ⎬ ⎥ ….. (A11.35)
ΔT , Δφ ⎢ ⎜ ∂ T ∂ φ ⎟ ⎥
⎢⎣ ⎪⎩ ⎝ ⎠ ⎪⎭ ⎥⎦
i =1 T * ,φ * T * ,φ *
Equating the partial derivative of the above objective function with respect to ΔT with
zero
⎡N ⎧ ⎛ ⎞ ⎪⎫ ∂s ( t , T , φ ) ⎤
⎢ ∑ −2 ⎪⎨ s0 ( i ) − ⎜ sc ( ti , T , φ ) * * + c ( i
∂s t , T , φ ) ∂sc ( ti , T , φ )
Δ T + Δ φ ⎟⎬ c i ⎥ = 0 (15)
⎢ i =1 ⎪ ⎜ T ,φ
∂T ∂φ ⎟ ∂T ⎥
⎣ ⎩ ⎝ T * ,φ * T * ,φ * ⎠ ⎪⎭ T * ,φ * ⎦
Simplifying
⎡ N ⎧⎛ ⎞ ⎫⎪⎤
2
⎢ ⎪ ∂s ( t , T , φ ) ⎟ ⎬ ⎥ ΔT
⎢∑ ⎨⎜
c i
i =1 ⎪⎜ ∂ T ⎟ ⎥
T * ,φ * ⎠ ⎪
⎣⎢ ⎩⎝ ⎭⎦⎥
⎡ N ⎧⎛ ∂s t , T , φ ⎞ ⎫⎪ ⎤
⎪ c( i ) ∂sc ( ti , T , φ )
+ ∑⎨
⎢ ⎜ ⎟ ⎬ ⎥ Δφ
⎢ i =1 ⎪⎜ ∂φ ∂ T ⎟ ⎥
T ,φ ⎠ ⎪
⎣ ⎩⎝ T ,φ ⎭⎦
* * * *
N ⎧ ∂s ( t , T , φ ) ⎫⎪ N ⎧⎪ ∂sc ( ti , T , φ ) ⎫⎪
⎪
= ∑ ⎨ s0 ( i ) c i ⎬ −∑ ⎨ sc ( ti , T , φ ) T * ,φ * ⎬ …….. (A11.37a)
i =1 ⎪ ∂T * * ⎪ i =1 ⎪ ∂T * * ⎪
⎩ T ,φ ⎭ ⎩ T ,φ ⎭
or
a (1,1)ΔT + a (1, 2)Δφ = c (1) …………………… (A11.37b)
where
N ⎧ ∂s ( t , T , φ ) ⎫⎪ N ⎧⎪ ∂sc ( ti , T , φ ) ⎫⎪
⎪
c(1) = ∑ ⎨ s0 ( i ) c i ⎬ −∑ ⎨ sc ( ti , T , φ ) T * ,φ * ⎬
i =1 ⎪ ∂T * * ⎪ i =1 ⎪ ∂T * * ⎪
⎩ T ,φ ⎭ ⎩ T ,φ ⎭
Similarly, equating the partial derivative of the objective function with respect to Δφ with
zero
⎡N ⎧ ⎛ ⎞ ⎫⎪ ∂s ( t , T , φ ) ⎤
⎢ ∑ −2 ⎪⎨ s0 ( i ) − ⎜ sc ( ti , T , φ ) * * + c ( i
∂s t , T , φ ) ∂sc ( ti , T , φ )
ΔT + Δφ ⎟ ⎬ c i ⎥ = 0 (17)
⎢ i =1 ⎪ ⎜ T ,φ
∂T ∂ φ ⎟ ∂ φ ⎥
⎣ ⎩ ⎝ T * ,φ * T * ,φ * ⎠ ⎪⎭ T * ,φ * ⎦
Simplifying
⎡ N ⎧⎛ ∂s t , T , φ ⎞ ⎫⎪⎤ ⎡ N ⎧⎛ ⎞ ⎫⎪⎤
2
⎪
⎢ ∑ ⎨⎜ c i ( ) ∂sc ( t i , T , φ ) ⎪
⎟ ⎬⎥ ΔT + ⎢ ∑ ⎨⎜ c i
∂ s ( t , T , φ ) ⎟ ⎬⎥ Δφ
⎢ i =1 ⎪⎜ ∂φ ∂ T ⎟ ⎥ ⎢ ⎜ ∂ φ ⎟ ⎥
T * ,φ * ⎠ ⎪
⎣ ⎩⎝ ⎢⎣ ⎪⎩⎝ T * ,φ * ⎠ ⎪
i =1
T * ,φ * ⎭⎦ ⎭⎥⎦
N ⎧ ∂s ( t , T , φ ) ⎫⎪ N ⎧⎪ ∂sc ( ti , T , φ ) ⎫⎪
⎪
= ∑ ⎨ s0 ( i ) c i ⎬ −∑ ⎨ sc ( ti , T , φ ) T * ,φ * ⎬ …. (A11.39a)
i =1 ⎪ ∂φ * * ⎪ i =1 ⎪ ∂φ * * ⎪
⎩ T ,φ ⎭ ⎩ T ,φ ⎭
or
a (2,1)ΔT + a (2, 2)Δφ = c ( 2 ) …………………. (A11.39b)
where
⎧⎪⎛ ∂s ( t , T , φ ) ⎞ ⎫⎪ N ⎧⎛ ⎞ ⎫⎪
2
N
a (2,1) = ∑ ⎨⎜ c i
∂sc ( t i , T , φ ) ⎪
⎟ ⎬; a (2, 2) = ∑ ⎨⎜ c i
∂s ( t , T , φ ) ⎟ ⎬ and
i =1 ⎪ ⎜ ∂ φ ∂ T ⎟⎪ ⎪⎜ ∂ φ ⎟
⎩⎝ T ,φ ⎠ ⎭ T ,φ ⎠ ⎪
⎩⎝
T ,φ i =1
⎭
* * * * * *
N ⎧ ∂s ( t , T , φ ) ⎫⎪ N ⎧⎪ ∂sc ( ti , T , φ ) ⎫⎪
⎪
c(2) = ∑ ⎨ s0 ( i ) c i ⎬ −∑ ⎨ sc ( ti , T , φ ) T * ,φ * ⎬
i =1 ⎪ ∂φ * * ⎪ i =1 ⎪ ∂φ * * ⎪
⎩ T ,φ ⎭ ⎩ T ,φ ⎭
The unknown ΔT and Δφ are solved from the algebraic equations (A11.37b) and
(A11.39b)
guess T * = 0.01m 2 / min and φ * = 0.2 . The drawdown is measured in the large diameter
C. Discrete Kernel, δ ( rw , Δt , N )
1 ⎡4 ∞ ⎤
U ( rw , N Δt ) = ⎢ ∫ {1 − exp(−τ x )} f1 ( x)dx ⎥
2
……………….… (A11.41)
4π T ⎣π 0 ⎦
in which,
Tt J ( x )Y0 ( ρx ) − J 0 ( ρx )Y1( x ) r
τ=
ϕ rw 2
; t = N Δt ; f1 ( x ) = 1
2
[
2 2
x J 1 ( x ) + Y1 ( x ) ] ; ρ=
rw
= 1 ; J0(x), J1(x) =
Bessel functions of first kind of zero and first order respectively; Y0(x) Y1(x)= Bessel
functions of second kind of zero and first order respectively; T= transmissivity (m2/day),
and φ=storativity of the aquifer; rw = radius of the well or shaft(m).
[ ]
∞
I= ∫ 1 − exp( −τx 2 ) f1( x )dx
0
[ ] [ ]
1 ∞
= ∫ 1 − exp( −τx 2 ) f1( x )dx + ∫ 1 − exp( −τx 2 ) f1( x )dx
0 1
=I1+I2
[ ]
1⎡
− τ ( 1 + v )2 ⎤
1
1+ v
I1= ∫ 1 − exp( −τx ) f1( x )dx = 0.5 ∫ ⎢1 − exp{
2
}⎥ f1 ( )dv
−1⎣⎢ ⎦⎥
4 2
0
[ ] [ ]
∞ 1
dv
I2= ∫ 1 − exp( −τx 2 ) f1( x )dx = ∫ 1 − exp( −τ / v 2 ) f1 ( 1 / v )
1 0 v2
1
⎡ ⎧ −4τ ⎫⎤ 2 4dy
= 0.5 ∫ ⎢1 − exp ⎨ 2 ⎬⎥ 1
f( )
⎩ (1 + y ) ⎭⎦ 1 + y (1 + y )
2
−1 ⎣
⎡ ⎧⎪ − 4τ ⎫⎪⎤
As y→ -1, ⎢1 − exp ⎨ ⎬⎥ → 1
⎢⎣ ⎪⎩( 1 + y ) 2 ⎪⎭⎥⎦
2 2 2 2
J1( )Y0 ( ρ ) − J0( ρ )Y1( )
4 2 ⎡ 4 ⎤ 1+ y 1+ y 1+ y 1+ y
[ ] f1 ( )= ⎢ ⎥
(1 + y )2 1 + y ⎢⎣ ( 1 + y ) 2 ⎥⎦ 2 2⎡ 2 2 2 ⎤
) ⎢ J1 ( ) + Y12 (
1 + y ⎥⎦
( )
1+ y ⎣ 1+ y
2 2 2 2
J1( )Y0 ( ρ ) − J0( ρ )Y1( )
1+ y 1+ y 1+ y 1+ y
=
⎡ 2 2 2 2 ⎤
⎢ J1 ( 1 + y ) + Y1 ( 1 + y )⎥
⎣ ⎦
2 2 2 2
J1( )Y0 ( ρ ) J0( ρ )Y1 ( )
1+ y 1+ y 1+ y 1+ y
= −
⎡ 2 2 2 2 ⎤ ⎡ 2 2 2 2 ⎤
⎢ J1 ( 1 + y ) + Y1 ( 1 + y )⎥ ⎢ J1 ( 1 + y ) + Y1 ( 1 + y )⎥
⎣ ⎦ ⎣ ⎦
2
As y→-1, Y1( ) →0 ; hence,
1+ y
2 2
J1( )Y0 ( ρ )
1+ y 1+ y
⎡ 2 2 2 2 ⎤
⎢ J1 ( 1 + y ) + Y1 ( 1 + y )⎥
⎣ ⎦
=1 (since ρ =1)
Similarly,
2 2
J0( ρ )Y1( )
1+ y 1+ y
→1
⎡ 2 2 2 2 ⎤
⎢ J1 ( 1 + y ) + Y1 ( 1 + y )⎥
⎣ ⎦
Location EC*106
S.NO. micro-
siemens/cm
TDS pH Na+ K+ Ca++ Mg++ Cl- SO42- Co32- HCo32- No3-
at 250c
Chand
1 Baori 440 240 7.9 14 5 62 12 28 0 0 207 14
2 Raghu nath 1220 645 7.7 62 16 108 44 85 34 0 500 46
Baori
3 Suraj kund 680 454 8.1 30 11 104 19 50 53 0 342 16
4 Ram Baori 450 252 8.1 17 30 46 12 28 0 0 220 8
Table A13.1 : Details of the small water bodies in and around the Jodhpur city
(Reference : Section 13.0).
S.No Waterbody Name Longitude Latitude Perimeter (km) Area (Sq. km)