Geology ReportsSearch

SEARCH · Geology Reports

Results for “Ground Water Modeling Investigation Report”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24Linked to original sources

Hydrochemical Regions of the Glacial Aquifer System, Northern United States, and Their Environmental and Water-Quality Characteristics

The glacial aquifer system in the United States is a large (953,000 square miles) regional aquifer system of heterogeneous composition. As described in this report, the glacial aquifer system includes all unconsolidated geologic material above bedrock that lies on or north of the line of maximum glacial advance within the United States. Examining ground-water quality on a regional scale indicates that variations in the concentrations of major and minor ions and some trace elements most likely are the result of natural variations in the geologic and physical environment. Study of the glacial aquifer system was designed around a regional framework based on the assumption that two primary characteristics of the aquifer system can affect water quality: intrinsic susceptibility (hydraulic properties) and vulnerability (geochemical properties). The hydrochemical regions described in this report were developed to identify and explain regional spatial variations in ground-water quality in the glacial aquifer system within the hypothetical framework context. Data analyzed for this study were collected from 1991 to 2003 at 1,716 wells open to the glacial aquifer system. Cluster analysis was used to group wells with similar ground-water concentrations of calcium, chloride, fluoride, magnesium, potassium, sodium, sulfate, and bicarbonate into five unique groups. Maximum Likelihood Classification was used to make the extrapolation from clustered groups of wells, defined by points, to areas of similar water quality (hydrochemical regions) defined in a geospatial model. Spatial data that represented average annual precipitation, average annual temperature, land use, land-surface slope, vertical soil permeability, average soil clay content, texture of surficial deposits, type of surficial deposit, and potential for ground-water recharge were used in the Maximum Likelihood Classification to classify the areas so the characteristics of the hydrochemical regions would resemble the characteristics of the clusters. The result of the Maximum Likelihood Classification is a map showing five hydrochemical regions of the glacial aquifer system. Statistical analysis of ion concentrations (calcium, chloride, fluoride, magnesium, sodium, potassium, sulfate, and bicarbonate) in samples collected from wells completed in the glacial aquifer system illustrates that variations in water quality can be explained, in part, by related environmental characteristics that control the movement of ground water through the aquifer system. A comparison of median concentrations of chemical constituents in ground water among the five hydrochemical regions indicates that ground water in the Midwestern Agricultural Region, the Urban-Influenced Region, and the Western Agriculture and Grassland Region has the highest concentrations of major and minor ions, whereas ground water in the Northern and Great Lakes Forested Region and the Mountain and Coastal Forested Region has the lowest concentrations of these ions. Median concentrations of barium, arsenic, lithium, boron, strontium, and nitrite plus nitrate as nitrogen also are significantly different among the hydrochemical regions.

Scientific Investigations Report

Numerical simulation of the alluvium and terrace aquifer along the North Canadian River from Canton Lake to Lake Overholser, central Oklahoma

Alluvium and terrace deposits of Quaternary age, which cover an area of about 400 square miles along the North Canadian River between Canton Lake and Lake Overholser, locally yield as much as 500 gallons per minute to wells. The deposits are as much as 100 feet thick and consist of clay, silt, sand, and gravel, with sand-sized material dominating. The underlying bedrock is Permian sandstone and shale. The amount of water stored in the aquifer during 1980 was estimated to be 4.00 x 1010 cubic feet. A digital model was used to estimate the ability of the aquifer to continue to supply water for irrigation, industry, and domestic use. A block-centered, finite-difference model with 1 mile node spacing was used to project the amount and distribution of water in the aquifer in the future. The model was calibrated using the aquifer discharge to the North Canadian River and the difference between computed and measured heads. The model was calibrated using recharge rate of 1 inch per year, a horizontal hydraulic conductivity of 4.5 x 10 -4 feet per second and a specific yield of 0.16. Model simulations using the 1979 pumping rate, a projected pumping rate (based on expected increases in water use), and double the projected increase in pumping rate were made to 1993. With the 1979 pumping rate, the volume of water in storage in 1993 is 3.88 x 10 10 cubic feet and ground-water discharge to the North Canadian River is 10.9 cubic feet per second. With the projected pumping increase, the volume of water in storage is 3.81 x 1010 cubic feet and ground-water discharge to the stream is 6.91 cubic feet per second. At double the projected increase in pumping rate, the volume of water in storage is 3.74 x 1010 cubic feet and the ground-water discharge is 2.93 cubic feet per second.

Water-Resources Investigations Report

Geohydrology and simulation of flow and water levels in the aquifer system in the Mud Lake area of the eastern Snake River plain, eastern Idaho

