USGS ScienceSearch

SEARCH · USGS Science

Results for “Water Surface Techniques”

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 307 records · Page 17Linked to original sources

Regional evaluation of evapotranspiration in the Everglades

Understanding the water budget of the Everglades system is crucial to the success of restoration and management strategies. Although the water budget is simple in concept, it is difficult to assess quantitatively. Models used to simulate changes in water levels and vegetation resulting from management strategies need to accurately simulate all components of the water budget. One of the most important components of the Everglades water budget is evapotranspiration (ET). ET is water removed from the surface and soils by direct evaporation and plant transpiration. In South Florida, ET rates may exceed 40 in/yr (inches per year) on the average; during dry years, the ET could exceed rainfall (about 50 in/yr). Thus, most of the water that falls on the land surface as rainfall is returned to the atmosphere by ET. Despite the importance of ET in the Everglades water budget, our knowledge of ET is, at present, only semi-quantitative. Recent advances in instrumentation and measurement techniques have made it possible to continuously measure ET, so that an accurate evaluation of ET in the Everglades can be made.

Florida

Pesticides data-collection activities of the U.S. Geological Survey in Texas

The water-resources investigations program of the Geological Survey has the objective of seeing that the Nation's water resources are appraised and that the necessary water data to develop and manage them efficiently are available when needed. In carrying out this objective, the Survey is the primary source of hydrologic data that are basic not only to the quantification of the water resource but also to its development and conservation. In recent years, the need has developed for the Survey to provide data and technology to predict the impact of man's activities upon the water resource and the water environment. To accomplish its mission, the Geological Survey determines the source, quantity, quality, distribution, movement, and availability of both surface and ground water. This work includes the assessment of the chemical, physical, and biological quality of water resources and its relation to various parts of the hydrologic cycle; research to improve the scientific basis of investigations and techniques; design and operation of a national water-data network; and publication of its investigations.

Texas

Development and Implementation of a Transport Method for the Transport and Reaction Simulation Engine (TaRSE) based on the Godunov-Mixed Finite Element Method

A model to simulate transport of materials in surface water and ground water has been developed to numerically approximate solutions to the advection-dispersion equation. This model, known as the Transport and Reaction Simulation Engine (TaRSE), uses an algorithm that incorporates a time-splitting technique where the advective part of the equation is solved separately from the dispersive part. An explicit finite-volume Godunov method is used to approximate the advective part, while a mixed-finite element technique is used to approximate the dispersive part. The dispersive part uses an implicit discretization, which allows it to run stably with a larger time step than the explicit advective step. The potential exists to develop algorithms that run several advective steps, and then one dispersive step that encompasses the time interval of the advective steps. Because the dispersive step is computationally most expensive, schemes can be implemented that are more computationally efficient than non-time-split algorithms. This technique enables scientists to solve problems with high grid Peclet numbers, such as transport problems with sharp solute fronts, without spurious oscillations in the numerical approximation to the solution and with virtually no artificial diffusion.

Scientific Investigations Report

Simulation of groundwater conditions and streamflow depletion to evaluate water availability in a Freeport, Maine, watershed

