Geology Reports⌕ Search

Geology topics

David L. Parkhurst

Publications and source records attributed to David L. Parkhurst.

At least 37 records · Page 2Linked to original sources
USGS2007

Geochemical effects of induced stream-water and artificial recharge on the Equus Beds Aquifer, South-Central Kansas, 1995-2004

Artificial recharge of the Equus Beds aquifer is part of a strategy implemented by the city of Wichita, Kansas, to preserve future water supply and address declining water levels in the aquifer of as much as 30 feet caused by withdrawals for water supply and irrigation since the 1940s. Water-level declines represent a diminished water supply and also may accelerate migration of saltwater from the Burrton oil field to the northwest and the Arkansas River to the southwest into the freshwater of the Equus Beds aquifer. Artificial recharge, as a part of the Equus Beds Ground-Water Recharge Project, involves capturing flows larger than base flow from the Little Arkansas River and recharging the water to the Equus Beds aquifer by means of infiltration or injection. The geochemical effects on the Equus Beds aquifer of induced stream-water and artificial recharge at the Halstead and Sedgwick sites were determined through collection and analysis of hydrologic and water-quality data and the application of statistical, mixing, flow and solute-transport, and geochemical model simulations. Chloride and atrazine concentrations in the Little Arkansas River and arsenic concentrations in ground water at the Halstead recharge site frequently exceeded regulatory criteria. During 30 percent of the time from 1999 through 2004, continuous estimated chloride concentrations in the Little Arkansas River at Highway 50 near Halstead exceeded the Secondary Drinking-Water Regulation of 250 milligrams per liter established by the U.S. Environmental Protection Agency. Chloride concentrations in shallow monitoring wells located adjacent to the stream exceeded the drinking-water criterion five times from 1995 through 2004. Atrazine concentrations in water sampled from the Little Arkansas River had large variability and were at or near the drinking-water Maximum Contaminant Level of 3.0 micrograms per liter as an annual average established by the U.S. Environmental Protection Agency. Atrazine concentrations were much smaller than the drinking-water criterion and were detected at much smaller concentrations in shallow monitoring wells and diversion well water located adjacent to the stream probably because of sorption on aquifer sediment. Before and after artificial recharge, large, naturally occurring arsenic concentrations in the recharge water for the Halstead diversion well and recharge site exceeded the Maximum Contaminant Level of 10 micrograms per liter established by the U.S. Environmental Protection Agency for drinking water. Arsenic and iron concentrations decreased when water was recharged through recharge basins or a trench; however, chemical precipitation and potential biofouling eventually may decrease the artificial recharge efficiency through basins and trenches. At the Sedgwick site, chloride concentrations infrequently exceeded regulatory criteria. Large concentrations of atrazine were treated to decrease concentrations to less than regulatory criteria. Recharge of treated stream water through recharge basins avoids potentially large concentrations of arsenic and iron that exist at the Halstead diversion site. Results from a simple mixing model using chloride as a tracer indicated that the water chemistry in shallow monitoring well located adjacent to the Little Arkansas River was 80 percent of stream water, demonstrating effective recharge of the alluvial aquifer by the stream. Results also indicated that about 25 percent of the water chemistry of the diversion well water was from the shallow part of the aquifer. Additionally, diverting water through a diversion well located adjacent to the stream removed about 75 percent of the atrazine, probably through sorption to aquifer sediment, and decreased the need for additional water treatment to remove atrazine.

Kansas↗
USGS2004

Hydrologic and geochemical evaluation of aquifer storage recovery in the Santee Limestone/Black Mingo Aquifer, Charleston, South Carolina, 1998-2002