Water users rely on surface water and ground water to irrigate crops and to maintain lakes on wildlife refuges in the 2,200-square-mile Mud Lake study area. Ground-water development between the late 1970's and 1989 increased withdrawals from about 240,000 acre-feet in 1983 to about 370,000 acre-feet in 1990. Concurrent with ground-water development, change from subirrigation to sprinkler irrigation was predicted to reduce recharge by 95,000 acre-feet, according to an independent study. Of the 660,000 acre-feet total estimated recharge from precipitation and irrigation in the study area in 1980, half was in the area in which irrigation methods were changed. Water managers need the ability to evaluate the effects of water-use changes on the future supply of surface water and ground water. Basalt and rhyolite predominate on the surface and in the subsurface of the study area. Total basalt thickness is less than 4,000 feet; total sediment thickness (clay, silt, sand, and gravel) is less than 1,000 feet. Basalt and sediment interbeds contribute to confined ground-water conditions and affect movement and supply of water in parts of the aquifer system. Estimated losses from and gains to perennial streams and lakes in 1980 were each about 110,000 acre-feet. Water-table altitudes ranged from about 4,500 to 6,200 feet above sea level, and water-table gradients were 3 to 120 feet per mile. Underflow from basins tributary to the study area was estimated to be about 450,000 acre-feet in 1980; measured discharge from flowing wells was about 10,000 acre-feet. A five-layer, three-dimensional, finite-difference, numerical ground-water flow model was calibrated by trial-and-error to assumed 1980 steady-state hydrologic conditions to obtain a better understanding of the geohydrology and provide a tool to evaluate water-use alternatives. Water-level gradients simulated by the model were similar to gradients measured in 1980. Simulated underflow across model boundaries for 1980 was 932,000 acre-feet. Simulated losses from and gains to most streams and lakes were within 2 percent of estimated values. Simulated discharge from flowing wells matched measurements for 1980. An attempt to calibrate the numerical model to transient hydrologic conditions in monthly increments from 1981 to 1990 was discontinued because available data did not justify changes that were indicated by model simulations.

Idaho

Hydrogeology, estimated impact, and regional well monitoring of effects of subsurface wastewater injection, Tampa Bay area, Florida

Six proposed injection sites are located in Pinellas County, Fla., and the city of St. Petersburg. Projected maximum injection rate, if all sites become operational, will be about 40 million gallons per day. The injection zone at the proposed sites is in a consistently dolomitized section of the Avon Park Limestone in the lower part of the Floridan aquifer. The injection zone contains saline ground water that has a chloride concentration of 19,000 to 20,000 milligrams per liter. Pressure and velocity changes were computed at selected regional locations in the upper and lower parts of the Floridan aquifer. Results of the model computations suggest that the regional impact after 20 years of injection will be small. Three locations are proposed for regional monitoring of subsurface injection. They are in the vicinity of the intersection of highways U.S. 19 and U.S. 60 in Pinellas County, Sun City in Hillsborough County, and the intersection of Sheldon Road and Gunn Highway in Hillsborough County. (USGS)

Water-Resources Investigations Report

Geohydrology and water quality of confined-drift aquifers in the Brooten-Belgrade area, west-central Minnesota

Confined-drift aquifers in six aquifer zones identified in a 1,300-square-mile area of west-central Minnesota near Brooten and Belgrade range in thickness from 5 to 110 feet. Transmissivities generally range from 500 to 10,000 feet squared per day, and theoretical well yields generally range from 100 to 900 gallons per minute. Regional ground-water flow in the confined-drift aquifers is to the southeast with local discharge to the East and Middle Branches of the Chippewa River, the North Fork Crow and Sauk Rivers, and to smaller streams, lakes, wetlands, and wells. Water levels near high-capacity pumped wells generally fluctuate 5 to 40 feet annually, compared to annual fluctuations of less than 5 feet in the unconfined aquifer. Water from confined-drift aquifers generally is suitable for most uses. The water is hard to very hard and contains locally elevated concentrations of iron, manganese, and dissolved solids. Results from a ground-water-flow model indicate that increased pumping from confined aquifers in the area would not adversely affect water levels. The addition of 10 to 20 hypothetical wells, pumping 123 to 246 million gallons per year, generally resulted in regional water-level declines of 0.1 to 1.0 feet. Simulations showed that the reduced recharge and increased pumping resulting from a 3-year drought probably would lower water levels between 5 and 10 feet regionally in the confined-drift aquifers and as much as 20 feet locally in the unconfined aquifer. Ground-water discharge to the East Branch Chippewa and North Fork Crow Rivers during the simulated drought would be reduced by 38 percent of 1984 conditions.

Minnesota

Simulated changes in ground-water levels resulting from proposed phosphate mining, west-central Florida; preliminary results

