Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Input”

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,009 records · Page 56Linked to original sources

Hydrogeology, Ground-Water-Age Dating, Water Quality, and Vulnerability of Ground Water to Contamination in a Part of the Whitewater Valley Aquifer System near Richmond, Indiana, 2002-2003

Assessments of the vulnerability to contamination of ground-water sources used by public-water systems, as mandated by the Federal Safe Drinking Water Act Amendments of 1996, commonly have involved qualitative evaluations based on existing information on the geologic and hydrologic setting. The U.S. Geological Survey National Water-Quality Assessment Program has identified ground-water-age dating; detailed water-quality analyses of nitrate, pesticides, trace elements, and wastewater-related organic compounds; and assessed natural processes that affect those constituents as potential, unique improvements to existing methods of qualitative vulnerability assessment. To evaluate the improvement from use of these methods, in 2002 and 2003, the U.S. Geological Survey, in cooperation with the City of Richmond, Indiana, compiled and interpreted hydrogeologic data and chemical analyses of water samples from seven wells in a part of the Whitewater Valley aquifer system in a former glacial valley near Richmond. This study investigated the application of ground-water-age dating, dissolved-gas analyses, and detailed water-quality analyses to quantitatively evaluate the vulnerability of ground water to contamination and to identify processes that affect the vulnerability to specific contaminants in an area of post-1972 greenfield development. The aquifer system in the study area includes an unconfined sand and gravel aquifer used for public-water supply (upper aquifer) and a confined sand and gravel aquifer (lower aquifer) separated by a till confining unit. Several hydrogeologic and cultural measures indicate that the upper aquifer is qualitatively vulnerable to contamination: the upper aquifer is unconfined and has a shallow depth to the water table (from about 4.75 to 14 feet below land surface), low-permeability sediments in the unsaturated zone are thin (less than 10 feet thick), estimated ground-water-flow rates through the upper aquifer are relatively rapid (the highest estimated rates ranged from 0.44 to about 5.0 feet per day), and potential contaminant sources were present. Ground-water-age dates indicate that ground-water samples represented recharge from about the time greenfield development began south of the ground-water-flow divide and that changes in water quality would lag changes in contaminant inputs. Estimates of ground-water age, computed with dichlorodifluoromethane (CFC-12) and trichlorotrifluoroethane (CFC-113) concentrations in water samples collected from seven observation wells in February and March 2003, indicated that water in the upper aquifer had recharged within about 13 to 30 years before sampling. Ground-water ages were youngest (from about 13 to 15 years since recharge) in water from the shallow wells along the glacial-valley margin and oldest (30 years) in water from a well at the base of the aquifer in the valley center. Ground-water ages determined for the shallow wells may be affected by mixing of recent recharge with older ground water from deeper in the aquifer, as indicated by upward hydraulic gradients between paired shallow and deep wells in the upper aquifer. Other parts of the Whitewater Valley aquifer system with similar hydrogeologic characteristics could be expected to have similarly young ground-water ages and residence times. Analyses of water samples collected from the seven observation wells in August and September 2002 indicated that concentrations of chloride, sodium, and nitrate generally were larger in ground water from the upper aquifer than in other parts of the Whitewater Valley aquifer system. Drinking-water-quality standards for Indiana were exceeded in water samples from one well for chloride concentrations, from four wells for dissolved-solids concentrations, and from one well for nitrate concentrations. Application of low-level methods for trace-element analyses determined that concentrations of aluminum, cobalt, iron, lithium, molybdenum, nickel, selenium, uranium, vanadium, and zinc were less than or equal to 8 micrograms per liter; concentrations of arsenic, cadmium, chromium, and copper were less than or equal to 1 microgram per liter. Application of low-level analytical methods to water samples enabled the detection of several pesticides and volatile, semivolatile, and wastewater-related organic compounds; concentrations of individual pesticides and volatile organic compounds were less than 0.1 microgram per liter and concentrations of individual wastewater organic compounds were less than 0.5 microgram per liter. The low-level analytical methods will provide useful data with which to compare future changes in water quality. Results of detailed water-quality analyses, ground-waterage dating, and dissolved-gas analyses indicated the vulnerability of ground water to specific types of contamination, the sequence of contaminant introduction to the aquifer relative to greenfield development, and processes that may mitigate the contamination. Concentrations of chloride and sodium and chloride/bromide weight ratios in sampled water from five wells indicated the vulnerability of the upper aquifer to roaddeicer contamination. Ground-water-age estimates from these wells indicated the onset of upgradient road-deicer use within the previous 25 years. Nitrate in the upper aquifer predates the post-1972 development, based on a ground-water-age date (30 years) and the nitrate concentration (5.12 milligrams per liter as nitrogen) in water from a deep well. Vulnerability of the aquifer to nitrate contamination is limited partially by denitrification. Detection of one to four atrazine transformation products in water samples from the upper aquifer indicated biological and hydrochemical processes that may limit the vulnerability of the ground water to atrazine contamination. Microbial processes also may limit the aquifer vulnerability to small inputs of halogenated aliphatic compounds, as indicated by microbial transformations of trichlorofluoromethane and trichlorotrifluoroethane relative to dichlorodifluoromethane. The vulnerability of ground water to contamination in other parts of the aquifer system also may be mitigated by hydrodynamic dispersion and biologically mediated transformations of nitrate, pesticides, and some organic compounds. Identification of the sequence of contamination and processes affecting the vulnerability of ground water to contamination would have been unlikely with conventional assessment methods.

Indiana, Ohio↗

Selenium and other elements in water and adjacent rock and sediment of Toll Gate Creek, Aurora, Arapahoe County, Colorado, December 2003 through March 2004

