CE 394k GIS in water resources

FINAL PROJECT REPORT

Beatriz Garcia-Fresca

12/5/02

 

 

POTENTIAL RECHARGE TO GROUNDWATER FROM MAINS LEAKAGE: AUSTIN, TEXAS (USA)

 

 

 

 

1. Introduction

2. Data collection and treatment

3. Spatial analysis and calculations

    leakage

    relationship with faulting

4. Results and discussion

5. Conclusions

6. Future work

7. References

 

Appendix: Building a Hydrologic Network

    elevation

    flow direction

    tracing

    discussion

 

 

 

 

 

 

 

 

 

 

 


 

1. INTRODUCTION

 

 

Humans are the number one geologic and environmental agent on the surface of earth. The effects are most notorious where population concentrates, i.e. in urban areas.

The urban underground is a complex network of buried structures, pipes, tunnels, etc., analogous to an epikarstic system (Krothe et al., 2002).

It is broadly recognized that recharge of rainfall to groundwater is inhibited in urban areas due to the presence of impervious cover. The present study discusses the fact that urban development can also produce the opposite effect, i.e. the increase of groundwater recharge. Water and wastewater mains leak, and lawns and parks are often over-irrigated, resulting in a potential increase on recharge to the groundwater, from strictly urban sources (e.g. Lerner, 1986, and Foster, 1996). Increases in recharge have been identified in cities all around the world, as shown in figure 1.

This study covers the Austin metropolitan area, which sits on the inactive Balcones fault zone, in central Texas (USA). The city occupies most of Travis County, and small parts of the surrounding counties.

Figure 2: location of the study area

 

ArcGIS is used to determine the leakage form the water distribution network to the Austin subsurface. The water distribution network will be superimposed to the geology of the Austin area in order to assess the potential recharge from the mains to the different hydrogeological units. The City of Austin recognizes a water loss in the distribution lines of 12% of the delivered water as shown in the newspaper clip in figure 3 (Austin American-Statesman, 1998).

 

Figure 3: unaccounted for water in Texas

 

The potential for the leakage to recharge to the groundwater depends to some extent on the geologic materials underlying the mains, as well as faulting.

 


 

2. DATA COLLECTION AND TREATMENT

 

Initially a large and varied amount of datasets were collected (see Future work section), the ones used in this study are the following:

Source: UT BEG (Trembley, Andrews and Dale, after Garner and Young, 1976)

Description: lithologic contacts and faults of the Austin Area

Figure 4: Geology dataset

Click here to see the footprint of the city over the geologic map

 

Source: City of Austin GIS datasets

Description: line and point shapefiles of the distribution lines. Attributes include year installed, diameter, material, etc.

Figure 5: water distribution network dataset

 

Source: City of Austin GIS datasets

Description: polygon shapefiles of the outlines of the different city jurisdictions. Only the “full purpose” jurisdiction was used.

Figure 6: Austin jurisdiction limits dataset

 

The principal data processing requirement consisted of representing all datasets in the same projection. The NAD 86 State Plane Texas Central (US ft) coordinate system was chosen, which was designed for this particular area and therefore minimizes distortion.

 

Minor processing consisted on assigning the standard colors to each of the lithologies displayed in the geologic map (figure 4).

 

 


 

3. SPATIAL ANALYSIS AND CALCULATIONS

 

Recharge to the groundwater depends, in part, on the permeability of the outcropping materials. For this reason, the first step is to reclassify the geologic formations into permeability units, by editing the attribute table of the shapefile containing the geology. The permeabilities assigned are listed in table 1, and correspond to the qualitative descriptions of the rocks by Duffin (1991), and Brune & Duffin (1983). The spatial distribution of the permeability units is shown in figure 7.

 

symbol

name

Age

description

relative permeability

ID

Qal

Alluvium

Recent

sand, silt, clay and gravel; tan to light gray

low to high

4

Qlcr

Lower Colorado River terrace deposits