A digital model of two-dimensional ground-water flow was used to simulate projected changes in the Floridan aquifer potentiometric surface in 1985 and 2000, resulting from proposed ground-water developments by the phosphate mining industry in west-central Florida. The .model was calibrated under steady-state conditions to simulate the September 1975 potentiometric surface. Under one development plan, existing phosphate mines in Polk County would continue to withdraw ground water at 1975 rates, until phased out as the ore is depleted; no new mines would be introduced. Preliminary results indicate that under this plan, maximum simulated recovery of the potentiometric surface is 11.9 feet by 1985 and 36.5 feet by 2000. Under an alternative plan, all proposed mines in Polk, Hardee, DeSoto, Hillsborough and Manatee Counties would begin operations: in addition to the continuation and phasing out of existing mines. Preliminary results indicate that the potentiometric surface would generally recover in Polk County and decline elsewhere in the modeled area. Maximum simulated recovery is 4.5 feet by 1985 and 29.6 feet by 2000; maximum simulated drawdown is 15.1 feet by 1985 and feet by 2000. All results are preliminary and subject to revision as the investigation continues.

Open-File Report

Processing, Analysis, and General Evaluation of Well-Driller Logs for Estimating Hydrogeologic Parameters of the Glacial Sediments in a Ground-Water Flow Model of the Lake Michigan Basin

In 2005, the U.S. Geological Survey began a pilot study for the National Assessment of Water Availability and Use Program to assess the availability of water and water use in the Great Lakes Basin. Part of the study involves constructing a ground-water flow model for the Lake Michigan part of the Basin. Most ground-water flow occurs in the glacial sediments above the bedrock formations; therefore, adequate representation by the model of the horizontal and vertical hydraulic conductivity of the glacial sediments is important to the accuracy of model simulations. This work processed and analyzed well records to provide the hydrogeologic parameters of horizontal and vertical hydraulic conductivity and ground-water levels for the model layers used to simulated ground-water flow in the glacial sediments. The methods used to convert (1) lithology descriptions into assumed values of horizontal and vertical hydraulic conductivity for entire model layers, (2) aquifer-test data into point values of horizontal hydraulic conductivity, and (3) static water levels into water-level calibration data are presented. A large data set of about 458,000 well driller well logs for monitoring, observation, and water wells was available from three statewide electronic data bases to characterize hydrogeologic parameters. More than 1.8 million records of lithology from the well logs were used to create a lithologic-based representation of horizontal and vertical hydraulic conductivity of the glacial sediments. Specific-capacity data from about 292,000 well logs were converted into horizontal hydraulic conductivity values to determine specific values of horizontal hydraulic conductivity and its aerial variation. About 396,000 well logs contained data on ground-water levels that were assembled into a water-level calibration data set. A lithology-based distribution of hydraulic conductivity was created by use of a computer program to convert well-log lithology descriptions into aquifer or nonaquifer categories and to calculate equivalent horizontal and vertical hydraulic conductivities (K and KZ, respectively) for each of the glacial layers of the model. The K was based on an assumed value of 100 ft/d (feet per day) for aquifer materials and 1 ft/d for nonaquifer materials, whereas the equivalent KZ was based on an assumed value of 10 ft/d for aquifer materials and 0.001 ft/d for nonaquifer materials. These values were assumed for convenience to determine a relative contrast between aquifer and nonaquifer materials. The point values of K and KZ from wells that penetrate at least 50 percent of a model layer were interpolated into a grid of values. The K distribution was based on an inverse distance weighting equation that used an exponent of 2. The KZ distribution used inverse distance weighting with an exponent of 4 to represent the abrupt change in KZ that commonly occurs between aquifer and nonaquifer materials. The values of equivalent hydraulic conductivity for aquifer sediments needed to be adjusted to actual values in the study area for the ground-water flow modeling. The specific-capacity data (discharge, drawdown, and time data) from the well logs were input to a modified version of the Theis equation to calculate specific capacity based horizontal hydraulic conductivity values (KSC). The KSC values were used as a guide for adjusting the assumed value of 100 ft/d for aquifer deposits to actual values used in the model. Water levels from well logs were processed to improve reliability of water levels for comparison to simulated water levels in a model layer during model calibration. Water levels were interpolated by kriging to determine a composite water-level surface. The difference between the kriged surface and individual water levels was used to identify outlier water levels. Examination of the well-log lithology data in map form revealed that the data were not only useful for model input, but also were useful for understanding th

Scientific Investigations Report

Changing storm conditions in response to projected 21st century climate change and the potential impact on an arctic barrier island–lagoon system—A pilot study for Arey Island and Lagoon, eastern Arctic Alaska

