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

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↗

Chloride in Groundwater and Surface Water in Areas Underlain by the Glacial Aquifer System, Northern United States

A study of chloride in groundwater and surface water was conducted for the glacial aquifer system of the northern United States in forested, agricultural, and urban areas by analyzing data collected for the National Water-Quality Assessment Program from 1991 to 2004. Groundwater-quality data from a sampling of 1,329 wells in 19 states were analyzed. Chloride concentrations were greater than the secondary maximum contaminant level established by the U.S. Environmental Protection Agency of 250 milligrams per liter in 2.5 percent of samples from 797 shallow monitoring wells and in 1.7 percent of samples from 532 drinking-water supply wells. Water samples from shallow monitoring wells in urban areas had the largest concentration of chloride, followed by water samples from agricultural and forested areas (medians of 46, 12, and 2.9 milligrams per liter, respectively). An analysis of chloride:bromide ratios, by mass, and chloride concentrations compared to binary mixing curves for dilute groundwater, halite, sewage and animal waste, potassium chloride fertilizer, basin brines, seawater, and landfill leachate in samples from monitoring wells indicated multiple sources of chloride in samples from wells in urban areas and agricultural areas. Water from shallow monitoring wells in urban areas had the largest chloride:bromide ratio, and samples with chloride:bromide ratios greater than 1,000 and chloride concentrations greater than 100 milligrams per liter were dominated by halite; however, the samples commonly contained mixtures that indicated input from sewage or animal waste. Chloride:bromide ratios were significantly larger in samples from public-supply drinking-water wells than from private drinking-water wells, and ratios were significantly larger in all drinking-water wells in eastern and central regions of the glacial aquifer system than in west-central and western regions of the glacial aquifer system. Surface-water-quality data collected regularly during varying time periods from 1991-2004 from 100 basins dominated by forested, agricultural, or urban land in 15 states were analyzed to determine maximum measured chloride concentrations. Samples from 15 sites in east, central, and west-central areas, collected primarily in winter, had chloride concentrations higher than the U.S. Environmental Protection Agency recommended chronic criterion concentration for aquatic life of 230 milligrams per liter. Concentrations of chloride in base-flow samples were predictive of maximum measured chloride concentrations, indicating that inputs of chloride from groundwater and (or) point-source wastewater discharges increase the likelihood of samples exceeding the recommended chronic aquatic criterion. Multiple linear regression analyses showed that the density of major roads, potential evapotranspiration, and the percentage of annual runoff from saturated overland flow were significant factors in describing the range of maximum measured chloride concentrations in the basins studied. Chloride loads and yields were determined at 95 surface-water-monitoring stations in basins dominated by forested, agricultural, or urban land. Annual chloride yield was largest in the urban basins (median of 88 tons per square mile) and smallest in the forested basins (median of 6.4 tons per square mile). The median chloride yield in the agricultural basins was 15.4 tons per square mile. Multiple linear regression analyses showed that the density of highways (roads in U.S. highway system), the number of major wastewater discharges in the basin, potential evapotranspiration, and urban minus agricultural land area were significant factors in describing the range of average annual chloride yields. Upward trends in chloride loads were apparent in several urban basins for which additional long-term data were available. Increases in chloride loads over time may be related to a variety of factors, including increases in road area and consequent deicing, incr

Scientific Investigations Report↗

Sediment Loading from Crab Creek and Other Sources to Moses Lake, Washington, 2007 and 2008