In order to evaluate water availability in the State of Maine, the U.S. Geological Survey (USGS) and the Maine Geological Survey began a cooperative investigation to provide the first rigorous evaluation of watersheds deemed "at risk" because of the combination of instream flow requirements and proportionally large water withdrawals. The study area for this investigation includes the Harvey and Merrill Brook watersheds and the Freeport aquifer in the towns of Freeport, Pownal, and Yarmouth, Maine. A numerical groundwater- flow model was used to evaluate groundwater withdrawals, groundwater-surface-water interactions, and the effect of water-management practices on streamflow. The water budget illustrates the effect that groundwater withdrawals have on streamflow and the movement of water within the system. Streamflow measurements were made following standard USGS techniques, from May through September 2009 at one site in the Merrill Brook watershed and four sites in the Harvey Brook watershed. A record-extension technique was applied to estimate long-term monthly streamflows at each of the five sites. The conceptual model of the groundwater system consists of a deep, confined aquifer (the Freeport aquifer) in a buried valley that trends through the middle of the study area, covered by a discontinuous confining unit, and topped by a thin upper saturated zone that is a mixture of sandy units, till, and weathered clay. Harvey and Merrill Brooks flow southward through the study area, and receive groundwater discharge from the upper saturated zone and from the deep aquifer through previously unknown discontinuities in the confining unit. The Freeport aquifer gets most of its recharge from local seepage around the edges of the confining unit, the remainder is received as inflow from the north within the buried valley. Groundwater withdrawals from the Freeport aquifer in the study area were obtained from the local water utility and estimated for other categories. Overall, the public-supply withdrawals (105.5 million gallons per year (Mgal/yr)) were much greater than those for any other category, being almost 7 times greater than all domestic well withdrawals (15.3 Mgal/yr). Industrial withdrawals in the study area (2.0 Mgal/yr) are mostly by a company that withdraws from an aquifer at the edge of the Merrill Brook watershed. Commercial withdrawals are very small (1.0 Mgal/yr), and no irrigation or other agricultural withdrawals were identified in this study area. A three-dimensional, steady-state groundwater-flow model was developed to evaluate stream-aquifer interactions and streamflow depletion from pumping, to help refine the conceptual model, and to predict changes in streamflow resulting from changes in pumping and recharge. Groundwater levels and flow in the Freeport aquifer study area were simulated with the three-dimensional, finite-difference groundwater-flow modeling code, MODFLOW-2005. Study area hydrology was simulated with a 3-layer model, under steady-state conditions. The groundwater model was used to evaluate changes that could occur in the water budgets of three parts of the local hydrologic system (the Harvey Brook watershed, the Merrill Brook watershed, and the buried aquifer from which pumping occurs) under several different climatic and pumping scenarios. The scenarios were (1) no pumping well withdrawals; (2) current (2009) pumping, but simulated drought conditions (20-percent reduction in recharge); (3) current (2009) recharge, but a 50-percent increase in pumping well withdrawals for public supply; and (4) drought conditions and increased pumping combined. In simulated drought situations, the overall recharge to the buried valley is about 15 percent less and the total amount of streamflow in the model area is reduced by about 19 percent. Without pumping, infiltration to the buried valley aquifer around the confining unit decreased by a small amount (0.05 million gallons per day (Mgal/d)), and discharge to the streams increased by about 8 percent (0.3 Mgal/d). A 50-percent increase in pumping resulted in a simulated decrease in streamflow discharge of about 4 percent (0.14 Mgal/d). Streamflow depletion in Harvey Brook was evaluated by use of the numerical groundwater-flow model and an analytical model. The analytical model estimated negligible depletion from Harvey Brook under current (2009) pumping conditions, whereas the numerical model estimated that flow to Harvey Brook decreased 0.38 cubic feet per second (ft 3 /s) because of the pumping well withdrawals. A sensitivity analysis of the analytical model method showed that conducting a cursory evaluation using an analytical model of streamflow depletion using available information may result in a very wide range in results, depending on how well the hydraulic conductivity variables and aquifer geometry of the system are known, and how well the aquifer fits the assumptions of the model. Using the analytical model to evaluate the streamflow depletion with an incomplete understanding of the hydrologic system gave results that seem unlikely to reflect actual streamflow depletion in the Freeport aquifer study area. In contrast, the groundwater-flow model was a more robust method of evaluating the amount of streamflow depletion that results from withdrawals in the Freeport aquifer, and could be used to evaluate streamflow depletion in both streams. Simulations of streamflow without pumping for each measurement site were compared to the calibratedmodel streamflow (with pumping), the difference in the total being streamflow depletion. Simulations without pumping resulted in a simulated increase in the steady-state flow rate of 0.38 ft 3 /s in Harvey Brook and 0.01 ft 3 /s in Merrill Brook. This translates into a streamflow-depletion amount equal to about 8.5 percent of the steady-state base flow in Harvey Brook, and an unmeasurable amount of depletion in Merrill Brook. If pumping was increased by 50 percent and recharge reduced by 20 percent, the amount of streamflow depletion in Harvey Brook could reach 1.41 ft 3 /s.

Maine

General introduction for the “National Field Manual for the Collection of Water-Quality Data”

Background As part of its mission, the U.S. Geological Survey (USGS) collects data to assess the quality of our Nation’s water resources. A high degree of reliability and standardization of these data are paramount to fulfilling this mission. Documentation of nationally accepted methods used by USGS personnel serves to maintain consistency and technical quality in data-collection activities. “The National Field Manual for the Collection of Water-Quality Data” (NFM) provides documented guidelines and protocols for USGS field personnel who collect water-quality data. The NFM provides detailed, comprehensive, and citable procedures for monitoring the quality of surface water and groundwater. Topics in the NFM include (1) methods and protocols for sampling water resources, (2) methods for processing samples for analysis of water quality, (3) methods for measuring field parameters, and (4) specialized procedures, such as sampling water for low levels of mercury and organic wastewater chemicals, measuring biological indicators, and sampling bottom sediment for chemistry. Personnel who collect water-quality data for national USGS programs and projects, including projects supported by USGS cooperative programs, are mandated to use protocols provided in the NFM per USGS Office of Water Quality Technical Memorandum 2002.13. Formal training, for example, as provided in the USGS class, “Field Water-Quality Methods for Groundwater and Surface Water,” and field apprenticeships supplement the guidance provided in the NFM and ensure that the data collected are high quality, accurate, and scientifically defensible.

Techniques and Methods

Development, Testing, and Sensitivity and Uncertainty Analyses of a Transport and Reaction Simulation Engine (TaRSE) for Spatially Distributed Modeling of Phosphorus in South Florida Peat Marsh Wetlands

