USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water Resource Report”

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

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

At least 1,693 records · Page 94Linked to original sources

Geohydrology of Storage Unit III and a combined flow model of the Santa Barbara and foothill ground-water basins, Santa Barbara County, California

The city of Santa Barbara pumps most of its ground water from the Santa Barbara and Foothill ground-water basins. The Santa Barbara basin is subdivided into two storage units: Storage Unit I and Storage Unit III. The Foothill basin and Storage Unit I of the Santa Barbara basin have been studied extensively and ground-water flow models have been developed for them. In this report, the geohydrology of the Santa Barbara ground- water basin is described with a special emphasis on Storage Unit III in the southwestern part of the basin. The purposes of this study were to summarize and evaluate the geohydrology of Storage Unit III and to develop an areawide model of the Santa Barbara and Foothill basins that includes the previously unmodeled Storage Unit III. Storage Unit III is in the southwestern part of the city of Santa Barbara. It is approximately 3.5 miles long and varies in width from about 2,000 feet in the southeast to 4,000 feet in the north-west. Storage Unit III is composed of the Santa Barbara Formation and overlying alluvium. The Santa Barbara Formation (the principal aquifer) consists of Pleistocene and Pliocene(?) unconsolidated marine sand, silt, and clay, and it has a maximum saturated thickness of about 160 feet. The alluvium that overlies the Santa Barbara Formation has a maximum saturated thickness of about 140 feet. The storage unit is bounded areally by faults and low-permeability deposits and is underlain by rocks of Tertiary age. The main sources of recharge to Storage Unit III are seepage from Arroyo Burro and infiltration of precipitation. Most of the recharge occurs in the northwest part of the storage unit, and ground water flows toward the southeast along the unit's long axis. Lesser amounts of recharge may occur as subsurface flow from the Hope Ranch subbasin and as upwelling from the underlying Tertiary rocks. Discharge from Storage Unit III occurs as pumpage, flow to underground drains, underflow through alluvium in the vicinity of Arroyo Burro across the Lavigia Fault, evapotranspiration, and underflow to the Pacific Ocean. The faults that bound Storage Unit III generally are considered to be effective barriers to the flow of ground water. Interbasin ground-water flow occurs where deposits of younger alluvium along stream channels cross faults. Ground-water quality in Storage Unit III deposits varies with location and depth. Upward leakage of poor-quality water from the underlying Tertiary rocks occurs in the storage unit, and such leakage can be influenced by poor well construction or by heavy localized pumping. The highest dissolved-solids concentration (4,710 milligrams per liter) in ground water resulting from this upward leakage is found in the coastal part of the storage unit. The ground-water system was modeled as two horizontal layers. In the Foothill basin and Storage Unit I the layers are separated by a confining bed. The upper layer represents the upper producing zone and the shallow zone near the coast. The lower layer represents the lower producing zone. In general, the faults in the study area were assumed to be no-flow boundaries, except for the offshore fault that forms the southeast boundary; the southeast boundary was simulated as a general-head boundary. The Storage Unit III model was combined with the preexisting Storage Unit I and Foothill basin models, using horizontal flow barriers, to form an areawide model. The areawide model was calibrated by simulating steady-state predevelopment conditions and transient conditions for 1978-92. The nonpumping steady- state simulation was used to verify that the calibrated model yielded physically reasonable results for predevelopment conditions. The calibrated areawide model calculates water levels in Storage Unit III that are within 10 feet of measured water levels at all sites of comparison. In addition, the model adequately simulates water levels in the Storage Unit I and Foothill basin areas. A total of 33,430 acre-feet of water was pumped from the study area during the simulation period. Model results indicate that 2,833 acre-feet came from storage and 5, 332 acre-feet crossed the general-head boundary from the ocean, thus indicating that seawater intrusion could occur. A sensitivity analysis indicates that, in general, the model is most sensitive to changes in transmissivity and total recharge.

California↗

Estimating design-flood discharges for streams in Iowa using drainage-basin and channel-geometry characteristics

Drainage-basin and channel-geometry multiple-regression equations are presented for estimating design-flood discharges having recurrence intervals of 2, 5, 10, 25, 50, and 100 years at stream sites on rural, unregulated streams in Iowa. Design-flood discharge estimates determined by Pearson Type-Ill analyses using data collected through the 1990 water year are reported for the 188 streamflow-gaging stations used in either the drainage-basin or channel-geometry regression analyses. Ordinary least-squares multiple-regression techniques were used to identify selected drainage-basin and channel-geometry characteristics and to delineate two channel-geometry regions. Weighted least-squares multiple-regression techniques, which account for differences in the variance of flows at different gaging stations and for variable lengths in station records, were used to estimate the regression parameters. Statewide drainage-basin equations were developed from analyses of 164 streamflow-gaging stations. Drainage-basin characteristics were quantified using a geographic-informationsystem procedure to process topographic maps and digital cartographic data. The significant characteristics identified for the drainage-basin ^equations included contributing drainage area, relative relief, drainage frequency, and 2-year, 24-hour precipitation intensity. The average standard errors of prediction for the drainagebasin equations ranged from 38.6 to 50.2 percent. The geographic-information-system procedure expanded the capability to quantitatively relate drainage-basin characteristics to the magnitude and frequency of floods for stream sites in Iowa and provides a flood-estimation method that is independent of hydrologic regionalization. Statewide and regional channel-geometry regression equations were developed from analyses of 157 streamflow-gaging stations. Channel-geometry characteristics were measured onsite and on topographic maps. Statewide and regional channel-geometry regression equations that are dependent on whether a stream has been channelized were developed on the basis of bankfull and active-channel characteristics. The significant channel-geometry characteristics identified for the statewide and regional regression equations included bankfull width and bankfull depth for natural channels unaffected by channelization, and active-channel width for stabilized channels affected by channelization. The average standard errors of prediction ranged from 41.0 to 68.4 percent for the statewide channel-geometry equations and from 30.3 to 70.0 percent for the regional channel-geometry equations. Procedures provided for applying the drainage-basin and channel-geometry regression equations depend on whether the design-flood discharge estimate is for a site on an ungaged stream, an ungaged site on a gaged stream, or a gaged site. When both a drainage-basin and a channel-geometry regression-equation estimate are available for a stream site, a procedure is presented for determining a weighted average of the two flood estimates. The drainage-basin regression equations are applicable to unregulated rural drainage areas less than 1,060 square miles, and the channel-geometry regression equations are applicable to unregulated rural streams in Iowa with stabilized channels.

Iowa↗

Lithology and fracture characterization from drilling investigations in the Mirror Lake area, Grafton County, New Hampshire

