USGS Science⌕ Search

SEARCH · USGS Science

Results for “Cells”

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 955 records · Page 53Linked to original sources

Simulation of ground-water flow in the Intermediate and Floridan aquifer systems in Peninsular Florida

A numerical model of the intermediate and Floridan aquifer systems in peninsular Florida was used to (1) test and refine the conceptual understanding of the regional ground-water flow system; (2) develop a data base to support subregional ground-water flow modeling; and (3) evaluate effects of projected 2020 ground-water withdrawals on ground-water levels. The four-layer model was based on the computer code MODFLOW-96, developed by the U.S. Geological Survey. The top layer consists of specified-head cells simulating the surficial aquifer system as a source-sink layer. The second layer simulates the intermediate aquifer system in southwest Florida and the intermediate confining unit where it is present. The third and fourth layers simulate the Upper and Lower Floridan aquifers, respectively. Steady-state ground-water flow conditions were approximated for time-averaged hydrologic conditions from August 1993 through July 1994 (1993-94). This period was selected based on data from Upper Floridan a quifer wells equipped with continuous water-level recorders. The grid used for the ground-water flow model was uniform and composed of square 5,000-foot cells, with 210 columns and 300 rows.

Water-Resources Investigations Report↗

Analytical and numerical simulation of the steady-state hydrologic effects of mining aggregate in hypothetical sand-and-gravel and fractured crystalline-rock aquifers

Analytical solutions and numerical models were used to predict the extent of steady-state drawdown caused by mining of aggregate below the water table in hypothetical sand-and-gravel and fractured crystalline-rock aquifers representative of hydrogeologic settings in the Front Range area of Colorado. Analytical solutions were used to predict the extent of drawdown under a wide range of hydrologic and mining conditions that assume aquifer homogeneity, isotropy, and infinite extent. Numerical ground-water flow models were used to estimate the extent of drawdown under conditions that consider heterogeneity, anisotropy, and hydrologic boundaries and to simulate complex or unusual conditions not readily simulated using analytical solutions. Analytical simulations indicated that the drawdown radius (or distance) of influence increased as horizontal hydraulic conductivity of the aquifer, mine penetration of the water table, and mine radius increased; radius of influence decreased as aquifer recharge increased. Sensitivity analysis of analytical simulations under intermediate conditions in sand-and-gravel and fractured crystalline-rock aquifers indicated that the drawdown radius of influence was most sensitive to mine penetration of the water table and least sensitive to mine radius. Radius of influence was equally sensitive to changes in horizontal hydraulic conductivity and recharge. Numerical simulations of pits in sand-and- gravel aquifers indicated that the area of influence in a vertically anisotropic sand-and-gravel aquifer of medium size was nearly identical to that in an isotropic aquifer of the same size. Simulated area of influence increased as aquifer size increased and aquifer boundaries were farther away from the pit, and simulated drawdown was greater near the pit when aquifer boundaries were close to the pit. Pits simulated as lined with slurry walls caused mounding to occur upgradient from the pits and drawdown to occur downgradient from the pits. Pits simulated as refilled with water and undergoing evaporative losses had little hydro- logic effect on the aquifer. Numerical sensitivity analyses for simulations of pits in sand-and-gravel aquifers indicated that simulated head was most sensitive to horizontal hydraulic conductivity and the hydraulic conductance of general-head boundaries in the models. Simulated head was less sensitive to riverbed conductance and recharge and relatively insensitive to vertical hydraulic conductivity. Numerical simulations of quarries in fractured crystalline-rock aquifers indicated that the area of influence in a horizontally anisotropic aquifer was elongated in the direction of higher horizontal hydraulic conductivity and shortened in the direction of lower horizontal hydraulic conductivity compared to area of influence in a homogeneous, isotropic aquifer. Area of influence was larger in an aquifer with ground-water flow in deep, low-permeability fractures than in a homogeneous, isotropic aquifer. Area of influence was larger for a quarry intersected by a hydraulically conductive fault zone and smaller for a quarry intersected by a low-conductivity fault zone. Numerical sensitivity analyses for simulations of quarries in fractured crystalline-rock aquifers indicated simulated head was most sensitive to variations in recharge and horizontal hydraulic conductivity, had little sensitivity to vertical hydraulic conductivity and drain cells used to simulate valleys, and was relatively insensitive to drain cells used to simulate the quarry.

Water-Resources Investigations Report↗

Evaluation of the Source and Transport of High Nitrate Concentrations in Ground Water, Warren Subbasin, California

