USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Monthly”

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 631 records · Page 35Linked to original sources

Groundwater conditions and trends, 2009–19, Saipan, Commonwealth of the Northern Mariana Islands

Groundwater is the primary source of municipal water for Saipan. Nearly all groundwater for the municipal water supply is withdrawn from a freshwater-lens system with a limited amount of freshwater that is susceptible to saltwater intrusion. The status of Saipan’s groundwater resources has not been thoroughly assessed since 2003. The U.S. Geological Survey—in cooperation with the Office of Grants Management, Commonwealth of the Northern Mariana Islands, and in collaboration with the Commonwealth Utilities Corporation—assessed the status and characteristics of Saipan’s groundwater resources by (1) evaluating groundwater withdrawals from municipal production wells during 2014–19, (2) evaluating chloride concentrations of municipal groundwater withdrawals during 2009–19, and (3) collecting salinity profiles at selected groundwater-monitoring wells during 2018–19. At the time of preparation of this report (2019), the periods of groundwater-withdrawal and chloride-concentration data represent the only periods of data available since 2003. During 2014–19, groundwater for the municipal water supply was withdrawn from about 143 production wells. Most of the wells are drilled into limestone formations in the southern plateau and the Kagman Peninsula and generally have withdrawal rates of about 40–60 gallons per minute. Records of monthly groundwater withdrawals from municipal production wells were available for May 2014–March 2019; during that period, monthly withdrawals ranged from 5.7 to 12.8 million gallons per day (Mgal/d) and averaged 9.3 Mgal/d, although records were unavailable for 9 months (May 2015–January 2016). Private wells, mainly located on the western coastal plain, currently are permitted to withdraw a total of about 7 Mgal/d of groundwater. Actual groundwater withdrawals from private wells, however, are uncertain because withdrawal records for private wells are not available. The Commonwealth Utilities Corporation measured the chloride concentration of groundwater pumped from each of its production wells about twice a year from 2009–19; during this period, 146 production wells were active and sampled. Only 32 of the 146 (22 percent) municipal production wells had median chloride concentrations less than or equal to 250 milligrams per liter (mg/L), the secondary drinking water standard set by the U.S. Environmental Protection Agency. Eighty-one wells (55 percent) pumped water with median chloride concentrations above 500 mg/L. The Mann-Kendall test was used to determine if chloride concentrations of groundwater withdrawals at 146 municipal production wells had statistically significant trends during December 2009–February 2019. Trends were considered statistically significant for probability values (p-values) less than or equal to 0.05. Test results indicate an upward trend at 9 wells, a downward trend at 52 wells, and no trend at 85 wells. Salinity profiles were measured in 12 selected monitor wells during July–August 2018 and were measured in six of the twelve selected monitor wells during March 2019. The salinity profiles were used to estimate the thickness of the freshwater lens at 10 monitor wells; freshwater-lens thickness was greatest (46 ft) in a monitor well in the Dan Dan well field near the northern part of the southern plateau. Freshwater-lens-thickness estimates elsewhere were (1) between 0 and 28 ft for the remaining monitor wells on the southern plateau, (2) between 19 and 21 ft for monitor wells on the Kagman Peninsula, (3) 2 ft for a monitor well in the Sablan Quarry well field on west-central Saipan, and (4) 8 ft for a monitor well in the Marpi Quarry well field on northern Saipan.

Commonwealth of the Norhtern Marianas Islands, Sai↗

Ground-Water Flow Model for the Spokane Valley-Rathdrum Prairie Aquifer, Spokane County, Washington, and Bonner and Kootenai Counties, Idaho

This report presents a computer model of ground-water flow in the Spokane Valley-Rathdrum Prairie (SVRP) aquifer in Spokane County, Washington, and Bonner and Kootenai Counties, Idaho. The aquifer is the sole source of drinking water for more than 500,000 residents in the area. In response to the concerns about the impacts of increased ground-water withdrawals resulting from recent and projected urban growth, a comprehensive study was initiated by the Idaho Department of Water Resources, the Washington Department of Ecology, and the U.S. Geological Survey to improve the understanding of ground-water flow in the aquifer and of the interaction between ground water and surface water. The ground-water flow model presented in this report is one component of this comprehensive study. The primary purpose of the model is to serve as a tool for analyzing aquifer inflows and outflows, simulating the effects of future changes in ground-water withdrawals from the aquifer, and evaluating aquifer management strategies. The scale of the model and the level of detail are intended for analysis of aquifer-wide water-supply issues. The SVRP aquifer model was developed by the Modeling Team formed within the comprehensive study. The Modeling Team consisted of staff and personnel working under contract with the Idaho Department of Water Resources, personnel working under contract with the Washington Department of Ecology, and staff of the U.S. Geological Survey. To arrive at a final model that has the endorsement of all team members, decisions on modeling approach, methodology, assumptions, and interpretations were reached by consensus. The ground-water flow model MODFLOW-2000 was used to simulate ground-water flow in the SVPR aquifer. The finite-difference model grid consists of 172 rows, 256 columns, and 3 layers. Ground-water flow was simulated from September 1990 through September 2005 using 181 stress periods of 1 month each. The areal extent of the model encompasses an area of approximately 326 square miles. For the most part, the model extent coincides with the 2005 revised extent of the Spokane Valley-Rathdrum Prairie aquifer as defined in a previous report. However, the model excludes Spirit and Hoodoo Valleys because of uncertainties about the ground-water flow directions in those valleys and the degree of hydraulic connection between the valleys and northern Rathdrum Prairie. The SVRP aquifer is considered to be a single hydrogeologic unit except in Hillyard Trough and the Little Spokane River Arm. In those areas, a continuous clay layer divides the aquifer into an upper, unconfined unit and a lower, confined unit. The model includes all known components of inflows to and outflows from the aquifer. Inflows to the SVRP aquifer include (1) recharge from precipitation, (2) inflows from tributary basins and adjacent uplands, (3) subsurface seepage and surface overflows from lakes that border the aquifer, (4) flow from losing segments of the Spokane River to the aquifer, (5) return percolation from irrigation, and (6) effluent from septic systems. Outflows from the SVRP aquifer include (1) ground-water withdrawals from wells, (2) flow from the aquifer to gaining segments of the Spokane River, (3) aquifer discharge to the Little Spokane River, and (4) subsurface outflow from the lower unit at the western limit of the model area near Long Lake. These inflow and outflow components are represented in the model by using MODFLOW-2000 packages. The parameter-estimation program PEST was used to calibrate the SVRP aquifer model. PEST implements a nonlinear least-squares regression method to estimate model parameters so that the differences between measured and simulated quantities are minimized with respect to an optimal criterion. Calibration data include 1,573 measurements of water levels and 313 measurements of streamflow gains and losses along segments of the Spokane and Little Spokane Rivers. Model parameters estimated during calib

Scientific Investigations Report↗

Nutrient Loading and Algal Response in West Thompson Lake, Thompson, Connecticut, 2003-2005