The hydrologic and geochemical effects of aquifer storage recovery were evaluated to determine the potential for supplying the city of Charleston, South Carolina, with large quantities of potable water during emergencies, such as earthquakes, hurricanes, or hard freezes. An aquifer storage recovery system, including a production well and three observation wells, was installed at a site located on the Charleston peninsula. The focus of this study was the 23.2-meter thick Tertiary-age carbonate and sand aquifer of the Santee Limestone and the Black Mingo Group, the northernmost equivalent of the Floridan aquifer system. Four cycles of injection, storage, and recovery were conducted between October 1999 and February 2002. Each cycle consisted of injecting between 6.90 and 7.19 million liters of water for storage periods of 1, 3, or 6 months. The volume of recovered water that did not exceed the U.S. Environmental Protection Agency secondary standard for chloride (250 milligrams per liter) varied from 1.48 to 2.46 million liters, which is equivalent to 21 and 34 percent of the total volume injected for the individual tests. Aquifer storage recovery testing occurred within two productive zones of the brackish Santee Limestone/Black Mingo aquifer. The individual productive zones were determined to be approximately 2 to 4 meters thick, based on borehole geophysical logs, electromagnetic flow-meter testing, and specific-conductance profiles collected within the observation wells. A transmissivity and storage coefficient of 37 meters squared per day and 3 x 10-5, respectively, were determined for the Santee Limestone/Black Mingo aquifer. Water-quality and sediment samples collected during this investigation documented baseline aquifer and injected water quality, aquifer matrix composition, and changes in injected/aquifer water quality during injection, storage, and recovery. A total of 193 water-quality samples were collected and analyzed for physical properties, major and minor ions, and nutrients. The aquifer and treated surface water were sodiumchloride and calcium/sodium-bicarbonate water types, respectively. Forty-five samples were collected and analyzed for total trihalomethane. Total trihalomethane data collected during aquifer storage recovery cycle 4 indicated that this constituent would not restrict the use of recovered water for drinking-water purposes. Analysis of six sediment samples collected from a cored well located near the aquifer storage recovery site showed that quartz and calcite were the dominant minerals in the Santee Limestone/Black Mingo aquifer. Estimated cation exchange capacity ranged from 12 to 36 milliequivalents per 100 grams in the lower section of the aquifer. A reactive transport model was developed that included two 2-meter thick layers to describe each of the production zones. The four layers composing the production zones were assigned porosities ranging from 0.1 to 0.3 and hydraulic conductivities ranging from 1 to 8.4 meters per day. Specific storage of the aquifer and confining units was estimated to be 1.5 x 10-5 meter-1. Longitudinal dispersivity of all layers was specified to be 0.5 meter. Leakage through the confining unit was estimated to be minimal and, therefore, not used in the reactive transport modeling. Inverse geochemical modeling indicates that mixing, cation exchange, and calcite dissolution are the dominant reactions that occur during aquifer storage recovery testing in the Santee Limestone/Black Mingo aquifer. Potable water injected into the Santee Limestone/Black Mingo aquifer evolved chemically by mixing with brackish background water and reaction with calcite and cation exchangers in the sediment. Reactive-transport model simulations indicated that the calcite and exchange reactions could be treated as equilibrium processes.

South Carolina↗
USGS2004

PHAST--a program for simulating ground-water flow, solute transport, and multicomponent geochemical reactions

The computer program PHAST simulates multi-component, reactive solute transport in three-dimensional saturated ground-water flow systems. PHAST is a versatile ground-water flow and solute-transport simulator with capabilities to model a wide range of equilibrium and kinetic geochemical reactions. The flow and transport calculations are based on a modified version of HST3D that is restricted to constant fluid density and constant temperature. The geochemical reactions are simulated with the geochemical model PHREEQC, which is embedded in PHAST. PHAST is applicable to the study of natural and contaminated ground-water systems at a variety of scales ranging from laboratory experiments to local and regional field scales. PHAST can be used in studies of migration of nutrients, inorganic and organic contaminants, and radionuclides; in projects such as aquifer storage and recovery or engineered remediation; and in investigations of the natural rock-water interactions in aquifers. PHAST is not appropriate for unsaturated-zone flow, multiphase flow, density-dependent flow, or waters with high ionic strengths. A variety of boundary conditions are available in PHAST to simulate flow and transport, including specified-head, flux, and leaky conditions, as well as the special cases of rivers and wells. Chemical reactions in PHAST include (1) homogeneous equilibria using an ion-association thermodynamic model; (2) heterogeneous equilibria between the aqueous solution and minerals, gases, surface complexation sites, ion exchange sites, and solid solutions; and (3) kinetic reactions with rates that are a function of solution composition. The aqueous model (elements, chemical reactions, and equilibrium constants), minerals, gases, exchangers, surfaces, and rate expressions may be defined or modified by the user. A number of options are available to save results of simulations to output files. The data may be saved in three formats: a format suitable for viewing with a text editor; a format suitable for exporting to spreadsheets and post-processing programs; or in Hierarchical Data Format (HDF), which is a compressed binary format. Data in the HDF file can be visualized on Windows computers with the program Model Viewer and extracted with the utility program PHASTHDF; both programs are distributed with PHAST. Operator splitting of the flow, transport, and geochemical equations is used to separate the three processes into three sequential calculations. No iterations between transport and reaction calculations are implemented. A three-dimensional Cartesian coordinate system and finite-difference techniques are used for the spatial and temporal discretization of the flow and transport equations. The non-linear chemical equilibrium equations are solved by a Newton-Raphson method, and the kinetic reaction equations are solved by a Runge-Kutta or an implicit method for integrating ordinary differential equations. The PHAST simulator may require large amounts of memory and long Central Processing Unit (CPU) times. To reduce the long CPU times, a parallel version of PHAST has been developed that runs on a multiprocessor computer or on a collection of computers that are networked. The parallel version requires Message Passing Interface, which is currently (2004) freely available. The parallel version is effective in reducing simulation times. This report documents the use of the PHAST simulator, including running the simulator, preparing the input files, selecting the output files, and visualizing the results. It also presents four examples that verify the numerical method and demonstrate the capabilities of the simulator. PHAST requires three input files. Only the flow and transport file is described in detail in this report. The other two files, the chemistry data file and the database file, are identical to PHREEQC files and the detailed description of these files is found in the PHREEQC documentation.