Pleistocene

sand, silt, clay and gravel; yellow to orange-brown

low to medium

3

Qucr

Upper Colorado River terrace deposits

Pleistocene

gravel, sand, silt and clay; orange-brown

low to medium

3

Qtt

Tributary terrace deposits

Pleistocene

gravel, sand, silt and clay; tan to light gray

low to medium

3

Qht

High terrace deposits

Pleistocene

gravel, sand, silt and clay; gray to tan

low to medium

3

Emi

Midway Group

Eocene

clay dark gray to brown-gray, sandy, micaceous, glauconitic; contains calcareous and ferruginous concretions

very low

1

Kna

Navarro Group

Cretaceous

clay, dark gray to brown, silty, calcareous, montmorillontic; contains locally interbedded sandy layers and calcareous concretions

very low

1

Kta

Taylor Group

Cretaceous

clay, dark gray to green gray, calcareous,  montmorillonitic; generally more calcareous in mid-portion of unit

very low

1

Kau

Austin Group

Cretaceous

chalk, marly limestone and limestone, light gray, soft to hard, thin to thick bedded, massive to slightly nodular

low

2

Kpt

Pilot Knob Tuff

Cretaceous

alterd tuff, green-brown to tan, nontronitic

very low

1

Kpb

Pilot Knob Basalt

Cretaceous

Basalt, black to dark green-gray, hard, fine grained

low

2

Kef

Eagle Ford Formation

Cretaceous

clay, dark gray, calcareous; contains sandy and silty flaggy limestone in mid-portion and a bentonite bed at base

impermeable

0

Kbu

Buda Formation

Cretaceous

limestone, gray to tan, hard, dense, slightly nodular, abundant fossil mollusks

impermeable

0

Kdr

Del Rio Formation

Cretaceous

clay, dark gray to olive brown, pyritic, gypsiferous, calcareous; contains abundant Exogyra arietina

impermeable

0

Kgt

Georgetown Formation

Cretaceous

limestone and marly limestone, gray to tan, hard to soft, bedded to nodular; contains abundant fossil mollusks

low to high

4

Ked

Edwards Formation

Cretaceous

limestone and dolomite, light gray to tan, hard to soft, thin to thick bedded, fine to medium grained; fossil rudist and nodular chert common; solution collapse zone near middle

low to high

4

Kcp

Comanche Peak Formation

Cretaceous

limestone, gray to tan, soft, marly, nodular, fine grained

low

2

Kwa

Walnut Formation

Cretaceous

limestone, marl, and marly limestone, gray to tan, soft to hard, thick to thin bedded, massive to nodular, fine to medium grained

low

2

Kgr

Glen Rose Formation

Cretaceous

limestone, dolomite, and marl, gray to tan, alternating hard and soft beds forming stair-step topography, thick to thin bedded, fine to medium grained

low to medium

3

Table 1: characteristics of the rocks outcropping in the study area

 

Figure 7: relative permeabilities of the geologic materials in the study area

 

Click here to see the footprint of the city over the permeability map.

 

The Edwards Group, Georgetown Formation, and the alluvium are the potentially most permeable materials. The former two form the Edwards Aquifer. The Glen Rose Limestone may not be as locally permeable as the Edwards and Georgetown limestones, but it contributes its drainage to the recharge of the Edwards Aquifer. River deposits can also present important permeability values. According to this interpretation, three different hydrogeologic environments can be distinguished:

 

 Cretaceous Karstic limestones: very permeable locally, and susceptible to pollution. The outcrop of the aquifer is shown in figure 8.

 

Figure 8: outcrop of the Edwards aquifer in the Austin area

 

Low permeability materials don’t allow the direct infiltration of water to the underground but could contribute drainage to adjacent aquifer units, and baseflow to creeks.

 

Quaternary deposits overly both limestone and impermeable materials (figure 9). In the study area they consist on fluvial alluvium and terrace deposits, which can reach significant permeability values. They are often closely related to surface waters, and can present perched water tables.

 