Water quality and nutrient loads were characterized for parts of the Quinebaug River and West Thompson Lake in northeastern Connecticut during 2003 to 2005. The West Thompson Lake watershed is a mainly forested watershed that receives treated municipal wastewater from several point sources in Massachusetts. The lake is a flood-control reservoir formed in 1966 by impoundment of the Quinebaug River. Median concentrations of total phosphorus in two inflow (upstream) and one outflow (downstream) sampling stations on the Quinebaug River were higher than the nutrient criteria recommended by the U.S. Environmental Protection Agency (USEPA) for rivers and streams in aggregate Ecoregion XIV. In general, concentrations of total phosphorus in West Thompson Lake also were above the nutrient criteria recommended by USEPA for lakes and impoundments in aggregate Ecoregion XIV. The trophic status of West Thompson Lake has changed since 1995 from a hypereutrophic lake to a eutrophic lake; however, the lake still has large algal blooms. These blooms are predominated by blue-green algae, with chlorophyll-a concentrations of more than 30 micrograms per liter and algal cell counts as high as 73,000 cells/mL. Water samples collected during the summer of 2005 identified phosphorus as the primary limiting nutrient early in the season, but algal growth is probably co-limited by phosphorus and nitrogen later in the season. Lake-bottom sediments were collected from several areas throughout the lake and ranged in thickness from less than 1 foot (ft) to more than 3 ft. Concentrations of phosphorus in sediments differed throughout the lake; the highest values were found in the middle of the lake. Concentrations of total phosphorus also increased from an average 1,800 milligrams per kilogram (mg/kg) in the upper layers of sediment to more than 6,000 mg/kg at depth in the sediment. Annual, seasonal, and monthly loads and yields of nutrients were calculated for the three sampling locations on the Quinebaug River to develop a nutrient mass-balance model (budget) for West Thompson Lake. The average annual yields of total phosphorus during 2000 to 2005 were 115 pounds per square mile per year (lb/mi2/yr) at Quinebaug (inflow station), 116 lb/mi2/yr at Red Bridge Road (inflow station), and 97.9 lb/mi2/yr at West Thompson (outflow station). The 18-percent decrease in the average annual yield of total phosphorus between the inflow station at Red Bridge Road and the outlet of West Thompson Lake at West Thompson indicates that a significant part of the phosphorus load is retained in the lake. Annual yields of total phosphorus at Quinebaug have decreased significantly since the 1980s, from 362 lb/mi2/yr (for 1981-1990) to 115 lb/mi2/yr (1996-2005). The annual net export of phosphorus in West Thompson Lake during water years 2000 to 2005 ranged from -36 percent (2005) to 1 percent (2002) of the incoming load. Seasonal mass-balance data for total phosphorus during the summers of 2000 to 2003, when streamflow was at or lower than normal, indicated a net export of phosphorus that ranged from 3.4 percent (2003) to 30.7 percent (2002) of the incoming load. During the summer of 2004, however, streamflows were much higher than normal, and there was a negative export of phosphorus in West Thompson Lake of -3.9 percent. The annual net export of nitrogen in West Thompson Lake during water years 2000 to 2005 ranged from -5 percent (2002) to 4 percent (2001) of the incoming load. No clear pattern was evident to relate total nitrogen export to seasonal variables or runoff. Removal of phosphorus during the summer by wastewater-treatment plants (WWTPs) in Massachusetts reduces the concentration and load of total phosphorus entering West Thompson Lake in the summer; however, the large amount of phosphorus retained in the lake during the other seasons, in addition to the phosphorus stored in the lake-bottom sediments, may become available to fuel algal blooms in the lake

Scientific Investigations Report↗

Water-quality trends for selected sites in the Boulder River and Tenmile Creek watersheds, Montana, based on data collected during water years 1997-2013

In the Boulder River and Tenmile Creek watersheds in southwestern Montana, there was intensive mining during a 40-year period after the discovery of gold in the early 1860s. Potential effects from the historic mining activities include acid-mine drainage and elevated concentrations of potentially toxic trace elements from mining remnants such as waste rock and tailing piles. In support of remediation efforts, water-quality monitoring by the U.S. Geological Survey began in 1997 in the Boulder River and Tenmile Creek watersheds and has continued to present (2014). The U.S. Geological Survey, in cooperation with the U.S. Forest Service, investigated temporal trends in water quality at 13 sites, including 2 adit (or mine entrance) sites and 11 stream sites. The primary purpose of this report is to present results of trend analysis of specific conductance, selected trace-elements (cadmium, copper, lead, zinc, and arsenic), and suspended sediment for the 13 sites. Trend results for most stream sites in the Boulder River watershed for water years 2000–13 (water year is the 12-month period from October 1 through September 30 and is designated by the year in which it ends) indicate decreasing trends in flow-adjusted specific conductance, in flow-adjusted concentrations (FACs) for most filtered and unfiltered-recoverable trace elements, and in suspended sediment. Overall, magnitudes of the decreasing trends in FACs of metallic contaminants are largest for Bullion Mine tributary at mouth (site 3), Jack Creek at mouth (site 4), and Cataract Creek at Basin (site 8). For sites 3, 4, and 8, magnitudes of decreasing trends generally ranged from about -5 to -10 percent per year. Notably, the watersheds upstream from sites 3, 4, and 8 have been targeted by substantial remediation activities. Consideration of trend patterns among all stream sites in the Boulder River watershed provides strong evidence that remediation activities are the primary cause of decreasing trends in metallic contaminants. Trend results for sites in the Tenmile Creek watershed generally are more variable and difficult to interpret than for sites in the Boulder River watershed. Trend results for Tenmile Creek above City Diversion (site 11) and Minnehaha Creek near Rimini (site 12) for water years 2000–13 indicate decreasing trends in FACs of cadmium, copper, and zinc. The magnitudes of the decreasing trends in FACs of copper generally are moderate and statistically significant for sites 11 and 12. The magnitudes of the decreasing trends in FACs of cadmium and zinc for site 11 are minor to small and not statistically significant; however, the magnitudes for site 12 are moderate and statistically significant. In general, patterns in FACs for Tenmile Creek near Rimini (site 13) are not well represented by fitted trends within the short data collection period, which might indicate that the trend-analysis structure of the study is not appropriate for describing trends in FACs for site 13. The large decreasing trend in FACs of suspended sediment is the strongest indication of change in water quality during the short period of record for site 13; however, this trend is not statistically significant.

Montana↗

Birds and Aircraft on Midway Islands, 1959-63 Investigations