The average sediment-accumulation rate on the bed of Moses Lake since 1980, based on the identification of Mount St. Helens ash in lakebed cores, was 0.24 inches per year. Summed over the lake surface area, the average sediment-accumulation rate on the lakebed is 190,000 tons per year. Based on USGS stream-gaging station data, the average annual sediment load to Moses Lake from Crab Creek was 32,000 tons per year between 1943 and 2008; the post Mount St. Helens eruption annual load from Crab Creek was calculated to be 13,000 tons per year. The total mass input from Crab Creek and other fluvially derived sediment sources since 1980 has been about 20,000 tons per year. Eolian sediment loading to Moses Lake was about 50,000 tons per year before irrigation and land-use development largely stabilized the Moses Lake dune field. Currently, eolian input to the lake is less than 2,000 tons per year. Considering all sediment sources to the lake, most (from 80 to 90 percent) of post-1980 lakebed-sediment accumulation is from autochthonous, or locally formed, mineral matter, including diatom frustuals and carbonate shells, derived from biogenic production in phytoplankton and zooplankton. Suspended-sediment samples collected from Crab Creek and similar nearby waterways in 2007 and 2008 combined with other USGS data from the region indicated that a proposed Bureau of Reclamation supplemental feed of as much as 650 cubic feet per second through Crab Creek might initially contain a sediment load of as much as 1,500 tons per day. With time, however, this sediment load would decrease to about 10 tons per day in the sediment-supply-limited creek as available sediment in the channel is depleted. Sediment loads in the supplemental feed ultimately would be similar to loads in other bypass canals near Moses Lake. Considering the hydrology and geomorphology of the creek over multiple years, there is little evidence that the proposed supplemental feed would substantially increase the overall sediment load from Crab Creek to Moses Lake relative to natural, background conditions. Because Moses Lake is relatively shallow and subject to significant wind-driven circulation currents, mixing also would redistribute some of the fluvial sediment load deposited from Crab Creek throughout Parker Horn and the rest of Moses Lake, further mitigating the local effect of Crab Creek sedimentation near the City of Moses Lake.

Scientific Investigations Report↗

Simulation of streamflow, evapotranspiration, and groundwater recharge in the lower San Antonio River Watershed, South-Central Texas, 2000-2007

The U.S. Geological Survey (USGS), in cooperation with the San Antonio River Authority, the Evergreen Underground Water Conservation District, and the Goliad County Groundwater Conservation District, configured, calibrated, and tested a watershed model for a study area consisting of about 2,150 square miles of the lower San Antonio River watershed in Bexar, Guadalupe, Wilson, Karnes, DeWitt, Goliad, Victoria, and Refugio Counties in south-central Texas. The model simulates streamflow, evapotranspiration (ET), and groundwater recharge using rainfall, potential ET, and upstream discharge data obtained from National Weather Service meteorological stations and USGS streamflow-gaging stations. Additional time-series inputs to the model include wastewater treatment-plant discharges, withdrawals for cropland irrigation, and estimated inflows from springs. Model simulations of streamflow, ET, and groundwater recharge were done for 2000-2007. Because of the complexity of the study area, the lower San Antonio River watershed was divided into four subwatersheds; separate HSPF models were developed for each subwatershed. Simulation of the overall study area involved running simulations of the three upstream models, then running the downstream model. The surficial geology was simplified as nine contiguous water-budget zones to meet model computational limitations and also to define zones for which ET, recharge, and other water-budget information would be output by the model. The model was calibrated and tested using streamflow data from 10 streamflow-gaging stations; additionally, simulated ET was compared with measured ET from a meteorological station west of the study area. The model calibration is considered very good; streamflow volumes were calibrated to within 10 percent of measured streamflow volumes. During 2000-2007, the estimated annual mean rainfall for the water-budget zones ranged from 33.7 to 38.5 inches per year; the estimated annual mean rainfall for the entire watershed was 34.3 inches. Using the HSPF model it was estimated that for 2000-2007, less than 10 percent of the annual mean rainfall on the study watershed exited the watershed as streamflow, whereas about 82 percent, or an average of 28.2 inches per year, exited the watershed as ET. Estimated annual mean groundwater recharge for the entire study area was 3.0 inches, or about 9 percent of annual mean rainfall. Estimated annual mean recharge was largest in water-budget zone 3, the zone where the Carrizo Sand outcrops. In water-budget zone 3, the estimated annual mean recharge was 5.1 inches or about 15 percent of annual mean rainfall. Estimated annual mean recharge was smallest in water-budget zone 6, about 1.1 inches or about 3 percent of annual mean rainfall. The Cibolo Creek subwatershed and the subwatershed of the San Antonio River upstream from Cibolo Creek had the largest and smallest basin yields, about 4.8 inches and 1.2 inches, respectively. Estimated annual ET and annual recharge generally increased with increasing annual rainfall. Also, ET was larger in zones 8 and 9, the most downstream zones in the watershed. Model limitations include possible errors related to model conceptualization and parameter variability, lack of data to quantify certain model inputs, and measurement errors. Uncertainty regarding the degree to which available rainfall data represent actual rainfall is potentially the most serious source of measurement error.

Texas↗

Suspended-Sediment Budget for the North Santiam River Basin, Oregon, Water Years 2005-08