The lithology and fracture network of the bedrock aquifer in the Mirror Lake area were characterized from hydrogeologic data collected from 1979-95 in Grafton County, N.H. The collection of these data is an integral part of an ongoing multidisciplinary study by the U.S. Geological Survey to characterize groundwater flow and solute transport in fractured rock. The data provide a physical framework and enable the characterization of the fractures and the rock types in the bedrock aquifer in the Mirror Lake study area. In addition, these data provide a detailed description of the subsurface intersected by boreholes that can be used to compare the results of other borehole testing. The Mirror Lake area is characterized by steep bedrock uplands that are mostly covered by colluvium, discontinuous stratified-drift deposits, and glacial till that varies locally in thickness from less than 10 meters to as much as 50 meters. The land-surface altitude ranges from 180 meters near the Pemigewasset River to 720 meters on the mountain top on the northwestern side of the study area. The bedrock in the area is predominantly sillimanite-grade pelitic schists that have been complexly folded and intruded by granitoids, pegmatites, and diabase dikes. The bedrock has been fractured in response to local and tectonic stress. The resulting interconnected network of fractures forms the bedrock aquifer. This report describes the lithologic units in the study area and provides a characterization of the lithology and fractures found in 40 boreholes, which range in depth from 60 to 305 meters, that were drilled for this study. Drilling logs and color video surveys were used to locate and characterize the fractures and rock types in the subsurface. Solid bedrock core was obtained from three of the boreholes. Petrographic thin-section, x-ray diffraction and scanning electron microscope with energy dispersive x-ray fluorescence spectrometry analyses were done on selected samples from boreholes and outcrops. Observations recorded at the time of drilling, descriptions of rock samples collected from the boreholes, interpretation of rock type and fractures based on boreholeimaging surveys, descriptions of rock core and petrographic analyses of selected rock samples are in tables and figures. Analysis of the data provided information on the distribution of fractures and lithology in the boreholes at Mirror Lake. The relative abundances of the rock types were computed for three groups of boreholes, including (1) the Forest Service Experimental (FSE) well field, (2) the Camp Osceola (CO) well field, and (3) the index boreholes, which are 15 boreholes distributed areally throughout the study area including the deepest borehole from each of the two well fields. The index boreholes are separated by hundreds of meters and are typically 100 meters deep. The FSE well field includes 13 boreholes that are separated by 10 to 40 meters. These 13 boreholes are approximately 100 meters deep, except for one borehole that is 230 meters deep. The rocks penetrated by the FSE wells are predominantly igneous. Approximately 70 percent of the rocks encountered in the boreholes in the FSE well field were granite, pegmatite, and aplite. The CO well field includes 9 boreholes that range from 60-70 meters deep and one borehole that is 175 meters deep. The rocks encountered in these boreholes were predominantly metamorphic. The distribution of rock types in the CO well field is similar to the distribution of rocks in highway roadcuts, that are approximately 90 to 150 meters east of the well field. Seventy percent of the roadcut exposures are schist. Collectively, in the 15 index boreholes, the metamorphic and igneous rocks are equally distributed. Analysis of the rock types in these boreholes indicates that the rock types tend to "change" every 5 to 9 meters. Although the metamorphic and igneous rocks each comprise approximately 50 percent of the rock types observed in the 15 index boreholes, 73 percent of the fractures were in the igneous rocks. This indicates that the granitoids tend to be more fractured than the metamorphic rocks. Pegmatite, diabase, quartzite, and gneissic rocks are relatively unfractured. Boreholes completed in bedrock generally have one or two water-bearing zones, which were identified during the drilling process. At the completion of drilling a borehole, the driller estimated the yield of the borehole with an air-lift test. Yields estimated by drillers ranged from less than 3 to 378 liters per minute. These yields are typical of the yields measured for domestic wells in Grafton County. Water levels measured in the open boreholes after the boreholes recovered from the hydraulic stresses of drilling were usually in the steel casing and were within 10 meters of the land surface. Water levels in eight of the boreholes were above the top of casing or above land surface.

New Hampshire↗

Evaluation of agricultural best-management practices in the Conestoga River headwaters, Pennsylvania: Characterization of surface-runoff and ground-water quantity and quality in a small carbonate basin near Churchtown, Pennsylvania, prior to terracing and implementation of nutrient management: Water-quality study of the Conestoga River headwaters, Pennsylvania

The U.S. Geological Survey, in cooperation with the Pennsylvania Department of Environmental Protection1 , conducted a study as part of the U.S. Department of Agriculture's Rural Clean Water Program to determine the effects of agricultural best-management practices on surface-water and ground-water quality in the Conestoga River headwaters basin. This report describes Field-Site 1 and characterizes the surface-runoff and ground-water quantity and quality at the site from January 1983 through September 1984, before the implementation of terracing and nutrient-management best-management practices. The 22.1-acre site, part of two dairy farms, was cropland used primarily for the production of corn and alfalfa, and is underlain by carbonate rock. During the 21-month study period, 91.2 inches of precipitation fell, of which 66 percent occurred during the 1984 water year. Of the 169 storms of 0.10 inch or larger, 97 produced measurable runoff. The average annual application of nutrients to the 14.4 acres of cornfields was 410 pounds per acre of nitrogen and 110 pounds per acre of phosphorus. About three times more nutrients were applied during 1984 water year than during the 1983 water year. The approximate soil-nitrate concentration as nitrogen for the study period was 8.2 milligrams per kilogram in the top 8 inches of soil. Runoff for the study period totalled 714,000 cubic feet, or 9.8 percent of total precipitation. Eighty-eight percent of the runoff occurred during the 1984 water year. Regression analyses indicate that total runoff was controlled primarily by total precipitation amounts during storms and by antecedent moisture conditions; mean event and maximum instantaneous discharges are controlled by precipitation intensity and by antecedent moisture conditions. Regression analyses also suggest that crop cover on the corn acreage affected runoff amounts and rates. Mean concentrations in runoff ranged from 24 to 73,000 milligrams per liter for suspended sediment, 0.45 to 24 milligrams per liter for total phosphorus, and 1.5 to 55 milligrams per liter for total nitrogen. Of the total nitrogen in instantaneous runoff samples, a median of 72 percent was total organic nitrogen. Results of multiple regression analyses suggest that mean nitrogen concentrations in runoff (1) increased with increased surface-nitrogen applications made prior to runoff, and (2) were diluted with increased precipitation duration. Runoff from storms on frozen ground produced the highest loads of nitrogen in runoff. Over the 21-month study period, an estimated 258 tons of suspended sediment, 314 pounds of nitrogen, and 176 pounds of phosphorus were transported with runoff from the site. Of this, 88 percent of the nitrogen and phosphorus loads and 94 percent of the suspended-sediment load were discharged during the wet 1984 water year. The loads for the study period represent 2.5 percent of the total nitrogen and 5.5 percent of the total phosphorus applied to the site as manure and commercial fertilizer. The ground-water basin at the site is slightly larger than the surface-water basin. Ground-water levels, except at one well, responded quickly to recharge; water levels peaked several hours to a day after precipitation. Surface-applied materials moved rapidly to the water table through macropores and slowly by transport through micropores in the unsaturated zone. Median concentrations of dissolved nitrate ranged from 9.2 to 13 milligrams per liter as nitrogen for ground-water samples collected monthly from five wells and a spring. Dissolved nitrate comprised 93 percent (median percentage) of the total nitrogen in ground-water samples. The effect of recharge on ground-water nitrate concentrations was affected by variations in nutrient availability at the land surface and in the unsaturated zone. A lag time of 1 to 3 months was observed between the time that nitrogen was applied to the land surface and local maximums in nitrate concentrations were detected in ground water unaffected by recharge events. Approximately 3,187,000 cubic feet of ground water and an associated 2,200 pounds of nitrate-nitrogen discharged from the site during the study period. For the study period, 42 percent of the precipitation recharged to ground water, 10 percent became runoff, and 48 percent evapotranspired. Inputs of nitrogen to the study area were estimated to be 93 percent from manure, 5 percent from commercial fertilizer, and 2 percent from precipitation. Nitrogen outputs from the system were estimated to be 38 percent to crop uptake, 39 percent to volatilization, 20 percent to ground-water discharge, and 3 percent to surface runoff.