Ground water historically has been the sole source of water supply for the Town of Yucca Valley in the Warren subbasin of the Morongo ground-water basin, California. An imbalance between ground-water recharge and pumpage caused ground-water levels in the subbasin to decline by as much as 300 feet from the late 1940s through 1994. In response, the local water district, Hi-Desert Water District, instituted an artificial recharge program in February 1995 using imported surface water to replenish the ground water. The artificial recharge program resulted in water-level recoveries of as much as 250 feet in the vicinity of the recharge ponds between February 1995 and December 2001; however, nitrate concentrations in some wells also increased from a background concentration of 10 milligrams per liter to more than the U.S. Environmental Protection Agency (USEPA) maximum contaminant level (MCL) of 44 milligrams per liter (10 milligrams per liter as nitrogen). The objectives of this study were to: (1) evaluate the sources of the high-nitrate concentrations that occurred after the start of the artificial-recharge program, (2) develop a ground-water flow and solute-transport model to better understand the source and transport of nitrates in the aquifer system, and (3) utilize the calibrated models to evaluate the possible effect of a proposed conjunctive-use project. These objectives were accomplished by collecting water-level and water-quality data for the subbasin and assessing changes that have occurred since artificial recharge began. Collected data were used to calibrate the ground-water flow and solute-transport models. Data collected for this study indicate that the areal extent of the water-bearing deposits is much smaller (about 5.5 square miles versus 19 square miles) than that of the subbasin. These water-bearing deposits are referred to in this report as the Warren ground-water basin. Faults separate the ground-water basin into five hydrogeologic units: the west, the midwest, the mideast, the east and the northeast hydrogeologic units. Water-quality analyses indicate that septage from septic tanks is the primary source of the high-nitrate concentrations measured in the Warren ground-water basin. Water-quality and stable-isotope data, collected after the start of the artificial recharge program, indicate that mixing occurs between imported water and native ground water, with the highest recorded nitrate concentrations in the midwest and the mideast hydrogeologic units. In general, the timing of the increase in measured nitrate concentrations in the midwest hydrogeologic unit is directly related to the distance of the monitoring well from a recharge site, indicating that the increase in nitrate concentrations is related to the artificial recharge program. Nitrate-to-chloride and nitrogen-isotope data indicate that septage is the source of the measured increase in nitrate concentrations in the midwest and the mideast hydrogeologic units. Samples from four wells in the Warren ground-water basin were analyzed for caffeine and selected human pharmaceutical products; these analyses suggest that septage is reaching the water table. There are two possible conceptual models that explain how high-nitrate septage reaches the water table: (1) the continued downward migration of septage through the unsaturated zone to the water table and (2) rising water levels, a result of the artificial recharge program, entraining septage in the unsaturated zone. The observations that nitrate concentrations increase in ground-water samples from wells soon after the start of the artificial recharge program in 1995 and that the largest increase in nitrate concentrations occur in the midwest and mideast hydrogeologic units where the largest increase in water levels occur indicate the validity of the second conceptual model (rising water levels). The potential nitrate concentration resulting from a water-level rise in the midwest and mideast hydrogeologic units was estimated using a simple mixing-cell model. The estimated value is within the range of concentrations measured in samples from wells, further indicating the validity of the second conceptual model. A ground-water flow model and a solute-transport model were developed for the Warren ground-water basin for the period 1956-2001. MODFLOW-96 was used for the ground-water flow model and MOC3D was used for the solute-transport model. The model cell size is about 500 feet by 500 feet and the models were discretized vertically into three layers. The models were calibrated using a trial-and-error approach using water-level and nitrate-concentration data collected between 1956-2001. In order to better match the measured data, low fault hydraulic characteristic values were required, thereby compartmentalizing the ground-water basin. In addition, it was necessary to parameterize the specific yield distribution for the top model layer where unconfined ground-water conditions occur into three homogeneous zones. Separate sets of specific- yield values were needed to simulate the drawdown and subsequent water-level recovery. In addition, the calibrated natural recharge was about 83 acre-feet per year. The entrainment of unsaturated-zone septage was simulated as recharge having an associated nitrate concentration. The volume of recharge was a function of the measured water-level rise between 1994-98 and the moisture content of the unsaturated zone. The nitrate concentration of the recharge water was a weighted function of the assumed nitrate concentration in the infiltrating water associated with the overlying land use. The simulated hydraulic head and nitrate concentration results were in good agreement with the measured data indicating that the mechanism for the increase in nitrate concentrations was rising water levels entraining high-nitrate septage in the unsaturated-zone. The calibrated models were used to simulate the possible effects of a planned conjunctive-use project in the western part of the ground-water basin. The simulated project included the addition of a new recharge pond and a new extraction well. In addition, recharge at two existing recharge ponds was increased and three existing production wells were pumped, treated in a nitrate-removal facility, and used for water supply. The simulated hydraulic heads increased in the west, the mideast, and parts of the east hydrogeologic units; however, the simulated hydraulic heads decreased in the midwest and northeast hydrogeologic units. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter (10 milligrams per liter as nitrogen) in parts of the west as a result of the increase in simulated hydraulic head. The simulated nitrate concentrations decreased in part of the midwest hydrogeologic unit as a result of the artificial recharge and pumping from the nitrate-removal wells. The simulated nitrate concentrations increased to above the MCL of 44 milligrams per liter in part of the mideast and parts of the east hydrogeologic units beneath commercial land-use areas.

California↗

Recalibration of a ground-water flow model of the Mississippi River Valley alluvial aquifer of northeastern Arkansas, 1918-1998, with simulations of water levels caused by projected ground-water withdrawals through 2049

A digital model of the Mississippi River Valley alluvial aquifer in eastern Arkansas was used to simulate ground-water flow for the period from 1918 to 2049. The model results were used to evaluate effects on water levels caused by demand for ground water from the alluvial aquifer, which has increased steadily for the last 40 years. The model results showed that water currently (1998) is being withdrawn from the aquifer at rates greater than what can be sustained for the long term. The saturated thickness of the alluvial aquifer has been reduced in some areas resulting in dry wells, degraded water quality, decreased water availability, increased pumping costs, and lower well yields. The model simulated the aquifer from a line just north of the Arkansas-Missouri border to south of the Arkansas River and on the east from the Mississippi River westward to the less permeable geologic units of Paleozoic age. The model consists of 2 layers, a grid of 184 rows by 156 columns, and comprises 14,118 active cells each measuring 1 mile on a side. It simulates time periods from 1918 to 1998 along with further time periods to 2049 testing different pumping scenarios. Model flux boundary conditions were specified for rivers, general head boundaries along parts of the western side of the model and parts of Crowleys Ridge, and a specified head boundary across the aquifer further north in Missouri. Model calibration was conducted for observed water levels for the years 1972, 1982, 1992, and 1998. The average absolute residual was 4.69 feet and the root-mean square error was 6.04 feet for the hydraulic head observations for 1998. Hydraulic-conductivity values obtained during the calibration process were 230 feet per day for the upper layer and ranged from 230 to 730 feet per day for the lower layer with the maximum mean for the combined aquifer of 480 feet per day. Specific yield values were 0.30 throughout the model and specific storage values were 0.000001 inverse-feet throughout the model. Areally specified recharge rates ranged from 0 to about 30 inches and total recharge increased from 1972 to 1998 by a factor of about four. Water levels caused by projected ground-water withdrawals were simulated using the calibrated model. Simulations represented a period of 50 years into the future in three scenarios with either unchanged pumpage, pumpage increased by historic trends, or pumpage increased by historic trends except in two areas of the Grand Prairie. If pumping remains at 1997 rates, this produces extreme water-level declines (areas where model cells have gone dry or where the water level in the aquifer is equal to or less than the original saturated thickness, assuming confined conditions in the aquifer everywhere in the formation in predevelopment times) in the aquifer in two areas of the aquifer (one in the Grand Prairie area between the Arkansas and White Rivers and the other west of Crowleys Ridge along the Cache River) with about 400 square miles going dry. Increasing the pumping rates to that which would be projected using historic data led to increased extreme water-level declines in both areas with about 1,300 square miles going dry. Declines in both scenarios generally occurred most rapidly between 2009 and 2019. Reducing the pumping rates to 90 percent of that used for projected historic rates in areas between the Arkansas and White Rivers relating to two diversion projects of the U.S. Army Corps of Engineers and other agencies did little to decrease the extreme water-level declines. However, these pumpage reductions are small (amounting to about 16 percent of the reductions that could result from implementation of these diversion projects).