Significant Findings An analysis of sediment transport in the North Santiam River basin during water years 2005-08 indicated that: Two-thirds of sediment input to Detroit Lake originated in the upper North Santiam River subbasin. Two-thirds of the sediment transported past Geren Island originated in the Little North Santiam River subbasin. The highest annual suspended-sediment load at any of the monitoring stations was the result of a debris flow on November 6, 2006, on Mount Jefferson. About 86 percent of the total sediment input to Detroit Lake was trapped in the lake, whereas 14 percent was transported farther downstream. More than 80 percent of the sediment transport in the basin was in November, December, and January. The variance in the annual suspended-sediment loads was better explained by the magnitude of the annual peak streamflow than by the annual mean streamflow.

Scientific Investigations Report↗

A Geochemical Mass-Balance Method for Base-Flow Separation, Upper Hillsborough River Watershed, West-Central Florida, 2003-2005 and 2009

Geochemical mass-balance (GMB) and conductivity mass-balance (CMB) methods for hydrograph separation were used to determine the contribution of base flow to total stormflow at two sites in the upper Hillsborough River watershed in west-central Florida from 2003-2005 and at one site in 2009. The chemical and isotopic composition of streamflow and precipitation was measured during selected local and frontal low- and high-intensity storm events and compared to the geochemical and isotopic composition of groundwater. Input for the GMB method included cation, anion, and stable isotope concentrations of surface water and groundwater, whereas input for the CMB method included continuous or point-sample measurement of specific conductance. The surface water is a calcium-bicarbonate type water, which closely resembles groundwater geochemically, indicating that much of the surface water in the upper Hillsborough River basin is derived from local groundwater discharge. This discharge into the Hillsborough River at State Road 39 and at Hillsborough River State Park becomes diluted by precipitation and runoff during the wet season, but retains the calcium-bicarbonate characteristics of Upper Floridan aquifer water. Field conditions limited the application of the GMB method to low-intensity storms but the CMB method was applied to both low-intensity and high-intensity storms. The average contribution of base flow to total discharge for all storms ranged from 31 to 100 percent, whereas the contribution of base flow to total discharge during peak discharge periods ranged from less than 10 percent to 100 percent. Although calcium, magnesium, and silica were consistent markers of Upper Floridan aquifer chemistry, their use in calculating base flow by the GMB method was limited because the frequency of point data collected in this study was not sufficient to capture the complete hydrograph from pre-event base-flow to post-event base-flow concentrations. In this study, pre-event water represented somewhat diluted groundwater. Streamflow conductivity integrates the concentrations of the major ions, and the logistics of acquiring specific conductance at frequent time intervals are less complicated than data collection, sample processing, shipment, and analysis of water samples in a laboratory. The acquisition of continuous specific conductance data reduces uncertainty associated with less-frequently collected geochemical point data.

Scientific Investigations Report↗

Arsenic, metals, and nutrients in runoff from two detention basins to Raccoon Creek, New Jersey Coastal Plain, 2008

Arsenic (As) concentrations in the waters of Raccoon Creek in southern New Jersey commonly exceed the State\'s Surface Water Quality Standard (SWQS) for freshwater of 0.017 microgram per liter (mu or ug/L). In order to assess contributions of As from residential runoff to the creek, samples of runoff water were collected from a detention basin in each of two residential developments underlain by different geologic formations and at the outlets of those basins. Samples of streamwater also were collected from Raccoon Creek adjacent to the developments. The samples were analyzed to determine concentrations of As, selected metals, organic carbon, and nutrients. Soil samples in and downgradient from the basins also were collected and analyzed. Concentrations of As in unfiltered water samples of runoff from the basin underlain by glauconitic clays generally were higher (up to 4.35 mu or ug/L) than in runoff from the basin underlain by predominantly quartz sands and silts (up to 2.68 mu or ug/L). Chromium (Cr) concentrations also were higher in runoff from the basin underlain by glauconitic clays than in runoff from the basin underlain by quartz sand and silt. In addition, Cr concentrations were higher in the glauconitic soils than in the quartz-rich soils. Metals such as aluminum (Al), iron (Fe), lead (Pb), and manganese (Mn) in the runoff and in the streamwater were mostly in particulate form. Arsenic, most metals, and phosphorus (P) however, were mostly in dissolved form in runoff but in particulate form in the streamwater. Total organic carbon concentrations in the runoff ranged from about 10 to nearly 16 milligrams per liter (mg/L). Given such levels of organic carbon and strong correlations between concentrations of some metals and organic carbon, it may be that many of the metals were complexed with dissolved organic carbon and transported in that form in the runoff. Although underlying geologic materials and soils appear to be major contributors of As to the streamwater, As also could have been contributed from lead arsenate pesticide residues. The residential development underlain by quartz-rich sediments formerly had been an orchard where such pesticides may have been used. The substantial inputs of As to runoff at this site may be attributable to this former land use, although Pb concentrations were about the same in runoff from both sites. The streamwater at both sites, however, contained Pb concentrations well above those in runoff, indicating that there are additional inputs of Pb, perhaps from roadside soils, upstream from the two sampling sites in this study. Positive relations between concentrations of As and some metals with dissolved organic carbon in runoff and streamwater indicate that complexation with organic carbon may provide a mechanism by which these constituents can be transported. Sorption of As, Pb, and P to Fe hydroxides may be indicated by the observed positive relation of particulate As, Pb, and P to particulate Fe, however, representing an additional mechanism for transport of these constituents.