Techniques and Methods↗
USGS2003

Reactive-transport simulation of phosphorus in the sewage plume at the Massachusetts Military Reservation, Cape Cod, Massachusetts

The subsurface transport of phosphorus introduced by the disposal of treated sewage effluent to ground-infiltration disposal beds at the Massachusetts Military Reservation on western Cape Cod was simulated with a three-dimensional reactive-transport model. The simulations were used to estimate the load of phosphorus transported to Ashumet Pond during operation of the sewage-treatment plant from 1936 to 1995 and for 60 years following cessation of sewage disposal. The model accounted for spatial and temporal changes in water discharge from the sewage-treatment plant, ground-water flow, transport of associated chemical constituents, and a set of chemical reactions, including phosphorus sorption on aquifer materials, dissolution and precipitation of iron- and manganese-oxyhydroxide and iron phosphate minerals, organic carbon sorption and decomposition, cation sorption, and irreversible denitrification. The flow and transport in the aquifer were simulated by using parameters consistent with those used in previous flow models of this area of Cape Cod, except that numerical dispersion was much larger than the physical dispersion estimated in previous studies. Sorption parameters were fit to data derived from phosphorus sorption and desorption laboratory column experiments. Rates of organic carbon decomposition were adjusted to match the location of iron concentrations in an anoxic iron zone within the sewage plume. The sensitivity of the simulated load of phosphorus transported to Ashumet Pond was calculated for a variety of processes and input parameters. Model limitations included large uncertainties associated with the loading of the sewage beds, the flow system, and the chemistry and sorption characteristics in the aquifer. The results of current model simulations indicate a small load of phosphorus transported to Ashumet Pond during 1965-85, but this small load was particularly sensitive to model parameters that specify flow conditions and the chemical process by which non-desorbable phosphorus is incorporated in the sediments. The uncertainties were large enough to make it difficult to determine whether loads of phosphorus transported to Ashumet Pond in the 1990s were greater or less than loads during the previous two decades. The model simulations indicate substantial discharge of phosphorus to Ashumet Pond after about 1965. After the period 2000-10 the simulations indicate that the load of phosphorus transported to Ashumet Pond decreases continuously, but the load of phosphorus remains substantial for many decades. The current simulations indicate a peak in phosphorus discharge to Ashumet Pond of about 1,000 kilograms per year during the 1990s; however, comparisons of simulated phosphorus concentrations with measured concentrations in 1993 indicate that the peak in phosphorus load transported to Ashumet Pond may be larger and moving more quickly in the model simulations than in the aquifer. The results of the three-dimensional reactive-transport simulations are consistent with the loading history, experimental laboratory data, and field measurements. The results of the simulations adequately reproduce the spatial distribution of phosphorus concentrations measured in 1993, the magnitude of changes in phosphorus concentration with time in a profile near the disposal beds following cessation of sewage disposal, the observed iron zone in the sewage plume, the approximate flow of treated sewage effluent into Ashumet Valley, and laboratory-column data for phosphorus sorption and desorption.

Massachusetts↗
USGS2002

Calculation of individual isotope equilibrium constants for implementation in geochemical models

Theory is derived from the work of Urey to calculate equilibrium constants commonly used in geochemical equilibrium and reaction-transport models for reactions of individual isotopic species. Urey showed that equilibrium constants of isotope exchange reactions for molecules that contain two or more atoms of the same element in equivalent positions are related to isotope fractionation factors by , where is n the number of atoms exchanged. This relation is extended to include species containing multiple isotopes, for example and , and to include the effects of nonideality. The equilibrium constants of the isotope exchange reactions provide a basis for calculating the individual isotope equilibrium constants for the geochemical modeling reactions. The temperature dependence of the individual isotope equilibrium constants can be calculated from the temperature dependence of the fractionation factors. Equilibrium constants are calculated for all species that can be formed from and selected species containing , in the molecules and the ion pairs with where the subscripts g, aq, l, and s refer to gas, aqueous, liquid, and solid, respectively. These equilibrium constants are used in the geochemical model PHREEQC to produce an equilibrium and reaction-transport model that includes these isotopic species. Methods are presented for calculation of the individual isotope equilibrium constants for the asymmetric bicarbonate ion. An example calculates the equilibrium of multiple isotopes among multiple species and phases.

