Search USGS⌕ Search

Geology topics

K. J. Halford

Publications and source records attributed to K. J. Halford.

At least 19 recordsLinked to original sources

Detecting drawdowns masked by environmental stresses with water-level models

Detecting and quantifying small drawdown at observation wells distant from the pumping well greatly expands the characterized aquifer volume. However, this detection is often obscured by water level fluctuations such as barometric and tidal effects. A reliable analytical approach for distinguishing drawdown from nonpumping water-level fluctuations is presented and tested here. Drawdown is distinguished by analytically simulating all pumping and nonpumping water-level stresses simultaneously during the period of record. Pumping signals are generated with Theis models, where the pumping schedule is translated into water-level change with the Theis solution. This approach closely matched drawdowns simulated with a complex three-dimensional, hypothetical model and reasonably estimated drawdowns from an aquifer test conducted in a complex hydrogeologic system. Pumping-induced changes generated with a numerical model and analytical Theis model agreed (RMS as low as 0.007 m) in cases where pumping signals traveled more than 1 km across confining units and fault structures. Maximum drawdowns of about 0.05 m were analytically estimated from field investigations where environmental fluctuations approached 0.2 m during the analysis period.

Ground Water↗

A computer program for flow-log analysis of single holes (FLASH)

A new computer program, FLASH (Flow-Log Analysis of Single Holes), is presented for the analysis of borehole vertical flow logs. The code is based on an analytical solution for steady-state multilayer radial flow to a borehole. The code includes options for (1) discrete fractures and (2) multilayer aquifers. Given vertical flow profiles collected under both ambient and stressed (pumping or injection) conditions, the user can estimate fracture (or layer) transmissivities and far-field hydraulic heads. FLASH is coded in Microsoft Excel with Visual Basic for Applications routines. The code supports manual and automated model calibration.

Ground Water↗

A refined characterization of the alluvial geology of yucca flat and its effect on bulk hydraulic conductivity

In Yucca Flat, on the Nevada National Security Site in southern Nevada, the migration of radionuclides from tests located in the alluvial deposits into the Paleozoic carbonate aquifer involves passage through a thick, heterogeneous section of late Tertiary and Quaternary alluvial sediments. An understanding of the lateral and vertical changes in the material properties of the alluvial sediments will aid in the further development of the hydrogeologic framework and the delineation of hydrostratigraphic units and hydraulic properties required for simulating groundwater flow in the Yucca Flat area. Previously published geologic models for the alluvial sediments within Yucca Flat are based on extensive examination and categorization of drill-hole data, combined with a simple, data-driven interpolation scheme. The U.S. Geological Survey, in collaboration with Stanford University, is researching improvements to the modeling of the alluvial section, incorporating prior knowledge of geologic structure into the interpolation method and estimating the uncertainty of the modeled hydrogeologic units.

Open-File Report↗

A unique approach to estimating lateral anisotropy in complex geohydrologic environments

Aquifers in fractured rock or karstic settings are likely to have anisotropic transmissivity distributions. Aquifer tests that are performed in these settings also are frequently affected by leakage from adjacent confining units. Finite-difference models such as MODFLOW are convenient tools for estimating the hydraulic characteristics of the stressed aquifer and adjacent confining units but are poor tools for the estimation of lateral anisotropy. This limitation of finite-difference methods can be overcome by application of the spin method, a technique whereby the positions of the observation wells are rotated about the production well to estimate anisotropy and orientation. Formal parameter estimation is necessary to analyze aquifer tests because of the number of parameters that are estimated. As a test, transmissivity, anisotropy, and orientation were successfully estimated for a simple hypothetical problem with known properties. The technique also was applied to estimate hydraulic properties of the Santee Limestone/Black Mingo (SL/BM) aquifer and a leaky confining unit beneath Charleston, South Carolina. A 9-day aquifer test with an average discharge of 6441/min was analyzed numerically. Drawdowns in the SL/BM aquifer and confining unit were simulated with a 12-layer MODFLOW model that was discretized into 81 rows of 81 columns. Simulated drawdowns at seven observation wells that ranged from 23 to 2700 m from the production well were matched to measured drawdowns. Transmissivity estimated along the minor axis ranged from 10 to 15 m 2 /day and along the major axis ranged from 80 to 10Om 2 /day. The major axis of transmissivity was oriented along compass heading 116° (degrees clockwise from north), which agrees with geologic interpretations. Vertical hydraulic conductivity and specific storage estimates for the overlying confining unit were 4 x 10 -5 m/day and 2 x 10 -4 1/m, respectively.