New Jersey↗

Sedimentation and occurrence and trends of selected nutrients, other chemical constituents, and cyanobacteria in bottom sediment, Clinton Lake, northeast Kansas, 1977-2009

A combination of available bathymetric-survey information and bottom-sediment coring was used to investigate sedimentation and the occurrence of selected nutrients (total nitrogen and total phosphorus), organic and total carbon, 25 trace elements, cyanobacterial akinetes, and the radionuclide cesium-137 in the bottom sediment of Clinton Lake, northeast Kansas. The total estimated volume and mass of bottom sediment deposited from 1977 through 2009 in the conservation (multi-purpose) pool of the reservoir was 438 million cubic feet and 18 billion pounds, respectively. The estimated sediment volume occupied about 8 percent of the conservation-pool, water-storage capacity of the reservoir. Sedimentation in the conservation pool has occurred about 70 percent faster than originally projected at the time the reservoir was completed. Water-storage capacity in the conservation pool has been lost to sedimentation at a rate of about 0.25 percent annually. Mean annual net sediment deposition since 1977 in the conservation pool of the reservoir was estimated to be 563 million pounds per year. Mean annual net sediment yield from the Clinton Lake Basin was estimated to be 1.5 million pounds per square mile per year. Typically, the bottom sediment sampled in Clinton Lake was at least 99 percent silt and clay. The mean annual net loads of total nitrogen and total phosphorus deposited in the bottom sediment of Clinton Lake were estimated to be 1.29 million pounds per year and 556,000 pounds per year, respectively. The estimated mean annual net yields of total nitrogen and total phosphorus from the Clinton Lake Basin were 3,510 pounds per square mile per year and 1,510 pounds per square mile per year, respectively. Throughout the history of Clinton Lake, total nitrogen concentrations in the deposited sediment generally were uniform and indicated consistent inputs to the reservoir over time. Likewise, total phosphorus concentrations in the deposited sediment generally were uniform. Although, for two of three coring sites, a possible positive trend in phosphorus deposition was indicated. The Wakarusa River possibly was a larger contributor of nitrogen and phosphorus to Clinton Lake than was Rock Creek. As a principal limiting factor for primary production in most freshwater environments, phosphorus is of particular importance because increased inputs can contribute to accelerated reservoir eutrophication and the production of algal toxins and taste-and-odor compounds. Trace-element concentrations in the bottom sediment of Clinton Lake generally were uniform over time. As is typical for eastern Kansas reservoirs, arsenic, chromium, and nickel concentrations typically exceeded the threshold-effects guidelines, which represent the concentrations above which toxic biological effects occasionally occur. Zinc concentrations frequently exceeded the threshold-effects guideline. Trace-element concentrations did not exceed the probable-effects guidelines (available for eight trace elements), which represent the concentrations above which toxic biological effects usually or frequently occur. Cyanobacterial akinetes (cyanobacteria resting stage) in the bottom sediment of Clinton Lake, combined with historical water-quality data on chlorophyll-a and total phosphorus concentrations, indicated that the reservoir likely has been eutrophic throughout most of its history. A statistically significant increase in cyanobacterial akinetes in the bottom sediment indicated that Clinton Lake may have become more eutrophic over the life of the reservoir. The increase in cyanobacterial akinetes may, in part, be related to a possible increase in total phosphorus concentrations.

Scientific Investigations Report↗