Water-Resources Investigations Report↗

Biogeochemical processes that produce dissolved organic matter from wheat straw

The chemical reactions that lead to the formation of dissolved organic matter (DOM) in natural waters are poorly understood. Studies on the formation of DOM generally are complicated because almost all DOM isolates have been derived from mixtures of plant species composed of a wide variety of different types of precursor compounds for DOM formation. This report describes a study of DOM derived mainly from bales of wheat straw that had been left in a field for several years. During this period of time, black water from the decomposing wheat straw accumulated in pools in the field. The nuclear magnetic resonance and infrared spectra of the black water DOM indicate that it is composed almost entirely of lignin and carbohydrate polymeric units. Analysis by high-performance size-exclusion chromatography with multi-angle laser-light scattering detection indicates that the number average molecular weight of the DOM is 124,000 daltons. The results presented in this report indicate that the black water DOM is composed of hemicellulose chains cross-linked to lignin oligomers. These types of structures have been shown to exist in the hemicellulose matrix of plant cell walls. The cross-linked lignin-hemicellulose complexes apparently were released from partially degraded wheat-straw cell walls with little alteration. In solution in the black water, these lignin-hemicellulose polymers fold into compact globular particles in which the nonpolar parts of the polymer form the interiors of the particles and the polar groups are on the exterior surfaces of the particles. The tightly folded, compact conformation of these particles probably renders them relatively resistant to microbial degradation. This should be especially the case for the aromatic lignin structures that will be buried in the interiors of the particles.

Water-Resources Investigations Report↗

Conjunctive-use optimization model of the Mississippi River Valley alluvial aquifer of Southeastern Arkansas

The Mississippi River Valley alluvial aquifer is a water-bearing assemblage of gravels and sands that underlies about 32,000 square miles of Missouri, Kentucky, Tennessee, Mississippi, Louisiana, and Arkansas. Because of the heavy demands placed on the aquifer, several large cones of depression have formed in the potentiometric surface, resulting in lower well yields and degraded water quality in some areas. A ground-water flow model of the alluvial aquifer was previously developed for an area covering 3,826 square miles, extending south from the Arkansas River into the southeastern corner of Arkansas, parts of northeastern Louisiana, and western Mississippi. The flow-model results indicated that continued ground-water withdrawals at rates commensurate with those of 1997 could not be sustained indefinitely without causing water levels to decline below half the original saturated thickness of the aquifer. Conjunctive-use optimization modeling was applied to the flow model of the alluvial aquifer to develop withdrawal rates that could be sustained relative to the constraints of critical ground-water area designation. These withdrawal rates form the basis for estimates of sustainable yield from the alluvial aquifer and from rivers specified within the alluvial aquifer model. A management problem was formulated as one of maximizing the sustainable yield from all ground-water and surface-water withdrawal cells within limits imposed by plausible withdrawal rates, and within specified constraints involving hydraulic head and streamflow. Steady-state conditions were selected because the maximized withdrawals are intended to represent sustainable yield of the system (a rate that can be maintained indefinitely).One point along the Arkansas River and one point along Bayou Bartholomew were specified for obtaining surface-water sustainable-yield estimates within the optimization model. Streamflow constraints were specified at two river cells based on average 7-day low flows with 10-year recurrence intervals. Sustainable-yield estimates were affected by the allowable upper limit on withdrawals from wells specified in the optimization model. Withdrawal rates were allowed to increase to 200 percent of the withdrawal rate in 1997. As the overall upper limit is increased, the sustainable yield generally increases. Tests with the optimization model show that without limits on pumping, wells adjacent to sources of water, such as large rivers, would have optimal withdrawal rates that were orders of magnitude larger than rates corresponding to those of 1997. Specifying an upper withdrawal limit of 100 percent of the 1997 withdrawal rate, the sustainable yield from ground water for the entire study area is 70.3 million cubic feet per day, which is about 96 percent of the amount withdrawn in 1997 (73.5 million cubic feet per day). If the upper withdrawal limit is increased to 150 percent of the 1997 withdrawal rate, the sustainable yield from ground water for the entire study area is 80.6 million cubic feet per day, which is about 110 percent of the amount withdrawn in 1997. If the upper withdrawal limit is increased to 200 percent of the 1997 withdrawal rate, the sustainable yield from ground water for the entire study area is 110.2 million cubic feet per day, which is about 150 percent of the amount withdrawn in 1997. Total sustainable yield from the Arkansas River and Bayou Bartholomew is about 4,900 million cubic feet per day, or about 6,700 percent of the amount of ground-water withdrawn in 1997. The large, sustainable yields from surface water represent a potential source of water that could supplement ground water and meet the total water demand. Unmet demand (defined as the difference between the optimized withdrawal rate or sustainable yield, and the anticipated demand) was calculated using different demand rates based on multiples of the 1997-withdrawal rate. Assuming that demand is the 1997 withdrawal rate, and that sustainable-