South Carolina↗

In situ arsenic removal in an alkaline clastic aquifer

In situ removal of As from ground water used for water supply has been accomplished elsewhere in circum-neutral ground water containing high dissolved Fe(II) concentrations. The objective of this study was to evaluate in situ As ground-water treatment approaches in alkaline ground-water (pH > 8) that contains low dissolved Fe (<a few tens of μg/L). The low dissolved Fe content limits development of significant Fe-oxide and the high-pH limits As adsorption onto Fe-oxide. The chemistries of ground water in the two aquifers studied are similar except for the inorganic As species. Although total inorganic As concentrations were similar, one aquifer has dominantly aqueous As(III) and the other has mostly As(V). Dissolved O 2 , Fe(II), and HCl were added to water and injected into the two aquifers to form Fe-oxide and lower the pH to remove As. Cycles of injection and withdrawal involved varying Fe(II) concentrations in the injectate. The As concentrations in water withdrawn from the two aquifers were as low as 1 and 6 μg/L, with greater As removal from the aquifer containing As(V). However, Fe and Mn concentrations increased to levels greater than US drinking water standards during some of the withdrawal periods. A balance between As removal and maintenance of low Fe and Mn concentrations may be a design consideration if this approach is used for public-supply systems. The ability to lower As concentrations in situ in high-pH ground water should have broad applicability because similar high-As ground water is present in many parts of the world.

Applied Geochemistry↗

Interpretation of transmissivity estimates from single-well pumping aquifer tests

Interpretation of single-well tests with the Cooper-Jacob method remains more reasonable than most alternatives. Drawdowns from 628 simulated single-well tests where transmissivity was specified were interpreted with the Cooper-Jacob straight-line method to estimate transmissivity. Error and bias as a function of vertical anisotropy, partial penetration, specific yield, and interpretive technique were investigated for transmissivities that ranged from 10 to 10,000 m2/d. Cooper-Jacob transmissivity estimates in confined aquifers were affected minimally by partial penetration, vertical anisotropy, or analyst. Cooper-Jacob transmissivity estimates of simulated unconfined aquifers averaged twice the known values. Transmissivity estimates of unconfined aquifers were not improved by interpreting results with an unconfined aquifer solution. Judicious interpretation of late-time data consistently improved estimates where transmissivity exceeded 250 m2/d in unconfined aquifers. ?? 2006 National Ground Water Association.

Ground Water↗

Estimating hydraulic properties using a moving-model approach and multiple aquifer tests

A new method was developed for characterizing geohydrologic columns that extended >600 m deep at sites with as many as six discrete aquifers. This method was applied at 12 sites within the Southwest Florida Water Management District. Sites typically were equipped with multiple production wells, one for each aquifer and one or more observation wells per aquifer. The average hydraulic properties of the aquifers and confining units within radii of 30 to >300 m were characterized at each site. Aquifers were pumped individually and water levels were monitored in stressed and adjacent aquifers during each pumping event. Drawdowns at a site were interpreted using a radial numerical model that extended from land surface to the base of the geohydrologic column and simulated all pumping events. Conceptually, the radial model moves between stress periods and recenters on the production well during each test. Hydraulic conductivity was assumed homogeneous and isotropic within each aquifer and confining unit. Hydraulic property estimates for all of the aquifers and confining units were consistent and reasonable because results from multiple aquifers and pumping events were analyzed simultaneously. Copyright ?? 2005 National Ground Water Association.