Streamwater and solid samples (rock, unconsolidated sediment, stream sediment, and efflorescent material) in the Toll Gate Creek watershed, Colorado, were collected and analyzed for major and trace elements to determine trace-element concentrations and stream loads from December 2003 through March 2004, a period of seasonally low flow. Special emphasis was given to selenium (Se) concentrations because historic Se concentrations exceeded current (2004) stream standards. The goal of the project was to assess the distribution of Se concentration and loads in Toll Gate Creek and to determine the potential for rock and unconsolidated sediment in the basin to be sources of Se to the streamwater. Streamwater samples and discharge measurements were collected during December 2003 and March 2004 along Toll Gate Creek and its two primary tributaries - West Toll Gate Creek and East Toll Gate Creek. During both sampling periods, discharge ranged from 2.5 liters per second to 138 liters per second in the watershed. Discharge was greater in March 2004 than December 2003, but both periods represent low flow in Toll Gate Creek, and results of this study should not be extended to periods of higher flow. Discharge decreased moving downstream in East Toll Gate Creek but increased moving downstream along West Toll Gate Creek and the main stem of Toll Gate Creek, indicating that these two streams gain flow from ground water. Se concentrations in streamwater samples ranged from 7 to 70 micrograms per liter, were elevated in the upstream-most samples, and were greater than the State stream standard of 4.6 micrograms per liter. Se loads ranged from 6 grams per day to 250 grams per day, decreased in a downstream direction along East Toll Gate Creek, and increased in a downstream direction along West Toll Gate Creek and Toll Gate Creek. The largest Se-load increases occurred between two sampling locations on West Toll Gate Creek during both sampling periods and between the two sampling locations on the main stem of Toll Gate Creek during the December 2003 sampling. These load increases may indicate that sources of Se exist between these two locations; however, Se loading along West Toll Gate Creek and Toll Gate Creek primarily was characterized by gradual downstream increases in load. Linear regressions between Se load and discharge for both sampling periods had large, significant values of r2 (r2 > 0.96, p < 0.0001) because increases in Se load (per unit of flow increase) were generally constant. This relation is evidence for a constant addition of water having a relatively constant Se concentration over much of the length of Toll Gate Creek, a result which is consistent with a ground-water source for the Se loads. Rock outcroppings along the stream were highly weathered, and Se concentrations in rock and other solid samples ranged from below detection (1 part per million) to 25 parts per million. One sample of efflorescence (a surface encrustation produced by evaporation) had the greatest selenium concentration of all solid samples, was composed of thenardite (sodium sulfate), gypsum (calcium sulfate) and minor halite (sodium chloride), and released all of its Se during a 30-minute water-leaching procedure. Calculations indicate there was an insufficient amount of this material present throughout the watershed to account for the observed Se load in the stream. However, this material likely indicates zones of ground-water discharge that contain Se. This report did not identify an unequivocal source of Se in Toll Gate Creek. However, multiple lines of evidence indicate that ground-water discharge supplies Se to Toll Gate Creek: (1) the occurrence of elevated Se concentrations in the stream throughout the watershed and in the headwater regions, upstream from industrial sources; (2) the progressive increase in Se loads moving downstream, which indicates a continuous input of Se along the stream rather than input from point sources; (3) the occurr

Colorado↗

Description and User Manual for a Web-Based Interface to a Transit-Loss Accounting Program for Monument and Fountain Creeks, El Paso and Pueblo Counties, Colorado

The U.S. Geological Survey, in cooperation with Colorado Springs Utilities, the Colorado Water Conservation Board, and the El Paso County Water Authority, began a study in 2004 with the following objectives: (1) Apply a stream-aquifer model to Monument Creek, (2) use the results of the modeling to develop a transit-loss accounting program for Monument Creek, (3) revise an existing accounting program for Fountain Creek to easily incorporate ongoing and future changes in management of return flows of reusable water, and (4) integrate the two accounting programs into a single program and develop a Web-based interface to the integrated program that incorporates simple and reliable data entry that is automated to the fullest extent possible. This report describes the results of completing objectives (2), (3), and (4) of that study. The accounting program for Monument Creek was developed first by (1) using the existing accounting program for Fountain Creek as a prototype, (2) incorporating the transit-loss results from a stream-aquifer modeling analysis of Monument Creek, and (3) developing new output reports. The capabilities of the existing accounting program for Fountain Creek then were incorporated into the program for Monument Creek and the output reports were expanded to include Fountain Creek. A Web-based interface to the new transit-loss accounting program then was developed that provided automated data entry. An integrated system of 34 nodes and 33 subreaches was integrated by combining the independent node and subreach systems used in the previously completed stream-aquifer modeling studies for the Monument and Fountain Creek reaches. Important operational criteria that were implemented in the new transit-loss accounting program for Monument and Fountain Creeks included the following: (1) Retain all the reusable water-management capabilities incorporated into the existing accounting program for Fountain Creek; (2) enable daily accounting and transit-loss computations for a variable number of reusable return flows discharged into Monument Creek at selected locations; (3) enable diversion of all or a part of a reusable return flow at any selected node for purposes of storage in off-stream reservoirs or other similar types of reusable water management; (4) and provide flexibility in the accounting program to change the number of return-flow entities, the locations at which the return flows discharge into Monument or Fountain Creeks, or the locations to which the return flows are delivered. The primary component of the Web-based interface is a data-entry form that displays data stored in the accounting program input file; the data-entry form allows for entry and modification of new data, which then is rewritten to the input file. When the data-entry form is displayed, up-to-date discharge data for each station are automatically computed and entered on the data-entry form. Data for native return flows, reusable return flows, reusable return flow diversions, and native diversions also are entered automatically or manually, if needed. In computing the estimated quantities of reusable return flow and the associated transit losses, the accounting program uses two sets of computations. The first set of computations is made between any two adjacent streamflow-gaging stations (termed 'stream-segment loop'); the primary purpose of the stream-segment loop is to estimate the loss or gain in native discharge between the two adjacent streamflow-gaging stations. The second set of computations is made between any two adjacent nodes (termed 'subreach loop'); the actual transit-loss computations are made in the subreach loop, using the result from the stream-segment loop. The stream-segment loop is completed for a stream segment, and then the subreach loop is completed for each subreach within the segment. When the subreach loop is completed for all subreaches within a stream segment, the stream-segment loop is initiated for the ne

Scientific Investigations Report↗

Water-balance model of a wetland on the Fort Berthold Reservation, North Dakota