A digital model for planning water management at Benton Lake National Wildlife Refuge, west-central Montana

Benton Lake National Wildlife Refuge is an important area for waterfowl production and migratory stopover in west-central Montana. Eight wetland units covering about 5,600 acres are the essential features of the refuge. Water availability for the wetland units can be uncertain owing to the large natural variations in precipitation and runoff and the high cost of pumping supplemental water. The U.S. Geological Survey, in cooperation with the U.S. Fish and Wildlife Service, has developed a digital model for planning water management. The model can simulate strategies for water transfers among the eight wetland units and account for variability in runoff and pumped water. This report describes this digital model, which uses a water-accounting spreadsheet to track inputs and outputs to each of the wetland units of Benton Lake National Wildlife Refuge. Inputs to the model include (1) monthly values for precipitation, pumped water, runoff, and evaporation; (2) water-level/capacity data for each wetland unit; and (3) the pan-evaporation coefficient. Outputs include monthly water volume and flooded surface area for each unit for as many as 5 consecutive years. The digital model was calibrated by comparing simulated and historical measured water volumes for specific test years.

Montana↗

Simulation of streamflows and basin-wide hydrologic variables over several climate-change scenarios, Methow River basin, Washington

The purpose of this project was to demonstrate the capabilities of an existing watershed model and downscaling procedures to provide simulated hydrological data over various greenhouse gas emission scenarios for use in the Methow River framework prototype. An existing watershed model was used to simulate daily time series of streamflow and basin-wide hydrologic variables for baseline conditions (1990–2000), and then for all combinations of three greenhouse gas emission scenarios and five general circulation models for future conditions (2008–2095). Input data for 18 precipitation and 17 temperature model input sites were generated using statistical techniques to downscale general circulation model data. The simulated results were averaged using an 11-year moving window to characterize the central year of the window to provide simulated data for water years 2008–2095. Simulation results indicate that substantial changes of monthly mean streamflows will occur. For all greenhouse gas emission scenarios, the future streamflows are greater in the winter than baseline conditions because a greater percentage of future precipitation is projected to fall as rain rather than as snow. The simulated future spring streamflows are less than baseline conditions because the spring snowpacks are smaller; therefore, flow contributions from snowmelt are less. A database was developed to automate model execution and to provide users with Internet access to voluminous data products ranging from summary figures to model output timeseries. Database-enabled Internet tools were developed to allow users to create interactive graphs of output results based on their analysis needs. For example, users were able to create graphs by selecting time intervals, greenhouse gas emission scenarios, general circulation models, and specific hydrologic variables.

Washington↗

Nutrient and suspended-sediment trends, loads, and yields and development of an indicator of streamwater quality at nontidal sites in the Chesapeake Bay watershed, 1985-2010

The U.S. Geological Survey (USGS) updates information on loads of, and trends in, nutrients and sediment annually to help the Chesapeake Bay Program (CBP) investigators assess progress toward improving water-quality conditions in the Chesapeake Bay and its watershed. CBP scientists and managers have worked since 1983 to improve water quality in the bay. In 2010, the U.S. Environmental Protection Agency (USEPA) established a Total Maximum Daily Load (TMDL) for the Chesapeake Bay. The TMDL specifies nutrient and sediment load allocations that need to be achieved in the watershed to improve dissolved oxygen, water-clarity, and chlorophyll conditions in the bay. The USEPA, USGS, and state and local jurisdictions in the watershed operate a CBP nontidal water-quality monitoring network and associated database that are used to update load and trend information to help assess progress toward reducing nutrient and sediment inputs to the bay. Data collected from the CBP nontidal network were used to estimate loads and trends for two time periods: a long-term period (1985-2010) at 31 "primary" sites (with storm sampling) and a 10-year period (2001-10) at 33 primary sites and 16 "secondary" sites (without storm sampling). In addition, loads at 64 primary sites were estimated for the period 2006 to 2010. Results indicate improving flow-adjusted trends for nitrogen and phosphorus for 1985 to 2010 at most of the sites in the network. For nitrogen, 21 of the 31 sites showed downward (improving) trends, whereas 2 sites showed upward (degrading) trends, and 8 sites showed no trends. The results for phosphorus were similar: 22 sites showed improving trends, 4 sites showed degrading trends, and 5 sites indicated no trends. For sediment, no trend was found at 40 percent of the sites, with 10 sites showing improving trends and 8 sites showing degrading trends. The USGS, working with CBP partners, developed a new water-quality indicator that combines the results of the 10-year trend analysis with results from a greater number of sites (64 primary sites) where loads and yields of total nitrogen and phosphorus and sediment could be calculated. The new indicator shows fewer significant trends for the 10-year time period than for the long-term time period (1985-2010). For 2001-10, total nitrogen trends were downward (improving) at 14 sites and upward (degrading) at 2 sites; no trend was found at 17 sites. For total phosphorus, 12 sites showed improving trends, 4 sites showed degrading trends, and 17 sites showed no trend. For total sediment, most sites (21) did not exhibit a significant trend; 3 sites showed improving trends, and 10 sites showed degrading trends. Few significant trends were seen at the 16 secondary sites: improving trends for total nitrogen at 4 sites, improving trends for total phosphorus at 2 sites, and a degrading trend for sediment at 1 site. Total streamflow to the Chesapeake Bay was 20 percent higher in 2010 than in 2009 and is considered to be within the normal range of flow, whereas annual streamflow at 28 sites was greater in 2010 than in 2009. No trends in daily streamflow were detected at the 31 long-term sites. Combined loads for the farthest downstream nontidal monitoring sites (called "River Input Monitoring sites") increased 33 percent for total nitrogen, 120 percent for total phosphorus, and 330 percent for total sediment from 2009 to 2010. The large increase in phosphorus and sediment loads in 2010 was caused in large part by two large storm events that occurred during the spring in the Potomac River Basin. Yields (load per watershed area) of total nitrogen in the Chesapeake Bay watershed decreased from north to south (New York to Virginia). No spatial patterns were discernible for total phosphorus or sediment.