Alterations to the predevelopment delivery of water and nutrients into the Everglades of southern Florida have been occurring for nearly a century. Major regional drainage projects, large-scale agricultural development, and changes to the hydrology of the Kissimmee River-Lake Okeechobee watershed have resulted in substantial phosphorus transport increases by surface waters. Excess phosphorus has accumulated in the soils of northern Everglades marshes to levels that have impaired the natural resources of the region. Regulations now limit the amount of phosphorous that enters the Everglades through an extensive network of water-control structures. This study involved the development and application of water-quality modeling components that may be applied to existing hydrologic models of southern Florida to evaluate the effects of different management scenarios. The result of this work is a spatially distributed water-quality model for phosphorus transport and cycling in wetlands. The model solves the advection-dispersion equation on an unstructured triangular mesh and incorporates a wide range of user-selectable mechanisms for phosphorus uptake and release parameters. In general, the phosphorus model contains transfers between stores; examples of stores that can be included are soil, water column (solutes), pore water, macrophytes, suspended solids (plankton), and biofilm. Examples of transfers are growth, senescence, settling, diffusion, and so forth, described with first order, second order, and Monod types of transformations. Local water depths and velocities are determined from an existing two-dimensional, overland-flow hydrologic model. The South Florida Water Management District Regional Simulation Model was used in this study. The model is applied to three case studies: intact cores of wetland soils with water, outdoor mesocosoms, and a large constructed wetland; namely, Cell 4 of Stormwater Treatment Area 1 West (STA-1W Cell 4). Different levels of complexity in the phosphorus cycling mechanisms were simulated in these case studies using different combinations of phosphorus reaction equations. Changes in water column phosphorus concentrations observed under the controlled conditions of laboratory incubations, and mesocosm studies were reproduced with model simulations. Short-term phosphorus flux rates and changes in phosphorus storages were within the range of values reported in the literature, whereas unknown rate constants were used to calibrate the model output. In STA-1W Cell 4, the dominant mechanism for phosphorus flow and transport is overland flow. Over many life cycles of the biological components, however, soils accrue and become enriched in phosphorus. Inflow total phosphorus concentrations and flow rates for the period between 1995 and 2000 were used to simulate Cell 4 phosphorus removal, outflow concentrations, and soil phosphorus enrichment over time. This full-scale application of the model successfully incorporated parameter values derived from the literature and short-term experiments, and reproduced the observed long-term outflow phosphorus concentrations and increased soil phosphorus storage within the system. A global sensitivity and uncertainty analysis of the model was performed using modern techniques such as a qualitative screening tool (Morris method) and the quantitative, variance-based, Fourier Amplitude Sensitivity Test (FAST) method. These techniques allowed an in-depth exploration of the effect of model complexity and flow velocity on model outputs. Three increasingly complex levels of possible application to southern Florida were studied corresponding to a simple soil pore-water and surface-water system (level 1), the addition of plankton (level 2), and of macrophytes (level 3). In the analysis for each complexity level, three surface-water velocities were considered that each correspond to residence times for the selected area (1-kilometer long) of 2, 10, and 20

Scientific Investigations Report

Hydrologic budget and conditions of Permian, Pennsylvanian, and Mississippian aquifers in the Appalachian Plateaus physiographic province

In response to challenges to groundwater availability posed by historic land-use practices, expanding development of hydrocarbon resources, and drought, the U.S. Geological Survey Groundwater Resources Program began a regional assessment of the Appalachian Plateaus aquifers in 2013 that incorporated a hydrologic landscape approach to estimate all components of the hydrologic system: surface runoff, base flow from groundwater, and interaction with atmospheric water (precipitation and evapotranspiration). This assessment was intended to complement other Federal and State investigations and provide foundational groundwater-related datasets in the Appalachian Plateaus. A regional Soil-Water-Balance model was constructed for a 160,000-square-mile study area that extended to the topographic divide of all streams originating outside but flowing into areas underlain by Appalachian Plateaus aquifers. The model incorporated soil, landscape, and climate variables to estimate an annual water budget for the 32-year period from 1980 to 2011 and was calibrated using base-flow data estimated by hydrograph separation techniques from 20 streamflow gaging stations across the study area. Over this period, an average of 47 inches per year (in/yr) of precipitation fell on Appalachian Plateaus aquifers. Simulations from the regional Soil-Water-Balance model indicate that only 19 percent of the precipitation or an average 9 in/yr recharged aquifers, and 19 percent resulted in surface runoff to streams. The remaining 62 percent, an average of 27 in/yr of water, was returned to the atmosphere via evapotranspiration. Because withdrawals from aquifers due to pumping equated to less than 1 percent of the water budget, differences in predevelopment and postdevelopment regional water budgets of the Appalachian Plateaus were minimal. Storage changes caused by filling of abandoned coal-mine aquifers and long-term differences in aquifer storage resulting from climate fluctuations constitute a small portion of the overall water budget. The percentage of precipitation that results in recharge, runoff, or evapotranspiration from the landscape varies annually by up to a factor of two depending on temporal changes in prevailing climate conditions and spatial changes in basin characteristics, precipitation patterns, and sources of atmospheric moisture over a large study area. A comparison of water-budget estimates from the regional Soil-Water-Balance model for a dry year (1988) and wet year (2004) showed that evapotranspiration accounts for most of the annual differences in precipitation. As a portion of annual precipitation, evapotranspiration ranged from 69 percent (dry year) to 52 percent (wet year), a range four times greater than the 15 percent (dry year) to 18 percent (wet year) range estimated for recharge. Evapotranspiration as a percentage of precipitation peaks during dry periods, whereas base flow and runoff tend to reach minimum values. During wet periods, this relationship is reversed and base flow and runoff as a percentage of precipitation generally peak while evapotranspiration percentages reach minimum values. Annual recharge in the Appalachian Plateaus reaches a maximum at near 20 percent of annual precipitation, regardless of the severity of wet conditions. Hydrograph separation data from 849 streamflow gaging stations in the study area were used to assess trends in streamflow, base flow, surface runoff, and base-flow index, or ratio of base flow to streamflow, in the Appalachian Plateaus for the period from 1930 to 2011. Annual data anomalies for each of the four variables were individually defined as the annual standard deviation from the mean at all 849 streamflow gaging stations. Annual data anomalies confirm the close relation of annual precipitation to both base flow and runoff components of streamflow, and both components increased during the period of analysis. Around 1970, conditions shifted streamflow from values generally below to above long-term means. At a regional scale, increases in base flow account for most of these observed increases in mean annual streamflow. The independence of the base-flow index to annual climate trends indicate that changes in the components of streamflow of the Appalachian Plateaus are probably in response to shifts in seasonal precipitation or widespread land-use practices. A subset of 77 index streamgages, defined as having 60 or more years of complete record between the years 1930 and 2011 with no more than 20 percent missing data, was selected to show spatial patterns of change in the water budget. Data from the index streamgages showed that the overall trends in base flow are dependent upon the period of evaluation. Long-term (1930–2011) increases in base flow were observed throughout the study area. For two shorter periods (1930–1969 and 1970–2011) trends in base flow were largely negative. In general, spatial patterns of change in streamflow, base flow, and runoff were mixed but generally consistent with prevailing climate patterns and land-use changes.