Ground Water↗

Effects of unsaturated zone on aquifer test analysis in a shallow-aquifer system

A comparison between two hypothetical flow models of an unconfined aquifer, one saturated and the other variably saturated, indicates that the variably saturated model which explicitly models drainage from the unsaturated zone provides a better conceptual framework for analyzing unconfined aquifer test data and better estimates of the lateral and vertical hydraulic conductivity in fine-grained sands. Explicitly accounting for multiple aquifers, well-bore storage, and the effects of delayed drainage from the unsaturated zone increases confidence in aquifer property estimates by removing some assumptions and allowing for the inclusion of early time data and water-table observations in an aquifer test analysis. The inclusion of the unsaturated zone expands the number of parameters to be estimated, but reasonable estimates of lateral and vertical hydraulic conductivity and specific storage of the unconfined aquifer can be obtained. For the cases examined, only the van Genuchten parameter a needed to be determined by the test, because the parameters n and 9r had a minimal effect on the estimates of hydraulic conductivities, and literature values could be used for these parameters. Estimates of lateral and vertical hydraulic conductivity using MODFLOW were not as good as the VS2DT based estimates and differed from the known values by as much as 30 percent.

Groundwater↗

Effects of steady-state assumption on hydraulic conductivity and recharge estimates in a surficial aquifer system

The ability of a calibrated flow model to predict the behavior of a surficial aquifer system is governed by the quality of the hydraulic conductivity and recharge estimates used. Reasonable lateral and vertical hydraulic conductivities can be estimated by steady-state simulations driven by effective recharge rates that approximate the net effects of evapotranspiration, and water released from storage during periods of recession. Results from a hypothetical, transient, cross-sectional model indicated that most of the water was contributed uniformly from storage from five to 25 days after a recharge event. Results also showed that a steady-state, snapshot calibration approach can be used on aquifers in a humid climate with diffusivities between 20 and 500 m 2 /d. Most estimates of the lateral and vertical hydraulic conductivities of the hypothetical aquifer system were within 30% of the actual values. Estimates of hydraulic conductivity from the transient cases were similar to those from the snapshot calibration cases. The long-term recharge rate could be identified by calibrating to multiple synoptic surveys that were sampled over the range of drier to wetter conditions. The effective recharge rates estimated for the driest and wettest conditions bracketed the long-term recharge rate. Results suggested that the effective recharge rate estimated for the synoptic survey with the lowest water level root-mean-square (RMS) error was the best estimate of the long-term recharge rate. A field application of the snapshot calibration approach simulated the surficial aquifer system beneath Cecil Field Naval Air Station well and provided reasonable estimates of the long-term recharge rate (0.4 mm/d) relative to the range of recharge rates that were independently estimated by the chloride concentration ratio method (0.2 to 0.6 mm/d).

Groundwater↗

Problems associated with estimating ground water discharge and recharge from stream-discharge records

Ground water discharge and recharge frequently have been estimated with hydrograph-separation techniques, but the critical assumptions of the techniques have not been investigated. The critical assumptions are that the hydraulic characteristics of the contributing aquifer (recession index) can be estimated from stream-discharge records; that periods of exclusively ground water discharge can be reliably identified; and that stream-discharge peaks approximate the magnitude and tinting of recharge events. The first assumption was tested by estimating the recession index from st earn-discharge hydrographs, ground water hydrographs, and hydraulic diffusivity estimates from aquifer tests in basins throughout the eastern United States and Montana. The recession index frequently could not be estimated reliably from stream-discharge records alone because many of the estimates of the recession index were greater than 1000 days. The ratio of stream discharge during baseflow periods was two to 36 times greater than the maximum expected range of ground water discharge at 12 of the 13 field sites. The identification of the ground water component of stream-discharge records was ambiguous because drainage from bank-storage, wetlands, surface water bodies, soils, and snowpacks frequently exceeded ground water discharge and also decreased exponentially during recession periods. The timing and magnitude of recharge events could not be ascertained from stream-discharge records at any of the sites investigated because recharge events were not directly correlated with stream peaks. When used alone, the recession-curve-displacement method and other hydrograph-separation techniques are poor tools for estimating ground water discharge or recharge because the major assumptions of the methods are commonly and grossly violated. Multiple, alternative methods of estimating ground water discharge and recharge should be used because of the uncertainty associated with any one technique.