Figure 9: Quaternary deposits in the Austin area

 

 

leakage

 

This section describes the calculations performed on ArcGIS to estimate the potential leakage from mains.

The water loss for the city of Austin is 12% of the water pumped into the distribution system. According to Thornton (2002), about 60% of the losses are “real losses”, and the remaining fraction “apparent losses”. The real losses correspond to the water lost to pipe leakage and breaks, and doesn’t include metering unbalances, thefts, etc. According to this, this study will adopt the value of 7% real losses, or a “leakage coefficient” (ℓ) of 0.07.

Water usage data from the City of Austin Water and Wastewater Service was used as the input to the calculations. Table 2 shows the historical water demand for Austin, for the period between 1993 and 2001.

 

Year

Pumpage

Peak

(millions of gallons)

Day

1993

39,824

189

1994

39,806

199

1995

39,542

192

1996

45,835

205

1997

42,812

195

1998

46,438

206

1999

46,422

211

2000

52,193

226

2001

50,140

243

Table 2: historical water demand, Austin, Texas

Source: City of Austin Water and Wastewater Utility

 

The input total annual pumpage chosen was that of 2001, i.e. 50,140 millions of gallons, which is an average of 364,442 m3/day. The volume of the mains the product of the diameter and length, both attributes of the dataset. The average daily average pumpage of water is then divided by the total volume of mains, and multiplied by the volume of each pipe, to obtain the daily volume of water delivered by each pipe. This value is then multiplied by the leakage coefficient to obtain the volumetric leakage per pipe. These calculations can be summarized by the following expression:

 

L = (xi · Di · PTA · ℓ) / (S [xi · Di] · 365)                                                                          (1)

 

where,

xi – length of each pipe

Di – diameter of each pipe

PTA – total annual pumpage

ℓ – leakage coefficient

 

A useful figure is the leakage per unit length of pipe (LLi), defined as:

 

LLi = L/xi                                                                                                                               (2)

 

Such simple calculations could be comfortably carried out in a spreadsheet, except for the fact that these particular datasets have too many records for the capacity of the Microsoft Excel spreadsheet.

 

 

relationship with faulting

 

Austin sits on the inactive Balcones Fault Zone, which is oriented NE-SW (see figure 4). Faults represent potential preferential pathways for water to reach the subsurface. Water leaking immediately above a conductive fracture is likely to recharge the aquifers below. ArcGIS allows the superposition of pipes and faults (figure 10), and the identification of the lines crossing above fractures (figure 11). Then the volume of water leaking above fractures, per unit length of pipe is computed.

 

 

Figure 10: geometric relationship between distribution pipes and faults

 

 

Figure 11: distribution pipes crossing over fault traces

 

 


 

4. RESULTS AND DISCUSSION

 

Table 3 summarizes the calculations and results of this study. The real losses are determined as 7% of the 520,631.78 m3/d pumped into the distribution network. The 36,434 m3 which leak on average every day, represent a potential groundwater recharge of 21.94 mm per year, over the area covered by the city.

 

average pumpage

520631.78

m3/d

City of Austin - Water and Wastewater Utility

total V of pipes

672465.55

m3

 

 

coeff.

0.774

 

 

 

 

 

 

 

 

area of Austin (full purpose jurisdiction)

6,523,836,222

sqft

 

 

606.1

Km2

 

 

 

 

 

 

 

annual precipitation

840

mm/a

NOOA

 

 

 

 

 

total mains loss

12

%

Austin American Statesman, July 16 1998

calculated              real loss         

7

%

real loss = 60% total loss (Thornton, 2002)

36434.18359

m3/d

 

 

13298477.01

m3/a

 

 

13.30

million m3/a

 

 

21.94

mm/a

 

 

 

 

 

 

 

leakage over faults

2918.07

l/d m

 

 

12422.25

m3/d

 

 

Table 3: summary of results