Delaware, Maryland, New York, Pennsylvania, Virgin↗

Comparison of two regression-based approaches for determining nutrient and sediment fluxes and trends in the Chesapeake Bay watershed

Nutrient and sediment fluxes and changes in fluxes over time are key indicators that water resource managers can use to assess the progress being made in improving the structure and function of the Chesapeake Bay ecosystem. The U.S. Geological Survey collects annual nutrient (nitrogen and phosphorus) and sediment flux data and computes trends that describe the extent to which water-quality conditions are changing within the major Chesapeake Bay tributaries. Two regression-based approaches were compared for estimating annual nutrient and sediment fluxes and for characterizing how these annual fluxes are changing over time. The two regression models compared are the traditionally used ESTIMATOR and the newly developed Weighted Regression on Time, Discharge, and Season (WRTDS). The model comparison focused on answering three questions: (1) What are the differences between the functional form and construction of each model? (2) Which model produces estimates of flux with the greatest accuracy and least amount of bias? (3) How different would the historical estimates of annual flux be if WRTDS had been used instead of ESTIMATOR? One additional point of comparison between the two models is how each model determines trends in annual flux once the year-to-year variations in discharge have been determined. All comparisons were made using total nitrogen, nitrate, total phosphorus, orthophosphorus, and suspended-sediment concentration data collected at the nine U.S. Geological Survey River Input Monitoring stations located on the Susquehanna, Potomac, James, Rappahannock, Appomattox, Pamunkey, Mattaponi, Patuxent, and Choptank Rivers in the Chesapeake Bay watershed. Two model characteristics that uniquely distinguish ESTIMATOR and WRTDS are the fundamental model form and the determination of model coefficients. ESTIMATOR and WRTDS both predict water-quality constituent concentration by developing a linear relation between the natural logarithm of observed constituent concentration and three explanatory variables—the natural log of discharge, time, and season. ESTIMATOR uses two additional explanatory variables—the square of the log of discharge and time-squared. Both models determine coefficients for variables for a series of estimation windows. ESTIMATOR establishes variable coefficients for a series of 9-year moving windows; all observed constituent concentration data within the 9-year window are used to establish each coefficient. Conversely, WRTDS establishes variable coefficients for each combination of discharge and time using only observed concentration data that are similar in time, season, and discharge to the day being estimated. As a result of these distinguishing characteristics, ESTIMATOR reproduces concentration-discharge relations that are closely approximated by a quadratic or linear function with respect to both the log of discharge and time. Conversely, the linear model form of WRTDS coupled with extensive model windowing for each combination of discharge and time allows WRTDS to reproduce observed concentration-discharge relations that are more sinuous in form. Another distinction between ESTIMATOR and WRTDS is the reporting of uncertainty associated with the model estimates of flux and trend. ESTIMATOR quantifies the standard error of prediction associated with the determination of flux and trends. The standard error of prediction enables the determination of the 95-percent confidence intervals for flux and trend as well as the ability to test whether the reported trend is significantly different from zero (where zero equals no trend). Conversely, WRTDS is unable to propagate error through the many (over 5,000) models for unique combinations of flow and time to determine a total standard error. As a result, WRTDS flux estimates are not reported with confidence intervals and a level of significance is not determined for flow-normalized fluxes. The differences between ESTIMATOR and WRTDS, with regard to model form and determination of model coefficients, have an influence on the determination of nutrient and sediment fluxes and associated changes in flux over time as a result of management activities. The comparison between the model estimates of flux and trend was made for combinations of five water-quality constituents at nine River Input Monitoring stations. The major findings with regard to nutrient and sediment fluxes are as follows: (1)WRTDS produced estimates of flux for all combinations that were more accurate, based on reduction in root mean squared error, than flux estimates from ESTIMATOR; (2) for 67 percent of the combinations, WRTDS and ESTIMATOR both produced estimates of flux that were minimally biased compared to observed fluxes(flux bias = tendency to over or underpredict flux observations); however, for 33 percent of the combinations, WRTDS produced estimates of flux that were considerably less biased (by at least 10 percent) than flux estimates from ESTIMATOR; (3) the average percent difference in annual fluxes generated by ESTIMATOR and WRTDS was less than 10 percent at 80 percent of the combinations; and (4) the greatest differences related to flux bias and annual fluxes all occurred for combinations where the pattern in observed concentration-discharge relation was sinuous (two points of inflection) rather than linear or quadratic (zero or one point of inflection). The major findings with regard to trends are as follows: (1) both models produce water-quality trends that have factored in the year-to-year variations in flow; (2) trends in water-quality condition are represented by ESTIMATOR as a trend in flow-adjusted concentration and by WRTDS as a flow normalized flux; (3) for 67 percent of the combinations with trend estimates, the WRTDS trends in flow-normalized flux are in the same direction and magnitude to the ESTIMATOR trends in flow-adjusted concentration, and at the remaining 33 percent the differences in trend magnitude and direction are related to fundamental differences between concentration and flux; and (4) the majority (85 percent) of the total nitrogen, nitrate, and orthophosphorus combinations exhibited long-term (1985 to 2010) trends in WRTDS flow-normalized flux that indicate improvement or reduction in associated flux and the majority (83 percent) of the total phosphorus (from 1985 to 2010) and suspended sediment (from 2001 to 2010) combinations exhibited trends in WRTDS flow-normalized flux that indicate degradation or increases in the flux delivered.