Executive Summary Arey Lagoon, located in eastern Arctic Alaska, supports a highly productive ecosystem, where soft substrate and coastal wet sedge fringing the shores are feeding grounds and nurseries for a variety of marine fish and waterfowl. The lagoon is partially protected from the direct onslaught of Arctic Ocean waves by a barrier island chain (Arey Island) which in itself provides important habitat for migratory shorebirds and waterfowl. In this work, numerically modeled waves and water levels are computed under the provision of sea-level rise and changing conditions brought about by 21st century climate variability. Model results, supported by observations, are used to assess the stability of the barrier chain and spatiotemporal changes in flood patterns across fringing coastal wet sedge areas. The results aim to support studies that investigate the possibility of new biological succession trajectories and loss or increase of habitat areas.

Alaska

Processes affecting the trihalomethane concentrations associated with the third injection, storage, and recovery rest at Lancaster, Antelope Valley, California, March 1998 through April 1999

The formation and fate of trihalomethanes (THM) during the third injection, storage, and recovery test at Lancaster, Antelope Valley, California, were investigated as part of a program to assess the long-term feasibility of using injection, storage, and recovery as a water-supply method and as a way to reduce water-level declines and land-subsidence in the Antelope Valley. The program was conducted by the U.S. Geological Survey in cooperation with the Los Angeles County Department of Public Works and the Antelope Valley-East Kern Water Agency. The water used for injection, storage, and recovery must be disinfected before injection and thus contains THMs and other disinfection by-products. THMs (chloroform, CHCl3, bromodichloromethane, CHCl2Br, dibromochloromethane, CHClBr2, and bromoform, CHBr3) are formed by reaction between natural dissolved organic carbon that is present in water and chlorine that is added during the disinfection step of the drinking water treatment process. THMs are carcinogenic compounds, and their concentrations in drinking water are regulated by the U.S. Environmental Protection Agency. During previous cycles of the Lancaster program, extracted water still contained measurable concentrations of THMs long after continuous pumping had extracted a greater volume of water than had been injected. This raised concerns about the potential long-term effect of injection, storage, and recovery cycles on ground-water quality in Antelope Valley aquifers. The primary objectives of this investigation were to determine (1) what controlled continued THM formation in the aquifer after injection, (2) what caused of the persistence of THMs in the extracted water, even after long periods of pumping, (3) what controlled the decrease of THM concentrations during the extraction period, and (4) the potential for natural attenuation of THMs in the aquifer. Laboratory experiments on biodegradation of THMs in microcosms of aquifer materials indicate that aquifer bacteria did not degrade CHCl3 or CHBr3 under aerobic conditions, but did degrade CHBr3 under anaerobic conditions. However, the aquifer is naturally aerobic and CHCl3 is the dominant THM species; therefore, biodegradation is not considered an important attenuation mechanism for THMs in this aquifer. The alluvial-fan sediments comprising the aquifer have very low contents of organic matter; therefore, sorption is not considered to be an important attenuation mechanism for THMs in this aquifer. Laboratory experiments on formation of THMs in the injection water indicate that continued THM formation in the injection water after injection into the aquifer was limited by the amount of residual chlorine in the injection water at the time of injection. After accounting for THMs formed by reaction of this residual chlorine, THMs behaved as conservative constituents in the aquifer, and the only process affecting the concentration of THMs was mixing of the injection water and the ground water. The mixing process was quantified using mass balances of injected constituents, the sulfur hexafluoride (SF6) tracer that was added to the injected water, and a simple descriptive mathematical mixing model. Mass balance calculations show that only 67 percent of the injected THMs and chloride were recovered by the time that a volume of water equivalent to 132 percent of the injection water volume was extracted. Pumping 250 percent of the injection water volume only increased recovery of injected THMs to 80 percent. THM and SF6 concentrations in the extracted water decreased concomitantly during the extraction period, and THM concentrations predicted from SF6 concentrations closely matched the measured THM concentrations. Because SF6 is a conservative tracer that was initially only present in the injection water, parallel decreases in SF6 and THM concentrations in the extracted water must be due to dilution of injection water with ground water. The simple descriptive mixing mode

Water-Resources Investigations Report

A Numerical model to evaluate proposed ground-water allocations in southwest Kansas

A computer model was developed to assist the Southwest Kansas Groundwater Management District No. 3 in the evaluation of applications to appropriate ground water. The model calculated the drawdown due from a proposed well at all existing wells in the section of the proposed well and at all wells in the adjacent eight sections. The depletion expected in the 9-square-mile area due to all existing wells and the proposed well is computed and compared with allowable limits defined by the management district. An optional program permits the evaluation of allowable depletion for one or more townships. All options are designed to run interactively, thus allowing for immediate evaluation of proposed ground-water withdrawals. (USGS)