Water-Resources Investigations Report↗
USGS1999

User's guide to PHREEQC (Version 2): A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations

PHREEQC version 2 is a computer program written in the C programming language that is designed to perform a wide variety of low-temperature aqueous geochemical calculations. PHREEQC is based on an ion-association aqueous model and has capabilities for (1) speciation and saturation-index calculations; (2) batch-reaction and one-dimensional (1D) transport calculations involving reversible reactions, which include aqueous, mineral, gas, solid-solution, surface-complexation, and ion-exchange equilibria, and irreversible reactions, which include specified mole transfers of reactants, kinetically controlled reactions, mixing of solutions, and temperature changes; and (3) inverse modeling, which finds sets of mineral and gas mole transfers that account for differences in composition between waters, within specified compositional uncertainty limits.New features in PHREEQC version 2 relative to version 1 include capabilities to simulate dispersion (or diffusion) and stagnant zones in 1D-transport calculations, to model kinetic reactions with user-defined rate expressions, to model the formation or dissolution of ideal, multicomponent or nonideal, binary solid solutions, to model fixed-volume gas phases in addition to fixed-pressure gas phases, to allow the number of surface or exchange sites to vary with the dissolution or precipitation of minerals or kinetic reactants, to include isotope mole balances in inverse modeling calculations, to automatically use multiple sets of convergence parameters, to print user-defined quantities to the primary output file and (or) to a file suitable for importation into a spreadsheet, and to define solution compositions in a format more compatible with spreadsheet programs. This report presents the equations that are the basis for chemical equilibrium, kinetic, transport, and inverse-modeling calculations in PHREEQC; describes the input for the program; and presents examples that demonstrate most of the program's capabilities.

Water-Resources Investigations Report↗
USGS1997

Geochemical mole-balance modeling with uncertain data

Geochemical mole-balance models are sets of chemical reactions that quantitatively account for changes in the chemical and isotopic composition of water along a flow path. A revised mole-balance formulation that includes an uncertainty term for each chemical and isotopic datum is derived. The revised formulation is comprised of mole-balance equations for each element or element redox state, alkalinity, electrons, solvent water, and each isotope; a charge-balance equation and an equation that relates the uncertainty terms for p H, alkalinity, and total dissolved inorganic carbon for each aqueous solution; inequality constraints on the size of the uncertainty terms; and inequality constraints on the sign of the mole transfer of reactants. The equations and inequality constraints are solved by a modification of the simplex algorithm combined with an exhaustive search for unique combinations of aqueous solutions and reactants for which the equations and inequality constraints can be solved and the uncertainty terms minimized. Additional algorithms find only the simplest mole-balance models and determine the ranges of mixing fractions for each solution and mole transfers for each reactant that are consistent with specified limits on the uncertainty terms. The revised formulation produces simpler and more robust mole-balance models and allows the significance of mixing fractions and mole transfers to be evaluated. In an example from the central Oklahoma aquifer, inclusion of up to 5% uncertainty in the chemical data can reduce the number of reactants in mole-balance models from seven or more to as few as three, these being cation exchange, dolomite dissolution, and silica precipitation. In another example from the Madison aquifer, inclusion of the charge-balance constraint requires significant increases in the mole transfers of calcite, dolomite, and organic matter, which reduce the estimated maximum carbon 14 age of the sample by about 10,000 years, from 22,700 years to 12,600 years.

Water Resources Research↗
USGS1997

PHREEQCI: A graphical user interface for the geochemical computer program PHREEQC

PhreeqcI is a Windows-based graphical user interface for the geochemical computer program PHREEQC. PhreeqcI provides the capability to generate and edit input data files, run simulations, and view text files containing simulation results, all within the framework of a single interface. PHREEQC is a multipurpose geochemical program that can perform speciation, inverse, reaction-path, and 1D advective reaction-transport modeling. Interactive access to all of the capabilities of PHREEQC is available with PhreeqcI. The interface is written in Visual Basic and will run on personal computers under the Windows(3.1), Windows95, and WindowsNT operating systems.

Water-Resources Investigations Report↗
USGS1996

Ground-water-quality assessment of the central Oklahoma Aquifer, Oklahoma — Geochemical and geohydrologic investigations