A numerical water-balance model was developed to simulate the responses of a wetland on the Fort Berthold Reservation, North Dakota, to historical and possible extreme hydrological inputs and to changes in hydrological inputs that might occur if a proposed refinery is built on the reservation. Results from model simulations indicated that the study wetland would likely contain water during most historical and extreme-precipitation events with the addition of maximum potential discharges of 0.6 acre-foot per day from proposed refinery holding ponds. Extended periods with little precipitation and above-normal temperatures may result in the wetland becoming nearly dry, especially if potential holding-pond discharges are near zero. Daily simulations based on the historical-enhanced climate data set for May and June 2005, which included holding-pond discharges of 0.6 acre-foot per day, indicated that the study-wetland maximum simulated water volume was about 16.2 acre-feet and the maximum simulated water level was about 1.2 feet at the outlet culvert. Daily simulations based on the extreme summer data set, created to represent an extreme event with excessive June precipitation and holding-pond discharges of 0.6 acre-foot per day, indicated that the study-wetland maximum simulated water volume was about 38.6 acre-feet and the maximum simulated water level was about 2.6 feet at the outlet culvert. A simulation performed using the extreme winter climate data set and an outlet culvert blocked with snow and ice resulted in the greatest simulated wetland water volume of about 132 acre-feet and the greatest simulated water level, which would have been about 6.2 feet at the outlet culvert, but water was not likely to overflow an adjacent highway.

North Dakota↗

Assessment of Areal Recharge to the Spokane Valley-Rathdrum Prairie Aquifer, Spokane County, Washington, and Bonner and Kootenai Counties, Idaho

A numerical flow model of the Spokane Valley-Rathdrum Prairie aquifer currently (2007) being developed requires the input of values for areally-distributed recharge, a parameter that is often the most uncertain component of water budgets and ground-water flow models because it is virtually impossible to measure over large areas. Data from six active weather stations in and near the study area were used in four recharge-calculation techniques or approaches; the Langbein method, in which recharge is estimated on the basis of empirical data from other basins; a method developed by the U.S. Department of Agriculture (USDA), in which crop consumptive use and effective precipitation are first calculated and then subtracted from actual precipitation to yield an estimate of recharge; an approach developed as part of the Eastern Snake Plain Aquifer Model (ESPAM) Enhancement Project in which recharge is calculated on the basis of precipitation-recharge relations from other basins; and an approach in which reference evapotranspiration is calculated by the Food and Agriculture Organization (FAO) Penman-Monteith equation, crop consumptive use is determined (using a single or dual coefficient approach), and recharge is calculated. Annual recharge calculated by the Langbein method for the six weather stations was 4 percent of annual mean precipitation, yielding the lowest values of the methods discussed in this report, however, the Langbein method can be only applied to annual time periods. Mean monthly recharge calculated by the USDA method ranged from 53 to 73 percent of mean monthly precipitation. Mean annual recharge ranged from 64 to 69 percent of mean annual precipitation. Separate mean monthly recharge calculations were made with the ESPAM method using initial input parameters to represent thin-soil, thick-soil, and lava-rock conditions. The lava-rock parameters yielded the highest recharge values and the thick-soil parameters the lowest. For thin-soil parameters, calculated monthly recharge ranged from 10 to 29 percent of mean monthly precipitation and annual recharge ranged from 16 to 23 percent of mean annual precipitation. For thick-soil parameters, calculated monthly recharge ranged from 1 to 5 percent of mean monthly precipitation and mean annual recharge ranged from 2 to 4 percent of mean annual precipitation. For lava-rock parameters, calculated mean monthly recharge ranged from 37 to 57 percent of mean monthly precipitation and mean annual recharge ranged from 45 to 52 percent of mean annual precipitation. Single-coefficient (crop coefficient) FAO Penman-Monteith mean monthly recharge values were calculated for Spokane Weather Service Office (WSO) Airport, the only station for which the necessary meteorological data were available. Grass-referenced values of mean monthly recharge ranged from 0 to 81 percent of mean monthly precipitation and mean annual recharge was 21 percent of mean annual precipitation; alfalfa-referenced values of mean monthly recharge ranged from 0 to 85 percent of mean monthly precipitation and mean annual recharge was 24 percent of mean annual precipitation. Single-coefficient FAO Penman-Monteith calculations yielded a mean monthly recharge of zero during the eight warmest and driest months of the year (March-October). In order to refine the mean monthly recharge estimates, dual-coefficient (basal crop and soil evaporation coefficients) FAO Penman-Monteith dual-crop evapotranspiration and deep-percolation calculations were applied to daily values from the Spokane WSO Airport for January 1990 through December 2005. The resultant monthly totals display a temporal variability that is absent from the mean monthly values and demonstrate that the daily amount and timing of precipitation dramatically affect calculated recharge. The dual-coefficient FAO Penman-Monteith calculations were made for the remaining five stations using wind-speed values for Spokane WSO Airport and other assumptions regarding

Scientific Investigations Report↗

Reconnaissance study of the hydrology of American Memorial Park, Island of Saipan, Commonwealth of the Northern Mariana Islands

American Memorial Park, a unit of the National Park Service on the Island of Saipan, includes among its features a 27-acre estuarine system that has become a rarity within the Commonwealth of the Northern Mariana Islands. The estuarine system's mosaic of marshy areas interspersed with emergent wetlands and mixed wet forests provides critical habitat for various migratory and resident waterfowl, including two Federally listed endangered species: the Marianas gallinule ( Gallinula chloropus guami ) and the nightingale reed warbler ( Acrocephalus luscinia ). With sensitivity to the park's ecologic assets and the uncertainty associated with locally rapid urbanization, a need to better understand the hydrology of American Memorial Park was recognized. To address that need, a reconnaissance study of the park was undertaken during August and September 2005. The goals of the study were (1) to describe the occurrence and salinity of surface and ground water within the park; (2) to develop a hydrologic model of the park area of the island, with emphasis on the 27-acre estuarine system; and (3) to identify additional data needed to further develop this model. With regard to surface water, three freshwater inputs to the park's natural wetland are possible: direct rainfall, seaward-flowing ground water, and overland flow. Direct rainfall, which is an important source of freshwater to the wetland, commonly exceeds evapotranspiration both seasonally and per storm. The seaward flow of ground water is likely to be a source of freshwater to the wetland because ground water generally has an upward vertical component in the nearshore environment. Overland flow upgradient of the park could potentially contribute a significant input of freshwater during periods of intense rainfall, but roads that flank the park's perimeter act as a barrier to surficial inflows. During the reconnaissance, four discrete bodies, or zones, of surface water were observed within the park's natural wetland. Conductivity within these surface-water zones typically ranged from 1,540 to 4,370 microsiemens per centimeter (µS/cm) at 25°C although values as low as 829 and as high as 8,750 µ S/cm were measured. As a result of these observations, the American Memorial Park wetland area meets the definition criteria of an estuarine system that is dominantly oligohaline. Conductivity was also measured in a constructed wetland that was built within the park to augment the storm-drainage infrastructure of the village of Garapan. Reverse-osmosis facilities, in operation at hotels adjacent to the park, have historically discharged highly saline wastewater into the storm-drainage system. This collective storm and wastewater flow is routed into the constructed wetland and from there into the ocean. The conductivity of water in the constructed wetland ranged from 45,000 to 62,500 µS/cm , exceeding nominal seawater values by as much as 25 percent, with the highest conductivities recorded near discharging storm drains. With regard to ground water, the reconnaissance included installation of a ground-water-monitoring network. Data collected from this network identified the presence of freshwater underlying the park and indicated that surface water is directly connected to ground water in the natural wetland because the water levels of both surface water and ground water directly varied with the tide. Conductivities of ground-water samples from wells in the monitoring network indicated that ground-water salinity was geographically related: conductivities were lower (801-2,490 (µS/cm) in surficially dry areas, intermediate (6,090-9,180 (µS/cm) in natural-wetland areas, and higher (18,250-27,700 (µS/cm) in areas adjacent to the constructed wetland and its associated ocean-discharge channel. Synoptic water-level surveys were made to enhance understanding of the spatial expression of the water table; they were scheduled to overlap with peak and trough tidal signals to enable limited characteri