Water-Resources Investigations Report

Development of the Hydroecological Integrity Assessment Process for Determining Environmental Flows for New Jersey Streams

The natural flow regime paradigm and parallel stream ecological concepts and theories have established the benefits of maintaining or restoring the full range of natural hydrologic variation for physiochemical processes, biodiversity, and the evolutionary potential of aquatic and riparian communities. A synthesis of recent advances in hydroecological research coupled with stream classification has resulted in a new process to determine environmental flows and assess hydrologic alteration. This process has national and international applicability. It allows classification of streams into hydrologic stream classes and identification of a set of non-redundant and ecologically relevant hydrologic indices for 10 critical sub-components of flow. Three computer programs have been developed for implementing the Hydroecological Integrity Assessment Process (HIP): (1) the Hydrologic Indices Tool (HIT), which calculates 171 ecologically relevant hydrologic indices on the basis of daily-flow and peak-flow stream-gage data; (2) the New Jersey Hydrologic Assessment Tool (NJHAT), which can be used to establish a hydrologic baseline period, provide options for setting baseline environmental-flow standards, and compare past and proposed streamflow alterations; and (3) the New Jersey Stream Classification Tool (NJSCT), designed for placing unclassified streams into pre-defined stream classes. Biological and multivariate response models including principal-component, cluster, and discriminant-function analyses aided in the development of software and implementation of the HIP for New Jersey. A pilot effort is currently underway by the New Jersey Department of Environmental Protection in which the HIP is being used to evaluate the effects of past and proposed surface-water use, ground-water extraction, and land-use changes on stream ecosystems while determining the most effective way to integrate the process into ongoing regulatory programs. Ultimately, this scientifically defensible process will help to quantify the effects of anthropogenic changes and development on hydrologic variability and help planners and resource managers balance current and future water requirements with ecological needs.

Scientific Investigations Report

Geohydrology and simulated ground-water flow in northwestern Elkhart County, Indiana

In 1994, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency and the City of Elkhart, developed a ground-water model of the Elkhart, Indiana, area to determine the avail-ability and source of water at potential new well fields. The modeled area covered 190 square miles of northwestern Elkhart County and a small part of southern Michigan. Three Superfund sites and several other sites in this area are undergoing environmental cleanup. The model would be used to guide the location of well fields so that Superfund sites and environmental cleanup areas would not be within recharge areas for the well fields. The City of Elkhart obtains its water supply from two aquifers separated by a generally continuous confining unit. The upper aquifer is composed primarily of sand and gravel of glacial origin. Thickness of the upper aquifer ranges from 0 to 116 feet and averages 47 feet. The lower aquifer is composed of sand and gravel with interbedded lenses of silt and clay. Thickness of the lower aquifer ranges from 1 to 335 feet and averages 35 feet. The intervening confining unit is composed of silt and clay with interbedded sand and gravel; the confining unit ranges from 0 to 177 feet, with an average thickness of 27 feet. Flow through the aquifers is generally horizontal vertically downward from the upper aquifer, through the confining unit, and into the lower aquifer, except where flow is vertically upward at the St. Joseph River and other large streams. The hydraulic characteristics of the aquifers and confining unit were estimated by analyzing aquifer-test data from well drillers? logs and by calibration of the model. The horizontal hydraulic conductivity of the upper aquifer is 170 feet per day within about 1 mile of the St. Joseph and Elkhart Rivers and 370 feet per day at distances greater than about 1 mile. The horizontal hydraulic conductivity of the lower aquifer is 370 feet per day throughout the modeled area, with the exception of an area near the center of the modeled area where the horizontal hydraulic conductivity is 170 feet per day. Transmissivity of the lower aquifer increases generally from southwest to northeast; transmissivity values range from near 0 where the lower aquifer is absent to 57,000 square feet per day and average about 8,100 square feet per day. The vertical hydraulic conductivity of the confining unit is 0.07 feet per day; the vertical conductivity of the streambeds commonly is 1.0 foot per day and ranges from 0.05 foot per day to 50 feet per day. The areal recharge rate to the outwash deposits was determined by a base-flow separation technique to be 16 inches per year, and the areal recharge rate to the till was assumed to be 4 inches per year. A two-layer digital model was used to simulate flow in the ground-water system. The model was calibrated on the basis of historical water-use data, water-level records, and gain/loss data for streams during May and June 1979. The model was recalibrated with water-use data and water-level records from 1988. For 1979 data, 49 percent of the inflow to the model area is from precipitation and 46 percent is ground-water inflow across the model boundaries. Most of the ground-water inflow across the model boundary is from the north and east, which corresponds to high values of transmissivity?as high as 57,000 feet squared per day?in the model layers in the northern and eastern areas. Eighty-two percent of the ground-water discharge is to the streams; 5 percent of the ground-water discharge is to wells. Source areas and flow paths to the City of Elkhart public well fields are affected by the location of streams and the geology in the area. Flow to the North Well Field originates north-west of the well field, forms relatively straight flow paths, and moves southeast toward the well field and the St. Joseph River. Flow to the South Well Field begins mostly in the out-wash along Yellow Creek south of the well field, moves northward, and t