Ground-water samples, core samples, and hydrologic measurements were obtained in the Central Oklahoma Aquifer as part of the pilot National Water-Quality Assessment Program. This report examines ground-water recharge and discharge, the potentiometric surface, the chemical and isotopic composition of ground water, and the abundances and textures of minerals in core materials to determine the rates and directions of ground-water flow and the geochemical reactions occurring within the aquifer.

Oklahoma↗
USGS1995

Gross-beta activity in ground water: natural sources and artifacts of sampling and laboratory analysis

Gross-beta activity has been used as an indicator of beta-emitting isotopes in water since at least the early 1950s. Originally designed for detection of radioactive releases from nuclear facilities and weapons tests, analysis of gross-beta activity is widely used in studies of naturally occurring radioactivity in ground water. Analyses of about 800 samples from 5 ground-water regions of the United States provide a basis for evaluating the utility of this measurement. The data suggest that measured gross-beta activities are due to (1) long-lived radionuclides in ground water, and (2) ingrowth of beta-emitting radionuclides during holding times between collection of samples and laboratory measurements. Although 40 K and 228 Ra appear to be the primary sources of beta activity in ground water, the sum of 40 K plus 228 Ra appears to be less than the measured gross-beta activity in most ground-water samples. The difference between the contribution from these radionuclides and gross-beta activity is most pronounced in ground water with gross-beta activities > 10 pCi/L, where these 2 radionuclides account for less than one-half the measured ross-beta activity. One exception is groundwater from the Coastal Plain of New Jersey, where 40 K plus 228 Ra generally contribute most of the gross-beta activity. In contrast, 40 K and 228 Ra generally contribute most of beta activity in ground water with gross-beta activities < 1 pCi/L. The gross-beta technique does not measure all beta activity in ground water. Although 3 H contributes beta activity to some ground water, it is driven from the sample before counting and therefore is not detected by gross-beta measurements. Beta-emitting radionuclides with half-lives shorter than a few days can decay to low values between sampling and counting. Although little is known about concentrations of most short-lived beta-emitting radionuclides in environmental ground water (water unaffected by direct releases from nuclear facilities and weapons tests), their activities are expected to be low. Ingrowth of beta-emitting radionuclides during sample holding times can contribute to gross-beta activity, particularly in ground water with gross-beta activities > 10 pCi/L. Ingrowth of beta-emitting progeny of 238 U, specifically 234 Pa and 234 Th, contributes much of the measured gross-beta activity in ground water from 4 of the 5 areas studied. Consequently, gross-beta activity measurements commonly overestimate the abundance of beta-emitting radionuclides actually present in ground water. Differing sample holding times before analysis lead to differing amounts of ingrowth of the two progeny. Therefore, holding times can affect observed gross-beta measurements, particularly in ground water with 238 U activities that are moderate to high compared with the activity of 40 K plus 228 Ra. Uncertainties associated with counting efficiencies for beta particles with different energies further complicate the interpretation of gross-beta measurements.

Applied Geochemistry↗
USGS1994

An interactive code (NETPATH) for modeling NET geochemical reactions along a flow PATH, version 2.0

NETPATH is an interactive Fortran 77 computer program used to interpret net geochemical mass-balance reactions between an initial and final water along a hydrologic flow path. Alternatively, NETPATH computes the mixing proportions of two to five initial waters and net geochemical reactions that can account for the observed composition of a final water. The program utilizes previously defined chemical and isotopic data for waters from a hydrochemical system. For a set of mineral and (or) gas phases hypothesized to be the reactive phases in the system, NETPATH calculates the mass transfers in every possible combination of the selected phases that accounts for the observed changes in the selected chemical and (or) isotopic compositions observed along the flow path. The calculations are of use in interpreting geochemical reactions, mixing proportions, evaporation and (or) dilution of waters, and mineral mass transfer in the chemical and isotopic evolution of natural and environmental waters. Rayleigh distillation calculations are applied to each mass-balance model that satisfies the constraints to predict carbon, sulfur, nitrogen, and strontium isotopic compositions at the end point, including radiocarbon dating. DB is an interactive Fortran 77 computer program used to enter analytical data into NETPATH, and calculate the distribution of species in aqueous solution. This report describes the types of problems that can be solved, the methods used to solve problems, and the features available in the program to facilitate these solutions. Examples are presented to demonstrate most of the applications and features of NETPATH. The codes DB and NETPATH can be executed in the UNIX or DOS1 environment. This report replaces U.S. Geological Survey Water-Resources Investigations Report 91-4078, by Plummer and others, which described the original release of NETPATH, version 1.0 (dated December, 1991), and documents revisions and enhancements that are included in version 2.0. 1 The use of trade, brand or product names in this report is for identification purposes only and does not constitute endorsement by the U.S. Geological Survey.