Scientific Investigations Report↗

Simulation of Surface-Water Conditions in the Nontidal Passaic River Basin, New Jersey

The Passaic River Basin, the third largest drainage basin in New Jersey, encompasses 950 mi2 (square miles) in the highly urbanized area outside New York City, with a population of 2 million. Water quality in the basin is affected by many natural and anthropogenic factors. Nutrient loading to the Wanaque Reservoir in the northern part of the basin is of particular concern and is caused partly by the diversion of water at two downstream intakes that is transferred back upstream to refill the reservoir. The larger of these diversions, Wanaque South intake, is on the lower Pompton River near Two Bridges, New Jersey. To support the development of a Total Maximum Daily Load (TMDL) for nutrients in the nontidal part of the basin (805 mi2), a water-quality transport model was needed. The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection and New Jersey EcoComplex, developed a flow-routing model to provide the hydraulic inputs to the water-quality model. The Diffusion Analogy Flow model (DAFLOW) described herein was designed for integration with the Water Quality Analysis Simulation Program (WASP) watershed water-quality model. The flow routing model was used to simulate flow in 108 miles of the Passaic River and major tributaries. Flow data from U.S. Geological Survey streamflow-gaging stations represent most of the model's upstream boundaries. Other model inputs include estimated flows for ungaged tributaries and unchanneled drainage along the mainstem, and reported flows for major point-source discharges and diversions. The former flows were calibrated using the drainage-area ratio method. The simulation extended over a 4+ year period representing a range in flow conditions. Simulated channel cross-sectional geometry in the DAFLOW model was calibrated using several different approaches by adjusting area and top width parameters. The model also was calibrated to observed flows for water year 2001 (low flow) at five mainstem gaging stations and one station at which flow was estimated. The model's target range was medium to low flows--the range of typical intake operations. Simulated flow mass balance, hydrographs (flood-wave speed, attenuation, and spread), flow-duration curves, and velocity and depth values were compared to observed counterparts. Mass balance and hydrograph fit were evaluated quantitatively. Simulation results generally were within the accuracy of the flow data at the measurement stations. The model was validated to observed flows for water years 2000 (average flow), 2002 (extreme low flow), and 2003 (high flow). Results for 19 of 20 comparisons indicate average mass-balance and model-fit errors of 6.6 and 15.7 percent, respectively, indicating that the model reasonably represents the time variation of streamflow in the nontidal Passaic River Basin. An algorithm (subroutine) also was developed for DAFLOW to simulate the hydraulic mixing that occurs near the Wanaque South intake upstream from the confluence of the Pompton and Passaic Rivers. The intake draws water from multiple sources, including effluent from a nearby wastewater-treatment plant, all of which have different phosphorus loads. The algorithm determines the proportion of flow from each source and operates within a narrow flow range. The equations used in the algorithm are based on the theory of diffusion and lateral mixing in rivers. Parameters used in the equations were estimated from limited available local flow and water-quality data. As expected, simulation results for water years 2000, 2001, and 2003 indicate that most of the water drawn to the intake comes from the Pompton River; however, during many short periods of low flow and high diversion, particularly in water year 2002, entrainment of the other flow sources compensated for the insufficient flow in the Pompton River. As additional verification of the flow model used in the water-quality model, a Branched Lagrangian Transport Model (B

Scientific Investigations Report↗

Analysis of salinity intrusion in the Waccamaw River and Atlantic Intracoastal Waterway near Myrtle Beach, South Carolina, 1995-2002

Six reservoirs in North Carolina discharge into the Pee Dee River, which flows 160 miles through South Carolina to the coastal communities near Myrtle Beach, South Carolina. During the Southeast's record-breaking drought from 1998 to 2003, salinity intrusions inundated a coastal municipal freshwater intake, limiting water supplies. To evaluate the effects of regulated flows of the Pee Dee River on salinity intrusion in the Waccamaw River and Atlantic Intracoastal Waterway, the South Carolina Department of Natural Resources and a consortium of stakeholders entered into a cooperative agreement with the U.S. Geological Survey to apply data-mining techniques to the long-term time series to analyze and simulate salinity dynamics near the freshwater intakes along the Grand Strand of South Carolina. Salinity intrusion in tidal rivers results from the interaction of three principal forces—streamflow, mean tidal water levels, and tidal range. To analyze, model, and simulate hydrodynamic behaviors at critical coastal gages, data-mining techniques were applied to over 20 years of hourly streamflow, coastal water-quality, and water-level data. Artificial neural network models were trained to learn the variable interactions that cause salinity intrusions. Streamflow data from the 18,300-square-mile basin were input to the model as time-delayed variables and accumulated tributary inflows. Tidal inputs to the models were obtained by decomposing tidal water-level data into a "periodic" signal of tidal range and a "chaotic" signal of mean water levels. The artificial neural network models were able to convincingly reproduce historical behaviors and generate alternative scenarios of interest. To make the models directly available to all stakeholders along the Pee Dee and Waccamaw Rivers and Atlantic Intracoastal Waterway, an easy-to-use decision support system (DSS) was developed as a spreadsheet application that integrates the historical database, artificial neural network models, model controls, streaming graphics, and model output. An additional feature is a built-in optimizer that dynamically calculates the amount of flow needed to suppress salinity intrusions as tidal ranges and water levels vary over days and months. This DSS greatly reduced the number of long-term simulations needed for stakeholders to determine the minimum flow required to adequately protect the freshwater intakes.

South Carolina↗

Nutrient Enrichment in Estuaries from Discharge of Shallow Ground Water, Mt. Desert Island, Maine

Nutrient enrichment from atmospheric deposition, agricultural activities, wildlife, and domestic sources is a concern at Acadia National Park because of the potential problem of water-quality degradation and eutrophication in its estuaries. Water-quality degradation has been observed at the Park?s Bass Harbor Marsh estuary but not in Northeast Creek estuary. Previous studies at Acadia National Park have estimated nutrient inputs to estuaries from atmospheric deposition and surface-water runoff, but the importance of shallow ground water that may contain nutrients derived from domestic or other sources is unknown. Northeast Creek and Bass Harbor Marsh estuaries were studied to (1) identify shallow ground-water seeps, (2) assess the chemistry of the water discharged from selected seeps, and (3) assess the chemistry of ground water in shallow ground-water hyporheic zones. The hyporheic zone is defined here as the region beneath and lateral to a stream bed, where there is mixing of shallow ground water and surface water. This study also provides baseline chemical data for ground water in selected bedrock monitoring wells and domestic wells on Mt. Desert Island. Water samples were analyzed for concentrations of nutrients, wastewater compounds, dissolved organic carbon, pH, dissolved oxygen, temperature and specific conductance. Samples from bedrock monitoring wells also were analyzed for alkalinity, major cations and anions, and trace metals. Shallow ground-water seeps to Northeast Creek and Bass Harbor Marsh estuaries at Acadia National Park were identified and georeferenced using aerial infrared digital imagery. Monitoring included the deployment of continuously recording temperature and specific conductance sensors in the seep discharge zone to access marine or freshwater signatures related to tidal flooding, gradient-driven shallow ground-water flow, or shallow subsurface flow related to precipitation events. Many potential shallow ground-water discharge zones were identified from aerial thermal imagery during flights in May and December 2003 in both estuaries. The occurrence of ground-water seeps was confirmed using continuous and discrete measurements of temperature and specific conductance in selected seeps and in the adjacent estuaries that showed salinity anomalies reflecting the input of freshwater in these complex tidal systems. Analysis of water samples from shallow ground water in the hyporheic zone and from ground-water seeps indicated the presence of elevated concentrations of dissolved nitrogen, compared to concentrations in the adjacent estuaries and surface-water tributaries draining into the estuaries. These findings indicate that shallow ground water is a source of dissolved nitrogen to the estuaries. Orthophosphate levels were low in ground water in the hyporheic zone in Bass Harbor Marsh, but somewhat higher in one hyporheic-zone well in Northeast Creek compared with the concentrations in both estuaries that were at or below detection limits. Household wastewater-related compounds were not detected in ground water in the hyporheic zone. Analysis of water samples from domestic and bedrock monitoring wells developed in fractured bedrock indicated that concentrations of dissolved nitrogen, phosphorus, and household wastewater-related compounds were typically at or below detection, suggesting that the aquifers sampled had not been contaminated from septic sources.

Scientific Investigations Report↗

Documentation of Computer Program INFIL3.0 - A Distributed-Parameter Watershed Model to Estimate Net Infiltration Below the Root Zone

This report documents the computer program INFIL3.0, which is a grid-based, distributed-parameter, deterministic water-balance watershed model that calculates the temporal and spatial distribution of daily net infiltration of water across the lower boundary of the root zone. The bottom of the root zone is the estimated maximum depth below ground surface affected by evapotranspiration. In many field applications, net infiltration below the bottom of the root zone can be assumed to equal net recharge to an underlying water-table aquifer. The daily water balance simulated by INFIL3.0 includes precipitation as either rain or snow; snowfall accumulation, sublimation, and snowmelt; infiltration into the root zone; evapotranspiration from the root zone; drainage and water-content redistribution within the root-zone profile; surface-water runoff from, and run-on to, adjacent grid cells; and net infiltration across the bottom of the root zone. The water-balance model uses daily climate records of precipitation and air temperature and a spatially distributed representation of drainage-basin characteristics defined by topography, geology, soils, and vegetation to simulate daily net infiltration at all locations, including stream channels with intermittent streamflow in response to runoff from rain and snowmelt. The model does not simulate streamflow originating as ground-water discharge. Drainage-basin characteristics are represented in the model by a set of spatially distributed input variables uniquely assigned to each grid cell of a model grid. The report provides a description of the conceptual model of net infiltration on which the INFIL3.0 computer code is based and a detailed discussion of the methods by which INFIL3.0 simulates the net-infiltration process. The report also includes instructions for preparing input files necessary for an INFIL3.0 simulation, a description of the output files that are created as part of an INFIL3.0 simulation, and a sample problem that illustrates application of the code to a field setting. Brief descriptions of the main program routine and of each of the modules and subroutines of the INFIL3.0 code, as well as definitions of the variables used in each subroutine, are provided in an appendix.

Scientific Investigations Report↗

Hydrologic and water-quality characterization and modeling of the Onondaga Lake Basin, Onondaga County, New York

Onondaga Lake in Onondaga County, New York, has been identified as one of the Nation’s most contaminated lakes as a result of industrial and sanitary-sewer discharges and stormwater nonpoint sources, and has received priority cleanup status under the national Water Resources Development Act of 1990. A basin-scale precipitation-runoff model of the Onondaga Lake basin was identified as a desirable water-resources management tool to better understand the processes responsible for the generation of loads of sediment and nutrients that are transported to Onondaga Lake. During 2003–07, the U.S. Geological Survey (USGS) developed a model based on the computer program, Hydrological Simulation Program–FORTRAN (HSPF), which simulated overland flow to, and streamflow in, the major tributaries of Onondaga Lake, and loads of sediment, phosphorus, and nitrogen transported to the lake. The simulation period extends from October 1997 through September 2003. The Onondaga Lake basin was divided into 107 subbasins and within these subbasins, the land area was apportioned among 19 pervious and impervious land types on the basis of land use and land cover, hydrologic soil group (HSG), and aspect. Precipitation data were available from three sources as input to the model. The model simulated streamflow, water temperature, concentrations of dissolved oxygen, and concentrations and loads of sediment, orthophosphate, total phosphorus, nitrate, ammonia, and organic nitrogen in the four major tributaries to Onondaga Lake–Onondaga Creek, Harbor Brook, Ley Creek, and Ninemile Creek. Simulated flows were calibrated to data from nine USGS streamflow-monitoring sites; simulated nutrient concentrations and loads were calibrated to data collected at six of the nine streamflow-monitoring sites. Water-quality samples were collected, processed, and analyzed by personnel from the Onondaga County Department of Water Environment Protection. Several time series of flow, and sediment and nutrient loads were generated for known sources of these constituents, including the Tully Valley mudboils (flow and sediment), Otisco Lake (flow and nutrients), the Marcellus wastewater-treatment plant (flow and nutrients), and springs from carbonate bedrock (flow). Runoff from the impervious sewered areas of the City of Syracuse was adjusted for the quantity that was treatable at the county wastewater-treatment plant; the excess flows were routed to nearby streams through combined-sanitary-and-storm-sewer overflows. The mitigative effects that the Onondaga Reservoir and Otisco Lake were presumed to have on loads of sediment and particulate constituents were simulated by adjustment of parameter values that controlled sediment settling rates, deposition, and scour in the reservoir and lake. Graphical representations of observed and simulated data, and relevant statistics, were compared to assess model performance. Simulated daily and monthly streamflows were rated “very good” (within 10 percent of observed flows) at all calibration sites, except Onondaga Creek at Cardiff, which was rated “fair” (10–15 percent difference). Simulations of monthly average water temperatures were rated “very good” (within 7 percent of observed temperatures) at all sites. No observed data were available by which to directly assess the model’s simulation of suspended sediment loads. Available measured total suspended solids data provided an indirect means of comparison but, not surprisingly, yielded only “fair” to “poor” ratings (greater than 30 percent difference) for simulated monthly sediment loads at half the water-quality calibration sites. Simulations of monthly orthophosphate loads ranged from “very good” (within 15 percent of measured loads) at three sites to “poor” (greater than 35 percent difference) at one site; simulations of ammonia nitrogen loads ranged from “very good” at one site to “fair” (25–35 percent difference) at two sites. Simulations of monthly total phosphorus, nitrate, and organic nitrogen loads were generally rated “very good” at all calibration sites. Sources of uncertainty in model results were identified, including (1) errors in precipitation data, (2) limitations in model structure, (3) nonuniqueness of values for highly sensitive parameters, (4) errors or bias in data used to calibrate the different components of the model, (5) misclassification of land-use and land-cover data, (6) changes in land use during the simulation period, (7) unidentified sources or sinks of chemical loads and water-quality processes that varied over time, and (8) differences in scale between large calibrated subbasins and small subbasins to which calibrated parameter values were transferred. Uncertainty in simulations of water-quality constituents was compounded by uncertainty in the processes on which the water-quality simulations were based. Therefore, sediment simulations were affected by uncertainty in the simulation of hydrology, and nutrient simulations were affected by uncertainty in both the hydrologic and sediment processes, as well as, in simulations of water temperature and dissolved oxygen concentrations. The calibrated model can be used to simulate scenarios that represent planned or hypothetical development and implementation of best-management practices in the Onondaga Lake basin and to assess the effects that these changes and practices are likely to have on rural and urban nonpoint sources of pollution to Onondaga Lake. Model results also can be used as input to a hydrodynamic model of Onondaga Lake that is being developed by Onondaga County and to prioritize areas of the basin where mitigative measures to decrease sediment and nutrient loads could provide the greatest benefits to Onondaga Lake.

New York↗

Simulations of ground-water flow and particle pathline analysis in the zone of contribution of a public-supply well in Modesto, eastern San Joaquin Valley, California

Shallow ground water in the eastern San Joaquin Valley is affected by high nitrate and uranium concentrations and frequent detections of pesticides and volatile organic compounds (VOC), as a result of ground-water development and intensive agricultural and urban land use. A single public-supply well was selected for intensive study to evaluate the dominant processes affecting the vulnerability of public-supply wells in the Modesto area. A network of 23 monitoring wells was installed, and water and sediment samples were collected within the approximate zone of contribution of the public-supply well, to support a detailed analysis of physical and chemical conditions and processes affecting the water chemistry in the well. A three-dimensional, steady-state local ground-water-flow and transport model was developed to evaluate the age of ground water reaching the well and to evaluate the vulnerability of the well to nonpoint source input of nitrate and uranium. Particle tracking was used to compute pathlines and advective travel times in the ground-water flow model. The simulated ages of particles reaching the public-supply well ranged from 9 to 30,000 years, with a median of 54 years. The age of the ground water contributed to the public-supply well increased with depth below the water table. Measured nitrate concentrations, derived primarily from agricultural fertilizer, were highest (17 milligrams per liter) in shallow ground water and decreased with depth to background concentrations of less than 2 milligrams per liter in the deepest wells. Because the movement of water is predominantly downward as a result of ground-water development, and because geochemical conditions are generally oxic, high nitrate concentrations in shallow ground water are expected to continue moving downward without significant attenuation. Simulated long-term nitrate concentrations indicate that concentrations have peaked and will decrease in the public-supply well during the next 100 years because of the low nitrate concentrations in recharge beneath the urban area and the increasing proportion of urban-derived ground water reaching the well. The apparent lag time between peak input concentrations and peak concentrations in the well is about 20 to 30 years. Measured uranium concentrations were also highest (45 micrograms per liter) in shallow ground water, and decreased with depth to background concentrations of about 0.5 microgram per liter. Naturally-occurring uranium adsorbed to aquifer sediments is mobilized by oxygen-rich, high-alkalinity water. Alkalinity increased in shallow ground water in response to agricultural development. As ground-water pumping increased in the 1940s and 1950s, this alkaline water moved downward through the ground-water flow system, mobilizing the uranium adsorbed to aquifer sediments. Ground water with high alkalinity and high uranium concentrations is expected to continue to move deeper in the system, resulting in increased uranium concentrations with depth in ground water. Because alkalinity (and correspondingly uranium) concentrations were high in shallow ground water beneath both the urban and the agricultural land, long-term uranium concentrations in the public-supply well are expected to increase as the proportion of uranium-affected water contributed to the well increases. Assuming that the alkalinity near the water table remains the same, the simulation of long-term alkalinity in the public-supply well indicates that uranium concentrations in the public-supply well will likely approach the maximum contaminant level; however, the time to reach this level is more than 100 years because of the significant proportion of old, unaffected water at depth that is contributed to the public-supply well.

California↗

Principal locations of major-ion, trace-element, nitrate, and Escherichia coli loading to Emigration Creek, Salt Lake County, Utah, October 2005

Housing development and recreational activity in Emigration Canyon have increased substantially since 1980, perhaps causing an observed decrease in water quality of this northern Utah stream located near Salt Lake City. To identify reaches of the stream that contribute to water-quality degradation, a tracer-injection and synoptic-sampling study was done to quantify mass loading of major ions, trace elements, nitrate, and Escherichia coli (E. coli) to the stream. The resulting mass-loading profiles for major ions and trace elements indicate both geologic and anthropogenic inputs to the stream, principally from tributary and spring inflows to the stream at Brigham Fork, Burr Fork, Wagner Spring, Emigration Tunnel Spring, Blacksmith Hollow, and Killyon Canyon. The pattern of nitrate loading does not correspond to the major-ion and trace-element loading patterns. Nitrate levels in the stream did not exceed water-quality standards at the time of synoptic sampling. The majority of nitrate mass loading can be considered related to anthropogenic input, based on the field settings and trends in stable isotope ratios of nitrogen. The pattern of E. coli loading does not correspond to the major-ion, trace-element, or nitrate loading patterns. The majority of E. coli loading was related to anthropogenic sources based on field setting, but a considerable part of the loading also comes from possible animal sources in Killyon Canyon, in Perkins Flat, and in Rotary Park. In this late summer sampling, E. coli concentrations only exceeded water-quality standards in limited sections of the study reach. The mass-loading approach used in this study provides a means to design future studies and to evaluate the loading on a catchment scale.

Utah↗

Water quality, hydrology, and simulated response to changes in phosphorus loading of Butternut Lake, Price and Ashland Counties, Wisconsin, with special emphasis on the effects of internal phosphorus loading in a polymictic lake

Butternut Lake is a 393-hectare, eutrophic to hypereutrophic lake in northcentral Wisconsin. After only minor improvements in water quality were observed following several actions taken to reduce the nutrient inputs to the lake, a detailed study was conducted from 2002 to 2007 by the U.S. Geological Survey to better understand how the lake functions. The goals of this study were to describe the water quality and hydrology of the lake, quantify external and internal sources of phosphorus, and determine the effects of past and future changes in phosphorus inputs on the water quality of the lake. Since the early 1970s, the water quality of Butternut Lake has changed little in response to nutrient reductions from the watershed. The largest changes were in near-surface total phosphorus concentrations: August concentrations decreased from about 0.09 milligrams per liter (mg/L) to about 0.05 mg/L, but average summer concentrations decreased only from about 0.055—0.060 mg/L to about 0.045 mg/L. Since the early 1970s, only small changes were observed in chlorophyll a concentrations and water clarity (Secchi depths). All major water and phosphorus sources, including the internal release of phosphorus from the sediments (internal loading), were measured directly, and minor sources were estimated to construct detailed water and phosphorus budgets for the lake during monitoring years (MY) 2003 and 2004. During these years, Butternut Creek, Spiller Creek, direct precipitation, small tributaries and near-lake drainage area, and ground water contributed about 62, 20, 8, 7, and 3 percent of the inflow, respectively. The average annual load of phosphorus to the lake was 2,540 kilograms (kg), of which 1,590 kg came from external sources (63 percent) and 945 kg came from the sediments in the lake (37 percent). Of the total external sources, Butternut Creek, Spiller Creek, small tributaries and near-lake drainage area, septic systems, precipitation, and ground water contributed about 63, 23, 9, 3, 1, and 1 percent, respectively. Because of the high internal phosphorus loading, the eutrophication models used in this study were unable to simulate the observed water-quality characteristics in the lake without incorporating this source of phosphorus. However, when internal loading of phosphorus was added to the BATHTUB model, it accurately simulated the average water-quality characteristics measured in MY 2003 and 2004. Model simulations demonstrated a relatively linear response between in-lake total phosphorus concentrations and external phosphorus loading; however, the changes in concentrations were smaller than the changes in external phosphorus loadings (about 25—40 percent of the change in phosphorus loading). Changes in chlorophyll a concentrations, the percentage of days with algal blooms, and Secchi depths were nonlinear and had a greater response to reductions in phosphorus loading than to increases in phosphorus loading. A 50-percent reduction in external phosphorus loading caused an 18-percent decrease in chlorophyll a concentrations, a 41-percent decrease in the percentage of days with algal blooms, and a 12-percent increase in Secchi depth. When the additional internal phosphorus loading was removed from model simulations, all of these constituents showed a much greater response to changes in external phosphorus loading. Because of Butternut Lake's morphometry, it is polymictic, which means it mixes frequently and does not develop stable thermal stratification throughout the summer. This characteristic makes it more vulnerable than dimictic lakes, which mix in spring and fall and develop stable thermal stratification during summer, to the high internal phosphorus loading that has resulted from historically high, nonnatural, external phosphorus loading. In polymictic lakes, the phosphorus released from the sediments is mixed into the upper part of the lake throughout summer. Once Butternut Lake became hypereutrophic (very productive), it became very difficult to alter its trophic state through reductions in external phosphorus loading because the high internal loading does not respond quickly to reductions in external nutrient loading. For Butternut Lake to become significantly less productive (change to a borderline mesotrophic/eutrophic state) a combined approach to reduce or eliminate internal phosphorus loading and reduce the external phosphorus loading by about 50 percent is needed.

Wisconsin↗

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

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

Scientific Investigations Report↗

Hydrologic and Water-Quality Conditions During Restoration of the Wood River Wetland, Upper Klamath River Basin, Oregon, 2003-05

Restoring previously drained wetlands is a strategy currently being used to improve water quality and decrease nutrient loading into Upper Klamath Lake, Oregon. In this 2003-05 study, ground- and surface-water quality and hydrologic conditions were characterized in the Wood River Wetland. Nitrogen and phosphorus levels, primarily as dissolved organic nitrogen and ammonium (NH4) and soluble reactive phosphorus (SRP), were high in surface waters. Dissolved organic carbon concentrations also were elevated in surface water, with median concentrations of 44 and 99 milligrams of carbon per liter (mg-C/L) in the North and South Units of the Wood River Wetland, respectively, reaching a maximum of 270 mg-C/L in the South Unit in late autumn. Artesian well water produced NH4 and SRP concentrations of about 6,000 micrograms per liter (ug/L), and concentrations of 36,500 ug-N/L NH4 and 4,110 ug-P/L SRP in one 26-28 ft deep piezometer well. Despite the high ammonium concentrations, the nitrate levels were moderate to low in wetland surface and ground waters. The surface-water concentrations of NH4 and SRP increased in spring and summer, outpacing those for chloride (a conservative tracer), indicative of evapoconcentration. In-situ chamber experiments conducted in June and August 2005 indicated a positive flux of NH4 and SRP from the wetland sediments. Potential sources of NH4 and SRP include diffusion of nutrients from decomposed peat, decomposing aquatic vegetation, or upwelling ground water. In addition to these inputs, evapoconcentration raised surface-water solute concentrations to exceedingly high values by the end of summer. The increase was most pronounced in the South Unit, where specific conductance reached 2,500 uS/cm and median concentrations of total nitrogen and total phosphorus reached 18,000-36,500 ug-N/L and about 18,000-26,000 ug-P/L, respectively. Water-column SRP and total phosphorus levels decreased during autumn and winter following inputs of irrigation water and precipitation, which have lower nutrient concentrations. The SRP concentrations, however, decreased faster than the dilution rate alone, possibly due to precipitation of phosphorus with iron, manganese, or calcium. The high concentrations of dissolved nitrogen and phosphorus during the growing season give rise to a rich plant community in the wetland consisting of emergent and submergent macrophytes and algae including phytoplankton and benthic and epiphytic algae that have pronounced effects on dissolved oxygen (DO) and pH. Midday readings of surface-water DO during summer often were supersaturated (as much as 310 percent saturation) with elevated pH (as much as 9.2 units), indicative of high rates of photosynthesis. Minimum DO concentrations in the shallow ground-water piezometer wells were 0.4 mg/L in the North Unit and 0.8 mg/L in the South Unit during summer, which is probably low enough to support microbial denitrification. Denitrification was confirmed during in-situ experiments conducted at the sediment-water interface, but rates were low due to low background nitrate (NO3). Nevertheless, denitrification (and plant uptake) likely contribute to low nitrate levels. Another possible cause of low nitrate levels is dissimilatory nitrate reduction to ammonia (DNRA), a microbial process that converts and decreases nitrate to ammonia. DNRA explains the excess ammonia production measured in the chambers treated with nitrate. Surface-water levels and standing surface-water volume in the Wood River Wetland reached a maximum in early spring, inundating 80-90 percent of the wetland. Surface-water levels and standing volume then declined reaching a minimum in August through November, when the South Unit was only 10 percent inundated and the North Unit was nearly dry. The shallow ground-water levels followed a trend similar to surface-water levels and indicated a strong upward gradient. A monthly water budget was developed individually for the North

Scientific Investigations Report↗

Primary Productivity in Meduxnekeag River, Maine, 2005

During August and September 2005, dissolved oxygen, temperature, pH, specific conductance, streamflow, and light intensity (LI) were determined continuously at six sites defining five reaches on Meduxnekeag River above and below Houlton, Maine. These data were collected as input for a dual-station whole-stream metabolism model to evaluate primary productivity in the river above and below Houlton. The river receives nutrients and organic matter from tributaries and the Houlton wastewater treatment plant (WWTP). Model output estimated gross and net primary productivity for each reach. Gross primary productivity (GPP) varied in each reach but was similar and positive among the reaches. GPP was correlated to LI in the four reaches above the WWTP but not in the reach below. Net primary productivity (NPP) decreased in each successive downstream reach and was negative in the lowest two reaches. NPP was weakly related to LI in the upper two reaches and either not correlated or negatively correlated in the lower three reaches. Relations among GPP, NPP, and LI indicate that the system is heterotrophic in the downstream reaches. The almost linear decrease in NPP (the increase in metabolism and respiration) indicates a cumulative effect of inputs of nutrients and organic matter from tributaries that drain agricultural land, the town of Houlton, and the discharges from the WWTP.

Scientific Investigations Report↗

Hydrogeologic Framework of Bedrock Units and Initial Salinity Distribution for a Simulation of Groundwater Flow for the Lake Michigan Basin

The U.S. Geological Survey is assessing groundwater availability in the Lake Michigan Basin. As part of the assessment, a variable-density groundwater-flow model is being developed to simulate the effects of groundwater use on water availability throughout the basin. The hydrogeologic framework for the Lake Michigan Basin model was developed by grouping the bedrock geology of the study area into hydrogeologic units on the basis of the functioning of each unit as an aquifer or confining layer within the basin. Available data were evaluated based on the areal extent of coverage within the study area, and procedures were established to characterize areas with sparse data coverage. Top and bottom altitudes for each hydrogeologic unit were interpolated in a geographic information system for input to the model and compared with existing maps of subsurface formations. Fourteen bedrock hydrogeologic units, making up 17 bedrock model layers, were defined, and they range in age from the Jurassic Period red beds of central Michigan to the Cambrian Period Mount Simon Sandstone. Information on groundwater salinity in the Lake Michigan Basin was compiled to create an input dataset for the variable-density groundwater-flow simulation. Data presented in this report are referred to as 'salinity data' and are reported in terms of total dissolved solids. Salinity data were not available for each hydrogeologic unit. Available datasets were assigned to a hydrogeologic unit, entered into a spatial database, and data quality was visually evaluated. A geographic information system was used to interpolate salinity distributions for each hydrogeologic unit with available data. Hydrogeologic units with no available data either were set equal to neighboring units or were vertically interpolated by use of values from units above and below.

Scientific Investigations Report↗