Ground Water↗

Plant-based plume-scale mapping of tritium contamination in desert soils

Plant-based techniques were tested for field-scale evaluation of tritium contamination adjacent to a low-level radioactive waste (LLRW) facility in the Amargosa Desert, Nevada. Objectives were to (i) characterize and map the spatial variability of tritium in plant water, (ii) develop empirical relations to predict and map subsurface contamination from plant-water concentrations, and (iii) gain insight into tritium migration pathways and processes. Plant sampling [creosote bush, Larrea tridentata (Sessé & Moc. ex DC.) Coville] required one-fifth the time of soil water vapor sampling. Plant concentrations were spatially correlated to a separation distance of 380 m; measurement uncertainty accounted for <0.1% of the total variability in the data. Regression equations based on plant tritium explained 96 and 90% of the variation in root-zone and sub-root-zone soil water vapor concentrations, respectively. The equations were combined with kriged plant-water concentrations to map subsurface contamination. Mapping showed preferential lateral movement of tritium through a dry, coarse-textured layer beneath the root zone, with concurrent upward movement through the root zone. Analysis of subsurface fluxes along a transect perpendicular to the LLRW facility showed that upward diffusive-vapor transport dominates other transport modes beneath native vegetation. Downward advective-liquid transport dominates at one endpoint of the transect, beneath a devegetated road immediately adjacent to the facility. To our knowledge, this study is the first to document large-scale subsurface vapor-phase tritium migration from a LLRW facility. Plant-based methods provide a noninvasive, cost-effective approach to mapping subsurface tritium migration in desert areas.

Vadose Zone Journal↗

Hydrogeology and simulation of ground-water flow in the aquifers underlying Belvidere, Illinois