Indiana

Water Budget for the Iao Area, Island of Maui, Hawaii

Ground-water recharge is estimated as the residual component of a monthly water budget calculated using soil characteristics and long-term average rainfall, streamflow, irrigation, and pan-evaporation data. The water-budget components of rainfall, direct runoff, evapotranspiration, and ground-water recharge are defined seasonally, through the use of monthly data, and spatially by land-use and geohydrologic areas, through the use of a geographic information system model. The long-term average ground-water recharge for the Iao area was estimated for four scenarios using natural land-use, and using 1926-79, 1980-85, and 1986-95 land-use and irrigation data. The recharge rate for natural conditions is 34 million gallons per day, which is 34 percent of rainfall. The average annual ground-water recharge rate for 1926-79 conditions is 51 million gallons per day, which is 41 percent of the sum of rainfall and irrigation. The recharge rates for 1980-85 and 1986-95 conditions are 40 and 36 million gallons per day, which are 37 and 35 percent of rainfall plus irrigation, respectively.

Water-Resources Investigations Report

Structure, temperature, pressure, and salinity of Cenozoic aquifers of south Texas

A study of the hydrogeology of deep sedimentary basins using the Neogene deposits of the northern Gulf of Mexico basin as a model was initiated by the U.S. Geological Survey in the 1960's (Jones, 1969). This study led to investigations of geopressured-geothermal resources in this basin (Papadopulos and others, 1975; Wallace and others, 1979) and of the waste storage capabilities of the Wilcox Group in Texas (Jones and others, 1976). Techniques for delineating pressure, temperature, and salinity variations within the geologic framework were developed in conjunction with these and related studies. The purpose of this hydrologic atlas is to delineate the hydrogeology of the onshore Cenozoic aquifers of south Texas, using these techniques. The location of the study area in Texas is shown in figure 1. Hydrogeologic cross section A-A' is shown in figure 2. The surface geology of the area and the line of section A–A' are shown in figure 3. Geophysical logs of oil tests provided the framework for this section. Interpretations presented are based upon data from these logs; information from published ground water and geological reports; and structural and geologic interpretations from maps prepared by Peppard, Souders and Associates for the U.S. Geological Survey under Purchase Order Numbers 67482 and 67.483.

Texas

Relation Between Solid-Phase and Dissolved Arsenic in the Ground-Water System Underlying Northern Preble County, Ohio

The U.S. Geological Survey (USGS), in cooperation with the Miami Conservancy District, collected and analyzed samples of the aquifer materials and ground water from multiple depths at two sites in northern Preble County, Ohio. The aquifer materials included glacial deposits and Silurian carbonate bedrock. In the study area, elevated arsenic concentrations have been detected in ground water from both types of aquifers. The aquifer materials were described in terms of the stratigraphy and the bulk elemental composition of 70 samples. In addition, six water-producing horizons were selected for more detailed study; ground-water quality was analyzed, microanalytical techniques were used to examine thin sections of the aquifer materials, and simplified geochemical modeling was done to identify plausible reactions between the ground water and aquifer materials. At both study sites, the highest solid-phase arsenic concentrations were from a roughly similar stratigraphic position - a transition zone that extends from just above the Wisconsinan/Illinoian contact to just below the Pleistocene/Silurian contact. For carbonate bedrock, the solid-phase arsenic concentrations were generally low (<1 to 4 mg/kg (milligrams per kilogram)). The one notable exception was a thin horizon about 10 feet below the top of bedrock, which had an arsenic concentration of 42 mg/kg. This horizon showed some textural and compositional evidence of alteration by geothermal fluids. Additional study might be warranted to investigate whether arsenic concentrations in ground water from carbonate bedrock could be decreased by excluding discrete horizons from the open intervals of wells. For glacial deposits, solid-phase arsenic concentrations were slightly higher in fine-grained deposits (2 to 20 mg/kg) than in coarse-grained deposits (2 to 9 mg/kg). In ground water, arsenic concentrations ranged from <1 to 51 ug/L (micrograms per liter); samples from two horizons had concentrations greater than the U.S. Environmental Protection Agency Maximum Contaminant Level (MCL) of 10 ug/L. Dissolved arsenic concentrations appear to be more closely related to redox conditions of the ground water than to the arsenic content of the aquifer materials. Geochemical modeling and thin-section analysis were generally consistent with the idea that arsenic was released to water from iron oxides under iron-reducing conditions. In addition, there was some evidence in support of the idea that arsenic can be removed from ground water by precipitation of sulfide minerals, which occurs under sulfate-reducing conditions. At one site, the dissolved arsenic concentrations in two water-bearing horizons increased from <1 to 51 ug/L over a depth of 15 feet. The large increase might be due to a shift from sulfate-reducing to methanogenic conditions; in the absence of sulfate reduction, arsenic is not sequestered in sulfide minerals and may accumulate in the ground water.