Delaware, Maryland, New York, Pennsylvania, Virgin↗

Estimates of gains and losses from unmeasured sources and sinks for streamflow and dissolved-solids load in selected reaches of the Arkansas River, southeastern Colorado, 2009-2010

The Arkansas River is an important municipal water supply and is the primary supply for about 400,000 acres of irrigated land in southeastern Colorado. The suitability of this water for domestic, agricultural, and industrial use is affected by high salinity in parts of the Arkansas River. There is a need to quantify mass loading of dissolved solids (DS) in the Arkansas River. In 2009, the U.S. Geological Survey, in cooperation with the Arkansas River Basin Regional Resource Planning Group and the Colorado Water Conservation Board, began a study to estimate gains and losses from unmeasured sources and sinks for streamflow and DS load in selected reaches of the Arkansas River in southeastern Colorado. Two study reaches were selected for investigation—Canon City to just upstream from Pueblo Reservoir (UARB) and Avondale to Las Animas (LARB). The results from the water-budget analyses indicated that potential areas of unmeasured sources and sinks of streamflow were identifiable in the two study reaches. In the UARB, a substantial volume of water in the subreach from Ark at Canon City to the seasonal gaging station 5 miles downstream (Ark nr Canon City) was unaccounted for by the methodology used in this analysis. The daily gain from unmeasured sources in this subreach was estimated to be about 100 cubic feet per second (ft 3 /s) or about 20 ft 3 /s per river mile. Water-budget estimates for the remaining 18 miles of the UARB study reach indicated that gains or losses from unmeasured sources or sinks were within the measurement error as defined for this report. In the LARB, gains and losses from unmeasured sources and sinks were identified in some of the subreaches but the magnitude of the flux generally was small. Unmeasured sources ranging from less than 2 to 3 ft 3 /s per mile were identified in the river subreaches from Ark at Catlin Dam downstream to Ark at Swink. A streamflow loss was indicated along the subreach from Ark at Nepesta to Ark at Catlin Dam, particularly in 2010. The mechanism and spatial extent of this sink was not identified, and further investigation would be required to better quantify the loss. The results from the analyses of unmeasured sources of DS load indicated that potential source areas were identifiable in the study areas. It might be expected that unmeasured DS load flux would be identified along the same reaches where unmeasured streamflow flux was identified. To that extent, some of the observed results from the analysis of daily DS loading did mirror the streamflow results. In some subreaches of the Arkansas River, however, unmeasured sources and sinks of DS load did not appear to be directly associated with unmeasured sources and sinks of streamflow. In the UARB from Ark at Canon City to Ark nr Canon City, unmeasured gains in DS load were estimated to range from 11 to 22 tons per day per mile in 2009 and from about 8 to 13 tons per day per mile in 2010; streamflow from unmeasured sources was estimated to be about 20 ft 3 /s per mile along this same reach. Downstream from this short reach, DS load to the river from unmeasured sources was estimated to range from 5.4 to 7.6 tons per day per mile in 2010 for Ark nr Canon City to Ark at Portland and from 11 to 16 tons per day per mile in 2009 for Ark at Portland to Ark nr Portland. Unmeasured gains in streamflow were not identified in either of these subreaches. Several small tributaries with DS concentrations ranging from 3,000 mg/L to as high as 6,000 mg/L enter the river along these subreaches. These inputs may indicate a potential source of groundwater that could affect DS loading in the river. Further investigation would be needed to identify the unmeasured source or sources of DS load to determine the nature and extent of unmeasured inputs. In the LARB, gains in DS load from unmeasured sources were identified for the subreach from Ark nr Avondale to Ark at Nepesta, although no substantial amounts of streamflow from unmeasured sources were identified for this subreach. In 2009, the estimated gain in DS load from unmeasured sources for this subreach was 4.7 tons per day per mile. An increase in DS load from unmeasured sources also was identified along the subreach of the river from Ark at Catlin to Swink; the DS load from unmeasured sources was estimated to range from 10 to 28 tons per day per mile. The only loss of DS load was identified for the subreach from Nepesta to Catlin Dam in 2010. The mechanism and spatial extent of the losses were not identified, and further investigation would be required to better understand the results.