Water-Resources Investigations Report↗

Evaluation of the Feasibility of Freshwater Injection Wells in Mitigating Ground-Water Quality Degradation at Selected Well Fields in Duval County, Florida

The Fernandina permeable zone contains brackish water in parts of Duval County, Florida. Upward flow from the Fernandina permeable zone to the upper zone of the Lower Floridan aquifer increases chloride concentrations in ground water in parts of Duval County. Numerical models of the ground-water flow system in parts of Duval, St. Johns, and Clay Counties, Florida, were used to (1) estimate the vertical flows between the low-quality water of the Fernandina permeable zone and the high-quality water of the upper zone of the Lower Floridan aquifer in the vicinity of Deerwood 3 and Brierwood well fields, based on 2000 ground-water withdrawal rates; (2) determine how such vertical flows change as several scenarios of injection, withdrawal, and intervening rest periods are simulated in the two well fields; and (3) evaluate the effects of changes in less certain hydraulic parameters on the vertical flows between the Fernandina permeable zone and the upper zone of the Lower Floridan aquifer. The ground-water flow system was simulated with a four-layer model using MODFLOW-2000, which was developed by the U.S. Geological Survey. The first layer consists of specified-head cells simulating the surficial aquifer system with prescribed water levels. The second layer simulates the Upper Floridan aquifer. The third and fourth layers simulate the upper zone of the Lower Floridan aquifer and the Fernandina permeable zone, respectively. Average flow conditions in 2000 were approximated with a steady-state simulation. The changes in upward flow from the Fernandina permeable zone due to periods of injections and withdrawals were analyzed with transient simulations. The grid used for the ground-water flow model was uniform and composed of square 250-foot cells, with 400 columns and 400 rows. The active model area encompasses about 360 square miles in parts of Duval, St. Johns, and Clay Counties, Florida. Ground-water flow simulation was limited vertically to the bottom of the Fernandina permeable zone. The steady-state ground-water flow model was calibrated using time-averaged 2000 heads at 20 control points. Environmental-water heads in the Fernandina permeable zone were calculated for wells with variable water density. Transmissivity of the Upper Floridan aquifer, the upper zone of the Lower Floridan aquifer, and the Fernandina permeable zone, and the leakance of the intermediate confining unit, the middle semiconfining unit, and the semiconfining unit were obtained from regional ground-water flow models and adjusted until a reasonable fit between simulated and computed heads was obtained. The calibrated hydraulic properties from the steady-state ground-water flow model, and the calibrated storage coefficient from the transient model, were used to simulate hypothetical transient scenarios of injection, withdrawal, and intervening rest periods to assess changes in flow between the Fernandina permeable zone and the upper zone of the Lower Floridan aquifer. Based on the simulated flows between the Fernandina permeable zone and the upper zone of the Lower Floridan aquifer and the 18 million gallons per day of water available for injection, the reversal of the prevailing upward flow from the Fernandina permeable zone was not achieved. However, steady-state and transient simulations indicate that the upward flow of water from the Fernandina permeable zone could be reduced by as much as 64 percent, from 0.11 to 0.04 cubic foot per second, if only injection periods are simulated.

Water-Resources Investigations Report↗

Limnological study of Lake Shastina, Siskiyou County, California

Lake Shastina provides water for irrigation in Shasta Valley, as well as recreation. Presently, its shoreline is being developed for summer homes. Surface water constituted more than 90 percent of the approximately 50,000 acre-foot (62-cubic hectometre) inflow to Lake Shastina in the 1972 water year. Controlled outflow is via the Montague Main Canal; however, leakage from the lake through volcanic rocks to the northwest was estimated to be greater than the measured outflow. Appreciable annual changes in the quantity of water in storage in the lake are related mainly to variations in annual inflow. From June through August the lake was thermally stratified. In the spring and summer the epilimnion was often supersaturated with oxygen, while at the same time the hypolimnion was undersaturated and 'often devoid of dissolved oxygen. Vertical stratification of carbon dioxide, carbonate, bicarbonate, hydrogen ion, nitrogen, and phosphorus was also recorded during the spring and summer. Orthophosphate, total phosphorus, and total nitrogen concentrations (organic, ammonium, and nitrate) were highest in the hypolimnion during the period of thermal stratification. Ten-inch (25-centimetre) core samples from the reservoir bottom were chemically analyzed at 0.8-inch (2-centimetre) intervals. The concentrations ranged from 6.3 to 28.9 milligrams per gram of iron, 0.07 to 0.43 milligrams per gram of manganese, 0.4 to 2.7 milligrams per gram of organic nitrogen plus ammonium, and 0.06 to 1.3 milligrams per gram of total phosphorus. Organic matter in the cores ranged from 4 to 14 percent. Green algae and diatoms were the dominant algal types, reaching maximum concentrations of about 7 and 30 million cells per litre, respectively. These phytoplankton occurred near the surface during thermally stratified periods, but were distributed at greater depths during nonthermally stratified periods. Blue-green algae were present only in the spring samples, and reached a maximum concentration of about 5 million cells per litre. Zooplankton numbers were greatest in March, July, and September, with lesser concentrations in June. Three major zooplankton groups, Cladocera, Copepoda, and Rotifera, were present. The major groups of benthic organisms were Oligochaeta, Chironomidae, and Chaoborus, with numbers ranging from 3350, 890, and 8450 per square metre, respectively. A discussion on algal control is included.