Pennsylvania↗

Distribution and variability of fecal-indicator bacteria in Scioto and Olentangy rivers in the Columbus, Ohio, area

This report presents the results of a study by the U.S. Geological Survey, in cooperation with the City of Columbus, Ohio, to determine the distribution and variability of fecal-indicator bacteria in Scioto and Olentangy Rivers. Fecal-indicator bacteria are among the contaminants of concern to recreational users of these rivers in the Columbus area. Samples were collected to be analyzed for fecal-coliform and Escherichia coli (E. coli) bacteria and selected water-quality constituents and physical properties at 10 sites-- 4 on the Olentangy River and 6 on the Scioto River during the recreational seasons in 1987, 1988, and 1989. Measurements of streamflow also were made at these sites at various frequencies during base flow and runoff. The concentrations of fecal-coliform and E. coli bacteria in the Scioto and Olentangy Rivers spanned a range of five orders of magnitude, from less than 20 to greater than 2,000,000 col/100 mL (colonies per 100 milliliters). In addition, the concentrations of fecal coliform and E. coli bacteria are well correlated (r=0.97) in the study area. At times, relatively high concentrations, for fecal-indicator bacteria (concentrations greater than 51,000 col/100 mL for fecal-coliform and E. coli ) were found in Olentangy River at Woody Hayes Drive and at Goodale Street, and in Scioto River at Greenlawn Avenue and at Columbus. Intermediate concentrations of fecal-indicator bacteria (from 5,100 to 50,000 col/100 mL for fecal coliform and (from 510 to 50,000 col/100 mL for E. coli ) were found in Scioto River at Town Street and below O'Shaughnessy Dam near Dublin, Ohio, and in Olentangy River at Henderson Road. The lowest (median) concentrations of fecal-indicator bacteria (from 20 to 5,000 col/100 mL for fecal coliform and from 20 to 500 col/100 mL for E. coli ) were found at Olentangy River near Worthington, Ohio, Scioto River at Dublin Road Water Treatment Plant and below Griggs Reservoir. Fecal-coliform concentrations exceeded the geometric mean and single-sample Ohio Water Quality Standards for recreation less frequently than E. coli concentrations. The E. coli numerical water-quality standards are more difficult to meet than the fecal coliform standards because they are as much as an order of magnitude lower in some instances. The geometric mean bathing-water and primary-contact standards for fecal-coliform and E. coli bacteria were exceeded in more samples for Olentangy River at Goodale Street than for any other site. The single-sample bathing-water standard for fecal-coliform bacteria was exceeded in 83 percent of all samples and for E. coli in 91 percent of samples for Olentangy River at Goodale Street. Compared to Olentangy River at Goodale Street, geometric means and single-samples exceeded the bathing-water standards somewhat less frequently for Scioto River at Town Street and far less frequently for Scioto River at Dublin Road Water Treatment Plant. In contrast to results for fecal-indicator bacteria, the differences between sites for pH and for concentrations for total alkalinity, total chloride, total nonfilterable residue, total nitrate plus nitrite as nitrogen, total phosphorus, and total organic carbon were small. The large contribution of streamflow and discharge of fecal-indicator bacteria from Olentangy River to Scioto River has a major effect on the Scioto River downstream from the confluence of Olentangy River during periods of rainfall and runoff. Fecal-indicator discharges were calculated for times before, during, and at 24-hour intervals for 48 to 72 hours after two runoff-producing storms. Fecal-coliform and E. coli concentrations were lower in samples collected before runoff and during receding streamflows at 24- to 48-hours after the storms than in samples collected during runoff. The fecal-indicator discharges entering Scioto River from Olentangy River ranged from 22.6 to nearly 100 percent of the total for two storms studied. Controlling nonpoint, unregulated,

Water-Resources Investigations Report↗

Water-quality assessment of the Rio Grande Valley, Colorado, New Mexico and Texas: Fish communities at selected sites, 1993-95

Fish communities at 10 sites in the Rio Grande Basin were sampled during low-flow periods between 1993 and 1995 as part of the U.S. Geological Survey National Water-Quality Assessment Program. The ecology of fish communities is one of several lines of evidence used to characterize water-quality conditions. This report describes the fish communities at selected sites in the Rio Grande Basin and relates the structure of these fish communities to the physical and chemical characteristics of the streams. Twenty-nine species of fish representing 10 families were identified in 25 samples collected during this study. Species richness ranged from 1 to 13. Cluster analysis of the 25 samples collected during this study delineated four groups of sites that were based on the similarity of the fish communities. The first two groups were individual sites with low species richness. The third group contained the most samples, and the fourth group consisted of samples from the Rio Grande at Isleta, New Mexico, and the Rio Grande at El Paso, Texas. The shift in community structure of samples from group 3 to group 4 reflects changes from predominantly coldwater fishes to warmwater fishes. Four metrics of biotic integrity (percentages of introduced individuals, omnivores, tolerant individuals, and anomalies) were used in this study to provide a broad overview of the community structure. The relative percentages of introduced species at the Rio Grande near Del Norte, Colorado; Saguache Creek near Saguache, Colorado; Rio Grande below Taos Junction Bridge, near Taos, New Mexico; and Rio Grande at Isleta are indicative of biological stress on the communities at these sites. The dominance of omnivores in samples from the Rio Grande below Taos Junction Bridge, near Taos; Rio Chama near Chamita, New Mexico; Rio Grande at Isleta; and Rio Grande at El Paso is an indication of environmental stress at these sites. In 1995, tolerant species accounted for the entire fish community at the Rio Grande at Isleta. In all samples the occurrence of anomalies was less than 2 percent of the individuals, with the exception of the sample from the Rio Grande at Isleta. On the basis of the relative percentages of introduced individuals, omnivores, tolerant individuals, and anomalies, the biotic integrity at the Rio Grande at Isleta appears to be the most impaired of all sites and shows indications of potential chemical and physical perturbations. Fish communities from three reaches at the Santa Fe River above Cochiti Lake, New Mexico, and the Rio Grande at Isleta were sampled in 1995 to assess small-scale spatial patterns in the structure of fish communities. The spatial pattern at these sites might be associated with natural variability of the fish communities or with the presence of habitat features such as pools. The total number of individuals and relative abundance in a sample varied at sites sampled yearly during this study. All sites, with the exception of the Rio Grande near Del Norte, had a decline in total number of individuals in a sample. The temporal decline in the total number of individuals at these site might be associated with the natural variability within the fish communities.

Colorado, New Mexico, Texas↗

Simulation of the interaction of karstic lakes Magnolia and Brooklyn with the upper Floridan Aquifer, southwestern Clay County, Florida

The stage of Lake Brooklyn, in southwestern Clay County, Florida, has varied over a range of 27 feet since measurements by the U.S. Geological Survey began in July 1957. The large stage changes have been attributed to the relation between highly transient surface-water inflow to the lake and subsurface conduits of karstic origin that permit a high rate of leakage from the lake to the Upper Floridan aquifer. After the most recent and severe stage decline (1990-1994), the U.S. Geological Survey began a study that entailed the use of numerical ground-water flow models to simulate the interaction of the lake with the Upper Floridan aquifer and the large fluctuations of stage that were a part of that process. A package (set of computer programs) designed to represent lake/aquifer interaction in the U.S. Geological Survey Modular Finite-Difference Ground-Water Flow Model (MODFLOW-96) and the Three-Dimensional Method-of-Characteristics Solute-Transport Model (MOC3D) simulators was prepared as part of this study, and a demonstration of its capability was a primary objective of the study. (Although the official names are Brooklyn Lake and Magnolia Lake (Florida Geographic Names), in this report the local names, Lake Brooklyn and Lake Magnolia, are used.) In the simulator of lake/aquifer interaction used in this investigation, the stage of each lake in a simulation is updated in successive time steps by a budget process that takes into account ground-water seepage, precipitation upon and evaporation from the lake surface, stream inflows and outflows, overland runoff inflows, and augmentation or depletion by artificial means. The simulator was given the capability to simulate both the division of a lake into separate pools as lake stage falls and the coalescence of several pools into a single lake as the stage rises. This representational capability was required to simulate Lake Brooklyn, which can divide into as many as 10 separate pools at sufficiently low stage. In the first of two calibrated models, recharge to the water table, specified as a monthly rate, was set equal to 40 percent of the monthly rainfall rate. The specified rate of inflow to the uppermost stream segment was set equal to outflows from Lake Lowry estimated from lake stage and the 1994-97 rating table. Leakage to the intermediate and Upper Floridan aquifers was assumed to occur from the surficial aquifer system through the confining layers directly beneath deeper parts of the lake bottom. A leakance coefficient value of 0.001 feet per day per foot of thickness was used beneath Lake Magnolia, and a value of 0.005 feet per day per foot of thickness was used beneath most of Lake Brooklyn. With these values, the conductance through the confining layers beneath Lake Brooklyn was about 19 times that beneath Lake Magnolia. The simulated stages of Lake Brooklyn matched the measured stages reasonably well in the early (1957-72) and later (1990-98) parts of the simulation time period, but the match was unsatisfactory in an intermediate time period (1973-89). To resolve this discrepancy, the hypothesis was proposed that undocumented losses of water from Alligator Creek upstream from Lake Brooklyn or from the lake itself occurred between 1973 and 1989 when there was sufficient streamflow. The resulting simulation of lake stages matched the measured lake stages accurately during the entire simulation time period. The model was then revised to incorporate the assumption that only 20 percent of precipitation recharged the water table (the second calibrated model). Recalibration of the model required that leakance values for the confining units under deeper parts of the lakes also be reduced by nearly 50 percent. The stages simulated with the new parameter assumptions, but retaining the assumption of surface-water losses, were an excellent match of the measured values. The stage of Lake Magnolia was also simulated accurately. The results of sensitivity analyses show that simulated streamflow between Lakes Magnolia and Brooklyn tends to be water-budget controlled, and is not appreciably affected by the specified outflow altitude or channel characteristics of the receiving stream. To match heads measured in observation wells of the surficial aquifer network, the assigned hydraulic conductivity values were zoned, and ranged from a minimum of 4 feet per day to a maximum of 400 feet per day in the first calibrated model. These values were reduced by about 50 percent in the second calibrated model. Differences between observation wells were noted in the abruptness of changes of measured head values, and in the relation of the timing of peak measured heads and simulated peak heads. These differences seemed to be correlated with the depth of the water table below land surface. Spatially uniform values of transmissivity were specified for the intermediate (10,000 feet squared per day) and Upper Floridan (100,000 feet squared per day) aquifers. Simulated heads in the Upper Floridan aquifer layer follow the trend of the heads measured in a long-term observation well with data beginning in 1960. This result suggests that the observed head decline could be explained entirely in terms of the stage decline in Lake Brooklyn and may not indicate a regional trend.

Florida↗

Evaluating evapotranspiration for six sites in Benton, Spokane, and Yakima counties, Washington, May 1990 to September 1992

This report evaluates evapotranspiration for six sites in Benton, Spokane, and Yakima Counties, Washington. Three sites were located on the Arid Lands Ecology Reserve in Benton County: one at a full-canopy grassland in Snively Basin (Snively Basin site), one at a sparse-canopy grassland adjacent to two weighing lysimeters (grass lysimeter site), and one at a sagebrush grassland adjacent to two weighing lysimeters (sage lysimeter site). Two sites were located on the Turnbull National Wildlife Refuge in Spokane County: one at a full-canopy grassland in a meadow (Turnbull meadow site), the other a full-canopy grassland near a marsh (Turnbull marsh site). The last site was located in a sagebrush grassland in the Black Rock Valley in Yakima County (Black Rock Valley site). The periods of study at the six sites varied, ranging from 5 months at the Black Rock Valley site to more than 2 years at the Snively Basin, grass lysimeter, and sage lysimeter sites. The periods of study were May 1990 to September 1992 for the Snively Basin, grass lysimeter, and sage lysimeter sites; May 1991 to September 1992 for the Turnbull meadow site; May 1991 to April 1992 for the Turnbull marsh site; and March to September 1992 for the Black Rock Valley site. Evapotranspiration and energy-budget fluxes were estimated for the Snively Basin site, the Turnbull meadow site, and the Black Rock Valley site using the Bowen-ratio and Penman-Monteith methods. Daily evapotranspiration for the Snively Basin site was also estimated using a deep-percolation model for the Columbia Basin. The Bowen-ratio method and weighing lysimeters were used at the grass and sage lysimeter sites. The Penman-Monteith method was used at the Turnbull marsh site. Daily evapotranspiration at the sites ranged from under 0.2 millimeter during very dry or cold periods to over 4\x11millimeters after heavy rainfall or during periods of peak transpiration. At all sites, peak evapotranspiration occurred in spring, coinciding with plant growth, and the lowest evapotranspiration occurred in late summer and winter, coinciding with plant dormancy and extremely hot or cold temperatures. Water budgets for the Snively Basin, grass lysimeter, and sage lysimeter sites were based on estimates of precipitation, evapotranspiration, and surface runoff. Surface runoff was estimated at zero for all sites. For the Snively Basin site, 1991 and 1992 water budgets using Penman-Monteith evapotranspiration estimates agreed within 1 percent of the annual budgets computed using deep-percolation model estimates; daily estimates of evapotranspiration by the two methods varied considerably, however. For the Snively Basin site, 100 percent of the precipitation became ET in 1991, and in 1992, about 91 percent of the precipitation became evapotranspiration. Water budgets based on weighing lysimeter data at the grass and sage lysimeter sites agreed within 1 percent of each other for 1991 and within 5 percent of each other for 1992. For 1991, 10 percent of the precipitation became ET at both lysimeter sites. For 1992, 94 to 98 percent of the precipitation became ET at the grass lysimeter site while 98 to 9 percent of the precipitation became ET at the sage lysimeter site. Though there were uncertainties in the methods used, recharge estimates for the Snively Basin and grass and sage lysimeter sites were of the same order of magnitude. The Penman-Monteith method (which incorporated Bowen-ratio measurements), deep- percolation model, and weighing lysimeters indicated that no recharge to subsurface systems (soil profile, unsaturated zone, and ground water) occurred in 1991 and that, in 1992, recharge to subsurface systems was probably less than 10 percent of the annual precipitation at the Snively Basin and grass and sage lysimeter sites.

Washington↗

Processes affecting dissolved-oxygen concentrations in the lower reaches of Middle Fork and South Fork Beargrass Creek, Jefferson County, Kentucky

This report provides data on dissolved-oxygen (DO) concentrations and identifies the environmental processes that most affect DO concentrations during base-flow periods in the lower reaches of Middle Fork and South Fork Beargrass Creek in Jefferson County, Kentucky. These reaches are affected by inputs from combined-sewer overflows. Sections of the lower reaches of the two streams run through single-family residential areas and public parks that are used extensively by local residents during the summer. Recreational fishing and wading also are common in the Middle Fork reach. Continuous-record data collected during the summer and early fall (July-September 1996 on the Middle Fork and July-October 1995 on the South Fork) at three monitoring sites along each reach indicate generally decreasing DO concentrations in the downstream direction except for the South Fork Beargrass Creek at Winter Avenue site where channel modifications have resulted in higher velocities along with shallower depths during low-flow conditions. The channel modifications at this site increased the reaeration-rate coefficient (a measure of the capacity of the stream to absorb oxygen through the air-water interface), increased the potential for algae to attach to the rough concrete surface, and increased algal exposure to sunlight. Synoptic data available for selected constituent concentrations were used to calibrate and verify a computer model (U.S. Environmental Protection Agency QUAL2E model) capable of simulating processes that affect DO concentrations in streams. The results of the study indicate that streamflow, reaeration, and sediment-oxygen demand (SOD) are the factors that most affect net production and depletion of DO in the lower reaches of Middle Fork and South Fork Beargrass Creek. For the QUAL2E model, streamflow is used in the determination of depth, which in tum is used to estimate the consumption of oxygen by SOD. Streamflow also is used in the determination of the reaeration-rate coefficient. From the QUAL2E simulations, DO concentrations (in the mass balance) attributed to reaeration and SOD were at least an order of magnitude greater than any of the other factors that can affect DO concentrations. Large diurnal variability in DO concentrations resulted at the monitoring sites located at upstream and downstream ends of the Middle Fork and South Fork reaches, but as indicated in model simulation, the net effect of photosynthesis and respiration on DO concentration was small. Nitrogen, ammonia, and carbonaceous biochemical-oxygen demand were present at low concentrations in each of the study reaches; the model results indicate these constituents did not have a substantial effect on DO concentrations. Model simulations indicated that lowering the SOD rate by 50 percent would result in a substantial improvement in DO concentrations in the Middle Fork Beargrass Creek reach for extremely low base-flow conditions but would result in only limited improvement in DO concentrations in the South Fork Beargrass Creek reach. However, no simulations for extremely low base-flow conditions were conducted for the South Fork Beargrass Creek reach. More information on SOD is needed for stream reaches affected by periodic inputs of effluent. In such stream systems, the temporal and spatial variability of SOD needs to be better defined.

Kentucky↗

Comparison of Bowen-ratio, eddy-correlation, and weighing-lysimeter evapotranspiration for two sparse-canopy sites in eastern Washington

This report compares evapotranspiration estimated with the Bowen-ratio and eddy-correlation methods with evapotranspiration measured by weighing lysimeters for two sparse-canopy sites in eastern Washington. The sites are located in a grassland area (grass lysimeter site) and a sagbrush- covered area (sage lysimeter site) on the Arid Lands Ecology Reserve in Benton County, Washington. Lysimeter data were collected at the sites from August 1990 to November 1994. Bowen-ratio data were collected for varying periods from May 1993 to November 1994. Additional Bowen-ratio data without interchanging air- temperature and vapor-pressure sensors to remove sensor bias (fixed-sensor system) were collected from October 1993 to June 1994. Eddy-correlation data were collected at the grass lysimeter site from March to April 1994, and at the sage lysimeter site from April to May 1994. The comparisons of evapotranspiration determined by the various methods differed considerably, depending on the periods of record being compared and the sites being analyzed. The year 1993 was very wet, with about 50 percent more precipitation than average; 1994 was a very dry year, with only about half the average precipitation. The study showed that on an annual basis, at least in 1994, Bowen-ratio evapotranspiration closely matched lysimeter evapotranspiration. In 1993, Bowen-ratio and lysimeter evapotranspiration comparisons were variable. Evapotranspiration estimated with the Bowen-ratio method averaged 5 percent more than evapotranspiration measured by lysimeters at the grass lysimeter site from October 1993 to November 1994, and 3 percent less than lysimeters at the sage lysimeter site from November 1993 to October 1994. From March 24 to April 5, 1994, at the grass lysimeter site, the Bowen-ratio method estimated 11 percent less, the Bowen-ratio method utilizing the fixed sensor system about 7 percent more, and the eddy-correlation method about 28 percent less evapotranspiration than the lysimeters measured. From May 7 to June 18, 1993, however, the Bowen-ratio method estimated only 54 percent of the evapotranspiration measured by lysimeters at the grass lysimeter site. This large difference possibly may be attributed to Bowen-ratio instrument variability or error, to the density of grasses in the lysimeters being greater than in the surrounding area, or to heating effects on the lysimeters. From September 1 to October 31, 1993, the Bowen-ratio method estimated more than 450 percent more evapotranspiration than was measured by lysimeters at the sage lysimeter site. This difference may have been due to conditions in the lysimeters at the sage lysimeter site that were unrepresentative of natural conditions. The Bowen-ratio instruments measured evapotrans- piration over sagebrush plants outside the lysimeters, which were blooming very heavily, possibly using supplemental ground water or spring water from nearby upslope areas. The sagebrush plants contained by the lysimeters showed very little evapotranspiration, possibly because they were root-bound and had already used all available water. Also, plants in the lysimeters would not have been able to access any supplemental water available to plants outside the confines of the lysimeters. Earlier in 1993, from June 17 to July 12, the Bowen-ratio method estimated only 1 percent less evapotranspiration than determined for the lysimeters at the sage lysimeter site. On the basis of lysimeter measurements from August 1990 to November 1994, cumulative evapotrans- piration ranged from about 97 to 103 percent of the annual precipitation each year. The evapotranspiration measurements made at the grass and sage lysimeter sites, which were based on weight changes in the lysimeters, showed that storage changes became nearly zero each year some time between August and November as average surface soil moisture decreased to about 2 percent and evapotranspiration rates decreased to less than 0.1 millimeter per day.

Water-Resources Investigations Report↗

Numerical simulation of streamflow distribution, sediment transport, and sediment deposition along Long Beach Creek in Northeast Missouri

This report presents the results of a study conducted by the U.S. Geological Survey in cooperation with the Missouri Department of Conservation to describe the hydrology, sediment transport, and sediment deposition along a selected reach of Long Branch Creek in Macon County, Missouri. The study was designed to investigate spatial and temporal characteristics of sediment deposition in a remnant forested riparian area and compare these factors by magnitude of discharge events both within and outside the measured range of flood magnitudes. The two-dimensional finite-element numerical models RMA2-WES and SED2D-WES were used in conjunction with measured data to simulate streamflow and sediment transport/deposition characteristics during 2-, 5-, 10-, and 25-year recurrence interval floods. Spatial analysis of simulated sediment deposition results indicated that mean deposition in oxbows and secondary channels exceeded that of the remaining floodplain areas during the 2-, 5-, 10-, and 25-year recurrence interval floods. The simulated mass deposition per area for oxbows and secondary channels was 1.1 to 1.4 centimeters per square meter compared with 0.1 to 0.60 centimeters per square meter for the remaining floodplain. The temporal variability of total incremental floodplain deposition during a flood was found to be strongly tied to sediment inflow concentrations. Most floodplain deposition, therefore, occurred at the beginning of the streamflow events and corresponded to peaks in sediment discharge. Simulated total sediment deposition in oxbows and secondary channels increased in the 2-year through 10-year floods and decreased in the 25- year flood while remaining floodplain deposition was highest for the 25-year flood. Despite increases in sediment inflows from the 2-year through 25-year floods, the retention ratio of sediments (the ratio of floodplain deposition to inflow load) was greatest for the 5-year flood and least for the 25-year flood. The decrease in retention ratio at greater flows is likely the result of higher velocities on the floodplain, resulting in higher bed shear stress, greater suspension time of deposited material, and greater sediment transport through the system. Simulated sediment deposition was most sensitive to sediment inflow concentrations and modification of floodplain roughness—factors that can be controlled through management practices. The increase in floodplain sediment deposition resulting from a simulated increase in vegetation density (increase in floodplain roughness from a Manning's n of 0.11 to 0.12) was 142,000 kilograms, or 6.5 percent for a 10-year recurrence interval flood. This increase was comparable to total oxbow and secondary channel deposition mass in the simulations, but would result in a mean increase in floodplain deposition thickness of only 0.025 centimeter. The hydrodynamic model results show the importance of the secondary channels and meander cutoff channels in this system because these areas quickly bring floodwaters and sediment to areas not close to the main channel. The meander cutoff channels in the simulation also effectively decrease flow and velocities in some main channel sections thereby affecting sediment deposition in the vicinity of these features.

Missouri↗

Hydrogeology and Ground-Water Quality of Brunswick County, North Carolina

Brunswick County is the southernmost coastal county in North Carolina and lies in the southeastern part of the Coastal Plain physiographic province. In this report, geologic, hydrologic, and chemical data were used to investigate and delineate the hydrogeologic framework and ground-water quality of Brunswick County. The major aquifers and their associated confining units delineated in the Brunswick County study area include, from youngest to oldest, the surficial, Castle Hayne, Peedee, Black Creek, upper Cape Fear, and lower Cape Fear aquifers. All of these aquifers, with the exception of the Castle Hayne aquifer, are located throughout Brunswick County. The Castle Hayne aquifer extends across only the southeastern part of the county. Based on available data, the Castle Hayne and Peedee confining units are missing in some areas of Brunswick County, which allows direct hydraulic contact between the surficial aquifer and underlying Castle Hayne or Peedee aquifers. The confining units for the Black Creek, upper Cape Fear, and lower Cape Fear aquifers appear to be continuous throughout Brunswick County. In examining the conceptual hydrologic system for Brunswick County, a generalized water budget was developed to better understand the natural processes, including precipitation, evapotranspiration, and stream runoff, that influence ground-water recharge to the shallow aquifer system in the county. In the generalized water budget, an estimated 11 inches per year of the average annual precipitation of 55 inches per year in Brunswick County is estimated to infiltrate and recharge the shallow aquifer system. Of the 11 inches per year that recharges the shallow system, about 1 inch per year is estimated to recharge the deeper aquifer system. The surficial aquifer in Brunswick County is an important source of water for domestic supply and irrigation. The Castle Hayne aquifer is the most productive aquifer and serves as the principal ground-water source of municipal supply for the county. The upper part of the Peedee aquifer is an important source of ground-water supply for domestic and commercial use. Ground water in the lower part of the Peedee aquifer and the underlying aquifers is brackish and is not known to be used as a source of supply in Brunswick County. Most of the precipitation that recharges the surficial aquifer is discharged to local streams that drain into the Waccamaw River, Cape Fear River, and Atlantic Ocean. Recharge to the Castle Hayne aquifer occurs primarily from the surficial aquifer. Recharge to the Peedee aquifer occurs primarily from the surficial and Castle Hayne aquifers, with some upward leakage of water also occurring from the underlying Black Creek aquifer. Discharge from the Castle Hayne and Peedee aquifers occurs to local streams, the Cape Fear River, and the Atlantic Ocean. Evaluation of water-level data for the period January 1970 through May 2002 indicated no apparent long-term temporal trends in water levels in the surficial and Castle Hayne aquifers and in the upper part of the Peedee aquifer. The most significant water-level trends were noted for wells tapping the lower part of the Peedee aquifer and tapping the Black Creek aquifer where water levels have declined as much as 41 and 37 feet, respectively. These ground-water-level declines are attributed to regional ground-water pumping in areas outside of Brunswick County. Water-level data for Brunswick County wells tapping the upper Cape Fear and lower Cape Fear aquifers tend to fluctuate within a fairly uniform range with no apparent temporal trend noted. Analysis of vertical hydraulic gradients during this same period primarily indicate downward flow of ground water within and among the surficial, Castle Hayne, and Peedee aquifers. The vertical flow of ground water in the Black Creek aquifer is upward into the overlying Peedee aquifer. Upward flow also is noted for the upper and lower Cape Fear aquifers. Historic and recent analytic data were evaluated to better understand the sources of water contained in Brunswick County aquifers and the suitability of the water for consumption. Based on analytical results obtained for recent samples collected during this study, ground water from the surficial aquifer, Castle Hayne aquifer, and upper part of the Peedee aquifer appears to be generally suitable for drinking water. Although concentrations of iron and manganese commonly exceeded the drinking-water standards, the concern generally associated with the occurrence of these analytes in a water supply is one of aesthetics. In all samples, nitrate, nitrite, and sulfate were detected at concentrations less than drinkingwater standards. Based on historic analytical data, the brackish water in the lower part of the Peedee aquifer and in the Black Creek, upper Cape Fear, and lower Cape Fear aquifers is classified as a sodium-chloride type water. The presence of brackish water in these deeper systems combined with upward vertical gradients presents the potential for upward migration of brackish water into overlying aquifers, or upconing beneath areas of pumping. The current (2001) location of the boundary between freshwater and brackish water in Brunswick County aquifers is unknown.

North Carolina↗

Arsenic loads in Spearfish Creek, western South Dakota, water years 1989-91

Numerous small tributaries on the eastern flank of Spearfish Creek originate within a mineralized area with a long history of gold-mining activity. Some streams draining this area are known to have elevated concentrations of arsenic. One such tributary is Annie Creek, where arsenic concentrations regularly approach the Maximum Contaminant Level of 50 mg/L (micrograms per liter) established by the U.S. Environmental Protection Agency. A site on Annie Creek was proposed for inclusion on the National Priorities List by the Environmental Protection Agency in 1991. This report presents information about arsenic loads and concentrations in Spearfish Creek and its tributaries, including Annie Creek. Stream types were classified according to geologic characteristics and in-stream arsenic concentrations. The first type includes streams that lack significant arsenic sources and have low in-stream arsenic concentrations. The second type has abundant arsenic sources and high in-stream concentrations. The third type has abundant arsenic sources but only moderate in-stream concentrations. The fourth type is a mixture of the first three types. Annual loads of dissolved arsenic were calculated for two reaches of Spearfish Creek to quantify arsenic loads at selected gaging stations during water years 1989-91. Mass-balance calculations also were performed to estimate arsenic concentrations for ungaged inflows to Spearfish Creek. The drainage area of the upstream reach includes significant mineralized areas, whereas the drainage area of the downstream reach generally is without known arsenic sources. The average load of dissolved arsenic transported from the upstream reach of Spearfish Creek, which is representative of a type 4 stream, was 158 kilograms per year, calculated for station 06430900, Spearfish Creek above Spearfish. Gaged headwater tributaries draining unmineralized areas (type 1) contributed only 16 percent of the arsenic load in 63 percent of the discharge. Annie Creek (type 2), which has the highest measured arsenic concentra- tions in the Spearfish Creek drainage, contributed about 15 percent of the arsenic load in about 2 percent of the discharge of the upstream reach. Squaw Creek, which drains another mineralized area, but has only moderate in-stream concentrations (type 3), contributed 4 percent of the arsenic load in 5 percent of the discharge. Ungaged inflows to the reach contributed the remaining 65 percent of the arsenic load in 30 percent of the discharge. The calculated loads from ungaged inflows include all arsenic contributed by surface- and ground-water sources, as well as any additions of arsenic from dissolution of arsenic-bearing solid phases, or from desorption of arsenic from solid surfaces, within the streambed of the upstream reach. Mass-balance calculations indicate that dissolved arsenic concentrations of the ungaged inflows in the upstream reach averaged about 9 mg/L. In-stream arsenic concentrations of ungaged inflows from the unmineralized western flank of Spearfish Creek probably are generally low (type 1). Thus, in-stream arsenic concentrations for ungaged inflows draining the mineralized eastern flank of Spearfish probably average almost twice that level, or about 18 mg/L. Some ungaged, eastern-flank inflows probably are derived from type 3 drainages, with only moderate arsenic concentrations. If so, other ungaged, eastern-flank inflows could have in-stream arsenic concentrations similar to those of Annie Creek. No significant arsenic sources were apparent in the downstream reach of Spearfish Creek. Over the course of the downstream reach, arsenic concentrations decreased somewhat, probably resulting from dilution, as well as from possible chemical adsorption to sediment surfaces or arsenic-phase precipitation. A decrease in arsenic loads resulted from various diversions from the creek and from the potential chemical removal processes. Because of a large margin of error associated with calculation of arsenic loads for the downstream reach, it is difficult to conclude decisively that chemical removal occurs. It can be concluded, however, that arsenic concentrations do not increase significantly in the downstream reach as a result of interactions with streambed and alluvial sediments. Thus, it also can be concluded that streambed interactions within the channel of Spearfish Creek probably are not a significant source of arsenic within the upstream reach

South Dakota↗

Twentieth century arroyo changes in Chaco Culture National Historical Park

Chaco Wash arroyo channel changes in the 20th century have become a major concern of the National Park Service. Several archeologic and cultural sites are located in the Chaco Wash corridor; thus, increased erosional activity of Chaco Wash, such as channel incision and increased meandering, may affect these sites. Through field surveys, photogrammetric analyses, and reviews of existing reports and maps, arroyo changes at Chaco Culture National Historic Park were documented. Arroyo changes were documented for the inner active channel and the entire arroyo cross section. The inner channel of Chaco Wash evolved from a wide, braided channel in the 1930's to a narrower channel with a well-developed flood plain by the 1970's. From 1934 to 1973 the active channel narrowed an average of 26 meters, and from the 1970's to 2000 the channel narrowed an average of 9 meters. Overall from 1934 to 2000, the inner channel narrowed an average of 30 meters. From 1934 to 2000, the top of Chaco Wash widened at four cross sections, narrowed at one, and remained the same at another. The top of Chaco Wash widened at a rate of 0.4 meter per year from the 1970's to 2000 compared with 0.2 meter per year from 1934 to 1973. At 50-percent depth or halfway down the arroyo channel, four cross sections widened and two cross sections narrowed from 1934 to 2000. Rates of widening at 50-percent depth decreased from 0.2 meter per year from 1934 to 1973 to 0.1 meter per year from the 1970's to 2000. From 1934 to 2000, arroyo depth decreased at five of six cross sections and increased at one cross section. Arroyo depth between 1934 and 1973 decreased an average 1.4 meters from aggradation and between the 1970's and 2000 increased an average 0.4 meter from channel scour. From 1934 to 2000, arroyo cross-sectional area decreased at all six cross sections. Cross-sectional areas in Chaco Wash decreased from 1934 to 1973 as a result of sediment deposition and both decreased and increased from the 1970's to 2000. The cross-sectional area decreased by the 1970's due to channel narrowing and flood-plain formation. Increases in cross-sectional area are from channel scour and channel widening. Photogrammetric analyses of volumetric changes for a 1.7-kilometer reach of Chaco Wash showed sediment deposition from 1934 to 1973 of 64 square meters per unit length of channel over 1.7 kilometers to erosion from 1973 to 2000 of 7 square meters per unit length of channel. Chaco Wash evolved from a braided channel in the 1930's to a narrow, sinuous inner channel by the 1970's. Chaco Wash was widening in the 1930's, leading to sediment deposition and formation of an inner flood plain. Channel narrowing resulted from increased sediment deposition on the flood plain. Sediment deposition may be related to a decrease in peak flows, an increase in flood-plain vegetation, or an increase in the transport of fine-grained sediment. Increases in bankfull depth of Chaco Wash between the 1970's and 2000 were due to aggradation of the flood plain and channel scour. Thus, rates of aggradation and cross-sectional filling were greater from 1934 to the 1970's than from the 1970's to 2000.

Water-Resources Investigations Report↗

A water-resources data-network evaluation for Monterey County, California, phase 3: Northern Salinas River drainage basin

This report evaluates existing data collection networks and possible additional data collection to monitor quantity and quality of precipitation, surface water, and groundwater in the northern Salinas River drainage basin, California. Of the 34 precipitation stations identified, 20 were active and are concentrated in the northwestern part of the study area. No precipitation quality networks were identified, but possible data collection efforts include monitoring for acid rain and pesticides. Six of ten stream-gaging stations are active. Two surface water quality sites are sampled for suspended sediment, specific conductance, and chloride; one U.S. Geological Survey NASOAN site and one site operated by California Department of Water Resources make up the four active sampling locations; reactivation of 45 inactive surface water quality sites might help to achieve objectives described in the report. Three local networks measure water levels in 318 wells monthly, during peak irrigation, and at the end of the irrigation season. Water quality conditions are monitored in 379 wells; samples are collected in summer to monitor saltwater intrusion near Castroville and are also collected annually throughout the study area for analysis of chloride, specific conductance, and nitrate. An ideal baseline network would be an evenly spaced grid of index wells with a density of one per section. When baseline conditions are established, representative wells within the network could be monitored periodically according to specific data needs. (USGS)

California↗

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↗

Methods for Estimating Peak Discharges and Unit Hydrographs for Streams in the City of Charlotte and Mecklenburg County, North Carolina

Procedures for estimating peak discharges and unit hydrographs were developed for streams in the city of Charlotte and Mecklenburg County in response to a need for better techniques for characterizing the flow of streams. The procedures presented in this report provide the means for estimating unit hydrographs as part of the process used in watershed modeling and(or) design of stormwater-management structures. The procedures include three statistical relations for use in estimating storm peak discharge, unit-hydrograph peak discharge, and unit-hydrograph lag time. A final component of the procedures is the development of a dimensionless unit hydrograph developed from streamflow and rainfall data collected during the 1995-2000 water years at 25 streamgaging stations and up to 60 raingages in the city and county. The statistical relation to estimate the storm peak discharge is based on analyses of observed peak discharges regressed against rainfall and basin characteristics using a database of 412 observations from 61 storm events among the 25 gaging stations. The rainfall characteristics included basin-average rainfall amounts as well as estimates of the maximum and minimum storm rainfall in the basin. The basin characteristics consisted of land-use information and other physical basin characteristics, such as drainage area, channel length, channel slope, percentage of impervious area, and percentage of the basin served by detention. The analyses resulted in a relation that can be used for estimating storm peak discharge based on drainage area, basin-average rainfall, and impervious area. Average unit hydrographs were developed for 24 of 25 streamgaging stations, using from three to nine storms at each site. The average unit hydrograph for each station was converted into four classes of unit hydrographs with durations corresponding to one-fourth, one-third, one-half, and three-fourths of the station-average lag time. For 23 sites, the lag-time-duration hydrographs were then translated into dimensionless unit hydrographs by dividing time ordinates by the lag time and discharge ordinates by peak discharge. For each lag-time-duration class, the dimensionless unit hydrographs for the sites were combined to create an average dimensionless unit hydrograph. The four average dimensionless unit hydrographs were later tested (with estimates of unit-hydrograph peak discharges and lag times) for selection of an overall dimensionless unit hydrograph to be used at ungaged sites in the study area. The two sites where the procedures did not produce unit hydrographs that could be included in the development of the overall dimensionless unit hydrograph had the smallest drainage areas among the sites used in the investigation. The statistical relations for estimating unit-hydrograph peak discharge and lag time were developed by regressing the dependent variables against explanatory variables that describe the basin characteristics. The statistical analyses resulted in a relation for use in estimating a unit-hydrograph peak discharge based on the drainage area. The estimation of the unit-hydrograph lag time is based on the drainage area and percentage of land use in the basin classified as 'woods/brush.' Both relations have coefficients of determination (R2 values) of 0.9 or better. The three components for estimating a unit hydrograph are the dimensionless unit hydrograph and two statistical relations for estimating the unit-hydrograph peak discharge and lag time. These components were applied by using each of the four lag-time-duration average dimensionless unit hydrographs to determine which would be selected as the final overall dimensionless unit hydrograph for streams in the city of Charlotte and Mecklenburg County. Comparisons of the simulated and observed hydrographs were based on the following: (1) hydrograph width at 50 percent of the peak discharge, (2) hydrograph width at 75 percent of the peak discharge, (3) peak discharge

North Carolina↗

Nutrients and organic compounds in Deer Creek and south branch Plum Creek in southwestern Pennsylvania, April 1996 through September 1998

This report presents results of an analysis of nutrient and pesticide data from two surface-water sites and volatile organic compound (VOC) data from one of the sites that are within the Allegheny and Monongahela River Basins study unit of the National Water-Quality Assessment Program of the U.S. Geological Survey. The Deer Creek site was located in a 27.0 square-mile basin within the Allegheny River Basin in Allegheny County. The primary land uses consist of small urban areas, large areas of residential housing, and some agricultural land in the upper part of the basin. The South Branch Plum Creek site was located in a 33.3 square-mile basin within the Allegheny River Basin in Indiana County. The primary land uses throughout this basin are mostly agriculture and forestland. Water samples for analysis of nutrients were collected monthly and during high-flow events from April 1996 through September 1998. Concentrations of dissolved nitrite, dissolved ammonia plus organic nitrogen, and dissolved phosphorus were less than the method detection limits in more than one-half of the samples collected. The median concentration of dissolved nitrite plus nitrate in South Branch Plum Creek was 0.937 mg/L and 0.597 mg/L in Deer Creek. The median concentration of dissolved orthophosphate was 0.01 mg/L in both streams. High loads of nitrate were measured in both streams from March to June. Concentrations of dissolved ammonia nitrogen, dissolved nitrate, and total phosphorus were lower during the summer months. Measured concentrations of nitrate nitrogen in both streams were well below the U.S. Environmental Protection Agency (USEPA) maximum contaminant level (MCL) of 10 mg/L. Water samples for analysis of pesticides were collected throughout 1997 in both streams and during a storm event on August 25-26, 1998, in Deer Creek. Samples were collected monthly at both sites and more frequently during the spring and early summer months to coincide with application of pesticides. Seventy-eight pesticides and 7 pesticide metabolites were analyzed in 31 samples collected in Deer Creek and in 18 samples collected in South Branch Plum Creek. Of the 85 pesticides and pesticide metabolites analyzed, 25 of the pesticides were detected at least once in Deer Creek, and 20 of the pesticides were detected at least once in South Branch Plum Creek. Atrazine was the most commonly detected pesticide in both streams. There was a distinct seasonal pattern of atrazine, simazine, and metolachlor concentrations measured at both sites. Prometon was detected in 3 of the 18 samples collected in South Branch Plum Creek in 1997 and in 28 of the 31 samples collected in Deer Creek in both 1997 and 1998. Prometon generally is applied in conjunction with asphalt paving projects and is commonly used in residential areas. The highest measured concentrations of prometon detected in Deer Creek were in the five storm samples collected on August 25-26, 1998. At the Deer Creek site, 9 of the 25 pesticides detected throughout the study were detected only in the sample collected on June 13, 1997. Those nine pesticides included acifluorfen, bentazon, bromoxynil, dicamba, dichlorprop, fenuron, linuron, MCPA, and neburon. Nine other pesticides also were detected in that sample. All concentrations of pesticides were well below established drinking-water guidelines. The maximum measured concentration of diazinon in Deer Creek (0.097 µg/L) and South Branch Plum Creek (0.974 µg/L) exceeded the aquatic life guideline of 0.009 µg/L established by the National Academy of Sciences/National Academy of Engineers. The maximum measured concentration of azinphos-methyl in South Branch Plum Creek (an estimated value of 0.033 µg/L) exceeded the chronic aquatic-life guideline of 0.01 µg/L established by the USEPA. Twenty-five samples were collected from Deer Creek and analyzed for volatile organic compounds (VOCs). Of 87 VOCs analyzed for, 22 were detected at least once, and 12 were gasoline-related compounds. Acetone, benzene, carbon disulfide, meta/paraxylene, methyl chloride, MTBE, p-isopropyl toluene, toluene, and 1,2,4-trimethylbenzene were each detected in five or more samples. VOCs generally were detected during the colder winter months and not frequently during the summer months. The maximum measured concentrations of benzene, ethylbenzene, o-dichlorobenzene, styrene, and toluene were two or more orders of magnitude lower than the MCLs established by the USEPA.

Pennsylvania↗