Colorado↗

Evaluation of the groundwater-flow model for the Ohio River alluvial aquifer near Carrollton, Kentucky, updated to conditions in September 2010

The Ohio River alluvial aquifer near Carrollton, Ky., is an important water resource for the cities of Carrollton and Ghent, as well as for several industries in the area. The groundwater of the aquifer is the primary source of drinking water in the region and a highly valued natural resource that attracts various water-dependent industries because of its quantity and quality. This report evaluates the performance of a numerical model of the groundwater-flow system in the Ohio River alluvial aquifer near Carrollton, Ky., published by the U.S. Geological Survey in 1999. The original model simulated conditions in November 1995 and was updated to simulate groundwater conditions estimated for September 2010. The files from the calibrated steady-state model of November 1995 conditions were imported into MODFLOW-2005 to update the model to conditions in September 2010. The model input files modified as part of this update were the well and recharge files. The design of the updated model and other input files are the same as the original model. The ability of the updated model to match hydrologic conditions for September 2010 was evaluated by comparing water levels measured in wells to those computed by the model. Water-level measurements were available for 48 wells in September 2010. Overall, the updated model underestimated the water levels at 36 of the 48 measured wells. The average difference between measured water levels and model-computed water levels was 3.4 feet and the maximum difference was 10.9 feet. The root-mean-square error of the simulation was 4.45 for all 48 measured water levels. The updated steady-state model could be improved by introducing more accurate and site-specific estimates of selected field parameters, refined model geometry, and additional numerical methods. Collection of field data to better estimate hydraulic parameters, together with continued review of available data and information from area well operators, could provide the model with revised estimates of conductance values for the riverbed and valley wall, hydraulic conductivities for the model layer, and target water levels for future simulations. Additional model layers, a redesigned model grid, and revised boundary conditions could provide a better framework for more accurate simulations. Additional numerical methods would identify possible parameter estimates and determine parameter sensitivities.

Kentucky↗