Scientific Investigations Report

Geochemistry of the Madison and Minnelusa aquifers in the Black Hills area, South Dakota

The Madison and Minnelusa aquifers are two of the most important aquifers in the Black Hills area because of utilization for water supplies and important influences on surface-water resources resulting from large springs and streamflow- loss zones. Examination of geochemical information provides a better understanding of the complex flow systems within these aquifers and interactions between the aquifers. Major-ion chemistry in both aquifers is dominated by calcium and bicarbonate near outcrop areas, with basinward evolution towards various other water types. The most notable differences in major-ion chemistry between the Madison and Minnelusa aquifers are in concentrations of sulfate within the Minnelusa aquifer. Sulfate concentrations increase dramatically near a transition zone where dissolution of anhydrite is actively occurring. Water chemistry for the Madison and Minnelusa aquifers is controlled by reactions among calcite, dolomite, and anhydrite. Saturation indices for gypsum, calcite, and dolomite for most samples in both the Madison and Minnelusa aquifers are indicative of the occurrence of dedolomitization. Because water in the Madison aquifer remains undersaturated with respect to gypsum, even at the highest sulfate concentrations, upward leakage into the overlying Minnelusa aquifer has potential to drive increased dissolution of anhydrite in the Minnelusa Formation. Isotopic information is used to evaluate ground-water flowpaths, ages, and mixing conditions for the Madison and Minnelusa aquifers. Distinctive patterns exist in the distribution of stable isotopes of oxygen and hydrogen in precipitation for the Black Hills area, with isotopically lighter precipitation generally occurring at higher elevations and latitudes. Distributions of 18O in ground water are consistent with spatial patterns in recharge areas, with isotopically lighter 18O values in the Madison aquifer resulting from generally higher elevation recharge sources, relative to the Minnelusa aquifer. Three conceptual models, which are simplifications of lumped-parameter models, are considered for evaluation of mixing conditions and general ground-water ages. For a simple slug-flow model, which assumes no mixing, measured tritium concentrations in ground water can be related through a first-order decay equation to estimated concentrations at the time of recharge. Two simplified mixing models that assume equal proportions of annual recharge over a range of years also are considered. An ?immediate-arrival? model is used to conceptually represent conditions in outcrop areas and a ?time-delay? model is used for locations removed from outcrops, where delay times for earliest arrival of ground water generally would be expected. Because of limitations associated with estimating tritium input and gross simplifying assumptions of equal annual recharge and thorough mixing conditions, the conceptual models are used only for general evaluation of mixing conditions and approximation of age ranges. Headwater springs, which are located in or near outcrop areas, have the highest tritium concentrations, which is consistent with the immediate-arrival mixing model. Tritium concentrations for many wells are very low, or nondetectable, indicating general applicability of the timedelay conceptual model for locations beyond outcrop areas, where artesian conditions generally occur. Concentrations for artesian springs generally are higher than for wells, which indicates generally shorter delay times resulting from preferential flowpaths that typically are associated with artesian springs. In the Rapid City area, a distinct division of isotopic values for the Madison aquifer corresponds with distinguishing 18O signatures for nearby streams, where large streamflow recharge occurs. Previous dye testing in this area documented rapid ground-water flow (timeframe of weeks) from a streamflow loss zone to sites located several miles away. These results are used to ill

Water-Resources Investigations Report

Hydrogeology and ground-water flow in the carbonate rocks of the Little Lehigh Creek basin, Lehigh County, Pennsylvania