California↗

The effects of using ground water to maintain water levels of Cedar Lake, Wisconsin

Cedar Lake, a kettle lake with no surface inlet or outlet, was studied to evaluate the feasibility of maintaining water levels of lakes in the glaciated kettle moraine area of eastern Wisconsin by pumping ground water into them. The general hydrogeology of the area around the lake was defined and the hydrology of the lake was studied in detail. There are two aquifer systems in the Cedar Lake area, and in general, throughout eastern Wisconsin. The shallow aquifer system is composed of unconsolidated glacial drift of Quaternary age and consolidated bedrock of Silurian age. Cedar Lake is located in a local recharge area within this system. The saturated thicknesses of the glacial drift and Silurian dolomite beneath the lake are approximately 100 feet and 500 feet, respectively. This aquifer system is separated from an underlying aquifer system by a layer of shale approximately '200 feet thick that forms an effective barrier to water movement. The deep aquifer system is composed of dolomites and sandstones of Ordovician and Cambrian age whose composite thickness is 600 to 650 feet. A water budget was prepared to quantify the components of the hydrologic system of Cedar Lake. Inflow to the lake averaged approximately ^5 inches per water year for the 2-year period October 197^ through September 1976. Inflow consisted of precipitation on the lake surface (28 inches), groundwater seepage to the lake (2 inches), and overland flow (15 inches) to the lake. Outflow averaged approximately U9 inches per year for the same 2- year period. Outflow consisted of evaporation (31 inches) from the lake surface and ground-water seepage from the lake (18 inches). A volume of water equivalent to that needed to raise the lake level 1*7 inches was pumped from the shallow aquifer system into Cedar Lake between February 1 and September 30, 1977. The water budget for the lake during this period indicated that approximately 90 percent of pumped water was either recycled from the lake to the well or otherwise lost as seepage from the lake. Selected physical and chemical characteristics of the lake water were measured periodically between October 197^ and September 1977- The water temperature ranged between 26°C in August to approximately 2°C in December during each year of the study. The dissolved-oxygen concentration of the water ranged from 13 milligrams per liter at the time ice cover first developed in December to h milligrams per liter at the time ice cover broke up the following April. The total nitrogen concentration of the water varied from 1.0 milligrams per liter during winter to 0.5 milligrams per liter during summer. The total phosphorus concentration of the water was fairly uniform throughout the year, averaging 0.01 milligrams per liter. The phytoplankton population was sampled periodically. The total phytoplankton count was small during winter, generally less than 250 cells per milliliter. These phytoplankton were composed of euglenoids, diatoms, green algae, and blue-green algae. The count increased to 50,000 to 100,000 cells per milliliter during summer. One or two species of blue-green algae dominated during these times. There were no identifiable changes in measured physical and chemical characteristics of lake water during sustained pumping of ground water into the lake, nor were there identifiable changes in the number or makeup of the phytoplankton community. Differences in physical and chemical characteristics of lake water and ground water added to the lake probably were not great enough to cause changes within the lake.

Wisconsin↗

Algal conditions in the Caloosahatchee River (1975-79), Lake Okeechobee to Franklin Lock, Florida

Maximum numbers of suspended algae occurred in late spring and early summer, in each of the years 1975-79, in the Caloosahatchee River. Numbers exceeded 100,000 cells per milliliter at all stations sometime during the study. Concentrations decreased during late summer and autumn and were low during winter, except in January 1979 when numbers at most sites exceeded 100,000 cells per milliliter. The January 1979 bloom coincided with large discharges from Lake Okeechobee. During previous winters, discharges and algal numbers were lower. During other seasons, algal blooms occurred most frequently under low-flow or stagnant conditions. The upstream site at Moore Haven, which had the least discharge and was most stagnant, had consistently higher algal concentrations than downstream sites. Blue-green algae were dominant in the river during the summer at the upstream site throughout the year. The percentage of blue-green algae decreased downstream. Concentrations of nitrite plus nitrate nitrogen were inversely correlated with concentrations of algae and decreased to near zero during algal blooms. The low concentrations of these forms of inorganic nitrogen relative to other major nutrients probably favor blue-green algae and limit growth of other algae. Contributions by the basin tributaries to the nutritive condition of the river were small because concentrations of nutrients, algal growth potential, and algae in the tributaries were generally less than those in the river. (USGS)

Florida↗

Ground-water flow in the Prairie du Chien-Jordan aquifer related to contamination by coal-tar derivatives, St. Louis Park, Minnesota

A three-dimensional, ground-water-flow model of the Prairie du Chien-Jordan aquifer and associated hydrogeologic units was developed to evaluate the movement of coal-tar derivatives from a coal-tar distillation and wood-preserving plant in St. Louis Park, Minnesota. A finite-difference grid was superimposed on the modeled area, which includes most of eastern Hennepin County. The individual cells are 400-foot squares in the center of the grid (St. Louis Park area); the cells increase in dimension toward the outside limits of the grid. Five geologic units are represented by four layers in the model. These units include the Jordan Sandstone, the Prairie du Chien Group (dolomite and sandy dolomite), the basal confining unit of the St. Peter Sandstone (silty and sandy shale), the St. Peter Sandstone, and glacial deposits in bedrock valleys. The model was calibrated for steady-state conditions for a period before significant ground-water development (1885-1930) and for a period of significant pumping stress (winter conditions, 1970 fs). A transient calibration was accomplished by simulation of a period during which seasonal changes in potentiometric head in the Prairie du Chien-Jordan aquifer were significant (1977-80). Sensitivity testing indicated that leakage to the upper model layer and the vertical hydraulic conductivity of the basal confining unit of the St. Peter Sandstone were the model hydrologic properties which, when changed, resulted in the greatest changes in model-calculated water levels. The calibrated model generally calculates water levels that are within 10 feet of measured values. Model simulations indicate that the potentiometric surface of the Prairie du Chien-Jordan aquifer would be raised by as much as 3 feet in the area of the plant site by water introduced into the aquifer through wells open to more than one aquifer system. The cones of impression created at these wells could have a significant impact on the transport of contaminants in the Prairie du Chien-Jordan. The presence of coal-tar derivatives in the aquifer has been difficult to explain in wells located upgradient from the plant site to the north, west, and southwest. Simulations suggest that, during periods of heavy withdrawal from certain of these wells (SLP10, SLP15, and SLP5), local hydraulic gradients may have been altered, resulting in the potential for the movement of contaminants from the area of the plant site to the wells. Cones of impression at multiaquifer wells near the plant site contributed to the alteration of local gradients. Simulation of a proposed gradient-control plan, in which lateral homogeneity and isotropy of individual hydro geologic units was assumed, indicates that the actions would be effective in limiting expansion of the contaminated volume in the Prairie du Chien-Jordan aquifer. The plan includes the control of withdrawal from five wells. The simulations also show, however, that modelcalculated potentiometric surfaces are sensitive to changes in withdrawal rates at wells not intended to be under the control of the plan. Management of discharge from these wells also will be important to overall effectiveness of the remedial-action plan.