At Midway Naval Station, 1.100 miles west-northwest of Honolulu, military aircraft collide with flying albatrosses at the rate of about 300 to 400 per year. One aircraft out of every five that hits an albatross on takeoff either aborts (stops before it is airborne), or dumps fuel and returns for appraisal of damage. About 70,000 pairs of Laysan albatrosses and 7,000 pairs of blackfooted albatrosses nest at Midway in any given year. The population is declining. Two-thirds or more of the birds of breeding age nest each year. The minimum breeding age recorded is 5 years (each species), but many individuals do not nest until at least 7 years of age. Young birds begin to return to Midway at 3 years of age and are found more frequently as breeding age approaches. They come ashore more frequently in March and April (the high bird strike months) than in midwinter. Even in midwinter the number of 'walkers' (birds not on nests) may comprise more than 40 percent of the albatrosses present on Sand Island, Midway. Maximum longevity of the Laysan albatross is believed to exceed 40 years; 6 out of 99 birds banded as breeding adults (7+ years old) were still alive 24 years after banding. Control methods tested experimentally include disturbance, gunfire, other sounds, radar beams, smoke, odors, destruction of nests, eggs, chicks, and adults, moving of birds, eggs, and chicks, erection of obstacles to flight, and habitat management. Habitat management (leveling and hardsurfacing of shoulders of runways) has been the most effective. Albatrosses were counted over the runways at 10 locations in 1957, 1958, and 1960 to determine the effects of wind direction, wind speed, and topography on the numbers of flying birds. Birds were most concentrated in areas where rising air currents were created as winds blew against dunes or tall trees. Soaring and strike rate both increased with greater wind speeds. There was a highly significant correlation between strike frequency and wind direction. Soaring and strike rate were greatest over runway 6-24 when winds were from a northerly quadrant, and greatest over runway 15-33 when winds were east to southeast. Counts (nearly 7,000 observations) before and after the leveling of dunes along part of the south side of the principal runway, 6-24, showed a 61 percent decrease in soaring over the principal runway as compared with a 4 percent increase over runway 15-33 where no major habitat management had been done. Studies of population dynamics of Laysan and black-footed albatrosses in 4 study plots (totaling 14 1/2 acres) are yielding information on nesting density, changes in nesting population from year to year, frequency of nesting of individual birds, closeness of return, reproductive success, rate of nest loss, age at which young birds begin nesting, age composition of the population, and life expectancy of adults. Recommendations are made for reducing the bird hazard and for protecting birds that are not involved in the hazard.

Special Scientific Report - Wildlife↗

Water Quality Conditions in Upper Klamath and Agency Lakes, Oregon, 2005

During June-October 2005, water quality data were collected from Upper Klamath and Agency Lakes in Oregon, and meteorological data were collected around and within Upper Klamath Lake. Data recorded at two continuous water quality monitors in Agency Lake showed similar temperature patterns throughout the field season, but data recorded at the northern site showed more day-to-day variability for dissolved oxygen concentration and saturation after late June and more day-to-day variability for pH and specific conductance values after mid-July. Data recorded from the northern and southern parts of Agency Lake showed more comparable day-to-day variability in dissolved oxygen concentrations and pH from September through the end of the monitoring period. For Upper Klamath Lake, seasonal (late July through early August) lows of dissolved oxygen concentrations and saturation were coincident with a seasonal low of pH values and seasonal highs of ammonia and orthophosphate concentrations, specific conductance values, and water temperatures. Patterns in these parameters, excluding water temperature, were associated with bloom dynamics of the cyanobacterium (blue-green alga) Aphanizomenon flos-aquae in Upper Klamath Lake. In Upper Klamath Lake, water temperature in excess of 28 degrees Celsius (a high stress threshold for Upper Klamath Lake suckers) was recorded only once at one site during the field season. Large areas of Upper Klamath Lake had periods of dissolved oxygen concentration of less than 4 milligrams per liter and pH value greater than 9.7, but these conditions were not persistent throughout days at most sites. Dissolved oxygen concentrations in Upper Klamath Lake on time scales of days and months appeared to be influenced, in part, by bathymetry and prevailing current flow patterns. Diel patterns of water column stratification were evident, even at the deepest sites. This diel pattern of stratification was attributable to diel wind speed patterns and the shallow nature of most of Upper Klamath Lake. Timing of the daily extreme values of dissolved oxygen concentration, pH, and water temperature was less distinct with increased water column depth. Chlorophyll a concentrations varied spatially and temporally throughout Upper Klamath Lake. Location greatly affected algal concentrations, in turn affecting nutrient and dissolved oxygen concentrations - some of the highest chlorophyll a concentrations were associated with the lowest dissolved oxygen concentrations and the highest un-ionized ammonia concentrations. The occurrence of the low dissolved oxygen and high un-ionized ammonia concentrations coincided with a decline in algae resulting from cell death, as measured by concentrations of chlorophyll a. Dissolved oxygen production rates in experiments were as high as 1.47 milligrams of oxygen per liter per hour, and consumption rates were as much as -0.73 milligrams of oxygen per liter per hour. Dissolved oxygen consumption rates measured in this study were comparable to those measured in a 2002 Upper Klamath Lake study, and a higher rate of dissolved oxygen consumption was recorded in dark bottles positioned higher in the water column. Data, though inconclusive, indicated that a decreasing trend of dissolved oxygen productivity through July could have contributed to the decreasing dissolved oxygen concentrations and percent saturation recorded in Upper Klamath Lake during this time. Phytoplankton self-shading was evident from a general inverse relation between depth of photic zone and chlorophyll a concentrations. This shading caused net dissolved oxygen consumption during daylight hours in lower parts of the water column that would otherwise have been in the photic zone. Meteorological data collected in and around Upper Klamath Lake showed that winds were likely to come from a broad range of westerly directions in the northern one-third of the lake, but tended to come from a narrow range of northwesterly directions

Scientific Investigations Report↗

Evaluation of land subsidence and ground failures at Bicycle Basin, Fort Irwin National Training Center, California, 1992–2017

Groundwater has been pumped in the Bicycle Basin at Fort Irwin National Training Center since the 1960s, and the amount pumped has generally increased since the 1990s. After a large crack (approximately 0.5-kilometer long) formed at the surface of Bicycle Lake playa during 2005–06 in the area used as an aircraft runway, a monitoring study was initiated by the U.S. Geological Survey, in cooperation with the U.S. Army Fort Irwin National Training Center, to help determine the cause. The extent and effect of groundwaterlevel declines and land-surface deformation in Bicycle Basin were evaluated using a number of approaches, including water-level measurements from December 2007 to June 2017, land surveys across the playa area, interferometric synthetic aperture radar (InSAR) analyses, geophysical surveys of the playa area, and numerical experiments to test hypotheses about soil mechanical processes. A specific objective of this study was to evaluate the recent development of ground failures in the form of large, surface cracks that pose a hazard to aircraft operations on the Bicycle Lake playa airstrip. Another large crack and set of interconnected cracks formed at the surface of Bicycle Lake playa in 2013 following a period of inundation similar to that of 2005–06; however, the 2013 cracks formed a network rather than a single feature. Groundwater pumping resulted in more than 27 meters (90 feet) of water-level declines in wells north of the Bicycle Lake playa from 1990 to 2017, and InSAR interferograms indicated more than 400 millimeters (16 inches) of subsidence during 1993–2015. Subsidence rates calculated from InSAR interferograms were variable, temporally and spatially. Results of leveling surveys indicated differential subsidence between 2009 and 2016 along a transect across the playa to the area of maximum subsidence; there was less subsidence south of the 2005–06 crack than north of it. The steepest subsidence gradient for this time was in the area of maximum subsidence. Repeat tape-extensometer measurements from April 2009 to November 2017 across the main-runway crack indicated slight opening along part of the crack. A baseline lidar (light detection and ranging) survey of the main-runway crack was done in January 2009 to track the development of this feature. During the 95-months from January 2009 to December 2016, 0.5 meters (19.7 inches) of subsidence was observed next to the subsidence pit in the crack along the western end of the scanned area, presumably due to erosion and backfilling when the lake was inundated. Standing water in the crack along the eastern end of the scanned area prevented determining change for that part of the crack. Time-series, shaded-relief images show the progression of change in the crack from January 2009 through December 2016, with the crack closing and opening, depending on the time of most recent inundation relative to when the scans were taken. Volumetric changes in the size of the crack, calculated between the sequential lidar surveys, show that from 2009 to 2016, the cumulative volumetric change was 1.5 cubic meters (54 cubic feet). This volume is a minimum estimate because the crack volume below the water surface in the pits was not included in the calculations. Subsurface imaging of the 2005–06 main-runway crack by a series of electrical resistivity tomography surveys in 2008 and 2017 indicated that the crack could extend 5 meters (16 feet) in the subsurface. Electromagnetic induction surveys in 2008, 2014, and 2015 evaluated the technique for ground-failure monitoring and showed that the technique was effective at revealing anomalies correlated with the features of concern in the playa area, such as the 2005–06 crack, numerous ”healed” macropolygon features, and scattered sink-like depressions. Results from numerical experiments simulating watertable decline at depth indicated that the material deposits in the desaturating capillary fringe zone might not be able to transmit large enough stresses up through the overlying soil to cause cracks at the land surface. Results from simulations of desiccation in the presence of a regional tectonically induced stress field, however, tended to support the hypothesis that the combined processes could control the formation of giant desiccation macropolygons in certain areas, such as Bicycle Lake playa.