Alabama, Kentucky, Maryland, Ohio, Pennslyvania, V

Surface-water-quality assessment of the Yakima River basin in Washington: major-and minor-element data for sediment, water, and aquatic biota, 1987-91

Major- and minor-element concentrations are presented for streambed and suspended sediment, filtered- and unfiltered-water, and aquatic-biota samples collected during 1987-91 from the Yakima River Basin in south-central Washington. The samples were collected as part of the U.S. Geological Survey's National Water-quality Assessment (NAWQA) program which is designed to provide results that are useful in understanding and managing the Nation's water resources. This report includes the sampling approach, field collection and processing techniques, and methods of chemical analysis, as well as a compilation of chemical data, statistical summaries, and quality- control data. These data may be used by scientists and resource managers to describe (1) spatial distribution of selected major and minor elements in sediment, water, and aquatic biota of the Yakima River Basin; (2) temporal variation for element concentrations in filtered water and in suspended sediment at selected sites; (3) suita- bility of surface water for preservation of aquatic life and protection of human health; and (4) major natural and anthropogenic sources of major and minor elements in the Yakima River Basin that affect observed water-quality conditions. Streambed-sediment samples were collected once from 27 sites in the basin during 1987-91. Suspended-sediment and filtered-water samples were collected monthly and during hydrologic events (including snowmelt and winter rainstorms) at seven sites, and filtered-water samples were collected at least once at an additional 37 sites during synoptic samplings. Unfiltered-water samples were collected at seven sites on a quarterly basis during 1987 only. Samples of aquatic plants were collected once in 1989, and aquatic insects, fish, and clams were collected from 34 sites three times during 1989-90.

Washington

Using full and partial unmixing algorithms to estimate the inundation extent of small, isolated stock ponds in an arid landscape

Many natural wetlands around the world have disappeared or been replaced, resulting in the dependence of many wildlife species on small, artificial earthen stock ponds. These ponds provide critical wildlife habitat, such that the accurate detection of water and assessment of inundation extent is required. We applied a full (linear spectral mixture analysis; LSMA) and partial (matched filtering; MF) spectral unmixing algorithm to a 2007 Landsat 5 and a 2014 Landsat 8 satellite image to determine the ability of a time-intensive (i.e., more spectral input; LSMA) vs. a more efficient (less spectral input; MF) spectral unmixing approach to detect and estimate surface water area of stock ponds in southern Arizona, USA and northern Sonora, Mexico. Spearman rank correlations ( r s ) between modeled and actual inundation areas less than a single Landsat pixel (< 900 m 2 ) were low for both techniques ( r s range = 0.22 to 0.62), but improved for inundation areas >900 m 2 ( r s range = 0.34 to 0.70). Our results demonstrate that the MF approach can model ranked inundation extent of known pond locations with results comparable to or better than LSMA, but further refinement is required for estimating absolute inundation areas and mapping wetlands <1 Landsat pixel.

Wetlands

Flood-plain delineation for Difficult Run Basin, Fairfax County, Virginia

Water-surface profiles of the 25-year and 100-year floods and maps on which the 25-, 50-, and 100-year flood boundaries are delineated for streams in the Difficult Run basin are presented in this report. The techniques used in the computation of the flood profiles and delineation of flood boundaries are presented, and specific hydraulic problems encountered within the study area are also included.

Virginia

Occurrence and distribution of pesticides and volatile organic compounds in ground water and surface water in Central Arizona Basins, 1996-98, and their relation to land use

Samples of ground water and surface water from the Sierra Vista subbasin, the Upper Santa Cruz Basin, and the West Salt River Valley were collected and analyzed to determine the occurrence and distribution of pesticides and volatile organic compounds in central Arizona. The study was done during 1996-98 within the Central Arizona Basins study unit of the National Water-Quality Assessment program. This study included 121 wells and 4 surface-water sites in the 3 basins and the analyses of samples from 4 sites along the Santa Cruz River that were part of a separate study. Samples were collected from 121 wells and 3 surface-water sites for pesticide analyses, and samples were collected from 109 wells and 3 surface-water sites for volatile organic compound analyses. Certain pesticides detected in ground water and surface water can be related specifically to agricultural or urban uses; others can be related to multiple land uses. Effects from historical agriculture are made evident by detections of DDE in ground-water and surface-water samples collected in the West Salt River Valley and detections of atrazine and deethylatrazine in the ground water in the Upper Santa Cruz Basin. Effects from present agriculture are evident in the seasonal variability in concentrations of pre-emergent pesticides in surface-water samples from the West Salt River Valley. Several detections of DDE and dieldrin in surface water were higher than established water-quality limits. Effects of urban land use are made evident by detections of volatile organic compounds in ground water and surface water from the West Salt River Valley. Detections of volatile organic compounds in surface water from the Santa Cruz River near Nogales, Arizona, also are indications of the effects of urban land use. One detection of tetrachloroethene in ground water was higher than established water-quality limits. Water reuse is an important conservation technique in the Southwest; however, the reuse of water provides a transport mechanism for pesticides and volatile organic compounds to reach areas that are not normally affected by manmade compounds from specific land-use activities. The most complex mixture of pesticides and volatile organic compounds is in the West Salt River Valley and is the result of water-management practices and the combination of land uses in this basin throughout history.