The U.S. Geological Survey investigated the ground-water-flow system and distribution of contaminants in the vicinity of Belvidere, Illinois, during 1992?2000. The study included the compilation, collection, and analyses of hydrogeologic and water-quality data and simulation of the ground-water-flow system. Hydrogeologic data include lithologic, stratigraphic, geophysical, hydraulic-property, water-level, ground-water withdrawal, and streamflow data. Water-quality data include analyses of water samples primarily for volatile organic compounds (VOC?s) and selectively for tritium and inorganic constituents. Data were collected from about 250 wells and 21 surfacewater sites. These data were used (1) to describe the hydrogeologic framework of the ground-waterflow system, preferential pathways and directions of ground-water movement and contaminant distribution, ground-water/surface-water relations, and the water budget and (2) to develop and calibrate the ground-water-flow model. The glacial drift (sand and gravel with some clay) and Galena-Platteville (fractured dolomite) aquifers and the sandstone aquifers of the Cambrian-Ordovician aquifer system compose the ground-water-flow system underlying Belvidere and vicinity. The Glenwood confining unit separates the Galena-Platteville aquifer from the underlying sandstone aquifers. The Galena- Platteville aquifer and confining unit may be absent in parts of the Troy Bedrock Valley, about 1.5 miles west of Belvidere. Throughout the study area, the Kishwaukee River and its tributaries seem to be gaining flow from shallow ground-water discharge. Potentiometric levels in the glacial drift and Galena-Platteville aquifers range from about 900 feet above sea level in the upland areas to 740 feet along the Kishwaukee River. Estimated horizontal hydraulic conductivity of the glacial drift aquifer ranges from about 0.13 to 280 feet per day. The Galena-Platteville aquifer is a dual-porosity unit with the greatest percentage of flow through fractures and bedding-plane partings. Estimated horizontal hydraulic conductivity ranges from about 0.005 to 2,500 feet per day. Estimated horizontal hydraulic conductivity of the St. Peter aquifer (the uppermost sandstone aquifer of the Cambrian-Ordovician aquifer system ranges from about 4.7 to 17.5 feet per day. Volatile organic compounds have been detected in all aquifers underlying Belvidere. Trichloroethene and tetrachloroethene are the principal VOC?s detected at concentrations above regulatory levels, with the largest number of detections and highest concentrations in the glacial drift aquifer. VOC?s generally are not detected in the glacial drift aquifer farther than 1,000 feet from known or potential source areas (industrial or disposal sites), because most source areas are near the Kishwaukee River, where shallow ground water discharges. Across most of the study area, the Glenwood confining unit seems to restrict downward movement of VOC?s into the underlying St. Peter aquifer; in the immediate vicinity of Belvidere, downward movement also seems restricted by lateral movement toward the municipal wells through permeable intervals in the 2 Hydrogeology and Simulation of Ground-Water Flow in the Aquifers Underlying Belvidere, Illinois Galena-Platteville aquifer. Fractures and (or) unused wells that may penetrate the confining unit seem to provide local pathways for limited movement of VOC?s to the sandstone aquifers. At least two municipal wells seem to intercept the bedding-plane partings at about 525 and 485 feet above sea level. Water levels in the lower one-third of the Galena-Platteville aquifer rapidly respond to withdrawals at these wells. The ground-water-flow system underlying Belvidere was simulated to test the conceptual model of the system. The three-dimensional, steady-state model represents the glacial drift, Galena-Platteville, and sandstone aquifers sep

Water-Resources Investigations Report↗

Ground-water flow in the surficial aquifer system and potential movement of contaminants from selected waste-disposal sites at Naval Station Mayport, Florida

Ground-water flow through the surficial aquifer system at Naval Station Mayport near Jacksonville, Florida, was simulated with a two-layer finite-difference model as part of an investigation conducted by the U.S. Geological Survey. The model was calibrated to 229 water-level measurements from 181 wells during three synoptic surveys (July 17, 1995; July 31, 1996; and October 24, 1996). A quantifiable understanding of ground-water flow through the surficial aquifer was needed to evaluate remedial-action alternatives under consideration by the Naval Station Mayport to control the possible movement of contaminants from sites on the station. Multi-well aquifer tests, single-well tests, and slug tests were conducted to estimate the hydraulic properties of the surficial aquifer system, which was divided into three geohydrologic units�an S-zone and an I-zone separated by a marsh-muck confining unit. The recharge rate was estimated to range from 4 to 15 inches per year (95 percent confidence limits), based on a chloride-ratio method. Most of the simulations following model calibration were based on a recharge rate of 8 inches per year to unirrigated pervious areas. The advective displacement of saline pore water during the last 200 years was simulated using a particle-tracking routine, MODPATH, applied to calibrated steady-state and transient models of the Mayport peninsula. The surficial aquifer system at Naval Station Mayport has been modified greatly by natural and anthropogenic forces so that the freshwater flow system is expanding and saltwater is being flushed from the system. A new MODFLOW package (VAR1) was written to simulate the temporal variation of hydraulic properties caused by construction activities at Naval Station Mayport. The transiently simulated saltwater distribution after 200 years of displacement described the chloride distribution in the I-zone (determined from measurements made during 1993 and 1996) better than the steady-state simulation. The advective movement of contaminants from selected sites within the solid waste management units to discharge points was simulated using MODPATH. Most of the particles were discharged to the nearest surface-water feature after traveling less than 1,000 feet in the ground-water system. Most areas within 1,000 feet of a surface-water feature or storm sewer had traveltimes of less than 50 years, based on an effective porosity of 40 percent. Contributing areas, traveltimes, and pathlines were identified for 224 wells at Naval Station Mayport under steady-state and transient conditions by back-tracking a particle from the midpoint of the wetted screen of each well. Traveltimes to contributing areas that ranged between 15 and 50 years, estimated by the steady-state model, differed most from the transient traveltime estimates. Estimates of traveltimes and pathlines based on steady-state model results typically were 10 to 20 years more and about twice as long as corresponding estimates from the transient model. The models differed because the steady-state model simulated 1996 conditions when Naval Station Mayport had more impervious surfaces than at any earlier time. The expansion of the impervious surfaces increased the average distance between contributing areas and observation wells.

Florida↗

Ground-water flow in the surficial aquifer system and potential movement of contaminants from selected waste-disposal sites at Cecil Field Naval Air Station, Jacksonville, Florida

As part of the Installation Restoration Program, Cecil Field Naval Air Station, Jacksonville, Florida, is considering remedialaction alternatives to control the possible movement of contaminants from sites that may discharge to the surface. This requires a quantifiable understanding of ground-water flow through the surficial aquifer system and how the system will respond to any future stresses. The geologic units of interest in the study area consist of sediments of Holocene to Miocene age that extend from land surface to the base of the Hawthorn Group. The hydrogeology within the study area was determined from gamma-ray and geologists' logs. Ground-water flow through the surficial aquifer system was simulated with a seven-layer, finite-difference model that extended vertically from the water table to the top of the Upper Floridan aquifer. Results from the calibrated model were based on a long-term recharge rate of 6 inches per year, which fell in the range of 4 to 10 inches per year, estimated using stream hydrograph separation methods. More than 80 percent of ground-water flow circulates within the surficial-sand aquifer, which indicates that most contaminant movement also can be expected to move through the surficial-sand aquifer alone. The surficial-sand aquifer is the uppermost unit of the surficial aquifer system. Particle-tracking results showed that the distances of most flow paths were 1,500 feet or less from a given site to its discharge point. For an assumed effective porosity of 20 percent, typical traveltimes are 40 years or less. At all of the sites investigated, particles released 10 feet below the water table had shorter traveltimes than those released 40 feet below the water table. Traveltimes from contaminated sites to their point of discharge ranged from 2 to 300 years. The contributing areas of the domestic supply wells are not very extensive. The shortest traveltimes for particles to reach the domestic supply wells from their respective contributing areas ranged from 70 to 200 years.

Florida↗

Hydrogeologic conditions and simulation of ground-water flow in the Greater Orlando Metropolitan Area, East-Central Florida

A finite-difference ground-water flow model was used to simulate the effects of both modern-day (1988) and projected 2010 ground-water withdrawals on the Floridan aquifer system in the greater Orlando metropolitan area. This area covers about 2,500 square miles and includes all of Orange and Seminole Counties and parts of Lake, Volusia, Brevard, Osceola, and Polk Counties. The hydrogeology of the area is characterized by a thin surficial aquifer underlain by the thick, highly productive rocks of the Floridan aquifer system. Water in the Upper Floridan aquifer is brackish (chloride concentrations greater than 1,000 milligrams per liter) in discharge areas beneath and near the St. Johns and Wekiva Rivers and is freshest (chloride concentrations less than 100 milligrams per liter) inrecharge areas. A slight trend toward increasing concentrations of dissolved solids, chloride, and sulfate has been observed at Upper Floridan aquifer springs. Chloride concentrations in the Upper Floridan aquifer measured between 1966 and 1993 at the Cocoa well field have increased from 50 milligrams per liter to 120 milligrams per liter; concentrations measured in the Lower Floridan aquifer between 1966 and 1993 have increasedfrom 600 milligrams per liter to 3,000 milligrams per liter. The flow model was calibrated by comparing (a) simulated and estimated Upper Floridan aquifer predevelopment (unstressed) potentiometric surfaces, (b) simulated and measured heads at 142 Upper Floridan aquifer monitoring wells in 1988 (averageabsolute error of 1.8 feet), (c) simulated and measured discharge rates at 15 Upper Floridan aquifer springs in 1988 (306 cubic feet per second), and (d) simulated and measured drawdowns at 134 Upper Floridan aquifer monitoring wells between 1988 and May 1990 (58 and 95 percent of simulated drawdowns were within plus or minus 25and 50 percent of measured drawdowns, respectively). Relative to predevelopment conditions, model simulations indicate that about half of the 305 million gallons per day of water pumped from the Floridan aquifer system in 1988 was accounted for by increased recharge from the surficial aquifer system. About 23 cubic feet persecond was derived from increased lateral inflow. A storage coefficient of 1x10-3 provided the best comparisons of measured-to-simulated data during the transient simulation from January to May 1990. This storativity probably is greater than the true storativity of the Upper Floridan aquifer because storage contributions from the intermediateconfining unit were not accounted for during model design and development. Calibrated transmissivity ranged from 10,000 to greater than 400,000 feet squared per day in the Upper Floridan aquifer, and from 5,000 to 600,000 feet squared per day in the Lower Floridan aquifer. Calibrated intermediate confining unit leakance ranged from 1x10-5 to 4x10-3 per day and was highest in areas where the unit is thin or has been breached by numerous sinkholes. In general,calibrated transmissivity and leakance values were higher than associated aquifer-test values. Simulated recharge rates to the Upper Floridan aquifer from the surficial aquifer system ranged from less than 3 to 21 inches per year. Recharge rates of greater than 10 inches per year were simulated in areas of west Seminole, west Orange, east Lake, and southwest Volusia Counties. Recharge rates of less than 3 inches per year were simulated in east Orange and northeast Osceola Counties. The calibrated model was used to simulate the effects of increased Floridan aquifer withdrawals in the year 2010 (542 million gallons per day) on water levels and spring flow. Projected effects were simulated for both "wet" conditions (using 1988 fixed-head arrays) and for "dry" conditions (using May 1990 fixed-head arrays), thus bracketing a potential range of effects. Relative to simulated 1988 conditions, simulated 2010 spring flow decreased by 43 cubic f

Water-Resources Investigations Report↗

Analysis of ground-water flow in the Catahoula aquifer system in the vicinity of Laurel and Hattiesburg, Mississippi

The upper, middle, and lower Catahoula aquifers in the vicinity of the cites of Laurel and Hattiesburg in southern Mississippi are made up of irregular, discontinuous sand zones in the Catahoula Formation of Miocene age. In places thee three aquifers may be hydraulically well connected, and are referred to as the Catahoula aquifer system. Withdrawal from the Catahoula aquifers increased from 28 million gallons per day (Mgal/d) to 41 Mgal/d during 1970 to 1985, and decreased to 38 Mgal/d during 1990. Most withdrawal in the Laurel area is from the lower and middle Catahoula, and most withdrawal in the Hattiesburg area is from the middle and upper Catahoula aquifers. In the Laurel area, water levels in selected wells in the lower Catahoula aquifer declined at rates ranging from about 1 to 3.6 feet/ year until the late 1980's in response to the increase in pumping. A three-dimensional model was developed to represent ground-water flow in the Catahoula aquifers. Simulated water levels in the lower Catahoula aquifer, the layer most affected by pumping, were lowered from predevelopment levels as much as 130 feet in the Laurel area and 100 feet in the Hattiesburg area, according to the model analysis of 1992 conditions. Three scenarios of increased pumpage, for the period 1992-2020, were simulated. Under the low-growth scenario, water- level declines would be 20 feet or less below 1992 water levels in the middle and upper Catahoula aquifer in the Hattiesburg area, and about 60 feet in the lower Catahoula aquifer in the Laurel area. Under the moderate-growth scenario, water-level declines would be 40 feet or less below 1992 water levels in the middle Catahoula aquifer in the Hattiesburg area. Water-level declines would be about 110 feet in the lower Catahoula aquifer in the Laurel area, and water levels would approach the top of the aquifer. Under the high-growth scenario, water-level declines would be 40 feet or less in the upper Catahoula aquifer and about 80 feet in the middle Catahoula, with the largest declines occurring in the Hattiesburg area. Water levels would decline about 130 feet and would be drawn down below the top of the lower Catahoula aquifer in the Laurel area under the high-growth scenario.

Water-Resources Investigations Report↗