The Little Lehigh Creek basin is underlain mainly by a complex assemblage of highly-deformed Cambrian and Ordovician carbonate rocks. The Leithsville Formation, Allentown Dolomite, Beekmantown Group, and Jacksonburg Limestone act as a single hydrologic unit. Ground water moves through fractures and other secondary openings and generally is under water-table conditions. Median annual ground-water discharge (base flow) to Little Lehigh Creek near Allentown (station 01451500) during 1946-86 was 12.97 inches or 82 percent of streamflow. Average annual recharge for 1975-83 was 21.75 inches. Groundwater and surface-water divides do not coincide in the basin. Ground-water underflow from the Little Lehigh Creek basin to the Cedar Creek basin in 1987 was 4 inches per year. A double-mass curve analysis of the relation of cumulative precipitation at Allentown to the flow of Schantz Spring for 1956-84 showed that cessation of quarry pumping and development of ground water for public supply in the Schantz Spring basin did not affect the flow of Schantz Spring. Ground-water flow in the Little Lehigh Creek basin was simulated using a finite-difference, two-dimensional computer model. The geologic units in the modeled area were simulated as a single water-table aquifer. The 134-squaremile area of carbonate rocks between the Lehigh River and Sacony Creek was modeled to include the natural hydrologic boundaries of the ground-water-flow system. The ground-water-flow model was calibrated under steady-state conditions using 1975-83 average recharge, evapotranspiration, and pumping rates. Each geologic unit was assigned a different hydraulic conductivity. Initial aquifer hydraulic conductivity was estimated from specific-capacity data. The average (1975-83) water budget for the Little Lehigh Creek basin was simulated. The simulated base flow from the carbonate rocks of the Little Lehigh Creek basin above gaging station 01451500 is 11.85 inches per year. The simulated ground-water underflow from the Little Lehigh Creek basin to the Cedar Creek basin is 4.04 inches per year. For steady-state calibration, the root-mean-squared difference between observed and simulated heads was 21.19 feet. The effects of increased ground-water development on base flow and underflow out of the Little Lehigh Creek basin for average and drought conditions were simulated by locating a hypothetical well field in different parts of the basin. Steady-state simulations were used to represent equilibrium conditions, which would be the maximum expected long-term effect. Increased ground-water development was simulated as hypothetical well fields pumping at the rate of 15, 25, and 45 million gallons per day in addition to existing ground-water withdrawals. Four hypothetical well fields were located near and away from Little Lehigh Creek in upstream and downstream areas. The effects of pumping a well field in different parts of the Little Lehigh Creek basin were compared. Pumping a well field located near the headwaters of Little Lehigh Creek and away from the stream would have greatest effect on inducing underflow from the Sacony Greek basin and the least effect on reducing base flow and underflow to the Ceda^r Creek basin. Pumping a well field located near the headwaters of Little Leh|igh Creek near the stream would have less impact on inducing underflow from|the Sacony Creek basin and a greater impact on reducing the base flow of Little Lehigh Creek because more of the pumpage would come from diverted base flow. Pumping a well field located in the downstream area of the Little Lehigh Creek basin away from the stream would have the greatest effect on the underflow to the Cedar Creek basin. Pumping a well field located in the downstream area of the Little Lehigh Creek basin near the stream would have the greatest effect on reducing the base flow of Little Lehigh Cteek. Model simulations show that groundwater withdrawals do not cause a proportional reduction in base flow. Under average conditions, ground-water withdrawals are equal to 48 to 70 percent of simulated base-flow reductions; under drought conditions, ground-water withdrawals are equal to 35 to 73 percent of simulated base-flow reductions. The hydraulic effects of pumping largely depend on well location. In the Little Lehigh basin, surface-water and ground-water divides do not coincide, and ground-water development, especially near surface-water divides, can cause ground-water divides to shift and induce ground-water underflow from adjacent basins. Large-scale ground-water pumping in a basin may not produce expected reductions of base flow in that basin because of shifts in the ground-water divide; however, such shifts can reduce base flow in adjacent surface-water basins.

Water-Resources Investigations Report

Hydrogeology of the Sarasota-Port Charlotte area, Florida

The surficial and intermediate aquifers are the major source of public water supplies in the Sarasota-Port Charlotte, Florida, area because of the relatively poor quality of Floridan aquifer water. The hydrogeologic framework consists of the surficial aquifer, intermediate aquifers (Tamiami-upper Hawthorn and lower Hawthorn-upper Tampa aquifers) and confining beds, Floridan aquifer, and lower confining bed (or base of the Floridan aquifer). The quality of ground water in the surficial and intermediate aquifers is generally good, except in the western (coastal) and southern parts where saltwater intrusion or incomplete flushing of connate water has occurred. The mineral content of ground water generally increases with depth and areally from the northeast towards the west and south. A water budget for the study area shows that an average annual rainfall of 51.0 inches minus an evapotranspiration of 38.0 inches per year and streamflow of 12.5 inches per year leaves 0.5 inch per year of recharge to the surficial aquifer. Combined pumpage from the aquifers is 1.06 inches per year. A preliminary quasi-time dimensional model has been applied to the study area to check the reasonableness of the hydrogeologic framework defined and of aquifer parameters. The model was considered calibrated when the final head matrix was within plus or minus 5 feet of the starting head.

Florida