Arizona

Surface-Water, Water-Quality, and Ground-Water Assessment of the Municipio of Mayaguez, Puerto Rico, 1999-2002

The surface-water assessment portion of this study focused on analysis of low-flow characteristics in local streams and rivers, because the supply of safe drinking water was a critical issue during recent dry periods. Low-flow characteristics were evaluated at one continuous-record gaging station based on graphical curve-fitting techniques and log-Pearson Type III frequency curves. Estimates of low-flow characteristics for 20 partial-record stations were generated using graphical-correlation techniques. Flow-duration characteristics for the continuous- and partial-record stations were estimated using the relation curves developed for the low-flow study. Stream low-flow statistics document the general hydrology under current land use, water-use, and climatic conditions. A survey of streams and rivers utilized 37 sampling stations to evaluate the sanitary quality of about 165 miles of stream channels. River and stream samples for fecal coliform and fecal streptococcus analyses were collected on two occasions at base-flow conditions. Bacteriological analyses indicate that a significant portion of the stream reaches within the municipio of Mayaguez may have fecal coliform bacteria concentrations above the water-quality goal (standard) established by the Puerto Rico Environmental Quality Board (Junta de Calidad Ambiental de Puerto Rico) for inland surface waters. Sources of fecal contamination may include: illegal discharge of sewage to storm-water drains, malfunctioning sanitary sewer ejectors, clogged and leaking sewage pipes, septic tank leakage, unfenced livestock, and runoff from livestock pens. Long-term fecal coliform data from five sampling stations located within or in the vicinity of the municipio of Mayaguez have been in compliance with the water-quality goal for fecal coliform concentration established in July 1990. Geologic, topographic, soil, hydrogeologic, and streamflow data were compiled into a database and used to divide the municipio of Mayaguez into five hydrogeologic terranes. This integrated database then was used to evaluate the ground-water potential of each hydrogeologic terrane. Lineament-trace analysis was used to help assess the ground-water development potential in the hydrogeologic terranes containing igneous rocks. Analyses suggest that areas with slopes greater than 15 degrees have relatively low ground-water development potential. The presence of fractures, independent of the topographic slope, may locally enhance the water-bearing properties in the hydrogeologic terranes containing igneous rocks. The results of this study indicate that induced streamflow generally is needed to sustain low to moderate ground-water withdrawal rates in the five hydrogeologic terranes. The ground-water flow systems in the hydrogeologic terranes are only able to sustain small withdrawal rates that rarely exceed 50 gallons per minute. Areas with a high density of fractures, as could be the case at the intersection of lineament traces in the upper parts of the Rio Ca?as and Rio Yaguez watersheds, are worthy of exploratory drilling for ground-water development.

Water-Resources Investigations Report

Effect of land-applied biosolids on surface-water nutrient yields and groundwater quality in Orange County, North Carolina

Land application of municipal wastewater biosolids is the most common method of biosolids management used in North Carolina and the United States. Biosolids have characteristics that may be beneficial to soil and plants. Land application can take advantage of these beneficial qualities, whereas disposal in landfills or incineration poses no beneficial use of the waste. Some independent studies and laboratory analysis, however, have shown that land-applied biosolids can pose a threat to human health and surface-water and groundwater quality. The effect of municipal biosolids applied to agriculture fields is largely unknown in relation to the delivery of nutrients, bacteria, metals, and contaminants of emerging concern to surface-water and groundwater resources. Therefore, the North Carolina Department of Environment and Natural Resources (NCDENR) collaborated with the U.S. Geological Survey (USGS) through the 319 Nonpoint Source Program to better understand the transport of nutrients and bacteria from biosolids application fields to groundwater and surface water and to provide a scientific basis for evaluating the effectiveness of the current regulations. The USGS conducted a paired agricultural watershed study in the Collins Creek and Cane Creek Reservoir watersheds in Orange County, North Carolina. Field activities were conducted from March 2011 through May 2013 at two field study sites, including biosolids field application sites owned by Orange County Water and Sewer Authority (OWASA) in the Collins Creek watershed and a background study site in the Cane Creek watershed that has no fields receiving biosolids applications. Samples of biosolids source material and soil were collected from the land-application fields for laboratory analyses. Soil samples were also collected from a background agricultural field in the Cane Creek watershed that has never received land-applied municipal biosolids. Shallow groundwater samples were collected quarterly from new monitoring wells installed by NCDENR along the edge of the biosolids land-application fields and a background agricultural field for laboratory analyses. Two surface-water monitoring sites were established on Collins Creek to compute continuous streamflow and collect discrete baseflow and stormwater runoff water-quality data upstream and downstream from the biosolids land-application fields. Surface water-quality samples were also collected for baseflow and stormwater runoff conditions at an existing USGS streamgage on Cane Creek to monitor water-quality conditions in the background study watershed. The study primarily focused on nutrients and bacteria; however, data for field properties and water-quality constituents, including metals, major ions, and contaminants of emerging concern (household-, industrial-, and agricultural-use compounds, pharmaceutical compounds, hormones, and antibiotics) also were collected and used in the analyses. There were no exceedances of the 10 elements with designated U.S. Environmental Protection Agency (EPA) ceiling concentrations for land-applied biosolids in any of the biosolids samples. Treatment processes and storage techniques used by OWASA are effective in eliminating Escherichia coli and fecal coliform bacteria from biosolids. Copper, molybdenum, total Kjeldahl nitrogen, and total phosphorus were elevated in the soil from biosolids land-application fields relative to the background field. The relative richness of these constituents in the biosolids land-application fields is consistent with biosolids being the source of the elevated concentrations given the relatively high concentrations of these constituents in the biosolids samples that were collected. Shallow groundwater in the transitional zone wells, which were located adjacent to and topographically downgradient from all the biosolids land-application fields, were found to be statistically different and had higher nitrate concentrations (medians greater than 12 milligrams per liter) than all the other wells sampled as part of the study. Surface-water nutrient concentrations and yields, primarily nitrate, were higher at the monitoring site on Collins Creek downstream from the biosolids land-application fields than the other study sites that drained watersheds without biosolids land application. The largest differences in concentrations between sites were measured at baseflow conditions, which indicate that the main cause of these differences, particularly between Cane Creek and the Collins Creek site downstream from the OWASA application fields, is related to nitrate contribution from the shallow groundwater. Contaminants of emerging concern were detected in approximately 40 percent of the laboratory analyses of the biosolids samples and more frequently in soil samples from the biosolids land-application fields (approximately 40 percent of laboratory analyses) relative to the soil samples from the background field (approximately 12 percent of laboratory analyses). However, contaminants of emerging concern detected in the laboratory analysis for this study do not appear to be good indicators of human-waste contaminants derived from land-applied biosolids in groundwater or surface-water because the number of detections and concentrations at the background wells and surface-water monitoring sites are similar to or higher than those at wells and monitoring sites adjacent to or downstream from the biosolids land-application fields. The data, analysis, and conclusions associated with this study can be used by regulatory agencies, resource managers, and wastewater-treatment operators to (1) better understand the quantity and characteristics of nutrients, bacteria, metals, and contaminants of emerging concern that are transported away from biosolids land-application fields to surface water and groundwater under current regulations for the purposes of establishing effective total maximum daily loads (TMDLs) and restoring impaired water resources, (2) assess how well existing regulations protect waters of the State and potentially recommend effective changes to regulations or land-application procedures, and (3) establish a framework for developing guidance on effective techniques for monitoring and regulatory enforcement of permitted biosolids land-application fields.