Minnesota↗

Computerized stratified random site-selection approaches for design of a ground-water-quality sampling network

Computer software was written to randomly select sites for a ground-water-quality sampling network. The software uses digital cartographic techniques and subroutines from a proprietary geographic information system. The report presents the approaches, computer software, and sample applications. It is often desirable to collect ground-water-quality samples from various areas in a study region that have different values of a spatial characteristic, such as land-use or hydrogeologic setting. A stratified network can be used for testing hypotheses about relations between spatial characteristics and water quality, or for calculating statistical descriptions of water-quality data that account for variations that correspond to the spatial characteristic. In the software described, a study region is subdivided into areal subsets that have a common spatial characteristic to stratify the population into several categories from which sampling sites are selected. Different numbers of sites may be selected from each category of areal subsets. A population of potential sampling sites may be defined by either specifying a fixed population of existing sites, or by preparing an equally spaced population of potential sites. In either case, each site is identified with a single category, depending on the value of the spatial characteristic of the areal subset in which the site is located. Sites are selected from one category at a time. One of two approaches may be used to select sites. Sites may be selected randomly, or the areal subsets in the category can be grouped into cells and sites selected randomly from each cell.

Water-Resources Investigations Report↗

A method to estimate canal leakage to the Biscayne Aquifer, Dade County, Florida

The leakage characteristics of channels that partially penetrate the Biscayne aquifer and have reduced bed permeability were studied. Leakage characteristics were described in terms of a reach transmissivity-defined as the volume flow rate out of the channel per unit length of the channel per unit drawdown, where drawdown is defined as the difference in altitude between the water surface in the canal and the water table in the adjacent aquifer. A theoretical expression was developed to relate the reach transmissivity to the transmissivity of the formation, mean channel width, distance of drawdown measurement from the channel centerline, ratio of drawdowns on both sides of the channel, and local reach transmissivity associated with reduced bed permeability. This theoretical expression was verified using a fine-scale numerical model, which gave accurate results when drawdowns were measured beyond 10 aquifer depths from the side of the channel. Using the theoretical formulation, it is shown that the reach transmissivity employed in regional ground-water models, which are based on average drawdowns within a cell, depends on the size of the cell as well as the transmissivity of the formation, channel width, and local reach transmissivity due to reduced bed permeability. The theoretical reach transmissivity function was compared with field measurements at L-31N Canal and Snapper Creek Extension Canal in Dade County, Florida. Analyses of the data for both canals showed good agreement between the estimated and measured reach transmissivities. At L- 31N Canal, field measurements indicated that the local reach transmissivity was relatively uniform over a 2-mile reach of the channel (averaging 630 cubic feet per second per mile per foot), and the formation transmissivity was 1.8 x106 feet squared per day. At Snapper Creek Extension Canal, an approximate analysis was necessary due to the inability of the acoustic velocity meter to measure very low water velocities in the channel. Assuming an aquifer transmissivity of 1 x 106 feet squared per day, drawdown measurements indicated that the local reach transmissivity was about 400 cubic feet per second per mile per foot. The theoretical relation, combined with the local reach transmissivity and formation transmissivity, was sufficient to predict the leakage out of L-31N Canal and Snapper Creek Extension Canal for any drawdown scenario.

Florida↗

Hydrology and water quality of Wind Lake in southeastern Wisconsin