Water-Resources Investigations Report↗
USGS1993

Ground-water-quality assessment of the Central Oklahoma Aquifer, Oklahoma: Geochemical and geohydrologic investigations

The National Water-Quality Assessment pilot project for the Central Oklahoma aquifer examined the chemical and isotopic composition of ground water, the abundances and textures of minerals in core samples, and water levels and hydraulic properties in the flow system to identify geochemical reactions occurring in the aquifer and rates and directions of ground-water flow. The aquifer underlies 3,000 square miles of central Oklahoma and consists of Permian red beds, including parts of the Permian Garber Sandstone, Wellington Formation, and Chase, Council Grove, and Admire Groups, and Quaternary alluvium and terrace deposits. In the part of the Garber Sandstone and Wellington Formation that is not confined by the Permian Hennessey Group, calcium, magnesium, and bicarbonate are the dominant ions in ground water; in the confined part of the Garber Sandstone and Wellington Formation and in the Chase, Council Grove, and Admire Groups, sodium and bicarbonate are the dominant ions in ground water. Nearly all of the Central Oklahoma aquifer has an oxic or post-oxic environment as indicated by the large dissolved concentrations of oxygen, nitrate, arsenic(V), chromium(VI), selenium(VI), vanadium, and uranium. Sulfidic and methanic environments are virtually absent. Petrographic textures indicate dolomite, calcite, sodic plagioclase, potassium feldspars, chlorite, rock fragments, and micas are dissolving, and iron oxides, manganese oxides, kaolinite, and quartz are precipitating. Variations in the quantity of exchangeable sodium in clays indicate that cation exchange is occurring within the aquifer. Gypsum may dissolve locally within the aquifer, as indicated by ground water with large concentra-tions of sulfate, but gypsum was not observed in core samples. Rainwater is not a major source for most elements in ground water, but evapotranspiration could cause rainwater to be a significant source of potassium, sulfate, phosphate and nitrogen species. Brines derived from seawater are the most likely source of bromide and chloride in the aquifer. The dominant reaction in recharge is the uptake of carbon dioxide gas from the unsaturated zone (about 2.0 to 4.0 millimoles per liter) and the dissolution of dolomite (about 0.3 to 1.0 millimoles per liter). This reaction generates calcium, magnesium, and bicarbonate water composition. If dolomite does not dissolve to equilibrium, pH values range from 6.0 to 7.3; if dolomite dissolves to equilibrium, pH values are about 7.5. By the time recharge enters the deeper flow system, all ground water is saturated or supersaturated with dolomite and calcite. After carbonate-mineral equilibration has occurred, cation exchange of calcium and magnesium for sodium is the dominant geochemical reaction, which occurs to a substantial extent only in parts of the aquifer. Mass transfers of cation exchange greater than 2.0 millimoles per liter occur in the confined part of the Garber Sandstone and Wellington Formation and in parts of the Chase, Council Grove, and Admire Groups. Associated with cation exchange is dissolution of small quantities of dolomite, calcite, biotite, chlorite, plagioclase, or potassium feldspar, which produces pH values that range from 8.6 to 9.1. Large tritium concentrations indicate ground-water ages of less than about 40 years for most samples of recharge. Carbon-14 ages for samples from the unconfined aquifer generally are less than 10,000 years. Carbon-14 ages of ground water in the confined part of the aquifer range from about 10,000 to 30,000 years or older. These ages produce a time trend in deuterium values that qualitatively is consistent with the timing of the transition from the last glacial maximum to the present interglacial period. The most transmissive geologic units in the Central Oklahoma aquifer are the Garber Sandstone and Wellington Formation and the alluvium and terrace deposits; the Chase, Council Grove, and Admire Groups are less transmissive on the basis of available specific-capacity data. The transmissivities of the Permian geologic units depend largely on the percentage of sandstone; the percentage is greatest in the central part of the aquifer and decreases in all directions from this central part. Because of large mudstone and siltstone contents, the Hennessey Group and the Vanoss Formation are assumed to be confining units above and below the aquifer. The Cimarron and Canadian Rivers are defined to be the northern and southern extent of the aquifer because of decreases in transmissivity beyond the rivers and because there is no indication of ground-water underflow at these rivers. The eastern boundary of the aquifer is the limit of the outcrop of the Chase, Council Grove, and Admire Groups. The presence of brines in the western part of the study unit and below the aquifer indicate the extent of the freshwater flow system in these directions. Regional ground-water flow is west to east; the Deep Fork is a major discharge area for the regional flow system. Local flow systems are present within the unconfined part of the study unit. Most streams are gaining streams, and very few losing streams are evident. Median values of aquifer properties were estimated as follows: recharge to the saturated zone, 1.6 inches per year; evapotranspiration of water that never reaches the saturated zone, 25 to 30 inches per year; porosity, 0.22; storage coefficient, 0.0002; transmissivity, 260 to 450 feet squared per day; horizontal hydraulic conductivity, 4.5 feet per day; and the ratio of horizontal to vertical hydraulic conductivity, 10,000. Reported ground-water withdrawals peaked in 1985 at 13,900 million gallons but had decreased to 7,860 million gallons by 1989. Unreported domestic withdrawals were estimated to be 1,685 million gallons in 1980. The flow system in the aquifer can be considered to have three major components: (1) A shallow, local flow system in the unconfined part of the aquifer, (2) a deep, regional flow system in the unconfined part of the aquifer, and (3) a deep, regional flow system in the confined part of the aquifer. In the shallow, local flow system, water flows relatively quickly along short flowlines from the point of recharge to the point of discharge at the nearest stream. Many water samples from shallow wells contain large concentrations of tritium, which indicate ground-water ages of less than 40 years. In the deep, regional flow system in the unconfined part of the aquifer, water takes more time to flow along longer flowlines than in the shallow, local flow system. Much of the water in this flow system is recharged along ridges that correspond to ground-water divides between drainage basins. Transit times for water recharging the aquifer along ridges is greater than 5,000 years, computed using a numerical flow model in conjunction with a particle-tracking model. The deep, regional flow system in the confined part of the Garber Sandstone and Wellington Formation is recharged from a small part of the outcrop area of the Garber Sandstone. From the recharge area, water flows west under the confining unit to discharge to streams as far away as the Cimarron River. Flowpaths are relatively long, as much as 50 miles. The transit times in this flow system range from thousands to tens of thousands of years. The long-term hydrogeochemical process occurring in the Central Oklahoma aquifer is removal of unstable minerals, including dolomite, calcite, biotite, chlorite, and feldspars, and the replacement of exchangeable sodium on clays with calcium and magnesium. Over geologic time, the flux of water through the rapidly moving, local flow system has been sufficient to remove most of the dolomite, calcite, and exchangeable sodium. In places, chlorite and feldspars have been removed. In the deep, regional flow system of the unconfined part of the Garber Sandstone and Wellington Formation, the flux of water has been sufficient to remove most of the exchangeable sodium, but carbonate minerals remain sufficiently abundant to maintain dolomite and calcite equilibrium. In the confined part of the Garber Sandstone and Wellington Formation and in the less transmissive parts of the unconfined aquifer, including the Chase, Council Grove, and Admire Groups, ground-water flow is slowest, and the flux of water and extent of reaction have been insufficient to remove either the carbonate minerals or the exchangeable sodium on clays.