North Carolina

Characterization of Interactions between Surface Water and Near-Stream Groundwater along Fish Creek, Teton County, Wyoming, by Using Heat as a Tracer

Fish Creek, a tributary of the Snake River, is about 25 river kilometers long and is located in Teton County in western Wyoming near the town of Wilson. Local residents began observing an increase in the growth of algae and aquatic plants in the stream during the last decade. Due to the known importance of groundwater to surface water in the area, the U.S. Geological Survey (USGS), in cooperation with the Teton Conservation District, conducted a study to characterize the interactions between surface water and near-stream groundwater along Fish Creek. The study has two main objectives: (1) develop an improved spatial and temporal understanding of water flow (fluxes) between surface water and groundwater, and (2) use a two-dimensional groundwater-flow and heat-transport model to interpret observed temperature and hydraulic-head distributions and to describe groundwater flow near Fish Creek. The study is intended to augment hydrologic information derived from previously published results of a seepage investigation on Fish Creek. Seepage measurements provide spatially averaged gains and losses over an entire reach for one point in time, whereas continuous temperature and water-level measurements provide continuous estimates of gain and loss at a specific location. Stage, water-level, and temperature data were collected from surface water and from piezometers completed in an alluvial aquifer at three cross sections on Fish Creek at Teton Village, Resor's Bridge, and Wilson from October 2004 to October 2006. The flow and energy (heat) transport model VS2DH was used to simulate flow through the streambed of Fish Creek at the Teton Village cross section from April 15 to October 14, 2006, (183 recharge periods) and at the Resor's Bridge and Wilson cross sections from June 6, 2005, to October 14, 2006 (496 recharge periods). A trial-and-error technique was used to determine the best match between simulated and measured data. These results were then used to calibrate the cross-sectional models and determine horizontal and vertical hydraulic conductivities. The fluxes of groundwater into the stream or fluxes of stream water into the alluvial aquifer were estimated by using the calibrated VS2DH model for each cross section. Results of the simulations indicated that surface water/groundwater interaction and hydraulic properties were different at the three cross sections. At the most upstream cross section, Teton Village, Fish Creek flowed intermittently and continually gained relatively large quantities of water from April through September. During other times of the year, the stream was dry near the cross section. Saturated hydraulic conductivity set at 1x10-4 m/s in both the horizontal and vertical directions resulted in the best match between simulated and measured temperatures. The Resor's Bridge cross section, about midway between the other two cross sections, was near the point where perennial flow begins. At this cross section, the stream gained water from groundwater during high flow in late spring and summer, was near equilibrium with groundwater during August and September, and lost water to groundwater during the remainder of the year. Horizontal hydraulic conductivity set at 5x10-5 m/s and vertical hydraulic conductivity set at 1x10-5 m/s resulted in the best match between simulated and measured temperatures. The Wilson cross section, the most downstream site, was at USGS streamflow-gaging station 13016450. This part of the stream is perennial and was almost always gaining a small volume of water from groundwater. Saturated hydraulic conductivity set at 1x10-4 m/s in the horizontal direction and at 5x10-6 m/s in the vertical direction resulted in the best match between simulated and measured temperatures. Quantitative values of the flux from groundwater into surface water were estimated by using VS2DH and ranged from 1.1 to 6.6 cubic meters per day (m3/d) at the Teton Village cross section, from -3.8 to 7.4 m3/d at t

Scientific Investigations Report

Assessment of water availability in the Osage Nation using an integrated hydrologic-flow model

