(.[I\ NZOO ,(r.s ~~ 110.11% A Digital Model Approach to Water-Supply Management of the Claiborne, Clayton, and Providence Aquifers in Southwestern Georgia Lee L. Gorday Jeny A. Lineback Anna F. Long and William H. McLemore GEORGIA DEPARTMENT OF NATURAL RESOURCES Lonice C. Barrett, Commissioner ENVIRONMENTAL PROTECfiON DMSION Harold F. Reheis, Director GEORGIA GEOLOGIC SURVEY William H. McLemore, State Geologist Atlanta 1997 BULLETIN 118 A Digital Model Approach to Water-Supply Management of the Claiborne, Clayton, and Providence Aquifers in Southwestern Georgia Lee L. Gorday Jerry A. Lineback Anna F. Long and William H. McLemore GEORGIA DEPARTMENT OF NATURAL RESOURCES Lonice C. Barrett, Commissioner ENVIRONMENTAL PROTECTION DIVISION Harold F. Reheis, Director GEORGIA GEOLOGIC SURVEY William H. McLemore, State Geologist Atlanta 1997 BULLETIN 118 CONTENTS ABSTRACT ............................................................... v INTRODUCTION .......................................................... 1 General Background . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 Climate and Runoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 Previous Investigations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 Acknowledgements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2 Scope of Work ....................................................... 2 CONCEP'fUAL MODEL ..................................................... 3 DIGITAL MODEL ANALYSIS ................................................ 4 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4 Model Description . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 Relation of Digital Model to Conceptual Model . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 PREDICTIVE MANAGEMENT SCENARIOS .................................... 6 General ............................................................ 6 Baseline Simulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 Drought Simulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7 Industrial Development Simulation ........................................ 7 The City of Albany Floridan Aquifer Usage Simulation ......................... 8 INTERPRETATIONS AND CONCLUSIONS ..................................... 8 REFERENCES ............................................................ 10 APPENDIX A: Baseflow Analysis ............................................ A-1 SUPPLEMENT 1: STRATIGRAPIDC AND HYDROGEOLOGIC SETIING STRATIGRAPHY ......................................................... I-1 General ........................................................... I-1 Ripley Formation .................................................... I-1 Providence Sand ..................................................... I-1 Clayton Formation ................................................... I-1 Wilcox Group ....................................................... I-1 Claiborne Group ..................................................... I-2 Ocala Limestone ..................................................... I-2 HYDROGEOLOGY ........................................................ I-2 Aquifer Framework ................................................... I-2 Upper Floridan Aquifer ............................................... I-2 Claiborne Aquifer .................................................... I-3 Clayton Aquifer ..................................................... I-3 Providence Aquifer ................................................... I-4 HYDROLOGIC BUDGET ................................................... I-4 SUPPLEMENT II: DIGITAL MODEL PACKAGE BOUNDARY CONDITIONS ................................................ II-1 i HYDROLOGIC PROPERTIES .............................................. 11-2 Transmissivity ...................................................... II-2 Storage Factor ..................................................... II-2 Leakance ......................................................... II-3 FLUX CONDIDONS ...................................................... II-3 Recharge .......................................................... II-3 Rivers ............................................................ II-3 Pumpage .......................................................... II-4 CALIBRATION .......................................................... II-5 Introduction ....................................................... II-5 Possible Sources of Error in Observed Head Data ........................... 11-6 Model Results ...................................................... II-7 Steady-State Simulation ............................................... 11-7 Transient Simulation ................................................. 11-9 Simulated Hydrographs ............................................... 11-9 Fluxes ........................................................... II-10 Summary ......................................................... II-11 SENSITIVITY ANALYSIS ................................................. 11-11 ll.LUSTRATIONS AND TABLES TEXT FIGURES 1. Study area . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14 2. Schematic diagram of regional ground-water flow in the Upper Floridan (A1), Claiborne (A2), Clayton (A3), and Providence (A4) aquifers ................................. 15 3. Local, intermediate, and regional ground-water flow . . . . . . . . . . . . . . . . . . . . . . . . . . 15 4. Stratigraphic column, hydrogeologic units, and model layers .................... 16 5. Predevelopment potentiometric surface of the Claiborne aquifer ................. 17 6. Predevelopment potentiometric surface of the Clayton aquifer . . . . . . . . . . . . . . . . . . 18 7. Finite-difference grid and boundary conditions .............................. 19 8. Simulated potentiometric surface in model layer A2 at the end of the five-year baseline simulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 9. Simulated potentiometric surface in model layer A3 at the end of the five-year baseline simulation ................................................................ 21 10. Simulated potentiometric surface in model layer A4 at the end of the five-year baseline simulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22 11. Change in simulated head in model layer A2 in response to a hypothetical two year drought .................................................................. 23 12. Change in simulated head in model layer A3 in response to a hypothetical two year drought ................................................................. 24 13. Change in simulated head in model layer A4 in response to a hypothetical two year drought ................................................................. 25 14. Change in simulated head in model layer A2 due to a simulated withdrawal of one MgaVd from both layers A2 and A3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26 15. Change in simulated head in model layer A3 due to a simulated withdrawal of one MgaVd from both layers A2 and A3 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 16. Change in simulated head in model layer A4 due to a simulated withdrawal of one MgaVd from both layers A2 and A3 ............................................. 28 17. Change in simulated head in model layer A2 due to a 20 percent reduction in pumpage from the Claiborne, Clayton, and Providence aquifers for the City of Albany . . . . . . . . . . . . . . 29 18. Change in simulated head in model layer A3 due to a 20 percent reduction in pumpage from the Claiborne, Clayton, and Providence aquifers for the City of Albany . . . . . . . . . . . . ... 30 19. Change in simulated head in model layer A4 due to a 20 percent reduction in pumpage ii from the Claiborne, Clayton, and Providence aquifers for the City of Albany ............. 31 TEXT TABLE 1. Stress periods and pumpage utilized in calibration of transient model . . . . . . . . . . . . . 13 SUPPLEMENT I FIGURES 1-1. Streamflow measurement, location, and general aquifer outcrop areas ............ 1-6 SUPPLEMENT II FIGURES II-1. Boundary for model layer A1 ......................................... 11-13 II-2. Outcrop, active, and specified-head cells in model layer A2 ................... 11-14 II-3. Outcrop, active, and specified-head cells in model layer A3 ................... 11-15 II-4. Outcrop, active, and specified-head cells in model layer A4 ................... 11-16 II-5. Transmissivity for model layer A2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . II-17 II-6. Transmissivity for model layer A3 ...................................... II-18 II-7. Transmissivity for model layer A4 ...................................... 11-19 11-8. Calibrated storage for model layer A2 ................................... 11-20 II-9. Calibrated storage for model layer A3 ................................... 11-21 II-10. Calibrated storage for model layer A4 ................................... 11-22 II-11. Calibrated leakance for model layer C2 ................................. 11-23 11-12. Calibrated leakance for model layer C3 ................................. '11-24 II-13. Simulated irrigation pumpage ......................................... 11-25 11-14. Simulated pumpage for all uses ........................................ 11-26 11-15. Simulated pumpage from model layer A2 in stress period 15 .................. 11-27 II-16. Simulated pumpage from model layer A3 in stress period 15 .................. 11-28 11-17. Simulated pumpage from model layer A4 in stress period 15 .................. 11-29 11-18. Steady-state simulated potentiometric surface for model layer A2 .............. 11-30 II-19. Steady-state simulated potentiometric surface for model layer A3 .............. 11-31 11-20. Steady-state simulated potentiometric surface for model layer A4 .............. 11-32 II-21. Estimated and simulated fluxes for model layer A2. model layer A3, and model layer A4 ............................................................... 11-33 II-22. Steady-state mode boundary fluxes for recharge, discharge, and leakage ......... 11-34 11-23. Simulated potentiometric surface for model layer A2 at the end of stress period 15 11-35 11-24. 1986 potentiometric surface for the Claiborne aquifer (model layer A2) ......... 11-36 11-25. Simulated potentiometric surface for model layer A3 at the end of stress period 15 11-37 II-26. 1986 potentiometric surface for the Clayton aquifer (model layer A3) ........... 11-38 II-27. Simulated potentiometric surface for model layer A4 at the end of stress period 15 11-39 11-28. 1986 potentiometric surface for the Providence aquifer (model layer A4) ........ 11-40 11-29. Simulated and observed heads at (a) USGS test well 4, Claiborne aquifer, Dougherty County and (b) USGS test well 2, Claiborne aquifer, Dougherty County ...... 11-41 11-30. Simulated and observed heads at (a) W. H. Fryer well, Claiborne aquifer, Lee County and (b) Cuthbert well, Clayton aquifer, Randolph County . . . . . . . . . . . . . . . . 11-42 II-31. Simulated and observed heads at (a) Don Foster well, Clayton aquifer, Terrell County and (b) H. T. McLendon #1 well, Clayton aquifer, Calhoun County ...... 11-43 11-32. Simulated and observed heads at (a) USGS test well 9, Clayton aquifer, Lee County and (b) Turner City well 2, Clayton aquifer, Dougherty County ............ 11-44 II-33. (a) Simulated and observed heads at USGS test well12, Clayton aquifer, Dougherty County and (b) simulated rate of release of water from storage through 15 stress periods ............................................................ 11-45 II-34. (a) Simulated ground-water discharge to rivers for selected stress periods and (b) net vertical flux across each confining unit ................................... 11-46 11-35. Simulated horizontal fluxes in (X) and out (0) for each aquifer and stress period. . 11-47 iii SUPPLEMENT II TABLES 11-1. Root mean square error and number of observations for calibrated model in steady-state model ........................................................ 11-12 11-2. Simulated and estimated ground-water discharge to rivers .................... 11-12 11-3. Root mean square error and number of observations for stress periods used in transient cahbration ....................................................... 11-12 iv ABSTRACT Management of ground-water resources can be effectively performed by state agencies responsible for issuing water-use permits that have operating numeric ground-water flow models for major aquifers. Proposed changes in pumpage can be simulated by such models to a sufficient degree of accuracy so as to provide the technical basis for decision making by management agencies. The Geologic Survey Branch of the Georgia Environmental Protection Division (EPD), with assistance from the United States Geological Survey (USGS), has developed a numeric ground-water flow model for the Claiborne, Clayton and Providence aquifers in southwestern Georgia to assist the Division with its water resource management responsibilities. The model utilizes the McDonald-Harbaugh modular quasi-three dimensional ground-water flow model code (MODFLOW). This code is capable of both steady-state and transient simulations. The model is calibrated to observed head and pumpage data through 1986. 1986 head and pumpage data are essentially unchanged through 1992. The model will be used by EPD to simulate potential changes in heads before issuing permits for any requested additional ground-water use in the affected aquifers. The model is used to predict heads in the aquifers for four hypothetical ground-water scenarios, developed in consultation with the Director of EPD. If present pumpage (the Baseline Simulation) continues unchanged, the heads in the Clayton and Providence aquifers will continue to decline. Simulated response to a major two-year drought (the Drought Simulation) indicates that heads in all three aquifers would decline during the duration of the drought, particularly in areas of heavy irrigation pumpage, and that the recharge to the aquifers would be replaced by the withdrawal of large volumes of ground-water from storage in the aquifers. The long term effect of withdrawing one million gallons per day for industrial use from the Claiborne and from the Clayton aquifers is simulated as the third scenario (the Industrial Development Simulation). Such additional pumpage would likely lower the heads in the Claiborne, Clayton, and Providence aquifers significantly near the industrial well(s). Lesser effects would be widely distnbuted in the aquifers. The Clayton aquifer, already the most stressed aquifer in the system, would be even more adversely affected. The fourth situation involves the City of Albany, the largest municipal water user in the region. In this scenario, pumpage from the Claiborne, Clayton, and Providence aquifers in the general vicinity of Albany is reduced by 20 percent, in response to a plan to utilize the Upper Floridan aquifer to meet some of the City's future water needs. Such a reduction in pumpage would cause the heads in the Claiborne, Clayton, and Providence aquifers to rise in the Albany area. The Upper Floridan, Claiborne, Clayton, and Providence aquifers in southwestern Georgia and southeastern Alabama are hydrologically interconnected, albeit not efficiently. Development of irrigation, industrial, and municipal wells in the region has resulted in declines in the potentiometric surfaces of the Claiborne, Clayton, and Providence aquifers when compared to predevelopment surfaces. Vertical leakage between aquifers is an important pathway of ground-water flow. The Claiborne aquifer, in particular, is well connected hydrologically with the overlying Floridan aquifer; and, as a result of this interconnection, additional ground-water resources may be developed at some locations. The Clayton aquifer, on the other hand, is poorly connected to the other aquifers, is heavily pumped in the region, and is not capable of supporting any significant additional pumpage in most parts of the region. The Clayton is somewhat better connected to the underlying Providence aquifer than it is to the overlying Claiborne aquifer; therefore, attempted additional development of the Clayton also would probably adversely affect the Providence and could adversely affect all three aquifers in places. The Claiborne, Clayton, and Providence aquifers are particularly heavily stressed in the Albany and Dawson areas where additional pumpage from these three aquifers is not recommended. v INTRODUCTION General Background Ground-water use in southwestern Georgia has increased significantly over the past few decades. Since the middle 1970's, combined increases in industrial, agricultural and municipal ground-water use have caused declines in ground-water levels in the deeper confined aquifers. The continued development of ground-water resources and the associated water-level declines, in conjunction with recent droughts, vividly demonstrate the need for proper and efficient management ofthese valuable resources. This report descnbes the development and application of a digital model of ground-water flow for the three major confined aquifers of southwestern Georgia: the Claiborne, Clayton, and Providence aquifers. Although not the focus of the study, the Upper Floridan is included in the model. The general purpose of this digital model, therefore, is to improve understanding of ground-water flow in the confined aquifers and to facilitate efficient and informed resource management. Specific elements of this study included: (1) development and calibration of a digital model that simulates ground-water flow in three vertically contiguous aquifers, and (2) application of this model to quantitatively descnbe pre-development, modem (1986), and hypothetical future ground-water flow conditions. Figure 1 shows the study area in southwest Georgia and southeast Alabama. It covers approximately 7,425 square miles and includes all or parts of 30 counties in the two states. The study area was defined by the distnbution and extent of the aquifer system of interest. The southwestern boundary of the study area is approximately defined by the western drainage divide of the Chattahoochee River in eastern Alabama. The southeastern boundary of the study area is approximately defined by the southern extent of the aquifer system. The northeastern boundary of the study area is the approximate drainage divide between the Flint and Ocmulgee Rivers in southwestern Georgia. The northwestern boundary of the study area is approximated by the inner margin of the Coastal Plain. The study area is drained by two of Georgia's largest rivers. The Chattahoochee River separates Alabama and Georgia and drains the western part of the study area. The Flint River drains the eastern part of the study area. Near the Georgia-Florida State line, the Flint River joins the Chattahoochee River to form the Appalachicola River, which flows southward across the Florida panhandle to the Gulf of Mexico. The Chattahoochee River has been dammed at Fort Gaines and forms the Walter F. George Reservoir. The Flint River has been dammed at the juncture of Crisp, Lee, and Worth Counties forming Lake Blackshear and at Albany forming a Georgia Power Company reservoir. Jim Woodruff Dam, creating Lake Seminole, is at the confluence of the Chattahoochee and Flint Rivers. Climate and Runoff Southwest Georgia and southeast Alabama generally have short, mild winters and hot, humid summers. Winter temperatures generally are above freezing, but do occasionally drop below 20 degrees F. Summer temperatures commonly are above 90 degrees F and temperatures above 100 degrees F are not rare. Precipitation in the study area occurs almost completely as rainfall, and ranges from about 46 to 54 inches per year (Carter and Stiles, 1983). Abundant rain occurs during winter months with a gradual increase to a maximum in February or March. Heaviest rains normally fall during July and August due to frequent summer thunderstorms. October and November are the driest months. A secondary period of diminished rainfall is in April and May. Annual runoff in southwestern Georgia is highly variable, and ranges from 12 to 24 inches per year. Previous Investigations The geology and hydrogeology of southwestern Georgia have been previously studied in either a localized or a general fashion. Stephenson and Veatch (1915) descnbed the geology of the Georgia Coastal Plain by formation, including the areal extent, lithology, stratigraphic position, strike and dip of beds, thickness of rock units, paleontology, and structure. Cooke (1943) descnbed the general geology of the Georgia Coastal Plain. Herrick (1961) advanced the knowledge of the geologic framework of the Coastal Plain of Georgia by descnbing detailed lithologic logs. Marsalis and Friddell (1975) gave an overview of the lithologic units exposed in the Chattahoochee River Valley area, including discussions of facies changes along strike and down dip. Gibson (1982) differentiated six Paleocene and Middle Eocene marine units in eastern Alabama and western Georgia, including composition, fossil assemblage, and descriptions of nonmarine and marine transitions. Numerous investigations of the ground-water hydrology were undertaken beginning in the early 1950's as the demand for ground-water increased. 1 As in most ground-water studies, investigation of the geology was commonly a substantial part of these efforts. In 1958, Wait described the stratigraphy and ground-water availability in Crisp County. Wait (1960 a, b, c,) also descnbed the geology and ground-water resources in Calhoun, Clay, and Terrell Counties, and discussed the geology and ground-water resources of Dougherty County (1963). Owen (1963) compiled existing data to extend the knowledge of ground-water conditions in Lee and Sumter Counties. Stewart (1973) reported Clayton aquifer hydraulic characteristics which were estimated from aquifer tests performed during the design and construction of the Walter F. George Lock and Dam in the Ft. Gaines area. Pollard and Vorhis (1980) descnbed the geohydrology of the Cretaceous aquifer system in Georgia. Hicks and others (1981) evaluated the development of ground-water resources in the Albany area. Ripy and others (1981) published an interim report on the hydrogeology of the Clayton and Oaiborne aquifers in southwest Georgia. McFadden and Perriello (1983) conducted a general study of water-level trends, ground-water quality, ground-water use, aquifer geometry, lithologic and hydrologic characteristics, and recharge and discharge mechanisms of the Clayton and Claiborne aquifers in southwest Georgia. Clarke and others (1983, 1984) descnbed and evaluated the effects of water use on the ground-water systems of the Providence and Clayton aquifers, respectively. In 1984, the Geologic Survey compiled an atlas (Arora, editor) descnbing aquifers in the Georgia Coastal Plain, including isopach, structure-contour, and potentiometric surface maps as well as cross-sections. Davis (1987) descnbed the stratigraphic and hydrogeologic framework of the Cretaceous, Tertiary, and Quaternary Systems in Alabama to aid in delineating aquifers and confining units within the Alabama Coastal Plain. Water-level, water-use, and water-quality information on the Clayton and Claiborne aquifers between 1982 and 1986 has been compiled (Long, 1989a). Ground-water flow and stream-aquifer relations in the outcrop areas of the Coastal Plain sediments were quantitatively descnbed by Faye and Mayer (1990). Hayes and others (1983) developed a digital finite-difference ground-water flow model of the Principal Artesian (Floridan) aquifer in the Dougherty Plain area of southwest Georgia. The United States Geological Survey (USGS), as part of their Regional Aquifer System Analysis (RASA) program, has two ground-water flow modeling studies pertinent to this study. Maslia and Hayes (1988) defined the regional flow system of the Floridan aquifer system in the Dougherty Plain. Faye and Mayer (1990) developed a ground-water flow model of regional ground-water flow in the Coastal Plain of Georgia. Although Faye and Mayer's study encompassed the entire hydrogeologic framework of the Coastal Plain, it primarily addressed the deep, regional component of the ground-water flow system, and was necessarily of rather large scale. A description of the stratigraphic and hydrogeologic framework of southwestern Georgia resulting from the previous studies mentioned in this section, and on which the conceptual and digital (or numeric) models presented in this report are based, are included in this report as Supplement I. Acknowledgments The development of a digital model for the Claiborne, Clayton, and Providence aquifer systems was begun by Anna Long of the Georgia Geologic Survey with the preparation ofHydrologicAtlas 19 (1989a). Long (1989b) proceeded with the development of a steady-state digital model. Lee Gorday expanded upon Long'swork and developed and calibrated a transient model and used the model to prepare predictive scenarios for various water-supplymanagementoptionsforsouthwestem Georgia (this report). Model development was carried out by the Georgia Geologic Survey to assist the Water Resources Management Branch of EPD in its assigned role of regional ground-water use management. The model was developed under the guidance ofthe USGS, Georgia District Office. The USGS provided technical assistance and guidance on a day-to-day basis through model calibration and the manuscript. Scope ofWork This report provides the results of an application of the McDonald-Harbaugh (1988) ground-water flow model (MODFLOW) of the Claiborne, Clayton, and Providence aquifers to four ground-water management simulations. For a description of the hydrogeological assumptions, mathematical assumptions, and computer codes, the reader is referred to their work. No field studies were performed for this study. Data entered into the model are direct measurements reported by others, estimates reported by others, calculated values made by the authors or by others, or estimates believed by the authors, after consultation with the USGS, to be reasonable. All modelling was done at the offices of the 2 Georgia District of the USGS using their computer facilities onto which MODFLOW had been installed. CONCEPTUAL MODEL Existing field data, previously published information on the aquifers, and theoretical concepts of ground-water flow were synthesized to develop a conceptual model of ground-water flow within the interconnected Floridan, Claiborne, Clayton, and Providence aquifers. The conceptual model addresses ground-water flow in the aquifer system from the point of initial recharge to the system to the point of ultimate discharge. The conceptual flow model developed for the present study closely follows the conceptual model presented by Faye and Mayer (1990) as a part of their digital model analysis that included additional aquifers and larger study area. The present study focuses on both the regional and intermediate flow systems, whereas the Faye and Mayer study focused strictly on the regional flow system. Figure 2 is a schematic representation of the aquifer system under investigation. The basic premise behind the conceptual model is that precipitation recharges the Claiborne, Clayton, and Providence aquifers in their outcrop areas in the upper Coastal Plain. Ground water then flows to the south-southeast down-gradient through the aquifers, which in tum become confined. Because the older aquifers crop out at higher elevations, the head in older aquifers is generally higher than the head in younger aquifers. This means that under predevelopment conditions there was an upward component of ground-water flow, across the confming (lower hydraulic conductivity) units (e.g. there was flow from the Providence to the Clayton, from the Clayton to the Claiborne, and from the Claiborne to the Floridan-lowering of head as a result of pumpage, however, could reverse such gradients). Ground-water flow in an aquifer recharge area is dynamic and complex and is controlled largely by topography and stream-aquifer relations. Toth (1963), in descnbing the flow of ground water in an unconfined area with local relief, introduced the concept of local, intermediate, and regional flow systems (Figure 3). The aquifers are unconfined and ground-water flow within each aquifer has significant vertical as well as horizontal components. The hilly topography of the aquifer recharge areas such as occurs in the western part of the upper Coastal Plain of Georgia, produces numerous subsystems within the major flow system. Most of the water that recharges the ground-water system flows to the closest stream and discharges; this is termed local flow. Local flow is characterized by short, shallow flow paths. Water that enters the system at the highest point and discharges to the stream or river at the lowest point in the area under investigation (or that flows downdip beneath younger semi-confining units) is known as regional flow. Regional ground-water flow follows the longest and deepest flow paths. Between these extremes is the intermediate flow system. Water in the intermediate flow system bypasses at least one local discharge site along its flow path (Toth, 1963). These terms are dependent to a large degree on scale, (e.g. Toth's concept of regional flow in unconfined outcrop area is not the same as regional flow throughout an entire aquifer system). Fluctuations in climatic conditions, such as droughts, affect the local flow system with its short, shallow flow paths, to a greater degree than the intermediate and regional flow systems. Because the purpose of this study was to develop a model for use in managing the ground-water resource (e.g. ground-water withdrawal permits), the conceptual model as well as the digital model focus on the intermediate and regional ground-water flow systems. The local flow system is beyond the scope of the management objective. For the purpose of this study, it was assumed that the drought of 1954, one of the severest on record, depleted flow in the local flow system. Precipitation that falls in the area of outcrop of the aquifers may run off or infiltrate the ground surface. A small amount may be held in puddles or ponds where it may evaporate. Much of the water that infiltrates the ground is transpired back to the atmosphere by plants. Water that percolates to the water table recharges the ground-water system. This is the ultimate source of water to the four aquifers under consideration in this study. The assumption is made that all streams that cross the outcrop areas either gain water from the aquifers by way of baseflow, or have no significant net loss of water to the aquifers. Although this assumption may not hold for the local flow system, it is probably valid for the intermediate and regional flow systems that are the focus of this study. Recharge to the aquifers is greatest in the interstream divides and lowest in the vicinity of streams. Water that recharges the aquifer may follow a number of flow paths. Most of the water that recharges the ground-water system is discharged to rivers and streams, as described above. Water that does not discharge to streams flows down the dip of the aquifer (in which it recharged) or moves 3 vertically to another aquifer as leakage. Vertical leakage between aquifers, prior to development, is conceptualized as being downward in areas of net recharge and upward in areas of net discharge. Between these areas, ground-water flow is essentially horizontal with little vertical movement. In the updip part of the study area, most vertical leakage is downward, the exception being in the vicinity of streams receiving discharge. This is corroborated by comparing the map of the potentiometric surface of the Claiborne aquifer (Figure 4) (McFadden and Perriello, 1983) with the potentiometric surface map of the Clayton aquifer (Figure 5) (McFadden and Perriello, 1983) of the same time period. Comparison of these maps indicates that the gradient between the Claiborne and Clayton aquifers generally is downward in the updip area of the study and generally is upward in the downdip area of study. Data are lacking to descnbe the head gradient in the entire downdip portions of the study area (e.g. south of Miller and Mitchell Counties). Head data that are solely from the Providence aquifer are sparse, which precludes a comparison between potentiometric surfaces of the Providence and Clayton aquifers. Because of the high hydraulic conductivity of the Providence, Clayton, and Claiborne aquifers relative to the confining units, flow is assumed to be chiefly horizontal in the aquifers and vertical through the confining units. The Flint and Chattahoochee Rivers influence ground-water flow in the downdip portion of the Claiborne and Clayton aquifers, as they do in the updip area. In Georgia, an approximate north-south ground-water flow divide has developed where ground water on one side flows west toward the Chattahoochee River; but on the other side, ground water flows east toward the Flint River. Because heads are lower on either side of this divide, there probably is no flow of water across the divide. The ground-water flow divide in the Clayton aquifer can be seen in Figure 5 extending from eastern Randolph County southward through central Calhoun County, to Miller County. Similar divides appear to exist between the Flint and the Ocmulgee Rivers and between the Chattahoochee and the Chochtawhatchee Rivers. The location of the ground-water flow divides can change with time in response to changes in hydrologic conditions such as recharge and pumpage. The eastern boundaries of the area of this study generally conform to the divide east of the Flint River for the Claiborne, Clayton, and Providence aquifers prior to development oflarge ground-water pumpage. The western boundary of the study area, in Alabama corresponds to the apparent divide west of the Chattahoochee River for the Providence and Clayton aquifers prior to development of large ground-water pumpage. Downdip facies changes that result in the decrease in the hydraulic conductivity of the Claiborne, Clayton, and Providence aquifers all occur in the same general area (see Supplement I). Although some water may continue to flow down dip, it is unlikely that much water flows across this boundary. Prior to the installation of large numbers of wells in the aquifers, the water flowing laterally in these aquifers would move downdip, then vertically upward (Faye and Mayer, 1990). After the water moved upward into the Floridan aquifer, it would be discharged to a stream or would flow laterally within the Floridan aquifer. The development of large quantities of pumpage (e.g. primarily irrigation) from the aquifers has led to a reduction in the upward flux in some down-dip areas. Locally, vertical gradients have been reversed. The increase in ground-water withdrawals probably has resulted in the lateral shifting of some of the ground water divides. DIGITAL MODEL ANALYSIS Introduction Water-level declines resulting from increasing ground-water withdrawal and a series of droughts have demonstrated a need for a better understanding of ground-water flow in the Claiborne, Clayton, and Providence aquifers in southwest Georgia. In response to this need, and in anticipation of future conflicting demands for the finite water resources of the region, a digital ground-water flow model was constructed to aid in the informed management of this vital resource. Details of model development and calibration and the results of sensitivity analysis are presented in this report as Supplement II. The primary focus of the digital model simulation is flow in the Claiborne and Clayton aquifers in the area of their greatest use. In order to adequately address these issues, the model had to include the Upper Floridan aquifer as a constant-head source/sink above the Claiborne aquifer, and the Providence aquifer as an active layer below the Clayton aquifer. In addition, the lateral boundaries of the model were extended, in most areas, well beyond the active use of the aquifers in Georgia to have stable conditions at the boundaries. A digital ground-water flow model computes the potentiometric head in an aquifer over space and, for transient simulations, time. Heads are 4 computed by solving ground-water flow equations, given the distribution of hydraulic parameters, boundary conditions, and initial conditions. Analytical solutions to these equations are available for a range of relatively simple boundary conditions. Digital models are useful in situations where the boundary and initial conditions are complex, and where hydraulic properties or characteristics vary through space. A model of ground-water flow conditions prior to development was constructed and calibrated against measured heads, estimated ground-water discharge to streams, and estimated flux at the boundaries of the model (Long, 1989b). The model was used in the steady-state mode to simulate conditions prior to the development of the aquifers having high-yielding wells. In the steady-state mode, conditions within the flow system do not change with time. The head distnbution from steady-state simulations of the model were used then as the initial head distnbution in transient simulations, where pumpage changed through time. Model parameters were adjusted to provide a closer match between the heads and fluxes simulated by the model in the steady-state mode and observed heads, estimated discharge to streams and estimated boundary fluxes. The simulated heads and model parameters (including boundary conditions) from the steady-state simulation were used along with additional parameters (storage and pumpage) in the transient simulation. Transient simulations were constructed for fifteen stress periods between 1900 and 1986 (Table 1). (Note: Heads were ftxed in the Upper Floridan aquifer (A1); this was deemed reasonable as the Upper Floridan aquifer in southwest Georgia is recharged every year and water level declines are only short term (ie. a few months).) The model heads computed in the transient simulations were compared to observed heads measured at seven times during the period of simulation. Observed heads sufficient for comparison were available for all aquifers at times corresponding to period 2 (1945-1959), period 8 (1978 and 1979), period 13 (1984) and period 15 (1986). Observed heads for the Providence aquifer were available for period 9 (1980) and for the Claiborne and Clayton aquifers for periods 10 and 11 (1981 and 1982), respectively. Details of the development, calibration, and sensitivity analysis of the transient model are in Supplement II of this report. Hydrographs of simulated head were prepared for various stress periods and compared to observed water-level fluctuations. Changes in any parameters, other than storage or pumpage, required the change to be implemented in the steady-state mode, and a complete evaluation was made of the match between simulated conditions and observed and estimated conditions prior to implementation in the transient mode. This loop approach was used until an acceptable match was achieved between simulated conditions and observed and estimated conditions for the model in both the steady-state and transient modes. Model Description The McDonald-Harbaughmodularquasi-three dimensional ground-water flow model code (MODFLOW), which is capable of both steady-state and transient simulations, was used in this study (McDonald and Harbaugh, 1988). This model is based on a fmite difference approach, which uses a rectangular grid, either uniform or variably-spaced. Within each cell of the grid, the hydraulic parameters are uniform. Fluxes simulated in each cell (such as recharge, well pumpage, and discharge to rivers) are distnbuted evenly over the cell. Potentiometric heads calculated by the model represent the head over the entire area of the cell. Three-dimensional flow is simulated by linking two-dimensional (lateral) flow in each layer with one-dimensional (vertical) flow between the layers as a representation of leakage. The model code constructs an equation describing ground-water flow for each node. These simultaneous equations are solved through an iterative process of matrix algebra known as the strongly implicit procedure (SIP). The procedure is descnbed in detail in the model documentation (McDonald and Harbaugh, 1988). Wang and Anderson (1982) present an overview of various solution techniques including SIP. Relation of Digital Model to Conceptual Model A finite-difference grid of 57 rows and 85 columns was used to subdivide the study area (Figure 6). The grid was oriented so that the predominant direction of ground-water flow would be parallel to the columns of the grid. The columns are aligned 30 degrees west of north. Grid spacing was either one or two miles for both rows and columns. Cell areas are 4 square miles for cells with both row and column spacings of 2 miles, 1 square mile for cells with both row and column spacings of 1 mile, and 2 square mile for 1x2 mile cells. The model area is bounded by either specified head and no-flow conditions both laterally and vertically. The specific application of these boundary conditions is discussed for each 5 layer in Supplement II of this report. The digital flow model is based upon the conceptual flow model described previously. The general inter-relationship between aquifers and semi-confining units and regional stratigraphy is shown in Figure 7. The Upper Floridan aquifer is represented by model layer A1, and is treated as a source-sink. Model layer A2 represents the Claiborne aquifer. The Clayton and Providence aquifers are represented by model layers A3 and A4, respectively. Flow across the confining units is simulated by one-dimensional flow based on the simulated heads in the adjacent layers and the leakance (vertical hydraulic conductivity divided by confining unit thickness). The confining unit between the Floridan and Claiborne aquifers is represented by C1 leakance. The Claiborne-Clayton and Clayton-Providence confining units are represented by C2leakance and C3 leakance, respectively. The digital flow model in this report is designed to simulate only the intermediate and regional ground-water flow systems. Simulation of local variation in head is not possible using the cell size in this model. Additionally, by addressing only the regional and intermediate flow systems, the independent estimate of flux to rivers and streams can be used to aid in the calibration of the model. Recharge as used in the digital model differs from the common concept of recharge. Only recharge to the intermediate and regional flow systems is considered in the model. Total recharge, therefore, is considered an upper limit to the recharge used in the digital model. Streams and rivers, which had significant flow during the 1954 drought, are considered in the digital model to estimate aquifer contnbutions in base flow. Streams that are simulated are identified in Appendix A. PREDICTIVE MANAGEMENT SCENARIOS General One of the purposes in developing the digital model of the Claiborne, Clayton, and Providence aquifers is to develop a tool for use in managing the ground-water resource. Four predictive scenarios or simulations were developed to demonstrate the usefulness of the model and to assess the response of the ground-water system to several hypothetical changes in pumpage. The four scenarios are: (1) The Baseline Simulation, (2) The Drought Simulation, (3) The Industrial Development Simulation, (4) The City of Albany Floridan Usage Simulation. For these four simulations, model parameters, including boundary conditions, were unchanged from transient simulations. The results of a simulation of a hypothetical change to the ground-water system are commonly assessed by examination of the trend and magnitude of changes in simulated head. In a system as dynamic as the Claiborne, Clayton, and Providence, changes in flux within the system and at the boundaries of the digital model also are of importance. Baseline Simulation The first predictive scenario simulates the response of continuing the model for five years at the pumpage rates similar used in the calibration of the model. This simulation represents a baseline or status-quo condition. Municipal and industrial pumpage are continued at the rates used in period 15 (1986) of the calibration simulations (see Table 1 and Supplement II). Irrigation pumpage is assigned a value that approximates the average pumpage for irrigation used in periods 10 through 15. The starting heads for the predictive simulations are the heads at the end of the calibration simulation. In the Baseline Simulation, simulated heads generally decline in areas having large withdrawals for municipal and industrial uses during the first year. Simulated heads rise in areas with large withdrawals for irrigation. This rise in heads is the result of the average irrigation withdrawal being smaller than the irrigation withdrawal for period 15, the last stress period of the calibration period. The changes in head are relatively small for layers A2 and A4. In layer A2, the simulated head ranges from six feet higher to three feet lower at the end of the first year of the Baseline Simulation compared to heads at the end of the cahbration period. Simulated heads decline across all of layer A4, with the maximum decline being less than six feet. Withdrawal for irrigation is much smaller in layer A4 than in layers A2 and A3; therefore, the reduced irrigation withdrawal in layer A4 has little impact on simulated heads. Simulated heads in layer A3 decline slightly over much of the area of the model. The maximum decline is less than six feet. A number of isolated areas are identified where simulated heads rose in response to the decreased withdrawals for irrigation. Where a number of irrigation systems are within a single cell, the rise in simulated head is quite large. The rise is as great as 23 feet in one cell. The area over which the simulated head rose was quite restricted in comparison to layer A2. Simulated heads at the end of the fifth year of the Baseline Simulation are similar to those at the 6 end of the first year of the simulation. Large areas in layer A2 have simulated heads that are higher at the end of the fifth year than at the end of the calibration period; however, the area is smaller than that at the end of the first year of the simulation. The maximum decline in simulated head in layer A2 is about seven feet. Figure 8 shows the simulated potentiometric surface in layer A2 at the end of the ftve-year Baseline Simulation. Simulated heads in layer A3 at the end of the fifth year of the simulation are significantly lower than at the end of the first year. Drawdown from the end of the calibration period is as great as 16 feet. Some areas of layer A3 continue to have simulated heads that are higher than at the end of the calibration period; but these are less numerous and of lower magnitude than at the end of the ftrst year of the simulation. The simulated potentiometric surface in layer A3 at the end of the Baseline Simulation is shown in Figure 9. Simulated heads in layer A4 are as much as 16 feet lower at the end ofthe simulation (Figure 10) than at the end of the first year or at the end of the calibration period. Simulated fluxes within the ground-water flow system and at the model boundaries had changed relatively little at the end of the Baseline Simulation from the simulated fluxes at the end of the calibration period. Net vertical flux across confining unit C3 increases from 15.9 cubic feet per second (cfs) at the end of the calibration period to 18.4 cfs at the end of the Baseline Simulation; this is within the range of values simulated at other stress periods in the calibration simulation. The simulated net vertical flux across confining units C1 and C2 increases to 6.7 and 9.7 cfs, respectively. These values are higher than the simulated flux at any stress period in the calibration simulation. Changes in horizontal constant-head fluxes and discharge to rivers and streams are minimal. The Baseline Simulation was developed not only to assess the effects of continuing pumpage at current values, but also to provide a basis for comparison with the other predictive scenarios. By comparing the results of the other scenarios to the Baseline Simulation, the effects of the modeled change can be isolated from the effect of continued pumpage at existing levels. Drought Simulation Drought conditions were simulated using reduced recharge rates and increased pumpage for irrigation. The recharge rate was reduced to 75 percent of the calibrated value, a reduction of 142 cfs. This reduction is similar in magnitude to the reduction in precipitation in a drought having a 10year recurrence interval. The pumpage for irrigation was arbitrarily assigned a value of 150 percent of the irrigation used in the Baseline Simulation. The pumpage in layer A2 represents an increase in pumpage of 25.5 cfs over the pumpage in the Baseline Simulation. Increases for layers A3 and A4 are 12.1 and 0.7 cfs, respectively. The change in simulated head in layer A2 at the end of the second year of the Drought Simulation (compared to simulated heads at the same time in the Baseline Simulation is shown in Figure 11. There are large areas having drawdowns of 5 feet or more. The maximum difference between the simulated head for the Drought Simulation and the Baseline Simulation is 24 feet. Large drawdowns are restricted to small areas. Head differences between the drought and baseline simulations in layer A3 are shown in Figure 12. Drawdowns in layer A3 are much greater than in layer A2 despite the fact that the amount of additional withdrawal compared to the Baseline Simulation is much smaller. The area having a simulated drawdown of 5 feet or more is significantly larger for layer A3. Drawdowns in layer A4, shown in Figure 13, cover a large area, similar to layer A3. The maximum drawdown is 18 feet, but this is produced by a very small additional withdrawal. Much of the difference in layer A4 is due to leakage of water to layer A3. The largest change in simulated flux compared to the Baseline Simulation is the movement of water from storage. The increase in water being removed from storage is 150 cfs at the end of two years of drought. The large change is due to the reduction in recharge. Upward vertical flow from layer A4 to A3 increases 4.8 cfs. Increased downward flow from layer A1 to layer A2, coupled with a decrease in upward flow from A2 to A1 results in a net change in flow across confining unit C1 of 1.9 cfs. Vertical flow, across confining unit C2, changes little between the Drought and Baseline Simulations. Fluxes from specified-head boundary cells change little in layers A3 and A4. Simulated fluxes from the specified-head boundary cells in layer A2 into the model area increase 2.7 cfs. Industrial Development Simulation The effects of a hypothetical industrial development were simulated. Withdrawals of 1.54 cfs (1 MgaVd) were taken from both layer A2 and layer A3 at row 37 and column 72 (2 Mawd total). This location corresponds to a site along the Flint River in Lee County northeast of Leesburg. A withdrawal of 3.1 cfs (2 MgaVd) is relatively typical 7 for many industrial purposes; thus the simulation is considered relatively conseiVative. The withdrawal in this simulation represents a 2 and 4 percent increase in total withdrawal for layers A2 and A3, respectively. Simulated heads in each layer decline in the area of the simulated industrial development compared to the Baseline Simulation. The drawdown due to the simulated industrial development in layer A2 after 5 years of pumpage is shown in Figure 14. The maximum drawdown is 22 feet at the cell in which the well is located. Drawdown in layer A2 decreases markedly within a few cells of the simulated development. The area having drawdowns of 5 feet or more is relatively small. Drawdown in layer A3, shown in Figure 15, is much greater than was noted in layer A2. The maximum drawdown is 46 feet. The area having a drawdown of 5 feet or more is large. Although there was no change in simulated pumpage in layer A4, drawdown did occur due to the pumpage and drawdown in layer A3. Drawdown in layer A4 (Figure 16) is as much as 11 feet. Although the maximum drawdown is small, the area having a drawdown of 5 feet or more is large. Changes in simulated flux between the Baseline Simulation and the simulation of the industrial development are relatively small, chiefly because of the small change in pumpage (3.1 cfs). Almost half of the additional withdrawal (1.4 cfs) is being removed from storage. Boundary flux from specified-head cells changes very little in layer A3. In layers A2 and A4, flux into the model from specified head cells increases and flux out of the model decreases. The net change is 05 cfs for layer A2 and 0.4 cfs for layer A4. Simulated vertical fluxes across the confining units also respond to the additional pumpage. An additional 0.4 cfs leak downward across Cl. A combination of decreasing flux from layer A3 to layer A2 and increasing flux from A2 to A3 results in a net change of 0.2 cfs in the flux across C2. The pumpage in layers A2 and A3 results in an increase in the upward flux across C3 of 0.8 cfs. The City of Albany Floridan Aquifer Usage Simulation Use of the Upper Floridan aquifer to reduce the demand on the Claiborne, Clayton, and Providence aquifers is an option available to the City of Albany. The effects of a shift in pumpage from the Claiborne, Clayton, and Providence aquifers to the Floridan aquifer was simulated by reducing the City of Albany's pumpage in layers A2, A3, and A4 by 20 percent. Pumpage from the Floridan aquifer is not simulated because the equivalent layer (A1) is simulated by specified heads. The reduction in pumpage from the baseline simulation is 3.3, 1.4, and 0.6 cfs for layers A2, A3, and A4 respectively. Simulated heads resulting from the reduction in pumpage are considerably higher than the simulated heads in the Baseline Simulation in the area of Albany. The difference in simulated head between the Baseline Simulation and the reduced pumpage simulation for layer A2 is shown in Figure 17. The positive numbers indicate a rise in simulated head. The maximum rise is 26 feet. Figure 17 also indicates that the simulated head rose at least 10 feet over a large area. The difference in simulated heads in layer A3 is shown in Figure 18. It is important to note that the maximum rise in simulated head for layer A3 is 32 feet, which is substantially greater than the rise in layer A2, even though the decrease in pumpage in A3 is less than half the decrease for layer A2. The area with a rise of 10 feet or more is slightly larger for layer A3 than for layer A2. The change in simulated head for layer A4 is shown in Figure 19. The maximum rise is 22 feet. The area with a rise in simulated heads of 10 feet or more is quite similar to the area for layer A3. Significant changes in flux occur as a result of the reduced pumpage in layers A2, A3, and A4. The rate of removal of water from storage drops 2.4 cfs from the baseline rate at the end of the five-year simulation. Little water is entering storage at the end of the simulation. Flux into storage at the end of the first year of the simulation is 2.0 cfs higher than the rate at the same time of the Baseline Simulation. Vertical flux from layer A1 to A2 decreases, whereas the flux from A2 to A1 increases, resulting in the net flux across confining unit C1 decreasing by 1.2 cfs. The simulated flux from layer A4 to A3 decreases by 05 cfs from the Baseline Simulation. Flux across confining unit C2 does not change significantly. A combination of decreasing rates of water entering layer A2 and increasing rates of water leaving the layer results in the simulated net flux from specified-head boundary cells decreasing 0.9 cfs. The net change in simulated flux from specified-head cells in layers A3 and A4 is very small. INTERPRETATIONS AND CONCLUSIONS The Upper Floridan, Claiborne, Clayton, and Providence aquifers comprise an interrelated aquifer system in southwest Georgia and in the adjacent area of Alabama. Development of wells 8 to supply water for municipal and industrial uses and for irrigation have resulted in declines in the potentiometric surface in these aquifers. A digital (or numeric) flow model was developed as a tool to assist in the management of this vital resource. The model used was a finite difference modular model published by McDonald-Harbaugh (1988). The main focus of the digital model was intermediate and regional ground-water flow in the Claiborne and Clayton aquifers. In order to model flow in these aquifers, it was necessary to model the interaction between these aquifers and the adjacent aquifers. The Floridan aquifer, which overlies the Claiborne aquifer, was simulated as a specified-head layer. The Providence aquifer, which underlies the Clayton aquifer, was included as an active layer (simulated aquifer). The model was used in the steady-state mode to simulate pre-development conditions and in the transient mode to simulate conditions between 1900 and 1986. Calibration of the model in the steady-state mode was conducted by comparison of the model results to observed heads, estimates of ground-water discharge to rivers, and estimated fluxes at the model boundaries. Cahbration of the model in the transient mode consisted of comparison of simulated heads with observed heads at seven different times, and comparisons of simulated hydrographs with observed hydrographs. The relatively small difference between simulated and observed heads along with the close match between estimated and simulated ground-water discharge to rivers, estimated and simulated boundary fluxes is deemed to be indicative that the model is well calibrated. Moreover, the close congruence between simulated and observed hydrographs indicates that the model is well validated. Indication that the model can be used for predictive purposes is provided by the fact that it is cahbrated and validated. The results of the model simulations indicate that vertical leakage is an important pathway of ground-water flow. The importance of vertical leakage in the understanding of the overall ground-water flow system is critical. Unfortunately, estimation ofthe rate of movement of ground water across confining units is very difficult, and has not been measured within the study area. Measurement of the leakance of the confining units between the aquifers in the area of this study could provide information useful in refining the hydrologic parameters used in the digital model. The results of the sensitivity analysis indicate that the model is most sensitive to changes in pumpage rates, recharge, and transmissivity. The rate of recharge to the aquifers is very difficult to measure. Furthermore, the focus of this model on the intermediate and regional flow systems makes the estimation of recharge for the purpose of this model difficult. The distnbution of transmissivity estimates from specific capacity values is sparse for the Claiborne and Clayton aquifers, and almost non-existent for the Providence aquifer. Additional measurements of transmissivity are needed for each of the aquifers included in this study. The pumpage data used in the calibration of the flow model are based on estimates and projections. The quality of the match between simulated and observed heads, especially with the independent verification provided by the hydrograph comparisons and comparison of simulated and estimated fluxes indicate that the estimates of these parameters are reasonable. Measurement of pumpage for irrigation (in contrast to estimation) and defining the contnbution of each aquifer to the flow from multi-aquifer wells would significantly improve the predictive capabilities of this model. The results of the cahbration and predictive simulations indicate that the Claiborne aquifer is well connected to the overlying Upper Floridan aquifer. This connection is indicated by the changes in leakage across confining unit C1 as a result of changes in stress in layer A2. Large changes in simulated pumpage from layer A2 resulted in only moderate changes in simulated head, thus suggesting that additional ground-water can be developed from the Claiborne aquifer. This probably represents a "real-world scenario" because layer A1 probably can supply more water in nature than in the model, in which A1 heads would be allowed to decline as a result of pumping in layer A2. Additional withdrawal in the Albany area, however, is not recommended. The amount of additional water that can be withdrawn depends on where the withdrawal is located and the allowable impact upon other ground-water users. The maximum acceptable increase in ground-water use from the Claiborne aquifer is likely to be small compared to the existing withdrawal. Large declines in the simulated head of the Clayton aquifer occurred over the period of the calibration simulation. Small increases in pumpage from layer A3 resulted in large declines in simulated heads in the predictive simulations. This response to changes in pumpage rates indicates that the Clayton aquifer is heavily stressed and is not capable of supporting significant additional withdrawal in most areas. Even small additional withdrawals are likely to produce unacceptable drawdowns at nearby users. 9 The simulation of an additional withdrawal in combined layers A2 and A3 resulted in significant changes in simulated head in layer A4. This indicates that the Oayton and Providence aquifers are reasonably well connected and reinforces the conclusion that further development ofthe Clayton aquifer would have an unacceptable effect on the entire flow system. The Albany and Dawson areas are particularly heavily stressed. Additional withdrawal from the Claiborne, Clayton, and Providence aquifers in these areas generally is not recommended. REFERENCES (Includes references appearing in Supplements I and II) Applin, P. L., and Applin, E. R., 1944, Regional subsurface stratigraphy and structure of Florida and southern Georgia: American Association of Petroleum Geologists Bulletin, v. 28, no. 12, pp. 1673-1753. Arora, R., editor, 1984, Hydrologic evaluation for underground injection control in Georgia: Georgia Geologic Survey Hydrologic Atlas 10, 41 pl. Carter, R. F., and Stiles, H. R., 1983, Average annual rainfall and runoff in Georgia 1941-1970: Georgia Geologic Survey Hydrologic Atlas 9, I pl. Clarke, J. S., Faye, R. E., and Brooks, R., 1983, Hydrogeology of the Providence aquifer of southwest Georgia: Georgia Geologic Survey Hydrologic Atlas 11, 5 pl. --, 1984, Hydrology of the Clayton aquifer of southwest Georgia: Georgia Geologic Survey Hydrologic Atlas 13, 6 pl. Clarke, J. S., Longsworth, S. A., Joiner, C. N., Peck, M. F., McFadden, K. W., and Milby, B. J., 1987, Ground-water data for Georgia, 1986: U. S. Geologic Survey Open File Report 87-376, 177 p. Cooke, C. W., 1943, Geology of the Coastal Plain of Georgia: U.S. Geological Survey Bulletin 941, 121 p. Davis, M. E., 1987, Stratigraphic and hydrologic framework of the Alabama Coastal Plain: U.S. Geological Survey Water Resources Investigations Report 87-4112, 39 p., 17 pl. Eargle, D. H., 1955, Stratigraphy of the outcropping Cretaceous rocks of Georgia: U. S. Geological Survey Bulletin 1014, 101 p. Faye, R. E., and Mayer, G. C., 1990, Digital model analysis of Coastal Plain aquifers in Georgia and adjacent parts of South Carolina and Alabama: U. S. Geological Survey Professional Paper 1410-F. Freeze, R. A., and Cherry, J. A., 1979, Groundwater: Prentice Hall, Inc., Englewood Cliffs, N. J., 604 p. Georgia Department of Natural Resources, 1976, Geologic Map of Georgia: Geologic and Water Resources Division, 1 pl. Gibson, T. G., 1982, Paleocene to Middle Eocene depositional cycles in eastern Alabama and western Georgia, in Arden, D. D., Beck, B. F., and Morrow, E. F., eds., Second Symposium on the geology ofthe southeastern Coastal Plain: Georgia Geologic Survey Information Circular 53, 219 p. Hayes, L. R., Maslia, M. L., and Meeks, W. C., 1983, Hydrology and model evaluation of the Principal Artesian Aquifer, Dougherty Plain, southwest Georgia: Georgia Geologic Survey Bulletin 97, 93 p. Herrick, S.M., 1961, Well logs of the Coastal Plain of Georgia: Georgia Geologic Survey Bulletin 70, 462 p. Hicks, D. W., Krause, R. E., and Clarke, J. S., 1981, Geohydrology of the Albany Area, Georgia: Georgia Geologic Survey Information Circular 57, 31 p. Hicks, D. W., Gill, H. E., and Longsworth, S. A., 1987, Hydrogeology, chemical quality, and availability of ground water in the Upper Floridan aquifer, Albany area, Georgia: U. S. Geological Survey Water-Resources Investigations Report 874145,52 p. Johnston, R. H., Krause, R. E., Meyer, F. W., Ryder, P. D., Tibbals, C. H., and Hunn, J. D., 1980, Estimated potentiometric surface for the Tertiary limestone aquifer system, southeastern United States, prior to development: U. S. Geological Survey Open-File Report 80-406, 1 pl. Long, A. F., 1989a, Hydrogeology of the Claiborne and Clayton aquifers in southwest Georgia, 1982-1986: Georgia Geologic Survey Hydrologic Atlas 19, 6 pl. 10 --, 1989b. Groundwater modeling of the Oaiborne, Clayton and Providence aquifer system in southwest Georgia and southeast Alabama: Unpublished Georgia State University Master's Thesis, 155 p. Ripy, B. J., McFadden, S. S., Perriello, P. D., and Gernazian, AM., 1981, An interim report on the hydrogeology of the Clayton and Claiborne aquifers in southwestern Georgia: Georgia Geologic Survey Open-File Report 82-2, 66 p. Marsalis, W. E., Jr., and Friddell, M.S., 1975, A guide to selected Upper Cretaceous and Lower Tertiary outcrops in the lower Chattahoochee River Valley of Georgia: Georgia Geologic Survey Guidebook 15, 79 p. Maslia, M. L., and Hayes, L. E., 1988, Hydrogeology and simulated effects of ground-water development of the Floridan aquifer system, southwest Georgia, northwest Florida and southernmost Alabama: U. S. Geological Survey Professional Paper 1403-H, 71 p., 24 pl. McDonald, M. G., and Harbaugh, A. W., 1988, A modular three-dimensional finite-difference ground-water flow model: U. S. Geological Survey Techniques of Water-resources Investigations, Book 6, Chapter A1, 528 p. McFadden, M. G., and Perriello, P. D., 1983, Hydrogeology of the Clayton and Claiborne aquifers in southwestern Georgia: Georgia Geologic Survey Information Circular 55, 59 p., 2 pl. Miller, J. A, 1986, Hydrogeologic framework of the Floridan aquifer system in Florida and in parts of Georgia, Alabama, and South Carolina: U. S. Geological Survey Professional Paper 1403-B. 91 p., 33 pl. Stephenson, L. W., and Veatch, J. 0., 1915, Underground waters of the Coastal Plain of Georgia: U.S. Geological Survey Water Supply Paper 341, 463 p. Stewart, J. W., 1973, Dewatering of the Clayton Formation during construction of the Walter F. George Lock and Dam, Fort Gaines, Clay County, Georgia: U. S. Geological Survey Water Resources Investigations 2-73, 22 p. Thomson, M. T., and Carter, R. F., 1963, Effect of a severe drought (1954) on streamflow in Georgia: Georgia Geologic Survey Bulletin 73, 97 p., 1 pl. Toth, J. A, 1963, A theoretical analysis of ground-water flow in small drainage basins: Journal of Geophysical Research, v. 68, no. 16, pp. 4795-4811. Toulmin, L. D., and LaMoreaux, P. E., 1963, Stratigraphy along Chattahoochee River, connecting link between Atlantic and Gulf Coastal Plains: The Bulletin of the American Association of Petroleum Geologists, v. 47, no. 3, pp. 385-404. Vorhis, R. E., 1972, Geohydrology of Sumter, Dooly, Pulaski, Lee, Crisp, and Wilcox Counties, Georgia: U. S. Geological Survey Hydrologic Investigations Atlas HA-435, 1 pl. --, 1990, Ground water atlas of the United States-Segment 6--Alabama, Florida, Georgia, and South Carolina: U. S. Geological Survey Hydrologic Investigations Atlas 0730-G, 28 p. Mitchell, G. D., 1981, Hydrogeologic data of the Dougherty Plain and adjacent areas, southwest Georgia: Georgia Geologic Survey Information Circular 58, 124 p., 3 pl. Owen, Vaux, 1963, Geology and ground-water resources of Lee and Sumter Counties, southwest Georgia: U. S. Geological Survey Water Supply Paper 1666, 70 p. Pollard, L. D., and Vorhis, R. C., 1980, The geohydrology of the Cretaceous aquifer system in Georgia: Georgia Geologic Survey Hydrologic Atlas 3, 5 pl. Wait, R. L., 1958, Summary of the ground-water resources of Crisp County, Georgia: Georgia Mineral Newsletter, v. 11, no. 1, pp. 44-47. --, 1960a, Summary of the ground-water resources of Calhoun County, Georgia: Georgia Mineral Newsletter, v. 13, no. 1, pp. 26-31. --, 1960b). Summary of the geology and ground-water resources of Clay County, Georgia: Georgia Mineral Newsletter, v. 13, no. 2, pp. 93-101. --, 1960c, Summary of the ground-water resources of Terrell County, Georgia: Georgia Mineral Newsletter, v. 13, pp. 117-122. --, 1963, Geology and ground-water resources of Dougherty County, Georgia: U. S. Geological 11 Sutvey Bulletin 1199-G, 37 p. Wang, H. F., and Anderson, M. P., 1982, Introduction to groundwater modeling: fmite-difference and finite-element methods: W. H. Freeman and Company, New York, 237 p. 12 Stress Period Time Steps SPI 1900-1955 1 SP2 1944-1959 4 SP3 1960-1963 4 SP4 1964-1967 4 SP5 1968-1971 4 SP6 1972-1975 4 SP7 1976-1977 2 SP8 1978-1979 2 SP9 1980 1 SP10 1981 1 SP11 1982 1 SP12 1983 1 SP13 1984 1 SP14 1985 1 SP15 1986 1 Pumpage (cfs) A2 A3 A4 3.08 7.56 12.9 9.81 12.0 12.6 26.5 25.5 27.0 87.2 54.9 80.5 62.8 51.9 94.6 2.68 8.72 8.88 9.94 11.2 12.8 18.8 17.8 18.8 55.5 26.2 43.2 30.4 23.2 48.0 1.97 6.05 2.96 3.18 3.88 4.06 5.78 6.21 6.36 8.19 5.35 5.86 5.34 5.16 6.81 Total 7.73 22.33 24.7 22.93 27.1 29.5 51.1 49.5 52.7 150.9 86.5 129.6 98.5 80.3 149.4 Table I. Stress periods and pumpage utilized in the calibration of transient model. 13 / Barbour -~!~ 'I ' ... I ..._ Dale I I I I 0;:-_-;1,=-0----=,2,=-0---=3r-0---;40 miles Figure 1. Study Area. 14 Coastal Plain - - - Figure 2. Schematic diagram of regional ground-water now in the Upper Floridan (Al), Claiborne (A2), Clayton (A3), and Providence (A4) aquifers. Semiconrming units or units of lower hydraulic conductivity are Cl, C2, and C3. Dark arrows indicate direction of ground-water now. Region of local system of groundwater flow. Region of intermediate system of groundwater flow. Region of regional system of ~oundwater flow. Figure 3. Local, intermediate, and regional ground-water now (modified from Toth, 1963). 15 [sYSTEIII EPOCH/SERIES GEOLOGIC FORMATIONS HYDROLOGIC FORMATIONS MODEL LAYERS Oligocene N/A N/A N/A Eocene Paleocene Gulfian Ocala Limestone Lisbon Formation Tallahatta Formation Hatchetigbee Formation Tuscahoma Formation Nanafalia Formation Upper Floridan Aquifer A1 Floridan -Claiborne Confining Unit C1 Claiborne Aquifer A2 Claiborne-Clayton C2 Confining Unit Clayton Formation Providence Sand Ripley Formation Clayton Aquifer A3 Clayton-Providence Confining Unit C3 Providence Aquifer A4 Providence N/A Confining Unit Figure 4. Stratigraphic column, hydrogeologic units, and model layers. 16 84W 85W 32N 0t===j10~===23:0~==~3I0~===~40 miles --400-Eievation of potentiometric surface in feet above mean sea level. Figure 5. Predevelopment potentiometric surface of the Claiborne aquifer. 17 I I \ \ T~vior l (;{-' ' ~ \ / Macon Marion I I Stewart / / I I I I I . I I I / I ' ' 0 ~ 0 10 20 30 40 miles -----200--Eievation of potentiometric surfacE: in feet above mean sea level. Figure 6. Predevelopment potentiometric surface of the Clayton aquifer. 18 .,--~ -.,."' (.. , N I - - . . . .~/' .... ._ I ,......,~ Turner t' I Tay I Tift I __,-...(.,._,...._,, ' ........... , I ', I .._ I .. , I I ..I. I I' I ( ' ... ...... I ',...... I I ' ... I ' I I I I Dale ''I Jackson 85 Ot:===j10~==~2~0~==~3I0~==~40. miles ~Upper ~ Floridan Aquifer A1 ~ Claiborne Aquifer A2 ~ Clayton Aquifer A3 Figure 7. Finite-difference boundary conditions. Providence Aquifer A4 19 I \ ..1,.. aylor \, Manon Stewart I 85W '\,, " I 'I \ Quitman ' '\ 32oN ',, I / I ' I , , ' I , , ,,._.I.'. ...., t .......... _ .. Barbour I I I I I Turner I / Tift -, ' -, I "' I I / -' ,Q / I I Dougherty F I I i I Mitchel! ? ' ) / I ,!-..,. f I I i I I ) I I 85W Ot:===j10~==~2~0~===3I0~==~40. miles - 2 0 0 - - Elevation of potentiometric surface in feet above mean sea level. Figure 8. Simulated potentiometric surface in model layer A2 at the end of the five-year baseline simulation. 20 \ \ Marion / / 32N / I I I I I / / / I / / ' 85W oc====1:Eo~==~2~o====:J3o~===~4o miles - - - 2 5 0 - - Elevation. of potentiometric surface in feet above mean sea level. Figure 9. Simulated potentiometric surface in model layer A3 at the end of the rwe-year baseline simulation. 21 0t===:j10E:==~2~0====3I0~==~40 miles - - - 300-- Elevation of potentiometric surface. in feet above mean sea level. Figure 10. Simulated potentiometric surface in model layer A4 at the end of the five-year baseline simulation. 22 I I I \ \ \ Taylor\ Macon .,' '-,II I I I I ' \ ' \ > I 'i" I ... ' 1,, (. I 1' ....... , , '' '...,' ' ,, I I I '......... ,'I ' ......I... ' I sp I I Turner I I I L I ,/ ,, ,... ', I _( Tift __, ........ ' I ..., , .... ... I ' I '-,I '.._, ' ... , I ,I ' ... , Schley / ... , / ,, I Sumter ', I I ' I I Marion '-...(, ,', / I ', /',I I ' .., 1 __.-../ f ........ ,--........ ... , .._ .. -..., ',<_.-.- / Worth ... '>~~ ',_ I -, I 11 I I I I ... , I ', I / ' ',, / I Colquitt I ', 'I: ( ... / ', Webster ', I ', , ', ' I 1 I I 1 Terrell Dougherty 1~I I I \ , I ,} ~, ',', , ', ' ', ' I I 1 ) 84W Stewart ..., , I ' ' ,, ... , , I I I I I' Mitchell I I I I , I I / ' - , ........1- . , _ ........_ ... , 1'-, I I ....... --... ,, ........... , / I I Randolph , I , ~ I I Baker I ), ' ' / ; I I I I 85W ,, I t'', I, /1 , , ',,, I " ' ... , , .. , I II Calhoun / I /1 /1 I , ( ' ,.............. I /Grady I ,I , I ' , .... \ Quitman ' 1' ' - I ', I ' ... I ...... ' )_ I ',, ..., ' '-1' ' I ', / '',"', ' .., I I 'I I I '.,' I ' ',,,I,..1...,.........--......._ __ , Clay I 1 1 I I I < I I ' .... I \ ' \, . , . . . . . ., . I 'I ,I ' ... I Decatur ,.......,. I I ~'-.. / ,,. I '--, Barbour ,I ~~ I ' ' ~ I \ ~ ' ... , ~' I I Sem1. nole ' ' ' \ ... , ' , Henry / '......._GE,O-.R-G-.I.Aw,.' \ I 1',, / -ALABAMA').. - I ,1 / I 1 ,' ' , '.... ' ... , ..... , , I / I Houston 1I .. __ '. I - -.......... FLORIDA ........ -\#~(, ~ 1 Dale I ', 1 ... , I / Jackson ; I I 0t===~10E:===32o====3]0~~~~40 miles m -5 feet ~ -10 feet -15 feet ~ -20 feet Figure 11. Change in simulated head in model layer A2 in response to a hypothetical twoyear drought. 23 I I I I I Macon I Dooly I I ,-- '/' ..... _,.. I I I \ \ (IIJ [ill / Crisp 1 Turner I > Taylor \ \ / I ,L----"-..._ __f_ / 1 Tift I /. < I / ---... ,/ ,, .. _,-...._,'-', ... ~' ... / ........." ,,.. ', I '..._, 1;...... I' (. 11 ..... ,1 ..... I ' ... ..., I' I I '-, I \ '' I ....... / ........ ..... I I ........ I ,I',' I I ........ ..... ........... , Schley I \ I 1 , -\ Worth 1 / I I Marion ..... I ''-< 1 ........ I I 1'- I ' ' ..... ,.._.....,I ' ', ..... I'-. I ..... I ', I '..._, Webster I ..... , ..... Stewart '..._, .....,I I I I I I / '\ - ''"~'~') ',, I ........ ',, I I' I I I ;' ,'1',, 1 Colquitt ...... ' ( I ..... , y I I I' ! l I \ . / ) I I I I I J ' "' \-.'I \ I I Baker } ' ........ '..._, ............. 84W' Mitchell I I I 85W 32N Barbour Jackson 85W Ot::===j10E:===23o~===3]0~==~40 miles a IIIIIliill -5 feet ~ -10 feet ~ -15 feet -20 feet -25 feet -45 feet Figure 12. Change in simulated head in model layer A3 in response to a hypothetical two- year drought. 24 __, 8,.4.-W"-/I I I \ I I Macon '1 ( ' \ Dooly I '..._ I I I 1 II I 1 Crisp I 1 Turner Taylor\l \ ,/ L ___ _,.-._ I I / I Tift I ( I.. 1~ ......... I' I'..!'/, I ,, ...... ,-... ,~'-'I' ... .., ' llli.1!~+L>o .. , Schley ... ~.... I ,t[ ~ IVI.AII m:a ------.... , /....... , , _ _,-......,-, .......... ,... ,i', t' '- I I I .... I '..... ' I I ' ......... .......... I ' .. , ,1 / \ '...._/ I I I -\ Worth I 1 ' .... I I F== ' ' < ...... I Marion /'..... "':'~ '> ' , , '- '- ... ~"...' I ...... I ................ 1 1 I I I I I 1 I ' ' ... / ', ' ~] ' I / 1''-, ....... I I Colquitt .... , ,--...1 / 1'--.' ' , , WvE:lb~S1,:e~rl I I ~~ugherty /\ ' .. , '',,, I '-, I . ' , , illt:: ' ,,,I Stewart 85W I I I / ,/ Riiiidl Jh ~II iI ~ ~!'~ ~ ~ Mini / I I ~ ' ') I I I' Mitchell I I ) 'y'' 'I', , / Baker \ ', 84W' I I 1 .. I ,/~rady ' \. I ' ' '-,1I ........ , _ ' I ~ .. , .., ,I ...... I , Quitman /'- I I , (,I ,I , I 1', I I ...... .... ' 1', I ... ' I I ',, I I ............ , , .. , ,' ' ... , -, , ' 32N ', ' ,, I I I ....... I I ,._._,,..._ ........... , . . . ' ,I 1 Cl ay l ----,... --..........-}', ,I I --- I I,I / .. , I ,"' .,, I I Miller I ... 'I ..... .,, I I I <. I I \ ...... I 1 1 Decatur \ '\. Early 1 /' \, I ... ............... I .......... / ,., '- I '--. Barbour / I \.. -~ ' ' ....... , / Seminole '-., I 1 ,_ I ( , .... , Henry I ... , / I I I I I I I ................ I ', '-, / / Dale ........... I '-( I I I \ I ' / '.... GEORGIA / \ "-....,A_L.A-.B...A...M-_A~)\ - I I Houston -- I - ---- ....... / / FLORIDA ......... , \ .. ~(, ' '. I \ I Jackson I 85W ot:====1:E0===~2r0====3]0~===~4Q. miles -5 -10 feet ~ feet -15 feet Figure 13. Change in simulated head in model layer A4 in response to a hypothetical two- year drought. 25 I I I .,"',.... ...I I I ' Macon I ( Dooly I I I I I ' I I -.... ... I < ' ... .... ... I I ...... I ' \ , I I \ 1 Crisp I Turner Taylor\ \ \ / L ___ _,.-.._ / I / (, / Tift I ' I ,'-,// ~... .I.. '-,,I ' ~-I.... I' 'v ', 1', I I) -, ' I '',"-,_:chley // sumter - - .. _ .;' - "':., --......... , .,...,/ '',, / __, ...... ..., /' ~lil~l~ ' ...' II ''......, I ' I I I // LeI'~l~ mI=~ '-~-' ' worth ' I ' '-' ~~ /,) -' / / '> ..... _ / Marion '....; ... , ' I I .... , t---, I , , .... I ,I ' ... , ... ........... __1_, ..-" / , - - '..I........... '"/<:-.1 '<~ ..,_ 11 I / I ,/ '... ' ......... ,I I I''-, / Colquitt I I ......... I / .......... / ' ... , Webster / / -,_ I ' ........ ' / Terrell ' ... , / Dougherty 1~ I I '. / I ~ ... , .. , 84 o'w':.: ' I I / ') ' .., / ' ..,, / / J Mitchell I II ' ... ,1- . . -. - - - 1/1' ......... - -. . . _____ / I ~I, / 'I1 Stewart / / I '"-8 5 ow I 1 -, 1... I I........ I '-,1 -........ ',,, / -... 'I 't Randolph / 1 / l ' c a lh oun I I l1 I 7' 1 } I 1 I / 1', I ........ ' , Baker , I ' t' i 1', , I ' ) / / I I 1 Grady II I ' . . I ', ;' ........ , '\ Quitman /'-,_ / ' .. , / ' ... , , / ' .. ,_ ' ' - / ' ... ', ,, 1 I 1 / ' I '..... ~--I-,.... '... / I ....... c lay ---~-~"'---' 1 I I 1 I 1 I ) I .,,I , M'1ller '/ I I ! I ( ----....... / \ .......... , I I ,I I I I Decatur ' I I I ,I \. ... '.. Early I .I.. .... ... .... ... 1 ' I / \. ......... I ---.. Barbour 1' I ' \. -...,,.... ' ... ' , I f 1 Seminole ' ',, / ( .. ' ... , Henry I / ' \ ' .... ,GE~!!_~~./ / ' '; I ... , 1 ', I .......... 1 , I ... / '......... 1 I / I I I I -ALABAMA"~ - ' I Houston I I - ...., _ _. FLORIDA .., ........... . .. ,\ ~111 _.. ,, / I I \ I 1 Dale ... , .../, ,I Jackson / : I 0 10 20 30 40 miles IDIIJiill -5 feet ~ -10 feet -20 feet Figure 14. Change in simulated head in model layer A2 due to a simulated withdrawal of one million gallons per day from both layers A2 and A3. 26 I I \ Taylor\ l I ';--... .I , , Macon I ( ' \ ' \ Dooly - ' I ' I ' I ' ' I ' I I I '1, I I ' I I I I Crisp / ,L----'-.._ __!_ --.m'r""m.,.,.,.TTT'M,., , / ....... ,;' - """ -,~N Turner ' ,J-..J I I I I I Tift I _ -_, .(.."..'.:........ I .. , I ', I I I ' ' ' - , J ... , I " ') ..''-......' I ' ... , Schley / ' I ' I I I I Marion ' '-...(I , ...... , I /',11 ' ... , ... , Sumter I I I I 1, ,--, _,.'.-/ I ', I ',1 Worth I 1 / I ....... I ',, I '- I ',I I ..I. I I' I I I I I I I Colquitt .... , ,--./ 1'-, ',< .. -- .......... I I ~ ... ,_ / I I 1' ' ...., Webster / ' , -... 1 1 Dougherty I I' 1 '-..., I I I ,, 1 Terrell ' ' I I I \ / I I ~ ' ') Stewart ... , , / ' I '-1 ......... I I ... , 1'- ' '......._1-..'-... .. _ I ... .....__ I -, I I ~ I I Randolph / I I I I ' I -., I I \' I ) \ I I ' ' Baker Mitchell I I I' I I I I I I } I I z 1 ' lcalhoun 1 , 1 I ', 1 Grady I ,, ;'',,t' ................. , ',, I 'I ' \ Q U.ltman ' 'I / I ' ,, ', I ', I /1 / /", I, I II ' ' ' I ' / , ' /'... '';~"'---'/ ' '-..., I 1 1 I ' ( I ) ' I I ' I ' , ... , I ............ ... .... , _ ' - / '/ ,, I '1 I ,'1 I ' ....... I I Cl ay I'I , ......._..._ ! "--.---......... ;' ,I ., I Miller I I (, I I I \ ........ I I ',, ..... ' ,I EarIy I I Decatur I -.. I I ' ... I .............. I I ', ......... Barbour I I / I I < ' ... , Henry ' ... , ' ' I I ... , \ ' ' , , 1 Seminole ', .. ,... 'I ' I \ I I '\ / -.. ......,GE~!!_~~./~ "'1 1 ' ' ... ,.... / 1 1/ ' .. ,, / I I I I Dale '.... / ,_ I ' I-, 1 -ALABAMA).....- .._--.. \ / Houston / I I I I FLORIDA ........... _ ... :::(, \ ' Jackson 0't===::J10E:===230====3I0~==~40 miles liliJI -5 feet ~ -10 feet ~ -15 feet ~ -20 feet lmi -25 feet -45 feet Figure 15. Change in simulated head in model layer A3 due to a simulated withdrawal of one million gallons per day from both layers A2 and A3. 27 I I I \ \ ' Taylor\ } Macon ,.I ' I Dooly ' I ' I I I I I - ...... ... ...... I I ......I... ..._ ( ' I I I \ ' I / Crisp \ / ,L ___ ...,-..__.!_ I / ... I ,._ '...._ I..._ I' 1'-._l I' ,1'..._1 .., 'v' ' '') I ---, ' E , ' ..., ... Schley / / I Marion Sumter / I I I I 1~1~ f ''~IIill .,-,~N -- , '-, ...... , I .... ._," / ' ....... }'._/ I Turner I Tift Worth I I' I I Colquitt Stewart Mitchell I Barbour Decatur 0t===:j10E===23o~===3l0~==~40 miles -10 ~ -5 feet feet Figure 16. Change in simulated head in model layer A4 due to a simulated withdrawal of one million gallons per day from both layers A2 and A3. 28 I I I ., I I Macon , ( Dooly . I I I I .... ... .... ... .... ... I I I < 1 ..... I -,~N -- --... '... ......, I I '/ ... ...... __,~ ,....,.~ I I \ \ I I Crisp I Turner I I Taylor\ l z " I ... , f ... , \ / -- ~L---J-.._ __/ / I , ........1 (, I~ I ' '-......' .. , .... ....,' ' ... ) I ' ... , Schley / ........ / , ........ ... I I Sumter I I I Tift I ( Marion ... ' ' ,..._.....,I I' ... II t' ... I ..... 1 ... , Webster I ... ,... I 1 I ', I Terrell ' I ..., , I I ' I ........ , .......... , I Stewart 85W ll I I I I I 1 I ' ..., r......-..__ / '...... Randolph ' ~I / I I 1 / Calhoun I I , ........ .... I .....(.. Colquitt ... ... .... .......... ......... ...... . 84W I I I I I I '-..., I 'I i I I \ ' Baker ' Mitchell I I 1 ... I I I I I / , /Grady .. , II', 1/',, ' ... '', I 'I ....... ' I I / 1 I II I 1', I ....... .... { I I'..._ I .... ) / I ..... I '-... ... I ............ \, Quitman ', /'......... I 1 '... / ' ... , / / ... , . .,....,___ / :I ' ... , / ' . ,..,/ '... , / ' 32oN ,, I ', I ....... I ,--...... 'I / Cl ay I I 1 l ----....... I \ ,,. I \, ,,{ Early ;> ,I 7 1 I /1 I Miller ! I I I /(, Decatur Barbour I I / / I \. ... ' , ' .. ,.._ ............. I \. .......... ', ' I I '~ . ' ... , 1 Sem1nole ', ',1 ' { I \ I ' , ........... Henry / '... GEORGIA / I ...... , . .............. _ '.I I', ' I I I I ... , ....... ' ... , I I I 1 I / -ALABAMA/' ~ .-.._. I \ I I F-LORI-D--A-......, _,:'::-1.' Houston I , '- , I I I Dale '..._, I ' ..( / \ / Jackson I I I n +5 feet Ot:===:J10~==~2~0====3I0~==~4.0 miles ~ ~ ~ +10 feet +15 feet +20 feet - +25 feet Figure 17. Change in simulated head in model layer A2 due to a 20 percent reduction in pumping from the Claiborne, Clayton, and Providence aquifers for the City of Albany. 29 ' ',1I .I , ' I I '' I I I I I I Macon , ( Dooly 1 II I \ ' \ I I Crisp I I Turner I ' Taylor \ , \ , I 1 I I I Tift } ( ,..... I .. , 1'..._ - ,L---..;--'- __/ , ....... , /,/ 1.. /!., ___-:I~.. ' ............ ,.., -, I ' 1'.. t' '....._ I ........... , I "-, ! ............... .... I II '' ,...,1 ' ') I \ ,/'......._' ' ... , , I Schley / I Sumter I \ I \ I I I ' I I I I I I I ' I I ... , I I ... , I I ' I I' I I Marion '-...( /'.... ,' I ' ' , ,.._ . -..., I I' ' I I' ',, I ,--, ',','<___ _':,.::./ ..I.. . . . . i'... 1 I '.... '... ~ Webster 1 I I I ', ', I I Terrell ' ............... ,'...I... II I I I I ', 1 Colquitt '-,t ' '...._ ... , , , . ' . 84W' ....... , ... , I I ... , 1'- ' ',1-...'-.. ... -... I '......__ I I I Mitchell I Stewart asow I I 1 / 1 l 1 Randolph ' ~I / I / I Calhoun -, I \1 / I ' Baker ' I / \ I I II / Grady 1 "\ .. I ' ' '-1I ' , ' ', 1 I I I ' ( , I II ' '', 1 ', / 1", I' I ',, 'I I /', / ', Quitman I/'- ' .. ,, / 1 , .. ,_1.,...., / 1 \' 1/ 'I / --- ,.... / Clay ,--..... l , I ----.... -...__ ',, ' I I / 1 I I I ..,1 \. ..... EarIy ',, / '' )_ I ... , 1 I'~ I .. , , ' - / ' 1-l ..,~~ 1 1 I Miller ', I I I I Decatur I 1 -. I ...,... I / \ . ......... ',, I Barbour I / I / ' .., \,. .. .,... I \ ', ' I I - .......' , ', / Seminole ', I ', I \ - , .... , Henry I .... , I ... , , 1/ ' .. , / / I / I I / I Dale ... , / '...._, I '{ I ' ' '-G-EA,O_LR-A.G.B...AI..A-.M~,. 'A'~ - \ I I I , - , _ _. , \ ,' Houston / FLORIDA ............ _ ..#~~, I ' / 1,. I Jackson I Ot:===j10~==~2l0====3I0~~~~40. miles 1111 mil +5 feet ~ +10 feet ~ +15 feet ~ +20 feet +25 feet +30 feet Figure 18. Change in simulated head in model layer A3 due to a 20 percent reduction in pumping from the Claiborne, Clayton, and Providence aquifers for the City of Albany. 30 / -........... / -,_,~ <' / -......( , I I Macon ' ( Dooly 1 / ' N /' I, '.. ~---.....; ' ......... I, \ \ I \ I 1 I crisp I 11 Turner ,.., 1 I > Taylor \ \ / L ___ _,-..._ / I / / 1 Tift I ,/ ';-... .... ,... (.,......../ ------.... ,/-',, ';;,_ 1,,.. ............ -.....,'-, ,1--- -............ _ I ' .., ' , ' 11 '_, , _1 ... ..._, I ' ...,... Schley 1 , 1 1 ...., _ I ..... / Sumter I I / 1 / ,, \ ' ~ I 1 Lee Marion '-..< ..... I ,........ I / I I ', I I ' ' ,..--, .. , _...,I ''' I . .,<____ ' ..._ .....-:./ .I. , .... , , , Worth - I ' .. , ,1 -,I I / I' I I I I , I I Colquitt ... 1 ..... 1'-. 1 I ' ' ,...... Webster I 1 I I '-, I , , ' , ..._ I I .... ,'....... Terrell ,' ......_, ',, I I } Mitchell 11 I '- .. ..._1- ..-. .. _ / '..._ 1 .. ..._ I ----.. 11 j , Stewart , 1 Randolph ~ / \1 \ 7 ' ' / I I Baker ' I I I II / I I 85W , I, 1' Calhoun / I , /Grady ' I .. , . I ,,',1I ' ....... - "'', I ', I ( I 11 ,I 1', I I I', II ' I ', I ,.. , Quitman '1 I, 1 ... , 11 ' . , -, I I ..... ........ -.. I .. ) I I , .. , / ' ... ,_ ,, , ' ....... ' .. ', I ....... I .. ,_.,...., I 'I -1 \., 1/ ............ / ,/ ---' ,; ~.., \~...~ / 1-......... Clay / 1 .:1 M"l ,I ,' 1ler <, \ ----........ 1 \ ,, I ',, ' ..,1 I I EarIy I I I Decatur '' I '.. (, I ', I . ... I -... ,I ' .. ,_ I .. , Barbour / I I ' \.. "~ I \ ~I ' .. .... 1I Seminole ' ', ' I \ <.. , Henry / ' ... GEORGIA / \ ', I ..... ,.......... I I', / -ALABAMA"')\ - I ,1 I / 1 I I ',, '... ' .................. -... I I / I Houston I ... '-----..... I I -- \ FLORIDA ...., -~\ I ":::-( 1 I \. / Dale ' ' ..( / Jackson I I I Ot:===j10~==~2~0====3I0~==~40. miles liiiiillil -5 feet ~ -10 feet -15 feet -20 feet Figure 19. Change in simulated head in model layer A4 due to a 20 percent reduction in pumping from the Claiborne, Clayton, and Providence aquifers for the City of Albany. 31 APPENDIX A: BASEFLOW ANALYSIS Streamflow data for the Chattahoochee River allow estimation of ground-water discharge from the combined Providence and Clayton aquifers, and from the Claiborne aquifer. Gaging station 02342960 at Eufala, Alabama is located very near the northern edge of the Providence aquifer outcrop area. Daily mean discharge at this station on October 5, 1954 was 877 cfs. Downstream, at Fort Gaines, daily mean discharge at station 02343260 on October 5, 1954 was 1,120 cfs. This station is located approximately at the southern edge of the outcrop of the Clayton aquifer. Daily mean discharge was used because minimum mean monthly streamflow was not available for these stations. Comparison with upstream and downstream stations that have daily and minimum mean monthly values indicates that the flow on October 5, 1954 approximates the minimum mean monthly flow. The reach delimited by these two gaging stations approximates the recharge area of the combined Providence and Clayton aquifers. The respective tnbutary streamflow into the Chattahoochee, however, must be subtracted to estimate baseflow discharge into the River. Streamflow was measured and minimum monthly streamflow computed for the Tobanee, Pataula, Holanna, and Cemochechobee Creeks in Georgia. Stream flow from the unmeasured lower reach of Pataula Creek as well as Barbour, Cheneyhatchee, and White Oaks Creeks in Alabama was estimated using the drainage area-ratio method using the Pataula Creek partial record station data. Addition of tnbutary streamflow and application of probable error produces a range of streamflow of 187-247 cfs with an average of 217 cfs. The minimum net gain along the reach of the Chattahoochee River between Eufala and Fort Gaines was calculated by applying the appropriate error and subtracting the maximum values of streamflow from the upstream station, (Chattahoochee River at Eufala [877 cfs x 1.05 cfs = 921 cfs]} and tnbutary streamflow (247 cfs) from the minimum value of streamflow for the downstream station (Chattahoochee River at Fort Gaines [1,120 X 0.95=1,064 cfs]). This yielded a minimum net gain of -104 cfs (1,064-247-921), which is interpreted as a minimum value of aquifer discharge from the Providence and Clayton aquifers. Average net gain along this reach was calculated by subtracting the value of average streamflow of the upstream station and the tnbutary streamflow from the average streamflow of the downstream station (1120-217-877). This, in tum, yielded 26 cfs. Maximum net gain was calculated by subtracting the minimum tnbutary streamflow (187) and minimum upstream streamflow from the maximum downstream streamflow {1176-187-833) yielding 156 cfs. The headwaters of lchawaynochaway Creek are in the extreme southeast comer of Stewart County and the southwest comer of Webster County. lchawaynochaway Creek dissects the Claiborne aquifer outcrop through Randolph and Terrell Counties as it flows to the Flint River in Baker County. Tnbutaries to Ichawaynochaway Creek are Nochaway Creek and Pachitla and Carter Creeks. Minimum mean monthly streamflow for Nochaway Creek just above the confluence with lchawaynochaway Creek is reported as 24 cfs (Thomson and Carter, 1963). Minimum mean monthly streamflow for Ichawaynochaway Creek just south of the outcrop downstream of the confluence of Pachitla and Carter Creeks upstream of the confluence with Ichawaynochaway Creek is reported as 51 cfs. The minimum net gain of lchawaynochaway Creek as aquifer discharge for the Claiborne aquifer was calculated by subtracting the maximum value of Pachitla and Carter Creeks (56 cfs) and the maximum value of Nochaway Creek (24 x 1.10 = 26 cfs) from the mmrmum value of Ichawaynochaway Creek (140 X 0.90=126 cfs) yielding 44 cfs. Average net gain was calculated by subtracting the mean value of Pachitla and Carter Creeks (51 cfs) and the mean value of Nochaway Creek (24 cfs) from the mean value of streamflow from Ichawaynochaway Creek (140 cfs) yielding 65 cfs. Maximum net gain of lchawaynochaway Creek as aquifer discharge from the Claiborne aquifer was calculated by subtracting the minimum value of Pachitla and Carter Creeks (51 x 0.90) and the minimum value of Nochaway Creek (24 x 0.90) from the maximum value of Ichawaynochaway Creek (140 X 1.10) yielding 86 cfs. A-1 SUPPLEMENT I: STRATIGRAPHIC AND HYDROGEOLOGIC SETTING STRATIGRAPHY General The Coastal Plain sediments form a seaward thickeningwedge ofgenerally unconsolidated units that range in age from Late Cretaceous to early Tertiary (Miller, 1986 and 1990). Successively younger sediments crop out seaward of, and overlie older sediments. Units in the Alabama part of the study area trend along an east-west line and dip to the south. Units in the Georgia part of the study area trend along a northeast-southwest line and dip to the southeast (see Figure 1-1). Figure 7 of the text is a generalized stratigraphic column of Upper Cretaceous and Tertiary sediments relevant to this investigation. From oldest to youngest, the lithologic units are descnbed below. Ripley Formation The Ripley Formation conformably underlies the Providence Sand. In the study area, the Ripley Formation is a fine-grained marine sand that is generally massive, clayey, calcareous, and highly fossiliferous (Eargle, 1955). The Ripley Formation undergoes a facies change to the east and becomes a clayey coarse-grained sand between the Flint and Ocmulgee Rivers (Clarke, et al, 1983). Providence Sand The Upper Cretaceous Providence Sand conformably overlies the Ripley Formation. It consists of medium- to very coarse-grained, cross-bedded sand (Marsalis and Friddell, 1975). The distinctive lower Perote Member of the Providence Sand consists of silt and fine-grained sand (Marsalis and Friddell, 1975). Farther west in Alabama, the upper part of the Providence Sand consists of cross-bedded, fine- to coarse-grained sand and clay with some lignite and kaolin. The basal Perote Member in Alabama consists chiefly of silty clay and very fine ro fine-grained sand. The Perote Member thins to the east and is not recognized east of the Flint River (McFadden and Perriello, 1983). Sediments of the Providence Sand become finer grained southward and downdip because they were deposited farther offshore, in the Suwannee Strait (P. Huddleston, 1991, personal communication). Clayton Formation The Clayton Formation unconformably overlies the Providence Sand. The lithology of the Clayton Formation varies considerably over the study area. Toulmin and LaMoreaux (1963) divided the Clayton Formation into three lithologic units. The lower unit is a basal conglomerate with overlying beds of calcareous sand and sandstone. The middle unit is a coquinoid limestone. The upper unit consists of soft marine limestone. In eastern Alabama, the Clayton Formation consists of a basal sand that is overlain by relatively pure limestone. This limestone is in tum overlain by sandy limestone and sandy clay that comprise the ~pper part of the formation (Davis, 1987). Very httle of the Clayton Formation is exposed at the land surface. The exposed limestone has undergone weathering that has produced a sandy clay residuum. The Clayton Formation was deposited on an irregular surface and later underwent a period of erosion resulting in variable formation thickness (McFadden and Perriello, 1983). Wilcox Group The Wilcox Group is composed of three formations: namely: the Nanafalia Formation, the Tuscahoma Formation, and the Hatchetigbee Formation. The Nanafalia Formation unconformably overlies the Clayton Formation. It is a massively bedded, fine-grained, glauconitic sand and sandy clay. Updip, north of Fort Gaines, the Nanafalia grades into a non-fossiliferous sand with some kaolinitic clay. Downdip, the Nanafalia is a silt and fine-grained sand with some siltstone nodules (Marsalis and Friddell, 1975). In eastern Alabama, the base of the Nanafalia Formation consists of cross-bedded, medium- to coarse-grained sand. The upper and middle parts grade laterally updip into a sequence of interbedded micaceous sand and kaolinitic and carbonaceous clay. The Tuscahoma Formation overlies the Nanafalia Formation. It has a basal quartz sand that is overlain by thinly bedded silt and clay that is interbedded with fine-grained sands (McFadden and Perriello, 1983). The basal part is more calcareous downdip. The Tuscahoma Formation does not change lithologically to the east in the study area and is readily distinguished from other formations (Marsalis and Friddell, 1975). In eastern Alabama, the basal part of the Tuscahoma Formation medium- to ccoonasrsiset-sgrapirnimedarsialynd.ofThcerousps-pbeerdpdead~ consists of silty clay. The Hatchetigbee Formation consists of sand and marl. It conformably overlies the Tuscahoma Formation. The marl is discontinuous and is interbedded with massively bedded sand. In eastern Alabama the updip part of the 1-1 Hatchetigbee Formation is sand; and the downdip part is the Bashi Marl Member. The Hatchetigbee is less sandy in Alabama and resembles the Tuscahoma lithologically. It is difficult to separate the Hatchetigbee from the Tuscahoma except where it interfingers with the Bashi Marl. As a general rule, the Wilcox Group sediments thicken and become fmer grained southward toward the Suwannee Channel axis (P. Huddlestun, 1991, personal communication, unpublished structure contour and isopach maps). Claiborne Group In the study area, the Claiborne Group is made up of the Tallahatta Formation and the Lisbon Formation. The Tallahatta Formation typically is a massive, quartz sand having increasing amounts of clay and shell fragments toward the top. The clay content of the Tallahatta increases in the eastern and northern parts of the study area (McFadden and Perriello, 1983). In eastern Alabama, the Tallahatta Formation consists chiefly of sand and sandy clay with beds of 'buhrstone' or massive siliceous siltstone present near the bottom and top of the formation. In northernmost outcrops, the basal part of the Tallahatta is in channels on the upper Hatchetigbee surface. The channel deposits consist of medium to coarse-grained sand that commonly is gravelly and cross-bedded and contains clay clasts at or near the base (Gibson, 1982). Downdip, the base is sandy and less glauconitic than Wilcox deposits. The upper part of the Tallahatta Formation is dominantly fine- to medium-grained sand, but consists of fossiliferous and sandy limestone in many places (Gibson, 1982). The Lisbon Formation is a dense, brownish-gray marl with thin beds of sandy limestone (McFadden and Perriello, 1983). Near its contact with the Tallahatta it is a calcareous sandstone. In eastern Alabama, the Lisbon Formation consists of fine- to medium-grained, quartzose sand, sandy limestone, and calcareous sandy clay. In updip areas, the Claiborne Group weathers to a moderately reddish-brown sand, which locally contains light-green waxy clay zones and silicified fossil fragments (Marsalis and Friddel,l1975). The Claiborne Group generally thickens to the south-southwest, ranging from 50 feet in the northeast to 200 feet in southwest Calhoun and Early Counties. The formations are very similar except that the Tallahatta is more clayey and siliceous, and the Lisbon is more calcareous and fossiliferous (Marsalis and Friddell, 1975). Ocala Limestone The Ocala Limestone unconformably overlies the Lisbon Formation in the study area. The Ocala consists of two rock types in most places (Applin and Applin, 1944). The upper part is a fossiliferous limestone; whereas the lower part is a fine-grained and granular, soft to semi-indurated, micritic limestone (Miller, 1986). The Ocala Limestone outcrop is locally irregular due to the dissolution of limestone. It ranges in thickness from a few feet at the updip limit to about 350 feet in the southeast part of the Dougherty Plain (Hayes, et al, 1983). Its residuum crops out along hilltops and uplands in the southern Fall Line Hills District. HYDROGEOLOGY Aquifer Framework This study is concerned with ground-water flow conditions in a series of interconnected aquifers contained in the wedge of sediments previously descnbed. The aquifers of primary interest are the Claiborne and Clayton aquifers. The Providence and Floridan aquifers were included in the study because oftheir interconnection with the Claiborne and Clayton aquifers. The relation between these aquifers and stratigraphic units is summarized in Figure 7. The brief descriptions of these aquifers in this report are intended to provide information necessary to understand the conceptual model of the ground-water flow system. A digital ground-water flow model was developed based upon the conceptual model. Detailed descriptions of the aquifers can be found in the references cited. Upper Floridan Aquifer The Upper Floridan aquifer is comprised primarily of the Ocala Limestone in the study area. Locally, permeable sands within the underlying Lisbon Formation also are considered part of the Floridan aquifer. The upper surface of the aquifer dips generally to the south and southeast. The outcrop of the Ocala Limestone and its residuum extends from Early County to Dooly County. The Floridan aquifer is simulated in the digital model only as a source/sink layer. Therefore, hydraulic properties of the Floridan aquifer are not relevant to this study. The reader is referred to Miller (1986), Mitchell (1981), Maslia and Hayes (1988), Hicks and others (1987), and Hayes and others (1983) for information on hydraulic properties of the Floridan aquifer in southwest Georgia. The Floridan aquifer is recharged by infiltrating rainfall and by vertical leakage from its 1-2 overlying residuum and from the underlying Claiborne aquifer, where head differences and the leakance of the confining unit permit. The Floridan aquifer discharges to streams that incise the aquifer and to adjacent aquifers by vertical leakage. Claiborne Aquifer The Tallahatta Formation is the main component of the Claiborne aquifer. Sands at the base of the overlying Lisbon Formation and at the top of the underlying Hatchetigbee Formation are locally present and are included as a part of the Claiborne aquifer. Clay layers of the Lisbon Formation separate the Claiborne aquifer from the overlying Floridan aquifer. The Hatchetigbee and Nanafalia Formations contain clay layers that separate the aquifer from the underlying Clayton aquifer (McFadden and Perriello, 1983; Long, 1989). Sediments that comprise the Claiborne aquifer crop out in a band extending from northern Early and southern Clay Counties northeastward to southern Macon and northern Dooly Counties (McFadden and Perriello, 1983). The Claiborne aquifer thickens toward the south and southeast. McFadden and Perriello (1983) suggest that the thickness and lithology of the sediments that comprise the Claiborne aquifer are relatively uniform. The sediments of the Claiborne aquifer become siltier in the vicinity of Decatur, Mitchell, and Worth Counties (P. Huddleston, 1995, personal communication). South of this area, these sediments are not considered to be a viable aquifer. Reported values of transmissivity range from 700 to 14,000 feet squared per day (fWd) (McFadden and Perriello, 1983). Many of the reported transmissivity values are based on specific capacity and should be considered as estimates only. The storativity of the Claiborne aquifer has been reported at only two locations in the study area, both near Albany. McFadden and Perriello (1983) report that the storativity at the two sites ranged between 0.001 and 0.0003. The primary source of recharge to the Claiborne aquifer is infiltration of precipitation where the sediments that comprise the aquifer crop out. Exclusive of the Floridan, the area of outcrop of the Claiborne Group is the broadest of the aquifers being studied. McFadden and Perriello (1983) estimate that the area of Claiborne outcrop receiving recharge is approximately 350 square miles. The leakance of the confining unit between the Floridan and Claiborne aquifers permits leakage of water to be an important source of recharge and discharge to the Claiborne aquifer where heads are sufficiently different in the two aquifers. Head relations between the Claiborne aquifer and the underlying Clayton aquifer are such that it is unlikely that the Claiborne aquifer receives significant leakage from the Clayton. Much of the water that recharges the Claiborne aquifer discharges to streams and rivers that dissect the area in which it crops out. The aquifer outcrop is incised by a number of creeks and rivers including the Chattahoochee and Flint Rivers and Kinchafoonee, Muckalee, Ichawaynochaway, Nochaway, Pachitla, and Carter Creeks. Smaller creeks and tnbutary streams also receive discharge from the Claiborne aquifer. Downdip, the Claiborne aquifer discharges water by vertical leakage to other aquifers, and by pumpage of wells. Clayton Aquifer The Clayton aquifer is composed mainly of limestone from the middle part of the Clayton Formation. Sands in the upper and lower parts of the Clayton Formation are locally present in hydraulic continuity with the limestone, and are included in the aquifer. Clay units within the upper part of the Clayton Formation and the Hatchetigbee and Nanafalia Formations separate the Clayton aquifer from the overlying Claiborne aquifer. Silts and clays in the lower part of the Clayton Formation and the upper part of the Providence Sand separate the Clayton and Providence aquifers at most locations. The outcrop of the sediments that comprise the Clayton aquifer stretches from southern Quitman and northern Clay Counties, eastward into west-central Macon County (McFadden and Perriello, 1983). The outcrop area is broadest in the west and is quite restricted in the east. The Clayton aquifer dips toward the south and southeast. The aquifer generally thickens to the south and southwest. The Clayton aquifer thins considerably east ofthe Flint River. The upper and lower contacts of the Clayton Formation are erosional in nature, therefore the thickness may vary considerably over relatively short distances (McFadden and Perriello, 1983). A facies change resulting in the silting and thinning of the sediments comprising the Clayton aquifer occurs along a line from southern Early County, across northern Miller County and into Mitchell County. This facies change occurs at the southern extent of both the Clayton aquifer and the Clayton Formation (P. Huddleston, 1995, unpublished isopach and structure-contour maps). Reported transmissivity values for the Clayton aquifer range from 300 to 13,000 ft/d (McFadden and Perriello 1-3 1983). Many of the reported transmissivity values are based on specific capacity, and should only be considered as estimates. The broad range in transmissivity of the Clayton aquifer is attnbuted to wide variation in the aquifer thickness due to the erosional nature of both the top and bottom of the aquifer. Calculated values of the storativity of the Clayton aquifer are very sparse, but range from 0.003 to 0.00003 (Stewart, 1973; McFadden and Perriello, 1983). Rainfall infiltrating the Clayton Formation outcrop area and vertical leakage from other aquifers recharge the Clayton aquifer. Rainfall infiltration to the Clayton aquifer is limited due to the aquifer's small outcrop area (estimated by McFadden and Perriello (1983) to be approximately 75 square miles), relatively low permeability of weathered Oayton Formation residuum, and steeper slope in the outcrop area relative to the other aquifers in the study. These factors combine to result in most of the rainfall running-off or evapotranspiring rather than percolating downward (McFadden and Perriello, 1983). The leakance of the confining units above and below the Clayton aquifer is large enough to allow vertical leakage into the aquifer from the overlying Claiborne aquifer and the underlying Providence aquifer, given a significant head difference between the aquifers. The vertical movement of water into the Clayton aquifer is an important source of water in areas where the Clayton aquifer is more heavily pumped. Water from the Clayton aquifer discharges to streams and rivers that cross the outcrop of the sediments that comprise the aquifer. The aquifer discharges water as vertical leakage to the overlying Claiborne aquifer, and to the underlying Providence aquifer where the head in the Clayton is greater than the head in the adjoining aquifer. Use of the Clayton aquifer for water for irrigation, public supply and industrial uses has resulted in pumpage being the most significant form of aquifer discharge. Providence Aquifer Sand units in the Providence Sand comprise the Providence aquifer. Overlying clays and silts within the lower part of the Clayton Formation and the uppermost part of the Providence Sand separate the aquifer from the overlying Clayton aquifer. Updip, this confining unit is frequently absent and the Providence aquifer merges with the Clayton aquifer to locally form the Clayton-Providence aquifer (Clarke, et al, 1983). The Providence aquifer is confined below by silts and fine sands within the lower Perote Member of the Providence Sand and the Ripley Formation. The sands that make up the Providence aquifer crop out in a band extending from Quitman County eastward to central and northern Macon County. The Providence aquifer dips to the southeast (Clarke, et al, 1983). Data on the thickness of the Providence aquifer, and other hydraulic parameters, are sparse. The greatest reported thicknesses within the area of this study are in western Dougherty County and southern Crisp County (Clarke, et al, 1983). The down-dip part of the Providence Sand becomes silty along a line extending from southern Early County, through northern Miller County and toward southern Worth County (Applin and Applin, 1944; P. Huddleston, 1995, personal communication). South of this area, the Providence Sand is not considered to be a viable aquifer. Clarke and others (1983) report transmissivity values ranging from 800 to 4,600 ft/d for wells that are screened in the Providence aquifer. Many of the reported transmissivity values were based on specific capacity and should be considered as estimates only. No estimates of the storativity are available. Studies of the Providence aquifer are hampered by the fact that few wells are constructed strictly in the Providence aquifer. The primary form of recharge to the Providence aquifer occurs as infiltration of precipitation in the outcrop of the Providence Sand. This is an important source of recharge in much of the outcrop of the Providence aquifer. Much of the water that recharges the Providence is discharged to streams and rivers that cross the outcrop of the Providence Sand. Other forms of discharge from the Providence are vertical flow to overlying or underlying aquifers and flow to wells. HYDROWGIC BUDGET The conceptual model, descnbed previously, qualitatively descnbes ground-water flow in the study area. Quantitative estimates of various components of the ground-water flow system are helpful in the development of a digital flow model. Stream-flow measurements taken in a period of drought can be used to estimate ground-water flow to streams. As noted earlier in the text, ground-water flow in the intermediate and regional flow systems is the focus of this study. Flow to streams during the drought of 1954 is assumed to be solely from the intermediate and regional flow systems. Data on the discharge to streams from the intermediate and regional flow systems provide a constraint to which the digital model can be calibrated. The stream discharge data are also 1-4 useful as a lower limit of recharge to the flow system. Minimum monthly streamflow at a number of sites for the 1954 drought was estimated by Thomson and Carter (1963). These data were utilized to estimate the discharge of water to streams from intermediate and regional ground-water flow for specific reaches of streams. The measurements upon which the estimates of minimum monthly streamflow are based generally do not coincide with the updip or downdip limits of the outcrop of the units comprising a specific aquifer (see Figure 1-1 for locations). There are a number of possible errors associated with the estimated minimum monthly flows and the discharge measurements upon which they are based (Thompson and Carter, 1963). However, for the basis of this analysis, the error associated with partial record stations (including all of the reaches in this analysis with the exceptions of those on the Flint and Chattahoochee Rivers) is assumed to be plus or minus 5 percent. Minimum monthly flows from continuous record stations are assumed to have an error of plus or minus five percent. Errors for estimates of monthly minimum flow based on basin drainage area are assumed to be plus or minus 10 percent. Minimum monthly flow rates estimated from drainage area comparisons were used for several tnbutaries to the Chattahoochee River, and for the reach of Pataula Creek between station number 788 and the Chattahoochee River (see Figure 1-1). The ground-water flow to Kinchafoonee Creek from the Providence aquifer was calculated by taking the minimum monthly flow at Preston, Georgia, at the southern margin of the Providence aquifer outcrop (Thomson and Carter map number 886, 40 cfs [1963]) and subtracting the minimum monthly flow at the upper margin of the Providence aquifer outcrop, southwest of Buena Vista (map number 885, not shown on Figure I-1, 12 cfs) and the contnbution of Slaughter Creek (map number 885A, 6.5 cfs). This results in an estimated discharge of 21.5 cfs from the intermediate and regional ground-water flow systems of the Providence aquifer. Application of the appropriate errors indicates that the discharge could be as great as 27.4 cfs (subtracting 90 percent of the upstream and tributary flows from 110 percent of the downstream flow), or as little as 15.7 cfs (subtracting 110 percent of the upstream and tnbutary flows from 90 percent of the downstream flow). Appendix A lists the estimated ground-water discharge for each stream reach simulated in the digital model along with the maximum and minimum estimates based upon application of potential errors. Data concerning drought flows in the Flint River and its tnbutaries were sparse. Most of the estimates of tnbutary flows were based on drainage areas. Therefore the drought flow estimates for the Flint River have a greater error range than the other values. The total ground-water discharge to streams from the intermediate and regional flow systems is estimated to be 509 cfs. Although some development of the system of aquifers had occurred at the time of these measurements, it is assumed that the development had little influence on ground-water discharge to streams. These estimates of ground-water discharge to streams were used as a tool in the calibration of the model in the steady-state mode, along with observed heads and boundary fluxes. Under the assumptions of the conceptual model of the flow system, this estimate represents a lower limit of the recharge to the intermediate and regional flow systems. Total recharge to the Claiborne, Clayton, and Providence aquifers will be greater, reflecting the inclusion of local flow systems. 1-5 I I I \ Taylor\ \ I I ' ;..., I ''I I' ...f ... , , , '" ' ... ... ,_ Marion ,, II ' ...... ,---..1I ' f'-~85A I I ' ... , I ... , , I Generalized Providence Aquifer Outcrop Area Stewart 85W Barbour I I I I I ( Dooly ,L ___ __,- ' I ' I I I I I I I < ' ' ' ' I I ' I I Crisp I I I -....... -"~,,N ,-- ', / ... .... __," ........ I J--J I Turner I I Tift ' / , I '/ ... 865 864') I I Sumter I I 897 I I /900 Lee I I I 8,8_7 .. .. ,< ___ _'..~/ ... I I' I I Colquitt ' I I I . I I .. ,,1II' .. , ... Terrell 911 I Mitchell I I ' ........../'"'-.. ... 912. Randolph I n --..., I I ' I.., I I I .916/ I Baker ' ' , I 1- I I I I /Grady I I I I I 1', { I II '.. I' I ' I '-, I '-, ) I ', I ' I ' ' I ... I I ' .. 1' ..,_ '-1' ' .. I 'i ) 794 I I I I I I I ,I I I I 7 I Miller < I , I I ... I 023"'4" 1 .. 2960 ', \. ' ''\ . ,... 02343500, I I I ... , , ....... , ,I, Decatur I I I I ' ' '-.._ ... , ' ... , I 1 Seminole ' ', -I 1 ' ' Henry I I "',~~!! Taylor\ \ , L ___ _,.-..__u._~.~o_~I.Lj I';.... ,,,1/ I I' !'~ ' ...., ,I ' .... 1 (. ...... .... ...., . . } -- lnlillllULIIIIIIU:IIIri~l~l I ' ... , Schley / Sumter ' ..... I I I ',, ,I I ~ I I Lee 1'.. .. ,--J Marion ',.....;_-,'..... , / / ~~ ' , ,.._.-.,., ~-~-.... , I 1'' .. 1 ', ', '<- .....' I ____ _.....,. .-.._/ I ..... .., 1 , I Webster / II .......... I ,, I I Terrell ,, ,' ', ' .... , I........... ....... , I I '' .... ,-.....-.. ... ._. - , I Stewart I 1 I Randolph 'I / I ,-- -...... -...\~ . " '1, ' , ' N ....... , /. . ' .... _,' ".1 Turner t' I Tift , ' , I <, ... _,- ........... .... ........ ' I ', I .... ... , I I I ort ..I. I I' I , I I I I Colquitt Dou hert I I I' Lt.' a er e~+~~ .......... Barbour I arly Mil r ecatur ............... \ \ , / ---- . ~.LORIDA --.., ... ':::( \ 1.. Jackson 85W, 0t====J10~==~2E:0====3l0~==~40 miles Figure 11-1. Boundary for model layer Al. All heads are specified. Patterned area represents the outcrop area of the Upper Floridan aquifer. II-13 --:~:r r ... m111:r..,..,.,.: ~ ~.t. 11 ..-..-.. .. ' ... , ... , ( (,' Macon :::: D~oo :; ~~t-:r-.- / , o o.o.~ ... . .. \ Taylor , \ I '' I! /"J't..., I =.:.:.:.:.:.:.:.:.:.:. ,'',, :.:.:.:.:.:,!.: I' ',, . ~.~~.~:~:.~::.-.:::...:::i'".-t':-:'"'~..;!:....~:.....:.t...~.:..:..:..~..:'.i~'...lI:'.lL.~.~I:,$..I.,....r.I"r..=I'r.'!Y,r..I~li~11 ' .... ... ..................-..&..............:.. ', Schey .:.,~::.Sumter..... r. ''~' I .............,t-' .... ' I .......e ...(.;. vrisp e . . . . ...... . . . ..... .._'.. ..... ....',_ .rn ... I ' ,..., I . . . . . . . . . . e . e . ~'- e e=~ ........ ..... ;(. Marion '... " I I 1'- t' ' '"'-! .....~.~....Y.~. I', I .... ... ' .. 1._..-..1 -:-...").-~- '.:~.:L.:..,::.-:.:.:.:.: .. 1I , ' I 1~-.... ...._ ' ' , , '.... , w e , b st..e.....r...:...:tJ..:,;...._~.:...:...:....:T.:e.:r.-r:e.:.ll:..:. ' ....z....:....%........................................ ...1~.......,..l.(\.,..,..:.s......................... .. ........................................... [ . .......~ . .,......~.......:..,..z. .................." .. Stewart ... J I :ftaJlPPIP!l:::::::.,.:::: ..............-,;..... II r-..,asow , I 1.1, : : : : e : e : : e : : e : : e : : e~_,. e e e e e ~ e : : Calhoun 4 ''~. .. r.r I Turner ~ , , I I Tift I _<~ ~ .. < ... _, ....... , , II ' ~MT~-J ',, I ' ... , / ..., I I I '' or~., I ~ I' b ,' .., ~ II ~I/ ~;,: 1-Colquitt L r ougherty ' ...' ..8. '4o.w..'.::, Mite el I I /' G~ \ I '""-,1' '.... ::::::::::,.:::::::: -., I ', .........4......... ' I ' ~~~~~~~~J~~~~~~~~~. .. ...............,...'.. ,, .ijilllay"' 'I ', Quitman ' ' v32N ' "'-l.. .. ,, ' I 1' ' 11 .....................r..t.4f...~.:.-..;..:...~..~.~...'.-..... _,''~.. 1 .... ',..1..-,.... .-......1... . - .... .... I' '--.. , .....;._'_... ..'... \'\~ ......... ........... ,I I J::::: ::::.:::t;.. EartV : : ::::: ::::: .~ :: ~ : : : : : : : : : ... , ... ecatur ,I, Barbour =.: ::.:.:.:.::=.:.:...:...:.-..~.;~.::.:.: ..'.... ::::::::::::Heriry .. / .... Ill (,I ' / 1 v:l I I 1 ' !IIIII ,,' Da.lE~:..!:::~:~..~::~~:~...::::..:.:::~:...::.::a...::'::~"...::"=.::....:::.=::.,..::.~=:..:..:-:.:.:=:.....::~:.::-=.:...:;:.:,J::.:..:.:~:.:.:.:.":.~:-.::..::j:..:.::.(::-:.::.;:.:~:~::...~:..:r.::~:~-,..:.~:-,J::...:~:H..:.':":1~:.1o:1.l.iuT.:s.!t:,.o;.._.~_.B~G~!rI~IV~.l.renl~~,~.8n-iAno"H1~'1.r1f,1',1>r1r'11" 1-..___ cc'F1 LORIDA ........., L Jackson ' \ ' \~'~// ..._.. '' \ " 1 1 :::::::::~"-: : : " 1> ot r- Jl 85W..._ Ot:===j10E:==~23o~===3J0~==~40 miles Figure ll-2. Outcrop, active, and specified head cells in model layer A2. Dot pattern represents the outcrop area of the Claiborne aquifer; the line pattern represents the additional active area; and the random "v" pattern represents the specified head cells. II-14 umter Marion ,, I ' t-, I ... ~ ''tI <:.........., . / ' '... , We ster I ' ;::::t I ...... , Lee Terrell Stewart ~ / 'asow / ... ,.... ,1.... \ I ',1 ' ... , ,, ' , , I ' I ' Quitman .'', ,, I / ,,,;..I..'., ..... \ J!i, 1:::1 ....... ~UTCR<)r Barbou ~AR 1r Randolph IT Henry Calhoun Early I Tift I 7 < ) "~ ' __,_<..."..."..:. ... ..... ...~ ..... r ' ... I I' I I I ' i;ol/ .L:--.Colquitt Dou herty Baker 'tch I ...................... ........ 84W~ I I /' ... Gr~ iller ~I I I I I ' Ot::==~1E:0===~2t0===:::3E0~~~~40 miles Figure ll-3. Outcrop, active, and specified head cells in model layer A3. Brick pattern represents the outcrop area of the Clayton aquifer; the line pattern represents the additional active area; and the random "v" pattern represents the specified head cells. II-15 Sumter Lee Terrell Dougherty Randol h aker Calhoun ... I I' I I . Colquitt < . ' . ....... 8'4W' Mitchell ,I..I Ia Early Decatur Henr I I I I ' I I , ' ' Dal / I Jackson I 85W O't:===1l0~==~20[:===3I0~==~40 miles Figure ll-4. Outcrop, active, and specified head cells in model layer A4. Diagonal pattern represents the outcrop area of the Providence aquifer; the line pattern represents the additional active area; and the random "v" pattern represents the specified head cells. IT-16 I ::,' =-~=~-~==;-=~=:-=2=7-o=o-=-==-3-=:-=:=:-3=::-:::=-===-====-====-====-=--=--=_-=-~_==I =-- ,, - ' Macon~' :::=.r-:::=:::=====~oooly::=:::=:::=::=-t::::::ssoo I ............/, I \ ,~: --~==f==1=4~2~Q~E=====*o~o==~~==--_ I : ~ ~,. i "' :_ ___:..g-!JW---::::-_-_ Turner -'"" :J:f /"::::;~:~-~~J~=~~~~~~~~~~~=~=~~~C:cE -, '-, Taylor\ , ~I;-.,, -.~"'.:,.. . . . . A ~66 ~q v ~: b ; : ot~ 9 4 b> .," =-=-=-=-=-=-=-=-=-=-=,=- fooo----_-_- ::!.' H '- I.,> v : .. ~-d.::o..~.'":".q.:;<.I.,.~ v... .,: ._.o... ..., --~---j-: -----------~+---------- 'v,:~ ~: ~==~~=====::::::~=~=~~===- r :<"';,".".. :..: ~---_'<-'_.'1_.'-.,, Tift I I -...I Mar::-~',V::}~i::);;rti:}y;=~~illiili~l~==--- '......... scti'le",,, su"iter. .,__-_:_-_-_-_"'1_-_-_-_-_-_-_-_-_-_-_-_-_- I I I I' I I '>-::/ ,/ ;?/!} ' tr::-~<::.:< ~ii~~~~~~i=~~;%=~~~~~~~~[fi~i~ ' / ,-..._~;3~-z~oob.-o,,b ;.,.tf;. .:4~1~0~0; " =-==-===-=~=7-o=o=-=--==--==--=-L===-~==-==:--=-2=8-o=o~~=--6=-o:-:o:0-::~:~:-=:=:=:=-=_= I 1 1 ,1 ' .... , ' ... , , 4 Wep,.s."'J-'~e.. vfe.r"r"t>e"v'll ., q " e.. :::-:::-:::-=-~-t:=~~-~~-!:.1f::::-e.--:::-:::--::::-:::-:::-:::-:::-:::- --=--=--=--=--.=r-.t=- --=--=--=--=--:-=---=--=-:=..r-_T3-i=--=--=--=---:--=- --=--=--=--- Stewart ,, I I ,, I ' ',,I ,,'' Quitman ' .... ',,,, I I ,....., > r' > b. ,J "" -...:!?..9Q:...-=---------t-_-_-....'"lflll~r_-_-_-~_-_-_ :-::.::~:?~~~{~: ~~i~~~l~~~~~~~~~I~~1~f-,. .l\,, "---- " ~~' I Pcv-/:t>:: ",. L l' 6 -=E="=-=-1-==--=-=-_==-_=-=_=-=/=_=-_~-=_=-=_=-=_=-=_=-_==-_=-=_=-=_==:-;-~-~~=~===~=~ecatur I I I I Barbour I I I I ( :., v '(' > < > ~ r., v 0 t> ~tc -_-_-_-_-_-_-_-_-_-_-_~-_-:_-~---=seminole \.A '-r->.1:"1&..,V1t_o,.,~'<'"'v>'C';O\ :~"J~A~ "q .. " > ., > < "" T" v r- ., " '- <1 "' , ' , , ' < r , , , ' ~ ~-_=T=_==- --1-Q--Q--( }=-=-=-==-==-==--=7=.~_*=_==_==_-- ---~,_,._-.y;;E-0-~~-A--L_-_-_-_-_ ', ',' , .;_; :,',<: -v"; ', , ,. \il...,.:.:::, ~"r"'>.,;H'(''e~.,n:'vr:'~y:;-,v.'"':v\~>,,~..~: " "I 4" i'~: =--=--~_=-=.=--~~~t -8~ -1i_-M_i--=-_~-_-_ ========j:"~===-=-=-=-=-=-=--#, ~ \ II 1 I 1 >t,..;...lo.lo"'"~'.".",.'(!-'".:.;,.),".":.,l..,. "'".".',,".'".,..,.,~..,."....~.,, ' " ) ,... 7 "., .,......:::. ., I. > AA V' Dale ; ,:- ,. J,.1 ~ ~ l. .l., ~ '., ..~" " ,. 1- f" "' v ") -I" ( v ..1 ..... ~ ~t~ v A qA q " Jackson I ~"1..'"'('"<..11..'-.,,.""-'.,<:.,1., VA A " ..1,.. ... ot:===11o~==~2ro~==:::3:Eo~==~4.o miles l:"/,J 3000 fWd [:=::::::3 3500 ft2/d l: :" ~ 4000 ft2/d Figure 11-S. Transmissivity for model layer A2. Black dots are estimated transmissivities from McFadden and Perriello, 1983. ll-17 I I I \ Taylor\ \ I I... 17-- ... I' ..1',' ,' ... ... ' Marion ... ', ,---./ i' .. I '-, I ' ... , I ' Stewart 85W 32oN Turner I I I Tift I _ - ' _, ..(..."...:. I ... , ... ', ' I I ' , ... , I I I I I I' I I I I I Colquitt 300 II M" ... ' Decatur "\ ' l I ...... .-...--...... FLORIDA ....., \ l .. ,~( \ '\. I I Jackson I I 85W ot:===j1o~==~2~o~==~3!o~~~~4o miles Fi,C!I 302 ft2/d IIIIID 605 ft2/d Figure ll-6. Transmissivity for model layer A3. Black dots are estimated transmissivities from McFadden and Perriello, 1983. 11-18 3400 Tift , I I ,I I I I I' I , I I olquitt Mi ,,.I Decatur .. , Seminole ',, '' ' \ I I FLOR--ID--A-......., \ ~(l I I I ................."'"\,.,.......UI / , I Jackson '"'- \ \. 85W Ot:===j10E:==~23:0~===3l0~~~~40 miles =:=a6o ft2/d 1380 ft2/d [[!]1560 ft2/d []1730 ff2/d ~1900 ff2/d ~2600 ft2/d Figure II-7. Transmissivity for model layer A4. Black dots are estimated transmissivities from Clarke and others, 1983. II-19 ow I - - I ... 8,4_' II (,'/ , I........... =:-::: =~;:.'<.., .. , .. ~- .... : .. - ...',",... I . .: { ' - ' ... , ... : ...... I ' , ... , ..... I - I I ... .. ;: ..',"". .. .. , ::::.:::.::-::':,,,Y..,. ~--',', ::-.:.,.~- I .. - ' ... I , ... - - I I .. ... I,_' ....... .. :!-'Dooly .. ',", .. l ..',"', .. :: ::~:.' "' f I ' ?!t.. }: --~~::.~-,~~:~:.:;::::_;,;~:::: ::.: :,::':;: / Macori:;':,,.,:~_.::.:",",-I I ., I 1 ,. I 1 1 "" ,.. I"' 1 1 "' (..., .., I I 1 I .. I 1 I ' ~...I ' :=:::=::::.. ,., ::-. -. ,- _'/ - Cnsp-._, , -, ~ """"" \, I N /"- ' , ' ...... , - - .... _, ' ':)'-...,/, Turner 1 Tift Marion 1 I ' ' ... ! :-j --<- . _,. .. ' ......',.,::.~;.:::.,::.:;,::,:.,=/...:.':.:,::::.:'.::..:::,.::..:.::,;..:.:.::..:::z:s::..::.:~.:~:.,}..:l.:.:....=...;;.,::::../:.:..~ ~- _1 , ~- .! I "".., , ,' ,.. 1 I ';: _'.1,.00028 . - :;.00112 "" I .-,I-'/,(':-~ I/~ ,. I t ., '511 f I ' 1 .., ... Colquitt Stewart a sow I .:::..;..;::::::::.:.;:.::::::;;::.;::)t::.Ca!houn /.-. - - I -- I .... - I- It' Grady' 32oN Barbour Jackson 85W Ot:===:J10~==~2~0====3I0~==~40 miles D l??Y] .1s .ooo11 Figure II-8. Calibrated storage for model layer A2. Black dots are reported values from McFadden and Perriello, 1983. 11-20 I I \ Taylor\ \ I I.. I ;-.. .... , I' .., 1',1 ..... ......... .... ...... ........ , ..... Marion I I Turner I I I Tift I - __, ..(....~, .......... ..... ',, I ', I ..... I ',I I ..I. I I' I I I Stewart I .. , I ''II \ Quitman ,,' ' ' ,, . , ... I ,......, I '.I I I Jackson I I miles t:===::J~===~====::E~===~40 Figure ll-9. Calibrated storage for model layer A3. Black dots are reported values from McFadden and Perriello, 1983. II-21 Tift ... ............ ~ ' I '.._ I .. , I I I I I' I I I '\ ' I ,,l \ I I Jackson I ot:===~1Eo=:=:=:~2ro~===a1o~===~4o. miles ~0.15 E):.wr1 o.ooo46 Figure ll-10. Calibrated storage for model layer A4. II-22 I \ \ Taylor\ \ I I... ;.... I ',,, I' .., ,1',' ..... ... '..., . ' I I ..... I Marion '"-(,....... , I ,, I .. ', ,---.1I ' l' .. I ..... I ',, I ' , ....... ..... I Turner ......... , }'../ , / I Tift I I' I , I I Stewart , / , I ' , ,, I ' ,1,I, .. , ,,I .. , ' Quitman 1 ' I ' ',, I ,I ,--- ',,.., I I ........ I\,, Barbour I I I I I I I (I.., . , --- I I I I I I I I , I I Jackson Ot:===j10~==~2~0====3I0~==~40 miles ' \ I I \ asw ~ ~ 2 3-12 sec-1 ~ 2.9-12 sec-1 Figure ll-11. Calibrated leakance for model layer C2. ~ ~5 a-13 sec-1 ll-23 ,__8r4-'W\/I I I I Macon \ \ Taylor\ \ I 1... ;-. I ', ,,I' ,; .... ... , ' ... , ... Schley '...,, ' Marion Webster I ................ ;' ' ....I... ' I I Turner I I Tift I ... I I' I I Stewart I as ow , ' ... I I I, /'',.../' ... ., ,, ' I ' I ' 32oN ' Quitman ' I \' I ' ',,, I I I ,---~ I \ ' Randol Barbour Henr , I I I I I I Dale Jackson Ot:===j10~==~2t0~==~3I0~==~40 miles ~ 5.8-11 sec -1 M\\\)J.i/i)1.2-11 sec -1 Figure ll-12. Calibrated leakance for model layer C3. II-24 \ I 1.. 85W 70 60 50 (i) ~ w Cl ~ 40 :::! :a::.>. 30 20 10 0 2 3 4 5 6 7 8 9 10 11 12 13 14 15 STRESS PERIOD Figure ll-13. Simulated irrigation pumpage. II-25 100 90 80 70 ! 60 i 50 ~ 40 30 20 10 0 A2 2 3 4 5 6 7 8 9 10 11 12 13 14 15 STRESS PERIOD Figure ll-14. Simulated pumpage for all uses. ll-26 I Macon Tift I Stewart 32N Jackson a sow Ot:===j10~==~2l0====3I0~==~40. miles R 0 to 0.5 ~ 0.51 to 1.0 1.~1 to 1.5 ~ Greater than 1.5 Figure 11-15. Simulated pumpage from model layer A2 in stress period 15. Pumpage is in cubic feet per second. TI-27 84W I . .__ ,-,1,/ ( ~~ ...,,/ lliill Macon I Dooly / ,, I / ',,(1 / . ' I I .,,,~ - \ N /.., ' ',,,.....,~,~-- ... '. \' flill IITill I 1 Crisp I 1 Turner 1 Taylor\ '. ,/ / / Tift } ,L-----'-..._ __!_ /.. - <-:. I / - ... -... ,./ ,, __, ........ ' ' 1 ' 1' ', ,,... ', I ............... ;..., I I ' ... , (' I ' 11 ,, '" 'fil ,1'..._1 .9 1'. / I ..._, I \ / ' , ' , 1I ' ... ' ... , ... , , 1 I 1 ' ... , Schley I I '~, /I ~ ... , I I Marion ' ... ,{ S~er ~ 1 ,/ I 1 \ ' - Worth \ Lee ' - , '~' ... ' mm IWJJ Jrrn ' , , , 1 I ... I 11 1 I 1 = I 1 ' '..._,..!.._.-.1 I' I I ... , ' , , I .... II I'- ' '-, "~o"~' "~"~' / - ' , '~d-...:_'-.-.-',~/' Webster ~~ ~ ~ ~ -, 1.:1.. 3'~ 1"=1 ~ ~ ,~ ... / ,, I m'mTrn ~/ ~ i l l Doughe r t y "'-' ' I'' I ... , / I ' /\ I / 1 Colquitt I '-.., '-, ' ... '-, I 1 '... ' ... ,J]]]l / I;ZJJ ~errIe'ZlZJl I / lll!l I / \ ow' , 84 ' ",/I' ..., , E3 lI1T1i1iJ. / , I' Mitchell Stewart -,___ / 1/ Rando'l ~~,],,...~-' mm - ... ~.// ; \ I / IIIID / Dlll 1 Baker ' I I 85W / I, 1 ~ ~/ Ca~ """"' 1 / I , /Grady I I '\ '' I ' ' II '-,1 .......... t::;;;;;1 I I II fiTTTT1i1l / ( I "', I '... 1 IWllUJJ I ' ... I' I I ' '' , , '1 ' Quitman ', I /'-... 1 '... I ', 1~ ~m, ,.7~' I I llJZI I , , -.. ...... ' ,, I '1 ~ I ,,,,,...I..,.., .... 5f""3!a"!y 1/ I .....,.,._., I mm '- ' ... ... I I ',II ', ) '- ' I .......... 'II I ') I .,,I I I Miller '/ I I I < I ' '-,..._ I \ ,' I I .. ~, I ' I.., \. I~ W.I ~I ar y I I I ....,, I Decatur I \ . ,..._ I .......... II / Barbour 1 ([]) ' .., ' -.. ,.. ', ' II s . , .......... ... , ..,I1 em1nole ',' I I \ I ' <,, .... , /''..,, Henry / '' GEORGIA/ I ..... , . .......... _ /1 -ALABAMA:~_.-... -...... \ I \ l / ', 1 ' I / Houston I .......... I I FLORIDA ....,_ ::::t, 1 - , I \. I / Dale I ' ',...'' I II I Jackson I Ot:===j10~==~23:0~===3I0~==~4.0 miles IIliJI 0 to 0.5 ~0.51 to 1.0 ~ 1.01 to 1.5 ~Greater than 1.5 Figure ll-16. Simulated pu"u.page from model layer A3 in stress period 15. Pumpage is in cubic feet per second. II-28 I I ' \ Taylor\ ; I ................ , ' I ... ...... ... I ' I I / Macon ., IJ]]J Dooly [IIl] 1 ' I ', I , <' llll I ' \ / I I Crisp I 1 Turner I I \ / llJ] L ___ ...... -__ / mL ,/' ------ .. , I / ,..,/J'', I I - <'...._ ~-- "-J ' , , Tift ;:-,, /, [I]]/WJ_/ ITIIll I' 'I ',/ '-, /'...._,/ I \ '.._, I ',1 ' 1 '-.., I ,,,, ',, ', / S~ley / 1 ' ..., I I ' ,.. I 1 Marion ' '-..!I :........... .~..Sumter I 1 / I I I I ,I \ : -\ Worth \ Lee \ - '---:--~ '; ' ~ I ' ... , , ....... ' , /1 / , I I , I I I I / ' .. , l:ttiiJ , - - ' t'... I ' ' .... I ' ...i. ',..,..--., IIII ' ', __ _..' / '<._: ..".. ~ ',', ~~ lllW~ ill.ll I ' / ,1'',..' , '1 Colquitt 11 ' I ... , / ' ...., Webster ,/ / liD'' .. , , / Te~ll '-../ Dougherty /\ / / l .......,, ' .. , I I / ) 84oW:::. llll Stewart ',, 1'', I ,, ......... I I ........ ,-...._ I' I I I 1 / ........ Randolph ..._ I ............ -, ' I 1 I / / I ' I I I ) --.., I '\ \I / Baker , Mitchell I I 1 .. I / I i'a5W ll 1 I } mill / Calhoun / I \ / Grady "\ I' 1 ' , / ',1 , 1 I / ( j / ,I .. , I ', I' 1', I, I ',_ '\/ \ Q u1'tman ' ' ' ,, 'I I ,II' ,_ ,II'.. , , I - I ''-,.,..., / ,I ....... I I 'I I --- I ' I .. ' , , I 'I ? I ' ) ' .. ' , ' I ' 'lI I - ' , ,, -_ ' ..,,, / Clay / (,-.......--mm.. / I \ ..... I \, .., , I Early ;1 , I 1 Miller / , I I ,, Decatur I I \. ... I ...,, I / ' , ......... _ / ,., ', I - ....... Barbour ,1 I I ' ' .. ,._ I \ ' .. ,_"",1, I Seminole ' ' ' I \ (-,, ' Henry / 1', / .... ,.... I I I ........ ' -,, I I 1 I I I Dale I ' I ', ..I, I / "'' GEORGIA/ -ALABAMA'~ .._ ,..--..-.~_ - '.I l / Houston I , I - --... -... FLORIDA ... ....... .... _ \ .. .. ::::~,I ' I I I ~ I ' Jackson /, 85W, ot:===j1o~==~23:o~==~a!o~~~~4o. miles llillJ ~ ~ 0 to 0.5 0.51 to 1.0 Greater than 1.0 Figure 11-17. Simulated pumpage from model layer A4 in stress period 15. Pumpage is in cubic feet per second. II-29 I T,ayK, lr \ , \ I 0 / ,' I' / ..... Marion i Stewart ~\ , , "' I 'I '\ Quitman ' ' '"\. i I '\., I ' __ .... ,\t ~...... _ , ........,. I I Barbour I "Turner I I / / I ,. I I Ot:===:J10~==~2~0====3I0~=~~40 miles ---250-- Elevation of potentiometric surface in feet above mean sea level. Figure ll-18. Steady-state simulated potentiometric surface for model layer A2. 11-30 \ I ;' Stevvart h / Quitrnan ,~--..,.__... / I /l I / / Ot:===~1E0~==~23:0====:J30~===~40 miles ---400-- Elevation of potentiometric surface in feet above mean sea level. Figure 11-19. Steady-state simulated potentiometric surface for model layer A3. 11-31 ot:===j10~==~23:0~==~3I0~~~~40. miles - 1 5 0 - - Elevation of potentiometric surface in feet above mean sea level. Figure ll-20. Steady-state simulated potentiometric surface for model layer A4. II-32 Net flux = 0 cfs A2 Note: where two values are present, estimated is on the left and simulated is on the right. Where one value is present, this is simulated. Figure 21b Net flux = 0 cfs Figure ll-21. Estimated and simulated fluxes for model layer A2, model layer A3, and model layer A4. ll-33 A2 2.3 1.6 "' I' + t 2.1 11.0 Dip- "'2.1 I' 3.3 Recharge = 339 River Flux = 320.1 A3 0.1 0.4 II' I' + t 6.5 6.9 Dip- II' I' 2.3 0 Recharge = 64 River Flux = 29.8 A4 2.2 0.5 II' I' + t 11.7 Dip- "5.9 I' 0.4 Recharge = 164 River Flux = 186.3 Figure ll-22. Steady-state mode boundary Ouxes for recharge, discharge, and leakance. All values are in cubic feet per second. ll-34 84W / Mitche!i ().t;!trnan 32N / I I ' ' 85W Ot:===j10~==~2~0====3l0~~~~40 miles - - 2 0 0 - - Elevation of potentiometric surface in feet above mean sea level. Figure ll-23. Simulated potentiometric surface for model layer A2 at the end of stress period lS. ll-35 ~ N Stewart 84 Macon /300 ......... Crisp- 250 ,(t; __ _ J------ - 'u-., ~ f ) I --, l I -----Miller - -t " )' '-...Worth ,150 ~----"j - - - j - -, ,...J - - - - 1 I .J'"" ,r-./ / ) Mitchell 100 0 10 20 30 40 50 kilometers ~==~==~~==~~====J~~~~ J. Potentiometric surface contour in feet above mean sea level. - - - 1 0 0 - - Dashed where approximated. Arrows indicate ground-water J. flow direction. Figure ll-24. 1986 potentiometric surface for the Claiborne aquifer (model layer A2). II-36 8 4-"W'~. \ \ ./ T:ft I Stevvart 85W 32N ! Baker / / / ' 85W OC:===:J10~==~2~0====3l0~==~40 miles ---450- Elevation of potentiometric surface in feet above mean sea level. Figure 11-25. Simulated potentiometric surface for model layer A3 at the end of stress period 15. II-37 ~ N Stewart Macon I r--- -..r.._~-v--j__ -: I Miller r I I I I j --1 I Worth Mitchell OE=:=:=:::13:0====2}0=:=:::::=-=::::--:J..,.~o========j~o miles 0 10 20 30 40 50 kilometers ~==~==~~==~~====J~~~~ ~ Potentiometric surface contour in feet above mean sea level. + - - - 1 0 0 - - Dashed where approximated. Arrows indicate ground-water flow direction. Figure 11-26. 1986 potentiometric surface for the Clayton aquifer (model layer A3). 11-38 84W ...... / Mitchel! ! / I / " Baher "' " .. ., ' ' ., I '! I ( ' ,, /'' ' I ..., ... ,_ / './ ' / I ? fv1i!k~r ' / , ___ i I I ' ' / ' I ' ' I \ Jackson \ 85W Ot:===j10E=:=:~2~0====3l0~~~~40 miles ---450--- Elevation of potentiometric surface in feet above mean sea level. Figure ll-27. Simulated potentiometric surface for model layer A4 at the end of stress period 15. ll-39 - - 250- - Water-Level Contour--Shows altitude at which water level , - ' \ would have stood in tightly cased wells. Dashed where approximately _______ __ located. Contour interval 50 feet. Datum is sea level ,,' , / \ ",~ _II I f' I t' V?, ,' Crawford r' '\_....._,\ I ./ . . ,,-..I .,___)I ,.,._ I I Taylor i I ' Bibb / ' /, ' "' "/' ' "',... I I I I !I'- _.,_- ,) .....1 - , ' -41- .-- '...( I ' I { t Wilcox I \ Crisp I I "If-DO c " I __.~.......... _ I ~ -/-,------ -'\ - - - - _ ;IJ2so~ '\1 .... _ 1/ -----7 I ,,/ \ \ l- ''"" ".", oo / ---r1 - - (---' Ben Hill 1 ' l 'r 1 Turner ' 1-' I --... ,- I ~----- J' -, Irwin ______ I I \ ..... _, I , r...1 I .r I__ -----., J- _}.----- ...r- I - I I \ I I \ 1 / Mitchell I ( ( .,., Colquitt / 1.. Miller /, \ I I (- I r 1 Area ' , - - - -I - - - - - - - - r - - -I -j - - - - - - - - - - - - - - - i I I _ - - - - - 1 I - - - - - - I 1 1 0 10 20 30 40 miles c:::===r:::===r::===::r=:==:J Figure ll-28. 1986 potentiometric surface for the Providence aquifer (model layer A4). II-40 230 --- ---- 220 --- --- --- --- 210 200 190 --------simulated - - - - - - observed 180 170 water level in 1960 feet above NGVD 1962 1964 1966 1968 1970 1972 1974 Year 1976 1978 1980 1982 1984 1986 a 170 160 150 --- --- - _______ -- ... ..... ..... .... .... ..... ..... ..... .._ ---- 140 ' '\ 130 120 110 --------simulated observed 100 ' 90 80 1960 1962 1964 1966 1968 1970 1972 1974 1976 Year 1978 1980 1982 1984 1986 b Figure 11-29. Simulated and observed heads at (a) USGS test well 4, Claiborne aquifer, Doughterty County and (b) USGS test well2, Claiborne aquifer, Dougherty County. 11-41 250-r-------------------------------------------------------------------------, 240 -------- -- --- ---- -- .... --- 230 220 \ 210 \ ' \ - \ /--, \ / \/ 200 --------simulated 190 - - - - - - observed 180 1960 Water level in feet above NGVD 1964 1966 1968 1970 1972 1974 Year 1976 1978 1980 1982 1984 1986 a 340-r----------------------------------------------------------------------------~ 330 320 310 300 ----- ---- ------ ----------- -- \ \ 290 \ --------simulated ' ' 280 - - - - - - observed ' ' 270 1960 1966 1968 1972 1974 Year 1976 1978 1980 1982 1984 1986 b Figure ll-30. Simulated and observed heads at (a) W.H. Fryer well, Claiborne aquifer, Lee County and (b) Cuthbert well, Clayton aquifer, Randolph County. II-42 330 320 310 300 290 280 270 260 250 240 230 220 210 200 1960 Water level in feet above NGVD 250 --- -- ---- ---- -- --- .... .... .... .... ' ' ' ' ' ... ... - - - - - - - - simulated observed 1962 1964 1966 1968 1970 1972 1974 1976 1978 1980 Year 1982 1984 1986 a 240 ... ... ..... ..... 230 \ 220 \ \ \ 210 '-"' \ \ ... .... ..... " "" \ \ 200 \ \ --------simulated 190 observed 180 170 1960 1962 1964 1966 1968 1970 1972 1974 Year 1976 1978 1980 1982 1984 1986 b Figure 11-31. Simulated and observed heads at (a) Don Foster well, Clayton aquifer, Terrell County and (b) H.T. McLendon #1 well, Clayton aquifer, Calhoun County. II-43 100 -r--------------------~,-----------------------------------------------, 180 ' ..... ..... 170 160 150 140 ' '''''\ \ 130 \ 120 110 100 90 80 --------simulated - - - - - - - observed 70 60 50 -+----~--~~---r----~----~---r----~----r---~----~----r---~----~~ 1960 1962 1964 1966 1968 1970 1972 1974 1976 1978 1980 1982 1984 1986 Water level in Year a feet above NGVD 170 -r--------------------------------------------------------------------------, 160 150 140 130 120 ------ ' ' ' ' ' ' ' ' ' ' ,_ ........... ...._ ~ 110 ' 100 90 80 --------simulated 70 - - - - - - - observed 60 50 40 1960 1962 1964 1966 1968 1970 1972 1974 Year 1976 1978 1980 1982 1984 1986 b Figure 11-32. Simulated and observed heads at (a) USGS test well9, Clayton aquifer, Lee County and (b) Turner City well 2, Clayton aquifer, Dougherty County. II-44 200 -r------------------------------------------------------------------~ ... .... .... 190 0180 > " z g! 170 0 .Q as ~ 160 .s: Gi ~ 150 .... ' ' ' ' ''''' ' ' ' '' ' - - - - - - Simulated observed 130 1960 1978 1980 1982 1984 1986 1600FS-r----------------------------------------------------------------------~a 0: total pumpage ,0. S: rate of release of 140 water from storage 120 100 80 60 40 20 0~------------------------------------------------------------------~ SP1 SP4 SP8 SP9 SP10 SP11 SP12 SP13 SP14 SP15 Figure ll-33. (a) Simulated and observed beads at USGS test well 12, Clayton aquifer, b Dougherty County and (b) simulated release of water from storage from 15 stress periods. ll-45 Figure ll-34. (a) Simulated ground-water discharge to rivers for selected stress periods and (b) net vertical flux across each confining unit. ll-46 'c0: 0 0 (]) Ul ~ 20 Q. Qj .! 0 :0 ::3 0 15 .s: x ~ iij "E 10 0 N " - J: 5 0 SP2 SP8 SP9 SP10 SP11 A2 in SP13 A4 in A4 out A3 out A3 in SP15 Figure ll-35. Simulated horizontal fluxes in (x) and out (o) for each aquifer and selected stress periods. II-47 Quantity: 250 Cost: $1,893.00 The Department ofNatural Resources is an equal opportunity employer and offers all persons the opportunity to compete and participate in each area of DNR employment regardless of race, color, religion, national origin, age, handicap, or other non-merit factors