California↗

Simulation of streamflow, evapotranspiration, and groundwater recharge in the Lower Frio River watershed, south Texas, 1961-2008

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, Fort Worth District; the City of Corpus Christi; the Guadalupe-Blanco River Authority; the San Antonio River Authority; and the San Antonio Water System, configured, calibrated, and tested a watershed model for a study area consisting of about 5,490 mi 2 of the Frio River watershed in south Texas. The purpose of the model is to contribute to the understanding of watershed processes and hydrologic conditions in the lower Frio River watershed. The model simulates streamflow, evapotranspiration (ET), and groundwater recharge by using a numerical representation of physical characteristics of the landscape, and meteorological and streamflow data. Additional time-series inputs to the model include wastewater-treatment-plant discharges, surface-water withdrawals, and estimated groundwater inflow from Leona Springs. Model simulations of streamflow, ET, and groundwater recharge were done for various periods of record depending upon available measured data for input and comparison, starting as early as 1961. Because of the large size of the study area, the lower Frio River watershed was divided into 12 subwatersheds; separate Hydrological Simulation Program-FORTRAN models were developed for each subwatershed. Simulation of the overall study area involved running simulations in downstream order. Output from the model was summarized by subwatershed, point locations, reservoir reaches, and the Carrizo-Wilcox aquifer outcrop. Four long-term U.S. Geological Survey streamflow-gaging stations and two short-term streamflow-gaging stations were used for streamflow model calibration and testing with data from 1991-2008. Calibration was based on data from 2000-08, and testing was based on data from 1991-99. Choke Canyon Reservoir stage data from 1992-2008 and monthly evaporation estimates from 1999-2008 also were used for model calibration. Additionally, 2006-08 ET data from a U.S. Geological Survey meteorological station in Medina County were used for calibration. Streamflow and ET calibration were considered good or very good. For the 2000-08 calibration period, total simulated flow volume and the flow volume of the highest 10 percent of simulated daily flows were calibrated to within about 10 percent of measured volumes at six U.S. Geological Survey streamflow-gaging stations. The flow volume of the lowest 50 percent of daily flows was not simulated as accurately but represented a small percent of the total flow volume. The model-fit efficiency for the weekly mean streamflow during the calibration periods ranged from 0.60 to 0.91, and the root mean square error ranged from 16 to 271 percent of the mean flow rate. The simulated total flow volumes during the testing periods at the long-term gaging stations exceeded the measured total flow volumes by approximately 22 to 50 percent at three stations and were within 7 percent of the measured total flow volumes at one station. For the longer 1961-2008 simulation period at the long-term stations, simulated total flow volumes were within about 3 to 18 percent of measured total flow volumes. The calibrations made by using Choke Canyon reservoir volume for 1992-2008, reservoir evaporation for 1999-2008, and ET in Medina County for 2006-08, are considered very good. Model limitations include possible errors related to model conceptualization and parameter variability, lack of data to better 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. A sensitivity analysis was performed for the Upper San Miguel subwatershed model to show the effect of changes to model parameters on the estimated mean recharge, ET, and surface runoff from that part of the Carrizo-Wilcox aquifer outcrop. Simulated recharge was most sensitive to the changes in the lower-zone ET (LZ

Texas↗

Characterization of the hydrologic resources of San Miguel County, New Mexico, and identification of hydrologic data gaps, 2011

The U.S. Geological Survey (USGS), in cooperation with San Miguel County, New Mexico, conducted a study to assess publicly available information regarding the hydrologic resources of San Miguel County and to identify data gaps in that information and hydrologic information that could aid in the management of available water resources. The USGS operates four continuous annual streamgages in San Miguel County. Monthly discharge at these streamgages is generally bimodally distributed, with most runoff corresponding to spring runoff and to summer monsoonal rains. Data compiled since 1951 on the geology and groundwater resources of San Miguel County are generally consistent with the original characterization of depth and availability of groundwater resources and of source aquifers. Subsequent exploratory drilling identified deep available groundwater in some locations. Most current (2011) development of groundwater resources is in western San Miguel County, particularly in the vicinity of El Creston hogback, the hogback ridge just west of Las Vegas, where USGS groundwater-monitoring wells indicate that groundwater levels are declining. Regarding future studies to address identified data gaps, the ability to evaluate and quantify surface-water resources, both as runoff and as potential groundwater recharge, could be enhanced by expanding the network of streamgages and groundwater-monitoring wells throughout the county. A series of seepage surveys along the lengths of the rivers could help to determine locations of surface-water losses to and gains from the local groundwater system and could help to quantify the component of streamflow attributable to irrigation return flow; associated synoptic water-quality sampling could help to identify potential effects to water quality attributable to irrigation return flow. Effects of groundwater withdrawals on streamflow could be assessed by constructing monitoring wells along transects between production wells and stream reaches of interest to monitor decline or recovery of the water table, to quantify the timing and extent of water-table response, and to identify the spatial extent of capture zones. Assessment of groundwater potential could be aided by a county-wide distribution of water-level information and by a series of maps of groundwater potential, compiled for each individual aquifer, including saline aquifers, for which the potential for municipal use through desalination could be explored. A county-wide geographic information system hydrologic geodatabase could provide a comprehensive picture of water use in San Miguel County and could be used by San Miguel County as a decision-support tool for future management decisions.

New Mexico↗

Numerical simulation of the groundwater-flow system in Chimacum Creek Basin and vicinity, Jefferson County, Washington

A groundwater-flow model was developed to evaluate potential future effects of growth and of water-management strategies on water resources in the Chimacum Creek Basin. The model covers an area of about 64 square miles (mi 2 ) on the Olympic Peninsula in northeastern Jefferson County, Washington. The Chimacum Creek Basin drains an area of about 53 mi 2 and consists of Chimacum Creek and its tributary East Fork Chimacum Creek, which converge near the town of Chimacum and discharge to Port Townsend Bay near the town of Irondale. The topography of the model area consists of north-south oriented, narrow, regularly spaced parallel ridges and valleys that are characteristic of fluted glaciated surfaces. Thick accumulations of peat occur along the axis of East Fork Chimacum Creek and provide rich soils for agricultural use. The study area is underlain by a north-thickening sequence of unconsolidated glacial (till and outwash) and interglacial (fluvial and lacustrine) deposits, and sedimentary and igneous bedrock units that crop out along the margins and the western interior of the model area. Six hydrogeologic units in the model area form the basis of the groundwater-flow model. They are represented by model layers UC (upper confining), UA (upper aquifer), MC (middle confining), LA (lower aquifer), LC (lower confining), and OE (bedrock). Groundwater flow in the Chimacum Creek Basin and vicinity was simulated using the groundwater-flow model, MODFLOW-2005. The finite-difference model grid comprises 245 columns, 313 rows, and 6 layers. Each model cell has a horizontal dimension of 200 × 200 feet (ft). The thickness of model layers varies throughout the model area and ranges from 5 ft in the non-bedrock units to more than 2,400 ft in the bedrock. Groundwater flow was simulated for steady-state conditions, which were simulated for calibration of the model using average recharge, discharge, and water levels for the 180-month period October 1994–September 2009. The model as calibrated has a mean residual of 4.5 ft and a standard error on the mean of 2.1 ft for heads, and 0.64±0.42 cubic feet per second for streamflows. After the model was calibrated, a Current Conditions simulation was developed to reflect current (October 2008–September 2009) hydrologic conditions, with representative pumping, return flows, and “normal” recharge (based on National Weather Service average precipitation for 1981 to 2010). The Current Conditions simulation was used to estimate current flow quantities, and as a basis to compare other simulations.Simulated steady-state inflow to the model area from precipitation and secondary recharge, or “return flow,” was 16,347 acre-feet per year (acre-ft/yr); groundwater inflow from other basins to the north of the model boundary was 1,518 acre-ft/yr (net, 3,114 acre-ft/yr in and 1,596 acre-ft/yr out) and simulated inflow from lake leakage was 613 acre-ft/yr (net, 684 acre-ft/yr in and 71 acre-ft/yr out). Simulated outflow from the model primarily was through discharge to Puget Sound (10,022 acre-ft/yr), streams (5,424 acre-ft/yr ), springs and seeps (1,521 acre-ft/yr), and through withdrawals from wells (1,506 acre-ft/yr). Four simulations were formulated using the calibrated model—one to represent current conditions (2009, the end of the period used for calibration) and three to provide representative examples of how the model can be used to evaluate the relative effects of potential changes in groundwater withdrawals and consumptive use on groundwater levels and stream base flows: Probable Future Use, based on population projections; Full Beneficial Use, based on Jefferson County Public Utility District #1 water rights; Sanitary Sewer, based on eliminating septic return flows in the Urban Growth Area. Particle tracking was used to assess flowpaths from sources and to sinks, and the effects of the presence of irrigation wells and their depths was assessed.

Washington↗

Hydrology and glacier-lake-outburst floods (1987-2004) and water quality (1998-2003) of the Taku River near Juneau, Alaska

The Taku River Basin originates in British Columbia, Canada, and drains an area of 6,600 square miles at the U.S. Geological Survey's Taku River gaging station. Several mines operated within the basin prior to 1957, and mineral exploration has resumed signaling potential for future mining developments. The U.S. Geological Survey in cooperation with the Douglas Indian Association, Alaska Department of Environmental Conservation, and the U.S. Environmental Protection Agency conducted a water-quality and flood-hydrology study of the Taku River. Water-quality sampling of the Taku River from 1998 through 2003 established a baseline for assessing potential effects of future mining operations on water quality. The annual mean discharge of the Taku River is 13,700 cubic feet per second. The monthly mean discharge ranges from a minimum of 1,940 cubic feet per second in February to a maximum of 34,400 cubic feet per second in June. Nearly 90 percent of the annual discharge is from May through November. The highest spring discharges are sourced primarily from snowmelt and moderate discharges are sustained throughout the summer by glacial meltwaters. An ice cover usually forms over the Taku River in December persisting through the winter into March and occasionally into April. Glacier-lake-outburst floods originating from two glacier-dammed lakes along the margin of the Tulsequah Glacier in British Columbia, Canada, are the source of the greatest peak discharges on the Taku River. The largest flood during the period of record was 128,000 cubic feet per second on June 25, 2004, resulting from an outburst of Lake No Lake. Lake No Lake is the larger of the two lakes. The outburst-flood contribution to peak discharge was 80,000 cubic feet per second. The volume discharged from Lake No Lake is relatively consistent indicating drainage may be triggered when the lake reaches a critical stage. This suggests prediction of the timing of these outburst floods might be possible if lake-stage data were available. Further increases in the volume of Lake No Lake are unlikely as all tributary glaciers have retreated out of the lake basin. Decreasing outburst-flood volumes from Tulsequah Lake suggests a continued decline in the volume of this lake. Physical and chemical parameters and concentrations of basic water-quality constituents indicate good water quality. Samples collected at the Taku River gaging station contained low concentrations of trace elements in the dissolved phase. Trace elements sampled were within acceptable limits when compared with the Alaska Department of Environmental Conservation aquatic-life criteria for fresh waters. The highest concentrations of total trace elements sampled were collected during glacial-outburst floods and likely are associated with suspended sediments. Total trace-element concentrations generally increase with increasing water discharge, although a high correlation for all constituents sampled does not always exist.

Alaska↗

Mercury methylation and bioaccumulation in Sinclair Inlet, Kitsap County, Washington

The U.S. Geological Survey evaluated the transformation of mercury to bioavailable methylmercury in Sinclair Inlet, Kitsap County, Washington, and assessed the effect of the transformation processes on the mercury burden in marine organisms and sediment. In August 2008, samples of sediment, water, and biota from six sites in Sinclair Inlet and three bays representative of Puget Sound embayments were collected. The extensive sediment sampling included analysis of methylmercury in sediment and porewater, estimates of methylation production potential, and analyses of ancillary constituents associated with organic carbon and reduction-oxidation (redox) conditions to assist in interpreting the mercury results. Analyses of methylmercury in water overlying incubated cores provided an estimate of the release of methylmercury to the water column. Collection of samples for mercury species in the aqueous, particulate (suspended solids), and biological phases, and for ancillary carbon and nitrogen constituents in surface water, continued, on about a monthly schedule, at four stations through August 2009. In February, June, and August 2009, seasonal sediment samples were collected at 20 stations distributed between greater Sinclair Inlet and Operable Unit B Marine of the Bremerton naval complex, Bremerton, Washington, to examine geographical and seasonal patterns of mercury biogeochemistry of sediment in Sinclair Inlet. At six of these seasonal sediment stations, porewater was collected and triplicate core incubation experiments were done. Median sediment-methylmercury concentrations were not statistically different between the representative bays and Sinclair Inlet. The percentage of sediment methylmercury (relative to total mercury) was actually lower in the Sinclair Inlet sites compared with the representative bays, reflecting the higher sediment total mercury concentration for the Sinclair Inlet stations compared with the representative bays. Likewise, median sediment methylmercury concentrations were not statistically different between the greater Sinclair Inlet stations and the Bremerton naval complex stations; whereas the percentage of sediment methylmercury to total mercury was lower in the Bremerton naval complex due to higher sediment total mercury concentrations than the greater Sinclair Inlet stations. The biogeochemical characteristics of each station, measured by redox, organic carbon, and the seasonal availability of nutrients controlled methylmercury biogeochemistry. Mercury methylation production potential was a function of temperature, concentration of total mercury in sediment, and the percentage of ferrous iron (relative to total measured iron) across all sites. Methylmercury porewater concentrations were best described by using concentrations of dissolved organic carbon and reduction-oxidation conditions. Likewise, the variable fluxes of methylmercury from incubated cores were best described using dissolved organic carbon and reduction-oxidation conditions. Sinclair Inlet exhibited the classic Puget Sound biological cycle, with spring and autumn phytoplankton blooms resulting in depletion of nitrate, orthophosphate, and silicate in the surface water. Although variable in timing between 2008 and 2009, a strong corresponding seasonal trend of increased availability, incorporation, and bioaccumulation of methylmercury into the food web of Sinclair Inlet occurred during the early spring and summer growing season.

Washington↗

Geology, water-quality, hydrology, and geomechanics of the Cuyama Valley groundwater basin, California, 2008--12

To assess the water resources of the Cuyama Valley groundwater basin in Santa Barbara County, California, a series of cooperative studies were undertaken by the U.S. Geological Survey and the Santa Barbara County Water Agency. Between 2008 and 2012, geologic, water-quality, hydrologic and geomechanical data were collected from selected sites throughout the Cuyama Valley groundwater basin. Geologic data were collected from three multiple-well groundwater monitoring sites and included lithologic descriptions of the drill cuttings, borehole geophysical logs, temperature logs, as well as bulk density and sonic velocity measurements of whole-core samples. Generalized lithologic characterization from the monitoring sites indicated the water-bearing units in the subsurface consist of unconsolidated to partly consolidated sand, gravel, silt, clay, and occasional cobbles within alluvial fan and stream deposits. Analysis of geophysical logs indicated alternating layers of finer- and coarser-grained material that range from less than 1 foot to more than 20 feet thick. On the basis of the geologic data collected, the principal water-bearing units beneath the monitoring-well sites were found to be composed of younger alluvium of Holocene age, older alluvium of Pleistocene age, and the Tertiary-Quaternary Morales Formation. At all three sites, the contact between the recent fill and younger alluvium is approximately 20 feet below land surface. Water-quality samples were collected from 12 monitoring wells, 27 domestic and supply wells, 2 springs, and 4 surface-water sites and were analyzed for a variety of constituents that differed by site, but, in general, included trace elements; nutrients; dissolved organic carbon; major and minor ions; silica; total dissolved solids; alkalinity; total arsenic and iron; arsenic, chromium, and iron species; and isotopic tracers, including the stable isotopes of hydrogen and oxygen, activities of tritium, and carbon-14 abundance. Of the 39 wells sampled, concentrations of total dissolved solids and sulfate from 38 and 37 well samples, respectively, were greater than the U.S. Environmental Protection Agency’s secondary maximum contaminant levels. Concentrations greater than the maximum contaminant levels for nitrate were observed in five wells and were observed for arsenic in four wells. Differences in the stable-isotopic values of hydrogen and oxygen among groundwater samples indicated that water does not move freely between different formations or between different zones within the Cuyama Valley. Variations in isotopic composition indicated that recharge is derived from several different sources. The age of the groundwater, expressed as time since recharge, was between 600 and 38,000 years before present. Detectable concentrations of tritium indicated that younger water, recharged since the early 1950s, is present in parts of the groundwater basin. Hydrologic data were collected from 12 monitoring wells, 56 domestic and supply wells, 3 surface-water sites, and 4 rainfall-gaging stations. Rainfall in the valley averaged about 8 inches annually, whereas the mountains to the south received between 12 and 19 inches. Stream discharge records showed seasonal variability in surface-water flows ranging from no-flow to over 1,500 cubic feet per second. During periods when inflow to the valley exceeds outflow, there is potential recharge from stream losses to the groundwater system Water-level records included manual quarterly depth-to-water measurements collected from 68 wells, time-series data collected from 20 of those wells, and historic water levels from 16 wells. Hydrographs of the manual measurements showed declining water levels in 16 wells, mostly in the South-Main zone, and rising water levels in 14 wells, mostly in the Southern Ventucopa Uplands. Time-series hydrographs showed daily, seasonal, and longer-term effects associated with local pumping. Water-level data from the multiple-well monitoring sites indicated seasonal fluctuations as great as 80 feet and water-level differences between aquifers as great as 40 feet during peak pumping season. Hydrographs from the multiple-well groundwater monitoring sites showed vertical hydraulic gradients were upward during the winter months and downward during the irrigation season. Historic hydrographs showed water-level declines in the Southern-Main, Western Basin, Caliente Northern-Main, and Southern Sierra Madre zone ranging from 1 to 7 feet per year. Hydrographs of wells in the Southern Ventucopa Uplands zone showed several years with marked increases in water levels that corresponded to increased precipitation in the Cuyama Valley. Investigation of hydraulic properties included hydraulic conductivity and transmissivity estimated from aquifer tests performed on 63 wells. Estimates of horizontal hydraulic conductivity ranged from about 1.5 to 28 feet per day and decreased with depth. The median estimated hydraulic conductivity for the older alluvium was about five times that estimated for the Morales Formation. Estimates of transmissivity ranged from 560 to 163,400 gallons per day per foot and decreased with depth. The median estimated transmissivity for the younger alluvium was about three times that estimated for the older alluvium. Geomechanical analysis included land-surface elevation changes at five continuously operating global positioning systems (GPS) and land-subsidence detection at five interferometric synthetic aperture radar (InSAR) reference points. Analysis of data collected from continuously operating GPS stations showed the mountains to the south and west moved upward about 1 millimeter (mm) annually, whereas the station in the center of the Southern-Main zone moved downward more than 7 mm annually, indicating subsidence. It is likely that this subsidence is inelastic (permanent) deformation and indicates reduced storage capacity in the aquifer sediments. Analysis of InSAR data showed local and regional changes that appeared to be dependent, in part, on the time span of the interferogram, seasonal variations in pumping, and tectonic uplift. Long-term InSAR time series showed a total maximum detected subsidence rate of approximately 12 mm per year at one location and approximately 8 mm per year at a second location, while short-term InSAR time series showed maximum subsidence of about 15 mm at one location and localized maximum uplift of about 10 mm at another location.

California↗

Pathogenic bacteria and microbial-source tracking markers in Brandywine Creek Basin, Pennsylvania and Delaware, 2009-10

The City of Wilmington, Delaware, is in the downstream part of the Brandywine Creek Basin, on the main stem of Brandywine Creek. Wilmington uses this stream, which drains a mixed-land-use area upstream, for its main drinking-water supply. Because the stream is used for drinking water, Wilmington is in need of information about the occurrence and distribution of specific fecally derived pathogenic bacteria (disease-causing bacteria) and their relations to commonly measured fecal-indicator bacteria (FIB), as well as information regarding the potential sources of the fecal pollution and pathogens in the basin. This study focused on five routinely sampled sites within the basin, one each on the West Branch and the East Branch of Brandywine Creek and at three on the main stem below the confluence of the West and East Branches. These sites were sampled monthly for 1 year. Targeted event samples were collected on two occasions during high flow and two occasions during normal flow. On the basis of this study, high flows in the Brandywine Creek Basin were related to increases in FIB densities, and in the frequency of selected pathogen and source markers, in the West Branch and main stem of Brandywine Creek, but not in the East Branch. Water exceeding the moderate fullbody-contact single-sample recreational water-quality criteria (RWQC) for Escherichia coli (E. coli) was more likely to contain selected markers for pathogenic E. coli ( eaeA , stx 1 , and rfb O157 gene markers) and bovine fecal sources ( E. hirae and LTIIa gene markers), whereas samples exceeding the enterococci RWQC were more likely to contain the same pathogenic markers but also were more likely to carry a marker indicative of human source ( esp gene marker). On four sample dates, during high flow between October and March, the West Branch was the only observed potential contributor of selected pathogen and bovine source markers to the main stem of Brandywine Creek. Indeed, the stx 2 marker, which indicates a highly virulent type of pathogenic E. coli , was found only in the West Branch and main stem at high flow but was not found in the East Branch under similar conditions. However, it must be noted that throughout the entire year of sampling there were occasions, during both high and normal flows, when both the East and West Branches were potential contributors of pathogen and microbial-source tracking markers to the main stem. Therefore, this study indicates that under selected conditions (high flow, October through March), West Branch Brandywine Creek Basin was the most likely source of elevated FIB densities in the main stem. These elevated densities are associated with more frequent detection of selected pathogenic E. coli markers ( rfb O157 stx 1 ) and are associated with MST markers of bovine source. However, during other times of the year, both the West Branch and East Branch Basins are acting as potential sources of FIB and fecally derived pathogens.

Pennsylvania, Delaware↗

Synthesis of data from high-frequency nutrient and associated biogeochemical monitoring for the Sacramento–San Joaquin Delta, northern California

Executive Summary This report is the second in a series of three reports that provide information about high-frequency (HF) nutrient and biogeochemical monitoring in the Sacramento–San Joaquin Delta of northern California (Delta). The purpose of this report is to synthesize the data available from a nutrient and water-quality HF (about every 15 minutes) monitoring network operated by the U.S. Geological Survey in the northern Delta. In this report, we describe the network and focus on the purpose of each station. We then present and discuss the available data, at various timescales—first at the monthly, seasonal, and inter-annual timescales, and second, for comparison, at the tidal and event timescales. As expected, we determined that there is substantial variability in nitrate-N concentrations at short timescales within hours, but also significant variability at longer timescales such as months or years. Resolving this variability is made possible by the HF data, with the largest variability caused by storms, tides, and diel biological processes. Given this large temporal variability, calculations of cumulative nutrient fluxes (for example, daily, monthly, or annual loads) is difficult without HF data. For example, in the Cache Slough, calculation of the annual load without the tidal variability resulted in a 30 percent underestimation of the true annual load value. We conclude that HF measurements are important for accurate determination of fluxes and loads in tidal environments, but, more importantly, provide important insights into processes and rates of nutrient cycling. This report, along with the other two reports of this series (Bergamaschi and others, 2017; Kraus, Bergamaschi, and others, 2017), was drafted in cooperation with the Delta Regional Monitoring Program to help scientists, managers, and planners understand how HF data improve our understanding of nutrient sources and sinks, drivers, and effects in the Delta. The first report in the series (Kraus, Bergamaschi, and others, 2017) provides an introduction to the reasons for and fundamental concepts behind using HF monitoring measurements, including a brief summary of nutrient status and trends in the Delta and an extensive literature review showing how and where other research and monitoring programs have used HF monitoring to improve our understanding of nutrient cycling. The report covers the various technologies available for HF nutrient monitoring and presents the different ways HF monitoring instrumentation may be used for fixed station and spatial assessments. Finally, it presents numerous examples of how HF measurements are currently (2017) being used in the Delta to examine how nutrients and nutrient cycling are related to aquatic habitat conditions. The third report in the series (Bergamaschi and others, 2017) provides the background, principles, and considerations for designing an HF nutrient-monitoring network for the Delta to address high-priority, nutrient-management questions. The report starts with discussion of the high‑priority management questions to be addressed, continues through discussion of the questions and considerations that place demands and constraints on network design, discusses the principles applicable to network design, and concludes with the presentation of three example nutrient-monitoring network designs for the Delta, proposed to address high-priority questions identified by the Delta Regional Monitoring Program (Delta Regional Monitoring Program Technical Advisory Committee, 2015).

California↗

Coeur d'Alene Lake, Idaho: Insights gained From limnological studies of 1991-92 and 2004-06

More than 100 years of mining and processing of metal-rich ores in northern Idaho’s Coeur d’Alene River basin have resulted in widespread metal contamination of the basin’s soil, sediment, water, and biota, including Coeur d’Alene Lake. Previous studies reported that about 85 percent of the bottom of Coeur d’Alene Lake is substantially enriched in antimony, arsenic, cadmium, copper, lead, mercury, silver, and zinc. Nutrients in the lake also are a major concern because they can change the lake’s trophic status—or level of biological productivity—which could result in secondary releases of metals from contaminated lakebed sediments. This report presents insights into the limnological functioning of Coeur d’Alene Lake based on information gathered during two large-scale limnological studies conducted during calendar years 1991–92 and water years 2004–06. Both limnological studies reported that longitudinal gradients exist from north to south for decreasing water column transparency, loss of dissolved oxygen, and increasing total phosphorus concentrations. Gradients also exist for total lead, total zinc, and hypolimnetic dissolved oxygen concentrations, ranging from high concentrations in the central part of the lake to lower concentrations at the northern and southern ends of the lake. In the southern end of the lake, seasonal anoxia serves as a mechanism to release dissolved constituents such as phosphorus, nitrogen, iron, and manganese from lakebed sediments and from detrital material within the water column. Nonparametric statistical hypothesis tests at a significance level of α=0.05 were used to compare analyte concentrations among stations, between lake zones, and between study periods. The highest dissolved oxygen concentrations were measured in winter in association with minimum water temperatures, and the lowest concentrations were measured in the Coeur d’Alene Lake hypolimnion during late summer or autumn as prolonged thermal stratification restricted mixing of the oxygenated upper water column and the hypolimnion, where oxygen was consumed. Large differences in median concentrations of dissolved inorganic nitrogen were measured between the euphotic zone and hypolimnion in the deep areas of the lake. These differences in nitrogen concentrations were attributable to several limnological processes, including seasonal inflow plume routing, isolation from wind-driven circulation and associated hypolimnetic enrichment, phytoplanktonic assimilation during summer months, and benthic flux. Increased chlorophyll- a and total phosphorus concentrations were measured throughout the lake in the 2004–06 study compared with results from the 1991–92 study. No significant change in hypolimnetic dissolved inorganic nitrogen concentration throughout the lake was noted even though total nitrogen loads into the lake decreased between study periods. Total zinc and total lead decreased throughout the lake from the 1991-92 study to the 2004-06 study except in the southern part of the lake, where concentrations were typically low. Median detected nitrogen-to-phosphorus ratios decreased from the 1991–92 study to the 2004–06 study. Whereas the lake was clearly phosphorus-limited in 1991–92, in 2004–06 the lake may have been much closer to the boundary value of 7.2 that separates nitrogen from phosphorus limitation. However, due to changes in analytical reporting limits in the period between the two studies, the data are insufficiently certain to draw reliable conclusions with regard to limiting nutrients. For both studies, the trophic state of the lake was classified as oligotrophic (less productive) or mesotrophic (moderately productive), depending on the constituent used for classification. Internal circulation from wind-generated waves and changes in the lake’s thermocline are important processes for distribution of water-quality constituents throughout Coeur d’Alene Lake. Surficial distribution of trace metals throughout most of the lake, including bays, is relatively uniform. Even south of the Coeur d’Alene River mouth, lakebed sediments are contaminated with trace metals. Inflow plume routing of the two primary inflow sources, the Coeur d’Alene and St. Joe Rivers, also significantly affects the fate and transport of contaminants. Most riverine inflows appear to move through the lake as overflow during summer, interflow during spring and autumn, and underflow during winter. Benthic flux is another key transport process for contaminants in Coeur d’Alene Lake. The results of in situ benthic flux chamber experiments indicated movement of dissolved metals, nutrients, and dissolved organic carbon out of the lakebed sediments. However, the lake is an overall sink for these constituents when they are associated with particulate material.

Idaho↗

Comparison of methods for estimating ground-water recharge and base flow at a small watershed underlain by fractured bedrock in the Eastern United States

This study by the U.S. Geological Survey (USGS), in cooperation with the Agricultural Research Service (ARS), U.S. Department of Agriculture, compared multiple methods for estimating ground-water recharge and base flow (as a proxy for recharge) at sites in east-central Pennsylvania underlain by fractured bedrock and representative of a humid-continental climate. This study was one of several within the USGS Ground-Water Resources Program designed to provide an improved understanding of methods for estimating recharge in the eastern United States. Recharge was estimated on a monthly and annual basis using four methods?(1) unsaturated-zone drainage collected in gravity lysimeters, (2) daily water balance, (3) water-table fluctuations in wells, and (4) equations of Rorabaugh. Base flow was estimated by streamflow-hydrograph separation using the computer programs PART and HYSEP. Estimates of recharge and base flow were compared for an 8-year period (1994-2001) coinciding with operation of the gravity lysimeters at an experimental recharge site (Masser Recharge Site) and a longer 34-year period (1968-2001), for which climate and streamflow data were available on a 2.8-square-mile watershed (WE-38 watershed). Estimates of mean-annual recharge at the Masser Recharge Site and WE-38 watershed for 1994-2001 ranged from 9.9 to 14.0 inches (24 to 33 percent of precipitation). Recharge, in inches, from the various methods was: unsaturated-zone drainage, 12.2; daily water balance, 12.3; Rorabaugh equations with PULSE, 10.2, or RORA, 14.0; and water-table fluctuations, 9.9. Mean-annual base flow from streamflow-hydrograph separation ranged from 9.0 to 11.6 inches (21-28 percent of precipitation). Base flow, in inches, from the various methods was: PART, 10.7; HYSEP Local Minimum, 9.0; HYSEP Sliding Interval, 11.5; and HYSEP Fixed Interval, 11.6. Estimating recharge from multiple methods is useful, but the inherent differences of the methods must be considered when comparing results. For example, although unsaturated-zone drainage from the gravity lysimeters provided the most direct measure of potential recharge, it does not incorporate spatial variability that is contained in watershed-wide estimates of net recharge from the Rorabaugh equations or base flow from streamflow-hydrograph separation. This study showed that water-level fluctuations, in particular, should be used with caution to estimate recharge in low-storage fractured-rock aquifers because of the variability of water-level response among wells and sensitivity of recharge to small errors in estimating specific yield. To bracket the largest range of plausible recharge, results from this study indicate that recharge derived from RORA should be compared with base flow from the Local-Minimum version of HYSEP.

Scientific Investigations Report↗

Methods for estimating selected low-flow frequency statistics and mean annual flow for ungaged locations on streams in North Georgia

The U.S. Geological Survey, in cooperation with the Georgia Department of Natural Resources, Environmental Protection Division, developed regional regression equations for estimating selected low-flow frequency and mean annual flow statistics for ungaged streams in north Georgia that are not substantially affected by regulation, diversions, or urbanization. Selected low-flow frequency statistics and basin characteristics for 56 streamgage locations within north Georgia and 75 miles beyond the State’s borders in Alabama, Tennessee, North Carolina, and South Carolina were combined to form the final dataset used in the regional regression analysis. Because some of the streamgages in the study recorded zero flow, the final regression equations were developed using weighted left-censored regression analysis to analyze the flow data in an unbiased manner, with weights based on the number of years of record. The set of equations includes the annual minimum 1- and 7-day average streamflow with the 10-year recurrence interval (referred to as 1Q10 and 7Q10), monthly 7Q10, and mean annual flow. The final regional regression equations are functions of drainage area, mean annual precipitation, and relief ratio for the selected low-flow frequency statistics and drainage area and mean annual precipitation for mean annual flow. The average standard error of estimate was 13.7 percent for the mean annual flow regression equation and ranged from 26.1 to 91.6 percent for the selected low-flow frequency equations. The equations, which are based on data from streams with little to no flow alterations, can be used to provide estimates of the natural flows for selected ungaged stream locations in the area of Georgia north of the Fall Line. The regression equations are not to be used to estimate flows for streams that have been altered by the effects of major dams, surface-water withdrawals, groundwater withdrawals (pumping wells), diversions, or wastewater discharges. The regression equations should be used only for ungaged sites with drainage areas between 1.67 and 576 square miles, mean annual precipitation between 47.6 and 81.6 inches, and relief ratios between 0.146 and 0.607; these are the ranges of the explanatory variables used to develop the equations. An attempt was made to develop regional regression equations for the area of Georgia south of the Fall Line by using the same approach used during this study for north Georgia; however, the equations resulted with high average standard errors of estimates and poorly predicted flows below 0.5 cubic foot per second, which may be attributed to the karst topography common in that area. The final regression equations developed from this study are planned to be incorporated into the U.S. Geological Survey StreamStats program. StreamStats is a Web-based geographic information system that provides users with access to an assortment of analytical tools useful for water-resources planning and management, and for engineering design applications, such as the design of bridges. The StreamStats program provides streamflow statistics and basin characteristics for U.S. Geological Survey streamgage locations and ungaged sites of interest. StreamStats also can compute basin characteristics and provide estimates of streamflow statistics for ungaged sites when users select the location of a site along any stream in Georgia.

Georgia↗