The hydrology and water quality of Wind Lake-a recreational lake in a densely populated area of southeastern Wisconsin was studied from October 1, 1987 through September 30,1989. A drought in 1988 affected the hydrologic budget of Wind Lake in water years 1988-89. Precipitation was 5.9 inches less than normal in water year 1988 but was 2.3 inches greater than normal in water year 1989. Streamflows were near normal in water year 1988 and 25 percent less than normal in water year 1989 as indicated by data from a nearby streamflow-gaging station. Surface runoff was the dominant source of water to the lake in water year 1988 and 75 percent of the inflow was from Big Muskego Lake. The water level in Big Muskego Lake was 1.1 feet below the dam crest at the start of the 1989 water year because of the 1988 drought. About 2,510 acre-feet of water had to fill Big Muskego Lake before water could discharge to Wind Lake. In water year 1989, surface runoff was still the dominant source of water to the lake, but Big Muskego Lake only contributed 52 percent of the water inflow. Streamflow dominated the outflow budget for both years. In water year 1988, 88 percent of the outflow budget left by way of Wind Lake outlet and 12 percent evaporated from the lake surface. In water year 1989, 81 percent of the outflow budget left by way of Wind Lake outlet and 19 percent evaporated from the lake surface. On the basis of outflow from Wind Lake for water year 1988, the hydraulic residence time was 0.46 year; in water year 1989 it was 1.05 years. The total phosphorus input to Wind Lake from external sources was the same for both years, 3,160 pounds. The largest percentage of the phosphorus load came from Big Muskego Lake-- 70 percent in water year 1988 and 65 percent in water year 1989. Analysis of data by use of Vollenweider's model indicates that the phosphorus loadings for each year would cause eutrophic conditions. Data from a nearby gaging station indicate that phosphorus loading to Wind Lake was less than normal. Phosphorus retention in the lake is small and averages 14 percent of the incoming load for both years. Oxygen depletion occurs in the bottom waters during winter and summer months. A maximum anoxic zone was reached on July 18, 1988, when depths greater than 15 feet (about 21 percent of the lake bottom area) were anoxic. Total phosphorus concentrations at the lake surface for both years ranged from 11 to 78 micrograms per liter. Mean total phosphorus concentrations in June, July, and August that had averaged 49 micrograms per liter in 1985 through 1987 declined to 20 micrograms per liter in water year 1988 and 22 micrograms per liter in water year 1989. This reduction was related to the drought and reduced phosphorus loadings. Phosphorus concentrations 1.5 feet above the lake bottom increase during summer anoxic periods. The phosphorus concentration increased at a rate of 5.2 and 4.8 micrograms per liter per day for total and dissolved orthophosphate phosphorus. A maximum concentration of 760 micrograms per liter of total phosphorus and 650 micrograms per liter of dissolved orthophosphate phosphorus occurred on September 21, 1988, just before autumn turnover. Internal loading of phosphorus for the period October 15, 1987 through October 14, 1988, was estimated to be 2,890 pounds. This represents 48 percent of the combined internal and external total-phosphorus input of 5,960 pounds. Algal populations in water year 1988 ranged from 28,200 to 1,610,000 cells per milliliter. A total of 143 species were identified. Blue-green algae dominated the algal population and ranged from 56 percent (February 16, 1988) to 99 percent (five other sampling dates). Aphanocapsa delicatissima caused the largest algal bloom, which reached a maximum concentration of 934,000 cells per milliliter (September 7, 1988). Zooplankton populations in water year 1988 ranged from 52.5 to 686 organisms per liter. Eighteen species were identified. The cladoceran, Daphnia, dominated 12 of the 18 samples.

Wisconsin↗

Identification of potential public water-supply areas of the Cape Cod aquifer, Massachusetts, using a geographic information system

Potential public water-supply areas of the Cape Cod aquifer, Massachusetts, were identified using a geographic information system (GIS) to aid regional and local ground-water resource management efforts. Criteria were selected to identify potential areas on the basis of data restrictions in addition to State requirements for siting new public water- supply wells, Federal or local restrictions on land use, and general hydrogeologic or water-quality concerns. Data layers were created for each criterion and overliad to eliminate areas from consideration as potential public water supplies. Remaining areas, those not included within the applied criteria, are the primary areas to consider for potential public water supplies. The areas identified in this analysis as potential public water supplies range from 0.5 to 7.9 percent of the individual flow cells, or 5.6 percent of the total flow cell area. The criteria were ranked so that criteria more limiting to potential public water supplies were given a higher rank than other criteria. The ranking scheme allows for the inclusion of areas with lower ranked criteria as potential public water supplies. Results can be viewed on a plat in this report, or accessed using the map-based, menu-driven GIS application, which provides interactive display and query of investigation results.

Water-Resources Investigations Report↗

Hydrologic and chemical interaction of the Arkansas River and the Equus Beds aquifer between Hutchinson and Wichita, south-central Kansas

Large chloride concentrations in Arkansas River water have the potential to degrade water quality in the adjacent Equus beds aquifer between Hutchinson and Wichita, Kansas. The aquifer is an important source of water for municipal, industrial, agricultural, and domestic uses. A three-dimensional, finite-difference, ground-water flow-model program (MODFLOW) was used with data from past studies and data collected during 1988-91 to simulate aquifer and stream conditions during the late 1930's, during 1940-89, and during 1990-2019. Results of ground-water flow-model simulations indicated that declining water levels in the Equus beds aquifer since the 1940's have caused base flow in the Arkansas and Little Arkansas Rivers to decrease. In 1940, the Arkansas and Little Arkansas Rivers had simulated net base-flow gains within the model area of about 21 and about 67 ft 3 /s (cubic feet per second), respectively. By the end of 1989, the Arkansas River had a simulated net base-flow loss of about 52 ft 3 /s, and the Little Arkansas River had a net base-flow gain of about 27 ft 3 /s. Simulations for 1990-2019 showed that the water-level changes in a selected model cell located in the central part of the Wichita well field could range from -0.2 to -78 feet. Waterlevel changes in a selected model cell located near the Arkansas River could range from +1.3 to -1.2 feet. In model simulations where only pumpage varied, net base-flow loss from the Arkansas River to the aquifer ranged from about 59 ft 3 /s (no increase in pumpage since 1989) to 117 ft 3 /s (a 3-percent per year increase in pumpage since 1989) by 2019. Assuming a chloride concentration of 630 milligrams per liter, the median concentration in Arkansas River water collected during 1988-91, the quantity of chloride discharged from the Arkansas River to the aquifer was estimated to have increased from about 21 tons per day in 1940 to about 100 tons per day in 1989. By 2019, chloride discharge was indicated to range from about 110 tons per day (associated with no increase in pumpage since 1989) to 200 tons per day (associated with a 3-percent per year increase in pumpage since 1989). A particle-tracking program (MODPATH), which used the results from the flow model, was used to simulate the distribution in the aquifer of chloride from the river during the same time periods. Particle-tracking simulations show that, during 1940-89, the simulated distribution of particles representing chloride from the Arkansas River expanded from relatively narrow bands near the river to a wider distribution within the aquifer and the Wichita well field. Particle-tracking simulations indicate that chloride discharge from the Arkansas River may have reached the edge of the Wichita well field as early as 1963.