The Osage Nation of northeastern Oklahoma, conterminous with Osage County, covers about 2,900 square miles. The area is primarily rural with 62 percent of the land being native prairie grass, and much of the area is used for cattle ranching and extraction of petroleum and natural gas. Protection of water rights are important to the Osage Nation because of its reliance on cattle ranching and the potential for impairment of water quality by petroleum extraction. Additionally, the potential for future population increases, demands for water from neighboring areas such as the Tulsa metropolitan area, and expansion of petroleum and natural-gas extraction on water resources of this area further the need for the Osage Nation to better understand its water availability. Therefore, the U.S. Geological Survey, in cooperation with the Osage Nation, completed a hydrologic investigation to assess the status and availability of surface-water and groundwater resources in the Osage Nation. A transient integrated hydrologic-flow model was constructed using the U.S. Geological Survey fully integrated hydrologic-flow model called the MODFLOW One-Water Hydrologic Model. The integrated hydrologic-flow model, called the Osage Nation Integrated Hydrologic Model (ONIHM), was constructed and uses an orthogonal grid of 276 rows and 289 columns, and each grid cell measures 1,312.34 feet (ft; 400 meters) per side, with eight variably thick vertical layers that represented the alluvial and bedrock aquifers within the study area, including the alluvial aquifer, the Vamoosa-Ada aquifer, and the minor Pennsylvanian bedrock aquifers, and the confining units. Landscape and groundwater-flow processes were simulated for two periods: (1) the 1950–2014 period from January 1950 through September 2014 and (2) the forecast period from October 2014 through December 2099. The 1950–2014 period ONIHM simulated past conditions using measured or estimated inputs, and the forecast-period ONIHM simulated three separate potential forecast conditions under constant dry, average, or wet climate conditions using calibrated input values from the 1950–2014 period ONIHM. The 1950–2014 period ONIHM was calibrated by linking the Parameter Estimation software (PEST) with the MODFLOW One-Water Hydrologic Model. PEST uses statistical parameter estimation techniques to identify the best set of parameter values to minimize the difference between measured or estimated calibration targets and their simulated equivalent values (residuals). Tikhonov regularization and singular-value decomposition-assist features of PEST were used during the calibration process. The 1950–2014 period ONIHM was calibrated to 713 measured groundwater levels at 195 wells; 95,636 estimated monthly mean groundwater levels at 124 wells; 5,307 measured streamflows at 13 streamgages; and 8,679 simulated mean monthly streamflows at 10 streamgages extracted from a surface-water model by adjusting 231 parameters. The estimated groundwater-level observations and streamflows were included as observations to improve the spatial and temporal density of observation targets during calibration. The best set of parameter values obtained during the calibration process of the 1950–2014 model was then used as the input parameter values for the forecast model simulations. A comparison of the calibration targets to their corresponding simulated values indicated that the model adequately reproduced streamflows and groundwater levels for some streamgages and wells and underestimated streamflows and groundwater levels at other locations. Measured and simulated streamflows correlated adequately with a coefficient of determination of 0.938, as did water levels with a coefficient of determination of 0.795. The 1950–2014 period ONIHM underestimated certain groundwater levels and streamflows, but generally measured or estimated calibration targets correlated well with simulated equivalents, which indicated that the model can adequately simulate the response of the hydrologic system to stresses in the 1950–2014 and forecast periods. In the 1950–2014 period ONIHM, the calibrated mean horizontal hydraulic conductivity for layer 1 alluvial aquifer was 30.7 feet per day, and the seven lower layers had a calibrated mean horizontal hydraulic conductivity of less than 3.3 feet per day. The mean calibrated groundwater-level residual was 16.6 ft, and the mean calibrated streamflow residual of the Arkansas River at Ralston, Oklahoma, streamgage (U.S. Geological Survey station 07152500) was within 6 percent (373 cubic feet per second) of mean measured streamflow for the 1950–2014 period ONIHM. The ONIHM simulated landscape fluxes of precipitation; groundwater applied by irrigation wells; evapotranspiration from precipitation, groundwater, and irrigation; runoff from precipitation; and deep percolation from precipitation. The largest loss of water from the landscape was evapotranspiration from precipitation with a calibrated mean annual outflow of 32 inches (in.): mean annual precipitation was about 36 in. Calibrated mean annual runoff and deep percolation (recharge to the water table) rates were 4.7 inches per year (in/yr) and 0.70 in/yr, respectively, for the 1950–2014 period ONIHM. The calibrated 1950–2014 period ONIHM groundwater fluxes included net farm net recharge (calculated as the difference between the inflow of recharge to the water table and the outflow of evapotranspiration from the water table such that negative values indicate that evapotranspiration from the water table was greater than deep percolation [recharge to the water table] and vice versa). Net farm net recharge was the largest flux from the groundwater system with a mean annual net outflow of 153.4 cubic feet per second. Stream leakage was the largest flux to the groundwater system with a mean annual net inflow of 152.5 cubic feet per second, indicating that, on average, the groundwater/surface-water interaction was a “losing” system where stream water leaked into the subsurface and recharged the water table. Simulated monthly trends demonstrated that net stream leakage was the largest inflow to the groundwater-flow system for 10 of the 12 months; for the other 2 months (January and March), farm net recharge (January) and net storage (March) were the largest inflow to the groundwater-flow system. A saline groundwater interface map was created for the study and compared to the water levels from the final stress period of the 1950–2014 model to identify the presence of fresh/marginal groundwater throughout the study area. Fresh/marginal groundwater was characterized as groundwater with less than 1,500 milligrams per liter of total dissolved solids. Fresh/marginal groundwater thickness ranged from 0 to 438.2 ft within the study area. The thickest regions of fresh/marginal groundwater were in the eastern part of the study area near Sand Creek, Bird Creek, and Hominy Creek and in the Arkansas River alluvial aquifer in the region downstream from the Arkansas River at Ralston, Okla. Like the 1950–2014 model, forecast model results for the landscape indicated that transpiration from precipitation was the largest flux out of the landscape for all three forecasts, constituting 77, 73, and 58 percent of precipitation for the dry, average, and wet forecasts, respectively. The dry and average forecast landscape fluxes demonstrated similar trends and magnitudes, whereas the wet forecast landscape fluxes indicated the largest changes compared to the average forecast fluxes. Most notably, runoff increased from a mean of 1.1 and 1.6 in/yr for the dry and average forecasts, respectively, to 10 in/yr for the wet forecast. Similar changes occurred for the other wet forecast landscape fluxes. The calibrated 1950–2014 period ONIHM simulated three forecasts to assess the effects of potential climatic changes on the hydrologic system from October 2014 to December 2099. The three forecasts simulated theoretical dry, average, and wet conditions using precipitation and potential evapotranspiration datasets from selected years in the calibrated 1950–2014 period ONIHM. Annual precipitation amounts were 26.89, 35.47, and 50.73 in. for the dry, average, and wet forecasts, respectively. Groundwater-flow component forecast results indicated that stream leakage is always a net inflow to the groundwater-flow system for dry, average, and wet conditions, meaning the study area stream network is always predominantly a “losing” regime where stream water infiltrates into the underlying aquifer. Storage was only a net outflow from the groundwater-flow system and indicated a replenishment to groundwater storage that resulted in an increase in groundwater levels only during the wet forecast. Further, these gains in groundwater storage for the wet forecast occurred only during February through June. Mean fresh/marginal groundwater saturated thicknesses were 125 and 126 ft for the dry and average forecast conditions, respectively, and wet forecast average thickness was 145 ft and ranged from 0 to 443 ft. The spatial extents of fresh/marginal groundwater at the end of the dry, average, and wet forecast model periods (December 2099) did not change substantially from the end of the 1950–2014 model period (September 2014).