Oklahoma↗
USGS1990

Geohydrology and water quality of the Roubidoux Aquifer, northeastern Oklahoma

The Roubidoux aquifer is an important source of freshwater for public supplies, commerce, industry, and rural water districts in northeastern Oklahoma. Ground-water withdrawals from the aquifer in 1981 were estimated to be 4.8 million gallons per day, of which about 90 percent was withdrawn in Ottawa County. Wells drilled at the beginning of the 20th century originally flowed at the land surface, but in 1981 water levels ranged from 22 to 471 feet below land surface. A large cone of depression has formed as a result of ground water withdrawals near Miami. Wells completed in the Roubidoux aquifer have yields that range from about 100 to more than 1,000 gallons per minute. An aquifer test and a digital ground-water flow model were used to estimate aquifer and confining-layer hydraulic characteristics. Using these methods, the transmissivity of the aquifer was estimated to be within a range of 400 to 700 square feet per day. The leakance of the confining layer was determined to be within a range from 0 to 0.13 per day, with a best estimate value in a range from 4.3 x 10-8 to 7.7 x 10-8 per day. Analyses of water samples collected as part of this study and of water-quality data from earlier work indicate that a large areal change in major-ion chemistry occurs in ground water in the Roubidoux aquifer in northeastern Oklahoma. The ground water in the easternmost part of the study unit has relatively small dissolved-solids concentrations (less than 200 milligrams per liter) with calcium, magnesium, and bicarbonate as the major ions. Ground water in the westernmost part of the study unit has relatively large dissolved-solids concentrations (greater than 800 milligrams per liter) with sodium and chloride as the major ions. A transition zone of intermediate sodium, chloride, and dissolved-solids concentrations exists between the easternmost and westernmost parts of the study unit. Three water-quality problems are apparent in the Roubidoux aquifer in northeast Oklahoma: (1) Contamination by mine water, (2) large concentrations of sodium and chloride, and (3) large radium-226 concentrations. Many wells in the mining area have been affected by mine-water contamination. At present (1990), all instances of ground-water contamination by mine water can be explained by faulty seals or leaky casings in wells that pass through the zone of mine workings and down to the Roubidoux aquifer. None of the data available to date demonstrate that mine water has migrated from the Boone Formation through the pores and fractures of the intervening geologic units to the Roubidoux aquifer. Ground water with large concentrations of sodium and chloride occurs at some depth throughout the study unit. In the eastern part of the study unit, chloride concentrations greater than 250 milligrams per liter are found at depths greater than approximately 1,200 to 1,500 feet. Data are too few to determine the depth to ground water with large concentrations of sodium and chloride in the southern and southwestern parts of the study unit. Large concentrations of gross-alpha radioactivity in ground water occur near the western edge of the transition zone. Generally, ground water with large concentrations of gross-alpha radioactivity was found to exceed the maximum contaminant level for radium-226. (available as photostat copy only)