Kansas↗

Benthic algae of benchmark streams in agricultural areas of eastern Wisconsin

Benthic algae were collected from 20 streams in the Western Lake Michigan Drainages by the U.S. Geological Survey in May and June of 1993 as part of the National Water-Quality Assessment program. These streams were selected to represent "benchmark" streams that were minimally affected by human activities, especially agriculture, for comparison to other streams in similar environmental settings. Streams were chosen from four relatively homogeneous units (RHU's) in agricultural areas with differing texture of surficial deposits and bedrock type. Blue-green algae were the dominant algal cells at all but 5 of the 20 stream sites, and the most abundant species at these sites was Calothrix parietina, a nitrogen-fixer typically found in pristine streams. Most of the taxa at all sites were diatoms. The dominant diatom guilds observed were the Achnanthes spp., erect forms, and Navicula spp. Except for three streams thought to have low productivity, the Shannon-Wiener diversity index for diatoms was high at all benchmark streams and indicated either minor stress or no stress on the diatom community. With regard to water quality, additional diatom indexes for 17 of 20 benchmark streams indicated no pollution effects and no significant siltation. All benchmark streams had good to excellent biological integrity and either minor or no impairment of aquatic life with regard to diatoms. A variety of algal metrics and relative abundances of diatom morphological guilds correlated with basin-, segment-and reach-level habitat characteristics, including drainage area, basin drainage density, basin soil permeability, Q/Q2 (instantaneous discharge measured at time of sampling divided by the estimated 2-year flood discharge), stream length, and average width of natural riparian vegetation. Algal taxa richness decreased with higher percentages of agricultural land and lower percentages of forested land. The relative abundance of pollution-tolerant diatoms was higher in streams where the basin land was primarily agricultural as compared to forested. The Shannon- Wiener diversity index for diatoms, the percentage of diatom taxa, and the percent relative abundances of diatom cells, pollution tolerant diatoms, Achnanthes spp., erect diatom forms, nitrogen-fixing algae, and blue-green algae differed significantly among either RHU's or ecoregions. Higher abundances of pollution-sensitive diatoms and a higher pollution index indicate that water quality in sampled streams in the North Central Hardwood Forests ecoregion may be less degraded than in streams in the Southeastern Wisconsin Till Plains ecoregion. Algal taxa richness decreased as specific conductance, dissolved nitrate plus nitrite, and suspended sediment increased. This relation may indicate a negative effect of agricultural activities on the algal taxa richness of the stream. Pollution-tolerant diatoms and the pollution index increased as these and additional factors correlated with agriculture increased. Multivariate analyses indicated multiple scales of environmental factors affect algae. Although two-way indicator species analysis (TWINSPAN), detrended correspondence analysis (DCA), and canonical correspondence analysis (CCA) generally separated sites according to RHU, only DCA ordination indicated a separation of sites according to ecoregion. Environmental variables con-elated with DCA axes 1 and 2 and therefore indicated as important explanatory factors for algal distribution and abundance were factors related to stream size, basin land use/cover, geomorphology, hydrogeology, and riparian disturbance. CCA analyses with a more limited set of environmental variables indicated that pH, average width of natural riparian vegetation (segment scale), basin land use/cover and Q/Q2 were the most important variables affecting the distribution and relative abundance of benthic algae at the 20 benchmark streams,

Wisconsin↗

Estimation of the recharge areas contributing water to the South Well Field, Columbus, Ohio

The city of Columbus, Ohio, operates four radial collector wells, designed to yield 42 Mgal/d (million gallons per day), in southern Franklin County, Ohio, as part of their municipal supply of water. The collector wells are adjacent to, and designed to induce infiltration from, Big Walnut Creek and Scioto River. A previously constructed, three-dimensional, steady-state and transient ground-water-flow model of this river-aquifer system was used to estimate contributing recharge areas (CRA's) and calculate particle flowpaths in southern Franklin County. The simulations were of two steady-state periods (October 1979 and March 1986) and one 5-year transient period (March 1986---June 1991). The first simulation (1979) was of conditions before construction of the collector wells. The second simulation (1986) was of conditions when the collector wells were producing 8 Mgal/d. During the 5 years covered in the transient simulation, production at the well field averaged 18.5 Mgal/d. Under the 1979 conditions, the largest ground-water contributing areas were of the quarries and Scioto River (41 and 47 percent of the study area, respectively). During 1986, when 8 Mgal/d was withdrawn, the primary contributing areas were of the quarries (40 percent), collector wells (34 percent), and rivers (8 percent). Travel times associated with simulated particles of water tracked from cells along Big Walnut Creek to their discharge points in cells along Scioto River were about 5 to 60 years in the 1979 simulation and about 7 to 41 years in the 1986 simulation. The endpoints of these particles varied as simulated pumping rates were increased to 22 Mgal/d. The 1986, 10-year CRA's of the collector wells under 8 Mgal/d-conditions totalled about 4.5 mi2. As the pumping rate was increased to 22 Mgal/d in a predictive simulation, 10-year CRA's of the collector wells increased to 6.7mi2. Because the transient simulation encompassed only 5 years, the 10-year CRA's could not be estimated from the transient simulation. However, the size of the 1- to 5-year CRA's for the transient simulation was similar to the size of the 1- to 5-year CRA's for a steady-state predictive simulation if well-field production were 16 Mgal/d. The transient simulations predicted discontinuous CRA's, especially adjacent to the rivers, due to changes in hydrologic stresses. Analyses of the steady-state and transient models showed that sizes of CRA's were most sensitive to changes in porosity, pumping rate, riverbed conductance, and horizontal hydraulic conductivity.

Ohio↗