Kansas, Oklahoma

User's guide to SAC, a computer program for computing discharge by slope-area method

This user's guide contains information on using the slope-area program, SAC. SAC can be used to compute peak flood discharges from measurements of high-water marks along a stream reach. The Slope-area method used by the program is the U.S. Geological Survey (USGS) procedure presented in Techniques of Water Resources Investigations of the U.S. Geological Survey, beok 3, chapter A2, "Measurement of Peak Discharge by the Slope-Area Method." The program uses input files that have formats compatible with those used by the water-surface profile program (WSPRO) described in the Federal Highways Administration publication FHWA-IP-89-027. The guide briefly describes the slope-area method documents the input requirements and the output produced, and demonstrates use of SAC.

Open-File Report

User's guide to the culvert analysis program

This user's guide contains information on using the culvert analysis program (CAP). The procedure used is based on that presented in Techniques of Water- Resources Investigations of the United States Geological Survey, book 3, chapter A3, "Measurement of Peak Discharge at Culverts by Indirect Methods." The program uses input files that have formats compatible with those used by the Water-Surface Profile (WSPRO) program. The program can be used to compute rating surfaces or curves that describe the behavior of flow through a culvert or to compute discharges from measurements of upstream and downstream water-surface elevations.

Open-File Report

Application of borehole geophysics in defining the wellhead protection area for a fractured crystalline bedrock aquifer

Wellbore geophysical techniques were used to characterize fractures and flow in a bedrock aquifer at a site near Blackwater Brook in Dover, New Hampshire. The primary focus ofthis study was the development of a model to assist in evaluating the area surrounding a planned water supply well where contaminants introduced at the land surface might be induced to flow towards a pumping well. Well logs and geophysical surveys used in this study included lithologic logs based on examination of cuttings obtained during drilling; conventional caliper and natural gamma logs; video camera and acoustic televiewer surveys; highresolution vertical flow measurements under ambient conditions and during pumping; and borehole fluid conductivity logs obtained after the borehole fluid was replaced with deionized water. These surveys were used for several applications: 1) to define a conceptual model of aquifer structure to be used in groundwater exploration; 2) to estimate optimum locations for test and observation wells; and 3) to delineate a wellhead protection area (WHPA) for a planned water supply well. Integration of borehole data with surface geophysical and geological mapping data indicated that the study site lies along a northeast-trending intensely fractured contact zone between surface exposures of quartz monzonite and metasedimentary rocks. Four of five bedrock boreholes at the site were estimated to produce more than 150 gallons per minute (gpm) (568 L/min) of water during drilling. Aquifer testing and other investigations indicated that water flowed to the test well along fractures parallel to the northeast-trending contact zone and along other northeast and north-northwest-trending fractures. Statistical plots of fracture strikes showed frequency maxima in the same northeast and north-northwest directions, although additional maxima occurred in other directions. Flowmeter surveys and borehole fluid conductivity logging after fluid replacement were used to identify water-producing zones in the boreholes; fractures associated with inflow into boreholes showed a dominant northeast orientation. Borehole fluid conductivity logging after fluid replacement also gave profiles of such water-quality parameters as fluid electrical conductivity (FEC), pH, temperature, and oxidation-reduction potential, strengthening the interpretation of crossconnection of boreholes by certain fracture zones. The results of this study showed that the application of these borehole geophysical techniques at the Blackwater Brook site led to an improved understanding of such parameters as fracture location, attitude, flow direction and velocity, and water quality; all of which are important in the determination of a WHPA.

Log Analyst