Oklahoma↗
USGS1989

Ground-water quality assessment of the central Oklahoma aquifer, Oklahoma - Analysis of available water-quality data through 1987

Beginning in 1986, the Congress annually has appropriated funds for the U.S. Geological Survey to test and refine concepts for a National Water-Quality Assessment (NAWQA) Program. The long-term goals of a full-scale program would be to: Provide a nationally consistent description of current water-quality conditions for a large part of the Nation's surface- and ground-water resources; Define long-term trends (or lack of trends) in water quality; and Identify, describe, and explain, as possible, the major factors that affect the observed water-quality conditions and trends. The results of the NAWQA Program will be made available to water managers, policy makers, and the public, and will provide an improved scientific basis for evaluating the effectiveness of water-quality management programs. At present (1988), the assessment program is in a pilot phase in seven project areas throughout the country that represent diverse hydrologic environments and water-quality conditions. The Central Oklahoma aquifer project is one of three pilot ground-water projects. One of the initial activities performed by each pilot project was to compile, screen, and interpret the large amount of water-quality data available within each study unit. The purpose of this report is to assess the water quality of the Central Oklahoma aquifer using the information available through 1987. The scope of the work includes compiling data from Federal, State, and local agencies; evaluating the suitability of the information for conducting a regional water-quality assessment; mapping regional variations in major-ion chemistry; calculating summary statistics of the available water-quality data; producing maps to show the location and number of samples that exceeded water-quality standards; and performing contingency-table analyses to determine the relation of geologic unit and depth to the occurrence of chemical constituents that exceed water-quality standards. This report provides an initial description of water-quality conditions in the Central Oklahoma aquifer study unit. No attempt was made in this report to determine the causes for regional variations in major-ion chemistry or to examine the reasons that some chemical constituents exceed water-quality standards.

Oklahoma↗
USGS1988

Chemical analyses of stream sediment in the Tar Creek basin of the Picher mining area, northeast Oklahoma

Chemical analyses are presented for 47 sediment samples from the Tar Creek drainage in the Picher mining area of northeast Oklahoma. The samples were taken in December 1983, June 1984, and June 1985. All of the samples were taken downstream from mine-water discharge points of abandoned lead and zinc mines. The 34 samples taken in December 1983 and June 1984 were analyzed semiquantitatively by emission spectrography for 64 elements and quantitatively for cadmium, copper, iron, manganese, nickel, lead, sulfur, zinc, and organic carbon. The 13 samples taken in June 1985 were analyzed quantitatively for aluminum, cadmium, cobalt, chromium, copper, iron, manganese, molybdenum, nickel, phosphorus, lead, sulfur, silicon, titanium, vanadium, zinc, and organic carbon.

Open-File Report↗
USGS1987

Chemical analyses of water samples from the Picher mining area, northeast Oklahoma and southeast Kansas

Chemical analyses are presented for 169 water samples from Tar Creek drainage and the Picher lead-zinc mining area of northeast Oklahoma and southeast Kansas. Water samples were taken from November 1983 through February 1986 from the abandoned mines, from points of mine-water discharge, and from surface-water locations upstream and downstream from mine discharge area. The pH, temperature, alkalinity, dissolved oxygen, and specific conductance were measured in the field. Laboratory analyses routinely included the major ions plus aluminum, cadmium, copper, iron, lead, manganese, nickel, and zinc. Non-routine analyses of dissolved gases and tritium are presented. Stable carbon-isotope ratios for 11 mine-water samples and three carbonate-rock samples are reported. Miscellaneous stream-discharge measurements made at the time of sampling or taken from gaging-station records are included in the report.

Open-File Report↗