ABSTRACT
Location of waste management systems is a factor affecting groundwater quality. Leakage from landfills and lagoons usually contain significant amounts of contaminants such as ammonia, nitrate, chloride and metals. These contaminants may reach aquifers and degrade the water resource such that it becomes hazardous to human health. This work is an example of how to determine acceptable areas for waste systems based on information of land use from Geographic Information Systems (GIS), hydraulic conductivity of soils and aquifer depth below the surface. The Equus Beds was the study site. The aquifer provides 50% of the water for municipal consumption for the city of Wichita, Hutchinson and towns in Reno, McPherson, Harvey and Sedgwick counties in Kansas. Also it provides water for crop developing and cattle raising as well as industrial water for private high water consumption companies. In this project we have developed maps assigning levels of contamination risk of aquifer for land zones based on characteristics of the soil such as hydraulic conductivity and aquifer depth. The assumptions, and initial and boundary conditions to calculate contaminant travel time through soil were: the system consisted of two homogeneous layers, the clay liner above and the natural soil underneath; the head above the system is constant; the seepage is vertical and governed by Darcy’s Law. The calculation of travel times was made using the data of hydraulic conductivity and aquifer depth as reference. The scheme of the problem was set up as two layers: compacted clay liner (hydraulic conductivity 10-7 cm/sec), and natural soil in which permeability values were determined from data provided by the Equus Aquifer Management District 2. After calculating the leachate travel time at hundreds of locations using the MS Excel database, risk levels of low, medium and high were assigned to those nodes. The remaining areas on the Land Use map satisfy the condition that they are at least 5 miles away from rivers, lakes and residential areas. To suggest final landfill locations, a distance of 20 miles from city centers was used to decrease transportation costs. Five percent of the studied area was determined acceptable for large waste systems.
KEYWORDS: Site characterization, GIS, waste management sites
INTRODUCTION
The use of Geographic Information Systems (GIS) in engineering situations have grown tremendously in the last 5 years. All major fields of CE can benefit by the appropriate use of GIS. Common benefits of GIS are its ability to (Miles and Ho, 1999; Parsons and Frost, 2000):
Perry et al., 1999, used GIS technology to study the effect of six different riparian-forest deforestation and reforestation management scenarios on surface water quality. A GIS based hydrologic model developed by Hellweger and Maidment, 1999, contained over 1000 inputs. Parsons and Frost, 2000, developed a geotechnical and geoenvironmental site investigation process using GIS. Complex aquifer geometry can also be characterized using GIS and water resources models (Gangopadhyay et al., 1999).
The Equus beds are an important source of high quality water for central Kansas. The aquifer provides 50% of the water for municipal consumption for the city of Wichita, Hutchinson and towns in Reno, McPherson, Harvey and Sedgwick counties. Also it provides water for crop developing and cattle raising as well as industrial water for private high water consumption companies.
Now, a project of groundwater recharge is being developed where recharge wells, lagoons and a recharge trench has been installed. Research about groundwater removal of oil field brine from the lower portion of the Equus Beds Aquifer and investigation on excessive nitrate-nitrogen in portions of the aquifer are taking place. The importance of avoiding contamination of any kind is crucial for the Equus Beds area and surroundings.
Local government policies for the Equus Beds are:
Our objective was to determine zones where waste containment is recommended for further analysis. Of course this analysis is just preliminary and should be accompanied of a carefully designed site investigation as well as a proper economic analysis.
METHODS
Site Selection and Specifications
A site investigation for a waste structure is important to ascertain the potential risk posed by the stored animal waste. Prior to onsite investigation, you should consult available geology or groundwater maps (GIS maps), published county soil surveys, previous designs in the same physiographic area and any other information that aids your assessment of the site. The intensity of investigation required depends on the experience in the given area, the types of soils and variability of the soil deposits, the size of the structure, the environmental sensitivity, and an assessment of the associated risks involved.
The subsurface investigation may employ auger holes, dozer pits, and backhoe pits. The site investigation should extent at least 2 feet below the planned bottom of the excavation. A field investigation can include site permeability testing or taking samples for laboratory testing, or may be limited to field classification of the soil. Permeability is an important factor in waste management structures and values of 1x10-6 cm/sec is a usual target value used in the design of these facilities. Records from site investigations are important, should be documented and included in the design documentation.
Compacted Clay liners
Liners are relatively impervious barriers used to reduce seepage losses to an acceptable level. Providing a liner for a waste storage structure is to improve the soil at the excavated grade by discing, watering, and compacting them to a thickness indicated by guidelines. Soils with suitable properties make excellent material for liners, but the liners must be designed and installed correctly.
Specific Discharge
The term specific discharge or unit seepage is the seepage rate for a unit cross-sectional area of a pond. This quantity is defined as follows from the Darcy’s Law:
![]() |
(1) |
where
K = Hydraulic conductivity
d = Thickness of soil liner
H = Maximum vertical distance measured from liquid level and the top of the liner
q = Q/A, and
<>Q = total seepage through area A
Modeling Work
First, it is important to state assumptions, and initial and boundary conditions to idealize the problem.
Software Used
ArcView was used in this project because it is the most popular software package. Some of its features include a well design and logically arranged user interface; wizards that simplify entering, querying manipulating and analyzing data; and tools that helps adapt the program to the needs of individual users. ArcView also includes comprehensive well-written and profusely illustrated documentation as well as an extensive online help file. ArcView’s two customization program allows you to quickly customize the interface the interface by adding or deleting buttons and attaching a script to automate a process.
Excel was used to calculate the travel times for the points in study. AutoCAD was used to resize maps and allocate points in the study.
SURFER (SURFER User's Guide. Golden Software, Inc., 1997) is one of the best map-generating software available. Its capability to interpolate values using different interpolation methods gives this software a preference to other software in the market. Its user-friendly interface allows the user to manipulate data and generated contour maps easily and effectively.
Site Selection Process
In order to calculate the contamination potential of the underground water, enough data should be provided to carry out necessary calculations such as landfill leachate travel time through the soil stratum. The data provided can be in the form of digitized data or maps. The digitized data is generally simpler and more valuable as the user can easily adjust it and edit it to complete the calculations.
There are three data sets necessary for this type of applications. The first is data regarding the soil type and soil permeability. Aquifer depth or the level at which the groundwater table exists is the 2nd set. These two data sets will be used to predict the leachate travel time. The third important data is the land-use or the land-cover.
In this study, four counties from the state of Kansas will be involved, Mc-Pherson, Reno, Harvey and Sedgwick counties. These counties are near Wichita City. Aquifer Depth, Soil Type and Land Use maps were all available at the State Data Access and Support Center (DASC) web site. The groundwater depth and hydraulic conductivity information was obtained from maps by the USGS in cooperation and the Kansas Geological Survey. The groundwater depth for the four counties is in the form of contour maps (Figure 1). The other map was the Hydraulic Conductivity map, which presents values in feet/day. For the Conductivity map, the area was divided into a set of adjacent square nodes hatched with different patterns (Figure 2). Each pattern is for a certain range of Conductivity.
Merging the DASC Files
Land Use maps were downloaded from the DASC site. Each map represents one of the aforementioned counties. Each map successfully displays and spots residential areas, grass lands, crop lands, wood lands, surface water, rivers, lakes and others. The GIS/ArcView tools allows the user to combine different maps and merge them together as long as they have same units and their location allows the merging to be carried out successfully.
Many of the new functions in ArcView are accessed through extensions that come with the software and can be turned on or off as needed. Choosing Extensions from the File menu opens the Extensions dialog box. Highlighting an extension provides information about it. Clicking on its check box loads it into ArcView. When the Geoprocessing extension is loaded, a menu choice called Geoprocessing Wizard appears on the View menu. This wizard organizes and performs six spatial data-processing tasks. Merging themes appends the features of one or more themes to those of the another theme of the same type. A new output theme is created that contains all the features of each input theme. Field names in the output theme table are taken from a selected input theme. Attributes from the other input themes are included in the output theme table as long as their field names match those of the selected input theme.
Figure 1. Contour map for groundwater depth (ft) for study region in central Kansas.
Figure 2. Hydraulic Conductivity (ft/day) map for study region in central Kansas.
Figure 3. Land use map for study region in central Kansas. City shown is Wichita.
Assigning Colors to Legends
In order to visualize the merged map effectively, colors were assigned to different land use types (Fig. 3). ArcView Legend Editor gives a variety of ways to classify and display feature attributes. Selecting the appropriate legend can choose a presentation style. The type of legend selected determines whether or not the features will be divided into classes and which feature attribute will be used for classification. Classification methods can be selected to determine how the attributes are assigned to classes and how many classes there will be. To easily understand the map, natural colors were assigned for each type by selecting the Unique Value from the Legend Type pull down menu. As an example, the color blue was assigned for water, green for grass and pink for urban areas. This made the look appealing and easy-to-understand as the theme containing the merged maps was a polygonal one and consequently, the assigned values filled the polygons with the appropriate color.
Land Use Pie Chart
Using the ArcView GIS, it is possible to create traditional presentation graphics as pie charts. Charts are a powerful way of presenting information to others. Figure 4 is a pie chart showing the distribution of Land Use areas. In order to create charts, a new Layout should be created first and then charts added into it. A pie chart for the Land-use merged view was created to graphically visualize how the area of the land was divided for different usage, (i.e. land use percentages). Data used in the pie chart were generated using the Summary feature in ArcView. This feature allows the user to summarize the data for a certain field in a certain Theme's Table and get statistics at any time about all records or selected records in a table. This feature was used to summarize the area field of the different records in the Land Use map. This chart was generated to demonstrate one of the many powerful capabilities in GIS.
Figure 4. Pie Chart showing Land Use Percentages for study region in central Kansas.
Digitizing the USGS Maps
At this stage, the other two maps, namely the Hydraulic Conductivity and the groundwater Depth maps are still useless for this work. This is due to being unable to edit and display them as well as applying ArcView features on them. In order to take full advantage of ArcView features, these maps needed to be digitized and imported into the GIS environment. Several steps were required to complete this task.
Maps shown in Figures 1 and 2 were divided into two-dimensional grid system. A coordinate system was chosen to represent each grid node by "x" and "y" values. These coordinates and grid points were drawn carefully as to have uniformity in both maps. Assigning values for each grid point in the hydraulic conductivity map was easy as the map itself was divided into a smaller grid system. Each grid point can have one of four different values. The first is for a permeability value of 15 feet/day. The second, third and fourth values are 150, 350 and 750, respectively. The north part of the map was recognized by its low values while the lower (southern) part was recognized by permeability values ranging between 350 and 750 feet/day. The work implemented for digitizing the groundwater contour map was more as the map was initially in the form of a contour map. Some grid points were not intersecting with the map contour lines and consequently, interpolation of the values was required. The generated data bank was stored in MS-Excel spreadsheet and prepared for editing. The total number of points at which values were recorded was 270.
Generating SURFER Maps
The easiest way to create digitized maps from the generated data bank is to use the SURFER program. Consequently, the data was imported into the SURFER program and two maps were generated. The first one was a "point" map that showed values of the permeability at different locations. This map was a replica of the initial permeability map. The second one was a generated contour map that showed contours for different water-depth values. This map accurately represented the original USGS contour map. This step will make the importing of these maps into ArcView as new themes possible. The only problem was in how to overlay the maps together as they have different units and how to locate a common point to overlay the map using its coordinates
Importing Maps to AutoCAD
To solve the problem of overlaying the maps, an engineering drawing software should be used as a medium between the SURFER and the ArcView software. The generated SURFER maps were imported into AutoCAD. The imported files were automatically transformed into entities and consequently lost every single assigned value to them, with the exception of their relative positions. Unfortunately, this made the process more complicated, but at the same time AutoCAD will allow to re-scale the maps as desired and to move them to new locations. From the Land Use map, a county boundary length was used to find the difference between the scales of the ArcView and AutoCAD maps. This difference was used to adjust the scale of the AutoCAD maps and have them in the same ArcView Land Use map scale. In order to position the AutoCAD maps in the correct location so that they overlay the ArcView Land Use map successfully, a boundaries-intersection point was chosen and spotted on all three maps for this reason. The AutoCAD drawings were moved so that the coordinates of this point are the same as the coordinates of that point in the Arcview display.
Importing Maps to ArcView
ArcView is capable of reading and displaying AutoCAD drawing files. The adjusted AutoCAD files were imported into ArcView as drawing files. At this stage of the work, there were two problems. First, these imported files are in drawing formats and editing them in not allowable. Second, the table of these added drawing themes contains no values at all. This means that all permeability and contour lines values were lost in the importing processes. To solve this problem, these drawing themes should be transferred into Shape files. This can be accomplished by selecting the "save as shape file" option from the ArcView menu bar.
Assigning Values for the RECORDS
All three maps are now displayed in ArcView. The hydraulic conductivity and the groundwater level maps are displayed, however, they are useless as they contain no information and consequently ArcView tasks can not be applied on them. Using the data bank available in Excel, for every node in the hydraulic conductivity map a value was assigned by editing the table of the theme. To edit a theme's table, the theme table should be active. From the Table menu, the Start Editing task should be selected. The field names become non-italic, confirming that the table is editable. New attribute fields were added as desired and values for each record were assigned. The work for the groundwater depth was harder as it is composed of contour lines instead of nodes. This is why the "Graphical Entity selection" was necessary for accomplishing this task. Also the "Promote" button made the job easier as it lists the selected entities one after another instead of having them scattered in the Theme's table.
Once again, the Unique Value color scheme was chosen for representing the legend of both the hydraulic conductivity and the groundwater depth maps. This allowed to graphically check (by eye inspection) whether the contamination potential results makes sense or not, as the two maps could be displayed and overlaid by each other easily.
Calculation of Travel Times
The calculation of travel times was made with the help of Excel using the data of hydraulic conductivity and aquifer depth as reference. The scheme of the problem was set up as two layers:
Risk Level Polygonal Themes
After calculating the leachate travel time at each node using the MS Excel database, risk levels of low, medium and high were assigned to those nodes. In order to display the risk levels as an ArcView theme, a new polygonal theme was created. To create a new empty polygonal theme and add shapes to it, the task "New Theme" should be selected from the View Menu. The New Theme dialog box displays allowing the user to select the feature type (polygon, line, point). After selecting the feature type, another dialog box allows the user to save the new shape file theme. Disc-shaped polygons were drawn on each grid node. To display the risk levels, a new polygonal theme was created (Fig. 5). The polygons were drawn to surround disk-shaped polygons that have the same risk level. Each drawn entity was given a certain value of low, medium and high. A color of red is for high-risk nodes. Yellow is for medium risk while green is for low risk nodes. This risk node-based map will be used later on for the purpose of selecting nodes within a certain distance from water and urban areas. These areas should not be chosen for placing the landfill.
Figure 5. Groundwater contamination risk polygon-based map for study region.
A color of red is for high-risk nodes. Yellow is for medium risk while green is for low risk nodes
BUILDING QUERIES
In the Land Use theme, Queries were implemented to filter out urban areas as well as surface water areas. After activating the merged maps theme, the Query Builder button is clicked (accessed from the theme properties dialog box). The name of the active theme is displayed at the top of the Query Builder. The dialog box contains a list of attribute fields and a set of operators, and a list of attribute values. When a field is clicked, all the unique values for that field display in the Values list, as long as the Update Values option is checked. As the query is being built, it is displayed in the query text box in the lower left corner of the dialog box. Once a query is entered in the text box, it is applied by selecting New Set button. New set creates a new set of selected features that match the query. This procedure was followed twice to create two different queries. The first one was used to display the urban areas in the four counties. After displaying the filtered urban areas, they were saved as a separate shape file. The same thing was applied for the second query for surface water area . Now the ArcView project contains two new themes, each one for a specific land use.
Finding the FEATURES NEARBY
The themes generated by using the query are used in the Select by Theme feature in ArcView. The purpose of this Select by Theme feature is to know which features are within a certain distance of other features or are adjacent to other features. This feature was used to locate select risk nodes that are within a certain distance from urban and water areas. The distance was selected to be 5 miles. All selected nodes was not considered as a reasonable place for a landfill. In order to carry out this task, Select By Theme feature should be selected from the Theme menu while having the risk nodes theme active. If the nodes are to be selected using the water areas theme, then all the records in the water areas theme should be selected. The same applies in the case of the urban areas theme. In the Select By Theme dialog box, the upper pull-down menu allows the user to select the "Are Within Distance of". The entry box allows the user to specify the distance (which is here 5 miles). This procedure selects some nodes in the active layer that will help marking lands unsuitable for landfill considerations.
New themes were generated to mark the lands in which nodes were selected from the Select By Theme feature. The two generated themes are assigned one color and no risk values were assigned for any of them, as risk level will not benefit the selection in this case. Instead, the two generated polygons are the generated restricted-areas polygon themes (Fig. 6 and 7).
Selecting Safe Landfill Sites
The remaining areas on the Land Use map satisfy the condition that they are at least 5 miles away from rivers, lakes and residential areas. But still there are some areas with high contamination risk. Having all themes activated at the same time (Figure 8), it is very easy to visualize the very few remaining areas and select the ones with low risk values. These selected areas meet the two conditions of having low contamination risk and safe distance from water and urban areas. The capabilities of Arcview and GIS-based maps made the task of selecting these areas very simple and easy. Figure 9 is a plot of the remaining available areas. To suggest final landfill locations, circles having a diameter of 20 miles were drawn in such a way that their center coincides with city centers (Fig. 10).
This will make the distance from the city center to the landfill site a convenient one. Areas circumscribed inside the circles are chosen as the proposed landfill locations.
Figure 6. Restricted areas around cities based on proximity.
Figure 7. Restricted areas based on proximity to surface waters.
Figure 8. All themes active. Red zone restricted area for landfill development.
Figure 9. Remaining land available for landfill development.
Figure 10. Suggested landfill sites based on groundwater contamination risk, proximity to surface waters and residential areas, and transportation costs.
CONCLUSIONS
The use of GIS for evaluation of future waste containment sites has shown to save time when there is need for fast evaluation. In this project we have developed maps assigning levels of contamination risk of aquifer for zones based on characteristics of the soil such as hydraulic conductivity and aquifer depth. The most time consuming task will likely be on gathering the information for hydraulic conductivity of representative soils for zones in study.
The main advantage found in the use of GIS is the globalization capability of the program as well as the easy access to the maps given by USGS which provides information about aquifers for all states in the United States. These maps may be used by local government agencies to anticipate suitable location for waste management systems and probable contaminant travel times. This also could be a very useful tool that could be used by private companies to recommend and bid waste containment structures with accuracy.
REFERENCES
Editor's note: The mention of commercial software ArcView and SURFER in this paper should not be construed as endorsement by EJGE.
![]() | |
© 2002 ejge |