Search USGSSearch

SEARCH · Search USGS

Results for “Irrigation and Drainage Systems”

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

Hydrogeology and simulation of the effects of reclaimed-water application in west Orange and southeast Lake counties, Florida

Wastewater reclamation and reuse has become increasingly popular as water agencies search for alternative water-supply and wastewater-disposal options. Several governmental agencies in central Florida currently use the land-based application of reclaimed water (wastewater that has been treated beyond secondary treatment) as a management alternative to surface-water disposal of wastewater. Water Conserv II, a water reuse project developed jointly by Orange County and the City of Orlando, began operation in December 1986. In 1995, the Water Conserv II facility distributed approximately 28 Mgal/d of reclaimed water for discharge to rapid-infiltration basins (RIBs) and for use as agricultural irrigation. The Reedy Creek Improvement District (RCID) began operation of RIBs in September 1990, and in 1995 these RIBs received approximately 6.7 Mgal/d of reclaimed water. Analyses of existing data and data collected during the course of this study were combined with ground-water flow modeling and particle-tracking analyses to develop a process-oriented evaluation of the regional effects of reclaimed water applied by Water Conserv II and the RCID RIBs on the hydrology of west Orange and southeast Lake Counties. The ground-water flow system beneath the study area is a multi-aquifer system that consists of a thick sequence of highly permeable carbonate rocks overlain by unconsolidated sediments. The hydrogeologic units are the unconfined surficial aquifer system, the intermediate confining unit, and the confined Floridan aquifer system, which consists of two major permeable zones, the Upper and Lower Floridan aquifers, separated by the less permeable middle semiconfining unit. Flow in the surficial aquifer system is dominated regionally by diffuse downward leakage to the Floridan aquifer system and is affected locally by lateral flow systems produced by streams, lakes, and spatial variations in recharge. Ground water generally flows laterally through the Upper Floridan aquifer aquifer to the north and east. Many of the lakes in the study area are landlocked because the mantled karst environment precludes a well developed network of surface-water drainage. The USGS three-dimensional ground-water flow model MODFLOW was used to simulate ground-water flow in the surficial and Floridan aquifer systems. A steady-state calibration to average 1995 conditions was performed by using a parameter estimation program to vary values of surficial aquifer system hydraulic conductivity, intermediate confining unit leakance, and Upper Floridan aquifer transmissivity. The calibrated model generally produced simulated water levels in close agreement with measured water levels and was used to simulate the hydrologic effects of reclaimed-water application under current (1995) and proposed future conditions. In 1995, increases of up to about 40 ft in the water table and less than 5 ft in the Upper Floridan aquifer potentiometric surface had occurred as a result of reclaimed-water application. The largest increases were under RIB sites. An average traveltime of 10 years at Water Conserv II and 7 years at the RCID RIBs was required for reclaimed water to move from the water table to the top of the Upper Floridan aquifer. Approximately 67 percent of the reclaimed water applied at the RCID RIB site recharged the Floridan aquifer system, whereas 33 percent discharged from the surficial aquifer system to surface-water features; 99 percent of the reclaimed water applied at Water Conserv II recharged the Floridan aquifer system, whereas only 1 percent discharged from the surficial aquifer system to surface-water features. The majority of reclaimed water applied at both facilities probably will ultimately discharge from the Floridan aquifer system outside the model boundaries. Proposed future conditions were assumed to consist of an additional 11.7 Mgal/d of reclaimed water distributed by the Water Conserv II and RCID facilities. Increases of up to about 20 ft in the water table and 2 ft in the potentiometric surface of the Upper Floridan aquifer were simulated. The directions of reclaimed water movement through the ground-water system generally were similar to those under 1995 conditions. However, the greater reclaimed-water application rate at the RCID RIBs caused approximately half of the RCID reclaimed water to discharge to surface-water features and half to recharge the Floridan aquifer system.

Florida

Ground water in the Thousand Oaks area, Ventura County, California

The ground-water basin beneath the city of Thousand Oaks, Calif. , corresponds closely in area with the surface-water drainage basin of Conejo Valley. Before World War II there was little ground-water development. After World War II, urban development put a stress on the ground-water basin; many wells were drilled and water levels in wells were drawn down as much as 300 feet in places. Beginning in 1963, imported water replaced domestic and municipal ground-water systems, and water levels rapidly recovered to predevelopment levels or nearly so. Most of the ground water in the Thousand Oaks area is stored in fractured basalt of the middle Miocene Conejo Volcanics. Depending on the degree of occurrence of open fractures and cavities in the basalt, recoverable ground water in the upper 300 to 500 feet of aquifer is estimated to be between 400,000 and 600,000 acre-feet. The yield of water from wells in the area ranges from 17 to 1,080 gallons per minute. Most of the ground-water in the eastern part of the valley is high insulfate and has a dissolved-solids concentration greater than 1,000 milligrams per liter. In the western part of the valley the ground-water is mostly of a bicarbonate type, and the dissolved-solids concentration is less than 800 milligrams per liter. In most areas of Conejo Valley, ground-water is a viable resource for irrigation of public lands and recreation areas. (USGS)

California

Possible changes in ground-water flow to the Pecos River caused by Santa Rosa Lake, Guadalupe County, New Mexico

In 1980 Santa Rosa Dam began impounding water on the Pecos River about 7 miles north of Santa Rosa, New Mexico, to provide flood control, sediment control, and storage for irrigation. Santa Rosa Lake has caused changes in the groundwater flow system, which may cause changes in the streamflow of the Pecos River that cannot be detected at the present streamflow gaging stations. Data collected at these stations are used to measure the amount of water available for downstream users. A three-dimensional groundwater flow model for a 950 sq mi area between Anton Chico and Puerto de Luna was used to simulate the effects of Santa Rosa Lake on groundwater flow to a gaining reach of the Pecos River for lake levels of 4,675, 4,715, 4,725, 4,750, 4,776, and 4,797 feet above sea level and durations of impoundment of 30, 90, 182, and 365 days for all levels except 4 ,797 feet. These simulations indicated that streamflow in the Pecos River could increase by as much as 2 cu ft/sec between the dam and Puerto de Luna if the lake level were maintained at 4 ,797 feet for 90 days or 4,776 feet for 1 year. About 90% of this increased streamflow would occur < 0.5 mi downstream from the dam, some of which would be measured at the streamflow gaging station located 0.2 mile downstream from the dam. Simulations also indicated that the lake will affect groundwater flow such that inflow to the study area may be decreased by as much as 1.9 cu ft/sec. This water may leave the Pecos River drainage basin or be diverted back to the Pecos River downstream from the gaging station near Puerto de Luna. In either case, this quantity represents a net loss of water upstream from Puerto de Luna. Most simulations indicated that the decrease in groundwater flow into the study area would be of about the same quantity as the simulated increase in streamflow downstream from the dam. Therefore, the net effect of the lake on the flow of the Pecos River in the study area appears to be negligible. Model simulations indicated that effect of lake levels below 4 ,750 feet on water levels in observation wells completed in the San Andres Limestone could not be distinguished from the effects of other hydrologic stresses. (Author 's abstract)

Water-Resources Investigations Report

Water resources of Rockland Basin, southeastern Idaho

Rockland basin comprises about 320 sq mi of the Snake River drainage in southeastern Idaho. Mountain ranges bordering the basin are composed predominantly of limestone and are complexly faulted. Major aquifers include Holocene alluvium, Quaternary-Tertiary volcanic rocks, and Tertiary sedimentary rocks. Groundwater occurs under water table conditions except where it is locally confined. Groundwater discharges to springs in the Deep Creek Mountains and maintains perennial streamflow. Near the mouth of Rock Creek, groundwater movement is northward toward the Snake River. Underflow is estimated to be 51,000 acre-ft/yr. Total water yield available to Rockland basin is estimated to be 5.0 in. (85,000 acre-ft) of the estimated 17.3 in. of annual precipitation. Evapotranspiration ranges from 9.9 to 17 in./yr, depending, in part, on altitude of the land surface. An estimated 12,000 acre-ft of surface water and 3,500 acre-ft of groundwater are used annually for irrigation. Less than 100 acre-ft of water is used for public supply, domestic, and stock supplies. East Fork Rock Creek supplies the most surface water for irrigation of agricultural lands. At the present (1980) state of groundwater development in Rockland basin, streams and aquifers are hydraulically connected. Pumping of groundwater in increased quantities from wells near streams will affect groundwater movement and may diminish streamflow. There are no long-term regional water table declines at present. Continued water level monitoring of selected wells may aid in documenting effects of future management practices on the groundwater system.

idaho

Hydrology of area 54, northern Great Plains, and Rocky Mountain coal provinces, Colorado and Wyoming

A nationwide need for information characterizing hydrologic conditions in mined and potential mine areas has become paramount with the enactment of the Surface Mining Control and Reclamation Act of 1977. This report, one in a series covering the coal provinces nationwide, presents information thematically by describing single hydrologic topics through the use of brief texts and accompanying maps, graphs, or other illustrations. The summation of the topical discussions provides a description of the hydrology of the area. Area 54, in north-central Colorado and south-central Wyoming, is 1 of 20 hydrologic reporting areas of the Northern Great Plains and Rocky Mountain coal provinces. Part of the Southern Rocky Mountains and Wyoming Basin physiographic provinces, the 8,380-square-mile area is one of contrasting geology, topography, and climate. This results in contrasting hydrologic characteristics. The major streams, the North Platte, Laramie, and Medicine Bow Rivers, and their principal tributaries, all head in granitic mountains and flow into and through sedimentary basins between the mountain ranges. Relief averages 2,000 to 3,000 feet. Precipitation in the mountains may exceed 40 inches annually, much of it during the winter, which produces deep snowpacks. Snowmelt in spring and summer provides most streamflow. Precipitation in the basins averages 10 to 16 inches annually, insufficient for sustained streamflow; thus, streams originating in the basins are ephemeral. Streamflow quality is best in the mountains where dissolved-solids concentrations generally are least. These concentrations increase as streams flow through sedimentary basins. The increases are mainly natural, but some may be due to irrigation in and adjacent to the flood plains. In the North Platte River, dissolved-solids concentrations are usually less than 300 milligrams per liter; in the Laramie and the Medicine Bow Rivers, the concentrations may average 500 to 850 milligrams per liter. However, water-quality stations on the Laramie and the Medicine Bow Rivers are farther removed from the mountain sources than the stations in the North Platte drainage. Because of the semiarid climate of the basins, soils are not adequately leached. Consequently, flow in ephemeral streams usually has a larger concentration of dissolved solids than that in perennial streams, averaging 1,000 to 1,600 milligrams per liter. Aquifers containing usable ground water are combined into three groups: (1) consolidated and unconsolidated non-coal-bearing Quaternary and Upper Tertiary deposits, (2) Mesozoic and Paleozoic sedimentary rocks, and (3) Lower Tertiary and Upper Cretaceous sedimentary rocks containing coal. These aquifers are used for municipal, domestic, irrigation, and stock supplies. Well yields range from about 5 to 1,000 gallons per minute, and depend on type of aquifer, saturated thickness, and degree of fracturing. The best quality ground water usually comes from the non-coal-bearing Quaternary and Upper Tertiary rocks or the Mesozoic and Paleozoic rocks; often it is dominated by calcium and bicarbonate ions. The coal-bearing formations have a large variability in water chemistry; dominant ions may be bicarbonate or sulfate and sodium, calcium, or magnesium. Dissolved-solids concentrations are generally larger than in the former two groups. The U.S. Geological Survey operates a network of hydrologic stations to observe the streamflow and groundwater conditions. This network currently includes 31 surface-water stations and 35 observation wells; information is available for many other sites observed in the past. Data available include rate of flow, water levels, and water quality; much of the data are available in published reports or from computer storage through the National Water Data Exchange (NAWDEX) or the National Water Data Storage and Retrieval System (WATSTORE). Five formations of Late Cretaceous and early Tertiary age contain coal.

Colorado, Wyoming

An update of hydrologic conditions and distribution of selected constituents in water, Snake River Plain aquifer and perched groundwater zones, Idaho National Laboratory, Idaho, emphasis 2006-08

Since 1952, radiochemical and chemical wastewater discharged to infiltration ponds (also called percolation ponds), evaporation ponds, and disposal wells at the Idaho National Laboratory (INL) has affected water quality in the eastern Snake River Plain aquifer and perched groundwater zones underlying the INL. The U.S. Geological Survey, in cooperation with the U.S. Department of Energy, maintains groundwater monitoring networks at the INL to determine hydrologic trends, and to delineate the movement of radiochemical and chemical wastes in the aquifer and in perched groundwater zones. This report presents an analysis of water-level and water-quality data collected from aquifer and perched groundwater wells in the USGS groundwater monitoring networks during 2006-08. Water in the Snake River Plain aquifer primarily moves through fractures and interflow zones in basalt, generally flows southwestward, and eventually discharges at springs along the Snake River. The aquifer primarily is recharged from infiltration of irrigation water, infiltration of streamflow, groundwater inflow from adjoining mountain drainage basins, and infiltration of precipitation. From March-May 2005 to March-May 2008, water levels in wells generally remained constant or rose slightly in the southwestern corner of the INL. Water levels declined in the central and northern parts of the INL. The declines ranged from about 1 to 3 feet in the central part of the INL, to as much as 9 feet in the northern part of the INL. Water levels in perched groundwater wells around the Advanced Test Reactor Complex (ATRC) also declined. Detectable concentrations of radiochemical constituents in water samples from wells in the Snake River Plain aquifer at the INL generally decreased or remained constant during 2006-08. Decreases in concentrations were attributed to decreased rates of radioactive-waste disposal, radioactive decay, changes in waste-disposal methods, and dilution from recharge and underflow. In April or October 2008, reportable concentrations of tritium in groundwater ranged from 810 ? 70 to 8,570 ? 190 picocuries per liter (pCi/L), and the tritium plume extended south-southwestward in the general direction of groundwater flow. Tritium concentrations in water from wells completed in shallow perched groundwater at the ATRC were less than the reporting levels. Tritium concentrations in deep perched groundwater exceeded the reporting level in 11 wells during at least one sampling event during 2006-08 at the ATRC. Tritium concentrations from one or more zones in each well were reportable in water samples collected at various depths in six wells equipped with multi-level WestbayTM packer sampling systems. Concentrations of strontium-90 in water from 24 of 52 aquifer wells sampled during April or October 2008 exceeded the reporting level. Concentrations ranged from 2.2 ? 0.7 to 32.7 ? 1.2 pCi/L. Strontium-90 has not been detected within the eastern Snake River Plain aquifer beneath the ATRC partly because of the exclusive use of waste-disposal ponds and lined evaporation ponds rather than using the disposal well for radioactive-wastewater disposal at ATRC. At the ATRC, the strontium-90 concentration in water from one well completed in shallow perched groundwater was less than the reporting level. During at least one sampling event during 2006-08, concentrations of strontium-90 in water from nine wells completed in deep perched groundwater at the ATRC were greater than reporting levels. Concentrations ranged from 2.1?0.7 to 70.5?1.8 pCi/L. At the Idaho Nuclear Technology and Engineering Center (INTEC), the reporting level was exceeded in water from two wells completed in deep perched groundwater. During 2006-08, concentrations of cesium-137, plutonium-238, and plutonium-239, -240 (undivided), and americium-241 were less than the reporting level in water samples from all wells and all zones in wells equipped with multi-level WestbayTM packer sampling systems

Scientific Investigations Report

Geology, ground-water hydrology, geochemistry, and ground-water simulation of the Beaumont and Banning Storage Units, San Gorgonio Pass area, Riverside County, California

Ground water has been the only source of potable water supply for residential, industrial, and agricultural users in the Beaumont and Banning storage units of the San Gorgonio Pass area, Riverside County, California. Ground-water levels in the Beaumont area have declined as much as 100 feet between the early 1920s and early 2000s, and numerous natural springs have stopped flowing. In 1961, the San Gorgonio Pass Water Agency (SGPWA) entered into a contract with the California State Department of Water Resources to receive 17,300 acre-feet per year of water to be delivered by the California State Water Project (SWP) to supplement natural recharge. Currently (2005), a pipeline is delivering SWP water into the area, and the SGPWA is artificially recharging the ground-water system using recharge ponds located along Little San Gorgonio Creek in Cherry Valley with the SWP water. In addition to artificial recharge, SGPWA is considering the direct delivery of SWP water for the irrigation of local golf courses and for agricultural supply in lieu of ground-water pumpage. To better understand the potential hydrologic effects of different water-management alternatives on ground-water levels and movement in the Beaumont and Banning storage units, existing geohydrologic and geochemical data were compiled, new data from a basin-wide ground-water level and water-quality monitoring network were collected, monitoring wells were installed near the Little San Gorgonio Creek recharge ponds, geohydrologic and geochemical analyses were completed, and a ground-water flow simulation model was developed. The San Gorgonio Pass area was divided into several storage units on the basis of mapped or inferred faults. This study addresses primarily the Beaumont and Banning storage units. The geologic units in the study area were generalized into crystalline basement rocks and sedimentary deposits. The younger sedimentary deposits and the surficial deposits are the main water-bearing deposits in the San Gorgonio Pass area. The water-bearing deposits were divided into three aquifers: (1) the perched aquifer, (2) the upper aquifer, and (3) the lower aquifer based on lithologic and downhole geophysical logs. Natural recharge in the San Gorgonio Pass area was estimated using INFILv3, a deterministic distributed- parameter precipitation-runoff model. The INFILv3 model simulated that the potential recharge of precipitation and runoff in the Beaumont and Banning storage units was about 3,710 acre-feet per year and that the potential recharge in 28 sub-drainage basins upstream of the storage units was about 6,180 acre-feet per year. The water supply for the Beaumont and Banning storage units is supplied by pumping ground water from wells in the Canyon (Edgar and Banning Canyons), Banning Bench, Beaumont, and Banning storage units. Total annual pumpage from the Beaumont and Banning storage units ranged from about 1,630 acre-feet in 1936 to about 20,000 acre-feet in 2003. Ground-water levels declined by as much as 100 feet in the Beaumont storage unit from 1926-2003 in response to ground-water pumping of about 450,160 acre-feet during this period. Since ground-water development began in the San Gorgonio Pass area, there have been several sources of artificial recharge to the basin including return flow from applied water on crops, golf courses, and landscape; septic-tank seepage; and infiltration of storm runoff diversions and imported water into recharge ponds. Return flow from applied water and septic-tank seepage was estimated to reach a maximum of about 8,100 acre-feet per year in 2003. Owing to the great depth of water in much of study area (in excess of 150 feet), the return flow and septic-tank seepage takes years to decades to reach the water table. Stable-isotope data indicate that the source of ground-water recharge was precipitation from storms passing through the San Gorgonio Pass as opposed to runoff from the higher altitudes of the San Bernar

Scientific Investigations Report

Forecasting selenium discharges to the San Francisco Bay-Delta Estuary: Ecological effects of a proposed San Luis Drain extension

Selenium discharges to the San Francisco Bay-Delta Estuary (Bay-Delta) could change significantly if federal and state agencies (1) approve an extension of the San Luis Drain to convey agricultural drainage from the western San Joaquin Valley to the North Bay (Suisun Bay, Carquinez Strait, and San Pablo Bay); (2) allow changes in flow patterns of the lower San Joaquin River and Bay-Delta while using an existing portion of the San Luis Drain to convey agricultural drainage to a tributary of the San Joaquin River; or (3) revise selenium criteria for the protection of aquatic life or issue criteria for the protection of wildlife. Understanding the biotransfer of selenium is essential to evaluating effects of selenium on Bay-Delta ecosystems. Confusion about selenium threats to fish and wildlife stem from (1) monitoring programs that do not address specific protocols necessary for an element that bioaccumulates; and (2) failure to consider the full complexity of the processes that result in selenium toxicity. Past studies show that predators are more at risk from selenium contamination than their prey, making it difficult to use traditional methods to predict risk from environmental concentrations alone. This report presents an approach to conceptualize and model the fate and effects of selenium under various load scenarios from the San Joaquin Valley. For each potential load, progressive forecasts show resulting (1) water-column concentration; (2) speciation; (3) transformation to particulate form; (4) particulate concentration; (5) bioaccumulation by invertebrates; (6) trophic transfer to predators; and (7) effects on those predators. Enough is known to establish a first-order understanding of relevant conditions, biological response, and ecological risks should selenium be discharged directly into the North Bay through a conveyance such as a proposed extension of the San Luis Drain. The approach presented here, the Bay-Delta selenium model, determines the mass, fate, and effects of selenium released to the Bay-Delta through use of (1) historical land-use, drainage, alluvial-fill, and runoff databases; (2) existing knowledge concerning biogeochemical reactions and physiological parameters of selenium (e.g., speciation, partitioning between dissolved and particulate forms, and bivalve assimilation efficiency); and (3) site-specific data mainly from 1986 to 1996 for clams and bottom-feeding fish and birds. Selenium load scenarios consider effluents from North Bay oil refineries and discharges of agricultural drainage from the San Joaquin Valley to enable calculation of (a) a composite freshwater endmember selenium concentration at the head of the estuary; and (b) a selenium concentration at a selected seawater location (Carquinez Strait) as a foundation for modeling. Analysis of selenium effects also takes into account the mode of conveyance for agricultural drainage (i.e., the San Luis Drain or San Joaquin River); and flows of the Sacramento River and San Joaquin River on a seasonal or monthly basis. Load scenarios for San Joaquin Valley mirror predictions made since 1955 of a worsening salt (and by inference, selenium) build-up exacerbated by an arid climate and massive irrigation. The reservoir of selenium in the San Joaquin Valley is sufficient to provide loading at an annual rate of approximately 42,500 pounds of selenium to a Bay-Delta disposal point for 63 to 304 years at the lower range of projections presented here, even if influx of selenium from the California Coast Ranges could be curtailed. Disposal of wastewaters on an annual basis outside of the San Joaquin Valley may slow the degradation of valley resources, but drainage alone cannot alleviate the salt and selenium build-up in the San Joaquin Valley, at least within a century. Load scenarios also show the different proportions of selenium loading to the Bay-Delta. Oil refinery loads from 1986 to 1992 ranged from 8.5 to 20 pounds of selenium per day; with treatment and cleanup, loads decreased to 3.0 pounds of selenium per day in 1999. In contrast, San Joaquin Valley agricultural drainage loads disposed of in a San Luis Drain extension could range from 45 to 117 pounds of selenium per day across a set of historical and future conditions. Components of this valley-wide load include five source subareas (i.e., Grassland, Westlands, Tulare, Kern, and Northern) defined by water and drainage management. Loads vary per subarea mainly because of proximity of the subarea to geologic sources of selenium and irrigation history. Loads from the Sacramento River, depending on flow conditions, range from 0.8 to 10 pounds of selenium per day. Loads from the San Joaquin River vary depending on restoration and flow conditions, which are considered. A consistent picture of ecological risk emerges under modeled selenium discharges from a proposed San Luis Drain extension. The threat to the estuary is greatest during low flow seasons and critically dry years. Where selenium undergoes reactions typical of low flow or longer residence time, highly problematic bioaccumulation in prey (food) is forecast. Surf scoter, greater and lesser scaup, and white sturgeon appear to be most at risk because these Bay-Delta predators feed on deposit and filter-feeding bivalves. Recent findings add Sacramento splittail and Dungeness crab to that list. During the low flow season of critically dry years, forecasted selenium concentrations in water, particulate matter, prey (diet), and predator tissue exceed guidelines with a high certainty of producing adverse effects under the most likely load scenario from a proposed San Luis Drain extension. High flows afford some protection under certain conditions in modeled San Joaquin River scenarios. However, meeting a combined goal of releasing a specific load during maximum flows and keeping selenium concentrations in the river below a certain objective to protect against bioaccumulation may not always be attainable. Management of the San Joaquin River on a constant concentration basis also could create problematic bioaccumulation during a wet year, especially during the low flow season, because high flows translate to high loads that are not always offset by seasonal river inflows. Prior to refinery cleanup, selenium contamination was sufficient to threaten reproduction in key species within the Bay-Delta ecosystems and human health advisories were posted based on selenium concentrations in tissues of diving ducks. During this time, selenium concentrations in the Bay-Delta were well below the most stringent recommended water quality criterion [1 microgram per liter (1 µg/L)]. Enhanced biogeochemical transformations to bioavailable particulate selenium and efficient bioaccumulation by bivalves characterized the system. If these biogeochemical conditions continue to prevail and agricultural selenium sources replace or exceed refinery sources, ecological forecasts suggest the risk of adverse effects will be difficult to eliminate under an out-of-valley resolution to the selenium problem. The Bay-Delta selenium model presented here is a systematic approach for conducting forecasts of the ecological effects from selenium on aquatic food webs. It is a new tool that links and models the major processes leading from loads through consumer organisms to predators. It also is a feasible approach for site-specific analysis and could provide a framework for developing new protective selenium foodweb guidelines and predator criteria. Model components that help ensure understanding ecosystems and the basis of environmental protection are (1) contaminant concentrations and speciation in sources, such as particulate material, that most influence bioavailability; (2) bioaccumulation models that calculate concentrations in diet, specifically in bivalves of the Bay-Delta that act as sensitive indicators of selenium contamination; (3) food-web type that determines what animals are threatened and when; and (4) multiple media concentrations (water, particulate material, and tissue of prey and predators) that, in-combination, determine risk or hazard.

California

Hydrogeology, land-surface subsidence, and documentation of the Gulf Coast Land Subsidence and Groundwater-Flow (GULF) model, southeast Texas, 1897–2018

Executive Summary As a part of the Texas Water Development Board groundwater availability modeling program, the U.S. Geological Survey developed the Gulf Coast Land Subsidence and Groundwater-Flow model (hereinafter, the “GULF model”) and ensemble to simulate groundwater flow and land-surface subsidence in the northern part of the Gulf Coast aquifer system (the study area) in Texas from predevelopment (1897) through 2018. Since the publication of a previous groundwater model for the greater Houston area in 2012, there have been changes to the distribution of groundwater withdrawals and advances in modeling tools. To reflect these changes and to simulate more recent conditions, the GULF model was developed in cooperation with the Harris-Galveston and Fort Bend Subsidence Districts to provide an updated Groundwater Availability Model. Since the early 1900s, most of the groundwater withdrawals in the study area have been from three of the hydrogeologic units that compose the Gulf Coast aquifer system—the Chicot, Evangeline, and Jasper aquifers and, more recently, from the Catahoula confining unit. Withdrawals from these hydrogeologic units are used for municipal supply, commercial and industrial use, and irrigation purposes. Withdrawals of large quantities of groundwater in the greater Houston area have caused widespread groundwater-level declines in the Chicot, Evangeline, and Jasper aquifers of more than 300 feet (ft). Early development of the aquifer system, which began before 1900, resulted in nearly 50 percent of the eventual historical groundwater-level minimums having been reached as early as 1946 in some areas. These groundwater-level declines led to more than 9 ft of land-surface subsidence—historically in central and southeastern Harris County and Galveston County, but more recently in northern, northwestern, and western Harris County, Montgomery County, and northern Fort Bend County—from depressurization and compaction of clay and silt layers interbedded in the aquifer sediments. In a generalized conceptual model of the Gulf Coast aquifer system, water enters the groundwater system in topographically high outcrops of the hydrogeologic units in the northwestern part of the aquifer system. Groundwater that does not discharge to streams flows to intermediate and deep zones of the aquifer system southeastward of the outcrop areas where it is discharged by wells and by upward leakage in topographically low areas near the coast. The uppermost parts of the aquifer system, which include outcrop areas, are under water-table (unconfined) conditions where the groundwater is not confined under pressure. As depth increases in the aquifer system and interbedded clay and silt layers accumulate, water-table conditions evolve into confined conditions where the groundwater is under pressure. Groundwater flow and land-surface subsidence in the GULF model and ensemble were simulated by using MODFLOW 6 with the Skeletal Storage, Compaction, and Subsidence package. The model consists of six layers, one for each of the five hydrogeologic units in the northern part of the Gulf Coast aquifer system and a surficial top layer that includes part of each hydrogeologic unit. Transient groundwater flow was simulated during 1897–2018 by using a combination of multiyear, annual, and monthly stress periods. An initial steady-state stress period was configured to represent predevelopment mean annual inflows and outflows. The subsidence package used in the GULF model and ensemble uses a head-based subsidence formulation that simulates the delayed drainage response from clay and silt sediment to changes in groundwater levels. The GULF model and ensemble were history matched to groundwater-level observations at selected wells, land-surface subsidence at benchmarks, aquifer compaction at borehole extensometers, and vertical displacement from Global Positioning System stations. A Bayesian framework was used to represent uncertainty in modeled parameters and simulated outputs of interest. History matching and uncertainty quantification were performed by using a Monte Carlo approach enabled through iterative ensemble smoother software to produce an ensemble of models fit to historical data. The iterative ensemble smoother substantially reduced the computational demand of parameter estimation by approximating the first-order relation between model inputs and outputs, thereby allowing 183,207 adjustable parameters to be used for history matching at a relatively low computational and time cost. The history-matched parameter values are within the ranges of previously published values and agree with the current understanding of the spatial and temporal patterns of parameter uncertainty for the Gulf Coast aquifer system. A good agreement between the observed (or estimated) and simulated groundwater levels, land-surface subsidence, compaction, and vertical displacement was obtained across the modeled area based on qualitative and quantitative comparisons. Ensemble mean annual groundwater-flow rates to the Chicot, Evangeline, Jasper aquifers and Catahoula confining unit were 0.0–0.49 inch (in.), 0.09–0.33 in., 0.01–0.07 in., and 0.01–0.05 in., respectively. GULF model mean annual groundwater-flow rates to the Chicot, Evangeline, and Jasper aquifers and Catahoula confining unit were 0.31 in., 0.19 in., 0.03 in., and 0.03 in., respectively. The GULF-model-simulated recharge to the outcrop area was the largest inflow (75 percent), and recharge to other areas was 25 percent of the model inflow. The simulated outflows included (1) net surface-water/groundwater exchange with study area streams (50 percent), (2) groundwater use (49 percent), and (3) net surface-water/groundwater exchange with the Gulf of Mexico (1 percent). The sum of the simulated values of the outflows (1,041,973 acre-feet per year [acre-ft/yr]) and the elastic expansion of the fine-grained sediment and numerical solver error (339 acre-ft/yr) minus the inflows (654,172 acre-ft/yr) represents the reduction of storage from the Gulf Coast aquifer system (388,140 acre-ft/yr). Most of the storage depletion is caused by the long-term groundwater-level declines that have resulted primarily in inelastic compaction. The GULF model was used to estimate Jasper aquifer compaction at selected benchmarks in Montgomery County and northern Harris County, which are the primary locations of Jasper aquifer groundwater use. Simulated Jasper aquifer compaction in northern Harris County was between 0.2 and 0.5 ft, or between about 5 and 16 percent of simulated subsidence at the benchmark locations. Simulated Jasper aquifer compaction in Montgomery County was between 0.8 and 1.2 ft, or between about 33 and 57 percent of simulated subsidence at the benchmark locations.

Texas

Occurrence and transport of nitrogen in the Big Sunflower River, northwestern Mississippi, October 2009-June 2011

The Big Sunflower River Basin, located within the Yazoo River Basin, is subject to large annual inputs of nitrogen from agriculture, atmospheric deposition, and point sources. Understanding how nutrients are transported in, and downstream from, the Big Sunflower River is key to quantifying their eutrophying effects on the Gulf. Recent results from two Spatially Referenced Regressions on Watershed attributes (SPARROW models), which include the Big Sunflower River, indicate minimal losses of nitrogen in stream reaches typical of the main channels of major river systems. If SPARROW assumptions of relatively conservative transport of nitrogen are correct and surface-water losses through the bed of the Big Sunflower River are negligible, then options for managing nutrient loads to the Gulf of Mexico may be limited. Simply put, if every pound of nitrogen entering the Delta is eventually delivered to the Gulf, then the only effective nutrient management option in the Delta is to reduce inputs. If, on the other hand, it can be shown that processes within river channels of the Mississippi Delta act to reduce the mass of nitrogen in transport, other hydrologic approaches may be designed to further limit nitrogen transport. Direct validation of existing SPARROW models for the Delta is a first step in assessing the assumptions underlying those models. In order to characterize spatial and temporal variability of nitrogen in the Big Sunflower River Basin, water samples were collected at four U.S. Geological Survey gaging stations located on the Big Sunflower River between October 1, 2009, and June 30, 2011. Nitrogen concentrations were generally highest at each site during the spring of the 2010 water year and the fall and winter of the 2011 water year. Additionally, the dominant form of nitrogen varied between sites. For example, in samples collected from the most upstream site (Clarksdale), the concentration of organic nitrogen was generally higher than the concentrations of ammonia and nitrate plus nitrite; conversely, at sites farther downstream (that is, at Sunflower and Anguilla), nitrate plus nitrite concentrations were generally higher than concentrations of organic nitrogen and ammonia. In addition to the routinely collected samples, water samples from the Big Sunflower River Basin were collected using a Lagrangian sampling scheme, which attempts to follow a single mass of water through time in order to determine how it changes through processing or other pathways as the water moves downstream. Lagrangian sampling was conducted five times during the study period: (1) April 8–21, 2010, (2) May 12–June 3, 2010, (3) June 15–July 1, 2010, (4) August 23–30, 2010, and (5) May 16–20, 2011. Streamflow conditions were variable for each sampling event because of input from local precipitation and irrigation return flow, and streamflow losses through the streambed. Streamflow and total nitrogen flux increased with drainage area, and the dominant form of nitrogen varied with drainage area size and temporally across sampling events. Results from each method indicate relatively conservative transport of nitrogen within the 160 miles between Clarksdale and Anguilla, providing further validation of the SPARROW models. Furthermore, these results suggest relatively conservative transport of nitrogen from the Big Sunflower River to the Gulf of Mexico and, therefore, imply a fairly close association of nutrient application and export from the Big Sunflower River Basin to the Mississippi River. However, within the Big Sunflower River Basin, two potential nitrogen sinks were identified and include the transport and potential transformation of nitrogen through the streambed and the sequestration and potential transformation of nitrogen above the drainage control structures downstream of Anguilla. By coupling these potential loss mechanisms with nitrogen transport dynamics, it may be possible to further reduce the amount of nitrogen leaving the Big Sunflower River Basin and ultimately arriving at the Gulf of Mexico.

Arkansas;Louisiana;Mississippi

Influences of environmental settings on aquatic ecosystems in the Apalachicola-Chattahoochee-Flint River basin

The watershed boundary of the Apalachicola-Chattahoochee-Flint (ACF) River basin defines an aquatic ecosystem whose water quality is the result of complex interactions of natural and human influences on land and water resources. Topics relating to the basin's environmental setting-its physical, biological, and cultural characteristics-are summarized to provide an understanding of factors that influence water quality and the health of aquatic ecosystems. The ACF River basin lies partly in southwestern Georgia, southeastern Alabama, and northwestern Florida and covers 19,800 square miles in the Blue Ridge, the Piedmont, and the Coastal Plain Provinces. The basin includes the drainages of the Chattahoochee River and the Flint River, which meet to form the Apalachicola River. The Apalachicola River flows into the Gulf of Mexico at Apalachicola Bay. Basin hydrology and water quality are influenced by 16 mainstem reservoirs, 13 of which are on the Chattahoochee River. Ground water in the basin is contained in six aquifers-the surficial aquifer system, the Floridan aquifer system, the Claiborne aquifer, the Clayton aquifer, the Providence aquifer, and the crystalline-rock aquifer. Physiography, climate, and hydrology of the ACF River basin provide natural conditions that support a rich and abundant diversity of plants and animals. Although most of the ACF River basin has been altered by human activities, the basin's environment is noteworthy for its remaining biological diversity and the role it plays in sustaining biological productivity in Apalachicola Bay. The Bay produces 90 percent of Florida's and 13 percent of the Nation's oyster harvest; and functions as a nursery for penaeid shrimp, blue crabs, and a variety of fin fish. The diversity of the basin's aquatic fauna is noteworthy because the basin is home to (1) the largest number of fish species among Gulf Coast drainages east of the Mississippi River, (2) the largest assemblage of freshwater fish in Florida, (3) the largest number of mollusc species among western Florida drainages, and (4) the highest species density of amphibians and reptiles on the continent north of Mexico. Population of the ACF River basin in 1990 was estimated at 2.6 million. Nearly 90 percent of the total population lived in Georgia, and nearly 60 percent lived in the Metropolitan Atlanta area. The 1990 basin population is projected to increase by 15 percent to 3.0 million by the year 2000, and by 30 percent to 3.4 million by 2010. The largest increases in populations are projected for the Metropolitan Atlanta area. In 1972-76, approximately 59 percent of the basin was covered by forest, 29 percent was agricultural, 5 percent was wetland, 4 percent was urban, and 3 percent was water or barren land. Most of the original land cover of the basin has been transformed by human activity. Timber is the basin's largest cash crop and most forests consist of second-growth stands or large acreages of planted pine. The dominant agricultural land use in the Piedmont Province is pasture and confined feeding for dairy, livestock, and poultry production. Row-crop agriculture, orchards, and silviculture are most common in the Coastal Plain Province. The top five crops in order from most to least acres harvested in 1990 were peanuts, corn, soybeans, wheat, and cotton. The water in the basin is used for public and industrial supply, irrigation, power generation, navigation, and recreation. Although most public-supply withdrawals in the Blue Ridge and Piedmont Provinces are from surface-water sources, with the exception of counties near or immediately below the Fall Line, all publicly supplied water in the Coastal Plain is withdrawn from ground-water sources. Ground water supplied 18 percent of the basin's population served by public supply. Total water withdrawn in the ACF River basin in 1990 was 2,098 million gallons per day (Mgal/d), of which Georgia withdrew 82 percent and Florida and Alabama each withdrew 9 percent. Power generation is the single largest water use. Sixteen of the basin's 22 power generating plants are located along the mainstem of the Chattahoochee River. The U.S. Army Corps of Engineers maintains a navigation channel from the mouth of the Apalachicola River to Columbus, Ga., on the Chattahoochee River and to Bainbridge, Ga., on the Flint River. Water quality in the basin is influenced by the operation of 137 municipal wastewater-treatment facilities. In 1990, 354 Mgal/d of municipal wastewater was discharged within the ACF River basin. Eighty-eight percent of the wastewater was discharged into the Chattahoochee River basin, 10.6 percent into the Flint River basin, and 1.4 percent into the Apalachicola River basin. Two-thirds of the 938 stream miles in the Georgia portion of the ACF River basin having water quality that does not meet or only partially meets the designated-use criteria in the Chattahoochee River basin. The Chattahoochee River is the most heavily-used water resource both in the ACF River basin and in Georgia. Urban runoff or unknown nonpoint sources are cited as the causes of water-quality regulations in 72 percent of violations. The remaining causes primarily are combined sewer overflows in the Atlanta area, and discharges from municipal or industrial treatment facilities with inadequate treatment capabilities or operational deficiencies.

Alabama, Florida, Georgia

Groundwater quality in the San Diego Drainages Hydrogeologic Province, California

More than 40 percent of California's drinking water is from groundwater. To protect this vital resource, the State of California created the Groundwater Ambient Monitoring and Assessment (GAMA) Program. The Priority Basin Project of the GAMA Program provides a comprehensive assessment of the State's groundwater quality and increases public access to groundwater-quality information. The San Diego Drainages Hydrogeologic Province (hereinafter referred to as San Diego) is one of the study units being evaluated. The San Diego study unit is approximately 3,900 square miles and consists of the Temecula Valley, Warner Valley, and 12 other alluvial basins (California Department of Water Resources, 2003). The study unit also consists of all areas outside defined groundwater basins that are within 3 kilometers of a public-supply well. The study unit was separated, based primarily on hydrogeologic settings, into four study areas: Temecula Valley, Warner Valley, Alluvial Basins, and Hard Rock (Wright and others, 2005). The sampling density for the Hard Rock study area, which consists of areas outside of groundwater basins, was much lower than for the other study areas. Consequently, aquifer proportions for the Hard Rock study area are not used to calculate the aquifer proportions shown by the pie charts. An assessment of groundwater quality for the Hard Rock study area can be found in Wright and Belitz, 2011. The temperatures in the coastal part of the study unit are mild with dry summers, moist winters, and an average annual rainfall of about 10 inches. The temperatures in the mountainous eastern part of the study unit are cooler than in the coastal part, with an annual precipitation of about 45 inches that occurs mostly in the winter. The primary aquifers consist of Quaternary-age alluvium and weathered bedrock in the Temecula Valley, Warner Valley, and Alluvial Basins study areas, whereas in the Hard Rock study area the primary aquifers consist mainly of fractured and decomposed granite of Mesozoic age. The primary aquifers are defined as those parts of the aquifers corresponding to the perforated intervals of wells listed in the California Department of Public Health (CDPH) database. Public-supply wells typically are drilled to depths between 200 and 700 feet, consist of solid casing from the land surface to a depth of about 60 to 170 feet, and are perforated, or consist of an open hole, below the solid casing. Water quality in the shallow and deep parts of the aquifer system may differ from water quality in the primary aquifers. Municipal water use accounts for approximately 70 percent of water used in the study unit; the majority of the remainder is used for agriculture, industry, and commerce. Groundwater accounts for approximately 8 percent of the municipal supply, and surface water, the majority of which is imported, accounts for the rest. Recharge to groundwater occurs through stream-channel infiltration from rivers and their tributaries, infiltration in engineered recharge basins, and infiltration of water from precipitation and irrigation. The primary source of discharge is water pumped from wells.

California

Reconnaissance investigation of water quality, bottom sediment, and biota associated with irrigation drainage in the Owyhee and Vale Projects, Oregon and Idaho, 1990-91

A reconnaissance investigation was conducted during 1990--91 in the Owyhee and Vale projects in eastern Oregon and southwestern Idaho, as well as at a number of sites in the Snake River and tributaries to the Snake River in the area of study. The objective of the study was to determine if agricultural drainwater entering the study area was causing, or had the potential to cause, significant harmful effects to human health, fish and wildlife, or may adversely affect the suitability of water for beneficial uses. Approximately 153,000 acres of land are irrigated annually within the areas of the Owyhee and Vale projects. Large quantities of water are required because of the semiarid climate and relatively high evaporation rates. Several reservoirs in the area are filled annually during the wet, nonirrigation season to sustain irrigation during the dry summer months. During the irrigation season, this impounded water, along with direct diversions from the Malheur, Owyhee, and Snake Rivers, is transported to the irrigated areas through a series of diversion tunnels, siphons, canals, aqueducts, ditches, and drains. Major crops grown in the area include sugar beets, alfalfa hay and other hay crops, onions, and winter wheat. Minor crops include corn, potatoes, mint, various seed crops, and fruit. In 1987, it was estimated that the following amounts of pesticides were used in the project areas: 2,4-D (21,000 lbs [pounds]), chlorpyrifos (1,000 lbs), dacthal (40,000 lbs), dicamba (320 lbs), endosulfan (2,500 lbs), ethion (11,000 lbs), malathion (24,000 lbs), parathion (5,000 lbs), and phorate (11,000 lbs). Median concentrations and values for total dissolved solids, alkalinity, sodium adsorption ratio, and hardness in the Vale project area were greater than 1.5 times those values observed in the Owyhee project area or at other Snake River locations. During irrigation (August 1990), total dissolved solids, alkalinity, sodium adsorption ratio, and hardness values increased in a downstream manner. Constituent values at drainwater sites generally were comparable to concentrations below the irrigated ureas in the Owyhee and Vale project areas. The trace elements arsenic, boron, copper, molybdenum, vanadium, and zinc were detected in most water samples; cadmium, chromium, lead, and selenium were detected in some samples at concentrations generally near the analytical reporting limit; mercury was not detected in any samples. In some water samples, concentrations of arsenic, boron, cadmium, copper, and lead exceeded State or Federal water-quality standards or criteria. Most trace elements in bottom sediment were detected at concentrations within the expected 95-percent baseline range for soils from the Western United States. Concentrations that exceeded the 95-percent baseline range for study area soils were: (1) arsenic and lead in one sample from a site in the Vale project area; (2) mercury, lead, and tin in one sample from a site in the Snake River system; (3) manganese in two samples from two sites in the Snake River system; and (4) manganese from one sample from a site in the Vale project area. Fifteen pesticides and metabolites were detected in whole-water samples collected from sites in the study area. DDT, plus its metabolites (DDE and DDD), dieldrin, endrin, 2,4-D, dicamba, and dacthal were detected in samples collected from seven or more sites. Other pesticides detected included chlorpyrifos, endosulfan, ethion, malathion, parathion, phorate, and lindane. Most of the detected pesticide concentrations generally were largest in drainwater and at the most downstream sampling locations in the Owyhee and Vale project areas. Concentrations exceeded water-quality criteria established for the protection of freshwater aquatic life in 86 percent of the whole-water samples analyzed for DDT plus its metabolites, 71 percent of the dieldrin samples, 14 percent of the endrin samples, and 10 percent of the parathion samples. Eight pesticides and

Water-Resources Investigations Report

Estimates of water use and trends in the Colorado River Basin, Southwestern United States, 1985–2010

The Colorado River Basin (CRB) drains 246,000 square miles and includes parts of California, Colorado, Nevada, New Mexico, Utah, and Wyoming, and all of Arizona (Basin States). This report contains water-use estimates by category of use for drainage basins (Hydrologic Unit Code 8; HUC‑8) within the CRB from 1985 to 2010, at 5-year intervals. Estimates for public supply, domestic, commercial, industrial, irrigation, livestock, mining, aquaculture, hydroelectric and thermoelectric power, and wastewater returns are tabulated as (1) water withdrawals from groundwater or surface‑water sources of fresh or saline quality, (2) water delivered for domestic use, (3) wastewater returns and instream use (hydroelectric), and (4) consumptive use, or water that is consumed (USGS definition) and not available for immediate reuse. Water transported outside of the CRB (interbasin transfers) is not included as part of withdrawals and are not accounted for in any category of use within the CRB. Total withdrawals in the CRB (excluding interbasin transfers) averaged about 17 million acre-feet (maf) from 1985 to 2010, peaked at about 17.76 maf in 2000, and reached their lowest levels of 16.43 maf in 1990. Interbasin transfers to serve mostly public-supply and irrigation needs outside of the CRB are reported for 2000, 2005, and 2010 only, and averaged 5.40 maf. More surface water was used in the CRB than groundwater, averaging about 78 percent of total withdrawals, and its use increased less than 2 percent from 1985 to 2010, while groundwater withdrawals decreased about 12 percent. From 1985 to 2010, surface water averaged 98 percent of withdrawals in the upper CRB, and about 59 percent in the lower CRB. Nearly all withdrawals were freshwater, but some saline groundwater was used for mining and self-supplied industrial. Interbasin transfers have a large effect on flows in the Colorado River and are listed in this report separately with no explanation of how the water is used outside of the CRB. There are 34 interbasin transfers that conveyed an estimated 5.83, 5.20, and 5.18 maf out of the CRB in 2000, 2005, and 2010, respectively. The largest interbasin transfers are in the lower CRB and convey surface water (Colorado River water) to southern California; these accounted for 80 to 84 percent of total interbasin transfers in the CRB from 2000 to 2010. Intrabasin transfers are conveyances of surface water that cross drainage basin or State boundaries in the CRB, but the water does not leave the CRB. There are many intrabasin transfers in the CRB, but this report lists 11 that are mostly in the State of Colorado. The largest is the Central Arizona Project (CAP), through which more than 1.00 maf of water was provided to irrigate nearly 1 million acres in Maricopa, Pinal, and Pima Counties, as well as provide municipal water for Phoenix and Tucson, Arizona, during 2000, 2005, and 2010. In 2010, interbasin and intrabasin transfers accounted for 24 and 11 percent of the total water withdrawals in CRB, respectively, with the larger volumes being conveyed out of the lower CRB. Total population in the CRB increased from 4.56 to 9.44 million people from 1985 to 2010. Most of those people were in the lower CRB, with 86 percent of the total in 1985, and 90 percent of the total in 2010. Total public-supply withdrawals in the CRB provided most people with their potable water, and averaged about 1.63 maf from 1985 to 2010, ranging from about 1.07 maf in 1985 to about 2.10 maf in 2010, when it peaked. Most of public-supply withdrawals occurred in the lower CRB, ranging from 87 to 91 percent of total public-supply withdrawals in the CRB over the 25 years. Total domestic use, comprised of public-supply deliveries and self-supply domestic withdrawals, increased more than 90 percent from 1985 to 2010, from about 0.80 maf to about 1.54 maf. Domestic daily per-capita use rates in the CRB ranged from about 144 (1985) to about 121 (2000) gallons (gal) per ‑ capita between 1985 and 2010. When comparing domestic daily per-capita rates for the upper and lower CRB, people in the lower CRB, on average, used less water for domestic purposes at 128 gal per-capita daily (1985–2010), while those in the upper CRB for the same time period averaged 133 gal per-capita daily. The trend in daily per-capita use rates for the entire CRB fluctuated between the reporting years, but decreased overall, indicating that more people used less water in 2010 than in 1985, likely due to improved infrastructure, conservation, and improvements to water using appliances in homes and businesses. Irrigation accounted for most total withdrawals in the CRB, excluding instream use for hydroelectric power and interbasin transfers, averaging 85 percent from 1985 to 2010. Far more surface water than groundwater was used for irrigation in both the upper and lower CRB, but in the upper CRB, it accounted for an average of more than 98 percent of the total withdrawals (1985–2010), whereas in the lower CRB, surface-water withdrawals for irrigation averaged 61 percent of total withdrawals. On average, the upper CRB accounted for 56 percent of total irrigated acres, and the irrigation systems in the upper CRB trended towards more efficient sprinkler systems from 1985 to 2010. Long-term drought in the CRB substantially decreased the amount of streamflow available for irrigation. Increases in micro-irrigation acres, which can have efficiencies that exceed 90 percent and require 20–50 percent less water than sprinkler systems, likely contributed to reduced withdrawals in the lower CRB. For thermoelectric power, total withdrawals, including the use of reclaimed wastewater, were greater in the upper CRB from 1985 through 2005. In 2010, the lower CRB exceeded the upper by only 11,000 acre-feet. On average, thermoelectric consumptive use accounted for about 80 percent of the total withdrawals; however, consumptive-use data in the upper CRB was incomplete. Surface water was the primary source in the upper CRB and groundwater was the primary source in the lower CRB. In the CRB overall, water withdrawals for thermoelectric generation has decreased since 2000, except for groundwater withdrawals in the lower CRB. Power generation at thermoelectric plants was greater in the upper CRB from 1985 to 2000, and after 2005 the difference in power generation was small; however, the upper CRB continued to have more power generation. In both the upper and lower CRB, power generation increased from 1985 to 2005.

Arizona, Colorado, Nevada, New Mexico, Utah, Wyomi

Geohydrology and water utilization in the Willcox Basin, Graham and Cochise Counties, Arizona

The Willcox basin is an area of interior drainage in the northern part of Sulphur Springs Valley, Cochise and Graham Counties, Ariz. The basin comprises about 1,500 square miles, of which the valley floor occupies about 950 square miles. The basin probably formed during middle and late Tertiary time, when the area was subjected to large-scale faulting accompanied by the uplift of the mountain ranges that presently border it. During and after faulting, large quantities of alluvium were deposited in the closed basin. The rocks in the basin are divided into two broad groups--the rocks of the mountain blocks, of Precambrian through Tertiary age, and the rocks of the basin, of Tertiary and Quaternary age. The mountain blocks consist of igneous, metamorphic, and sedimentary rocks; the water-bearing characteristics of these rocks depend primarily on their degree of weathering and fracturing. Even in areas where these rocks are fractured and jointed, only small amounts of water have been developed. The rocks of the basin consist of moderately consolidated alluvium, poorly consolidated alluvium, and unconsolidated alluvium. The water-bearing characteristics of the moderately and poorly consolidated alluvium are not well known. The unconsolidated alluvium underlies most of the valley floor and consists of two facies, stream deposits and lake beds associated with the old playa. The lenticular sand and gravel layers interbedded in silt- and clay-size material of the unconsolidated alluvium constitute the principal aquifer in the basin. The other aquifers, which yield less water, consist of beds of poorly to moderately consolidated sand- and gravel-size material; these beds occur in both the poorly consolidated and moderately consolidated alluvium. In the Stewart area the median specific capacity of wells per 100 feet of saturated unconsolidated alluvium was 20 gallons per minute, and in the Kansas Settlement area the specific capacity of wells penetrating the poorly and moderately consolidated alluvium, undifferentiated, was only 7.4 gallons per minute per 100 feet of saturated material penetrated. The aquifer in the Kansas Settlement area is much less permeable but more homogeneous than the aquifer in the Stewart area. The coefficient of transmissibility of the aquifers, which was estimated from the specific-capacity data, ranged from 58,000 to 160,000 gal. tons per day per foot. Prior to extensive ground-water pumpage, the ground-water system probably was in equilibrium, with discharge equaling recharge. At that time, ground water moved toward the playa, where it was discharged by transpiration and evaporation. The estimate of the evapotranspiration in the playa area before large-scale development was about 75,000 acre-feet per year. On the basis of estimates of coefficients of transmissibility of the aquifer and on the basis of the water-table configuration, underflow toward the playa was computed to be about 54,000 acre-feet per year. By 1963, large-scale pumping had caused marked changes in the shape of the piezometric surface; large cones of depression had developed, and ground-water movement was toward the centers of pumping. The cones of depression caused by large-scale pumping have since expanded, and water-level declines have been measured in the recharge areas along the mountain fronts. Ground water has been used for irrigation since 1910. In 1928, about 4,000 acre-feet of ground water was pumped, and by 1963 180,000 acre-feet per year was being pumped. An estimated 1,860,000 acre-feet of water has been pumped for irrigation in the Willcox basin through 1963; 680,000 acre-feet from the Stewart area, 990,000 acre-feet from the Kansas Settlement area, and 190,000 acre-feet from the Pearce-Cochise area. In the Sierra Bonita Ranch area and the north playa area, ground-water withdrawal for irrigation through 1963 was small. From the spring of 1952 to the spring of 1964 water-level declines resulting from the

Water Supply Paper

Selenium loading through the Blackfoot River watershed--linking sources to ecosystem

The upper Blackfoot River watershed in southeast Idaho receives drainage from 11 of 16 phosphate mines that have extracted ore from the Phosphoria Formation, three of which are presently active. Toxic effects from selenium (Se), including death of livestock and deformity in aquatic birds, were documented locally in areas where phosphatic shales are exposed ( Piper et al., 2000 ; Presser et al., Chapter 11). Current drainage conditions are leading to Se bioaccumulation at concentrations that pose a risk to fish in the Blackfoot River and its tributaries (Hamilton et al., Chapter 18). A gaging station on the Blackfoot River was re-activated in April 2001 to assess hydrologic conditions and concentration, load, and speciation for Se discharges on a watershed scale. The gaging-station data are considered to represent regional drainage conditions in the upper Blackfoot River water- shed because of its location near the outlet of the watershed and directly upstream of the Blackfoot Reservoir. Watershed discharges for 2001 and 2002 were below minimum hydrologic conditions for the gage as documented by the historical record. Drought emergencies were declared in the area in both 2001 and 2002. Unmonitored diversions for irrigation that routinely take place during the snowmelt season also affected conditions downstream. Annual cycles in Se concentration, load, and selenate (Se 6+ ) reached maxima in the spring during the period of maximum flow at the gaging station. Thirty-seven to 44% of annual flow occurred dur- ing the three-month high-flow season (April through June) in 2001 and 56% of annual flow occurred during that time period in 2002. Extrapolation from historical hydrographs for average and wet years and a limited data set of regional Se concentrations for 2001 and 2002 indicated potential for a 3.6- to 7.4-fold increase in Se loading because of increased seasonal flows in the Blackfoot River watershed. Supplementation data indicate that: (a) the difference between total Se and dissolved Se, as a measure of the contribution of particulate Se, was < 10% except at the peak of con- centration when total Se was 18% more than dissolved Se; (b) selenite (Se 4+ ) represented less than 10% of the dissolved species during all months of 2001; and (c) dissolved Se was approximately a 50:50 mixture of selenate and organic selenide (operationally defined Se 2- ) during summer 2001 (June through August). Ecological risk based on regional Se drainage occurred during both the high- and low-flow seasons. Seventy to 83% of the Se load occurred during the high-flow season. During early May of both years, dissolved-Se concentrations exceeded the criterion for the protection of aquatic life and the ecological threshold of 5 gL 1 Se at which sub- stantive risk occurs. During the majority of the three-month high-flow season, dissolved- Se concentrations exceeded the 2 gL 1 Se concern level for aquatic biota. The Se concentration in suspended material during high flow in 2002 was within the range of marginal risk to aquatic life (2-4 gg 1 Se, dry weight). Selenate was the major species during peak flows, with both selenate and organic selenide being major species during relatively low-flow periods in summer. A change in speciation to reduced Se may indicate elevated biotic productivity during summer months and could result in enhanced Se uptake in food webs. In addition to the magnitude of regional Se release in the Blackfoot River watershed, Se concentrations in individual source drains and waste-rock seeps, and those predicted by experimental column leaching of proposed mining overburden materials, also indicate that drainage options that currently meet existing demands for phosphate mining cause eco- logical risk thresholds to be exceeded. At times, the drinking-water Se standard (50 g L 1 Se) and the criterion for hazardous Se waste (1000 L -1 Se) ( US Department of the Interior, 1998 ; US Environmental Protection Agency, 1987 ) are also exceeded. For water-years 2001 and 2002, seasonal increased input of water in the mining area resulted in increased Se transport, suggesting a mechanism of contamination that involves a significant Se reservoir. Hence, recognition and monitoring of Se loading to the envi- ronment on a mass balance basis (i.e. inputs, fluxes and storage within environmental media, and outputs) are essential to evaluating how to control Se concentrations within environmentally protective ranges ( Presser and Piper, 1998 ). In areas where release of Se to aquatic systems is anticipated as a product of future expansion of phosphate mining, continuous monitoring of flow and development of seasonal Se loading patterns would help to model watersheds in terms of sources, flow periods, and environmental-Se con- centrations that most influence bioavailability. These data, in turn, could be linked to Se- bioaccumulation models specific to food webs and vulnerable species of the impacted areas to accurately project ecological effects. Gaging at this site on the Blackfoot River is planned to continue in order to establish a long-term (>10 year) record of hydrologic conditions.

Idaho

Estimated dissolved-solids loads and trends at selected streams in and near the Uinta Basin, Utah, Water Years 1989–2013

The U.S. Geological Survey (USGS), in cooperation with the Colorado River Basin Salinity Control Forum, studied trends in dissolved-solids loads at selected sites in and near the Uinta Basin, Utah. The Uinta Basin study area includes the Duchesne River Basin and the Middle Green River Basin in Utah from below Flaming Gorge Reservoir to the town of Green River. Annual dissolved-solids loads for water years (WY) 1989 through 2013 were estimated for 16 gaging stations in the study area using streamflow and water-quality data from the USGS National Water Information System database. Eight gaging stations that monitored catchments with limited or no agricultural land use (natural subbasins) were used to assess loads from natural sources. Four gaging stations that monitored catchments with agricultural land in the Duchesne River Basin were used to assess loads from agricultural sources. Four other gaging stations were included in the dissolved-solids load and trend analysis to help assess the effects of agricultural areas that drain to the Green River in the Uinta Basin, but outside of the Duchesne River Basin. Estimated mean annual dissolved-solids loads for WY 1989–2013 ranged from 1,520 tons at Lake Fork River above Moon Lake, near Mountain Home, Utah (UT), to 1,760,000 tons at Green River near Green River, UT. The flow-normalized loads at gaging stations upstream of agricultural activities showed no trend or a relatively small change. The largest net change in modeled flow-normalized load was -352,000 tons (a 17.8-percent decrease) at Green River near Green River, UT. Annual streamflow and modeled dissolved-solids loads at the gaging stations were balanced between upstream and downstream sites to determine how much water and dissolved solids were transported to the Duchesne River and a section of the Green River, and how much was picked up in each drainage area. Mass-balance calculations of WY 1989–2013 mean annual dissolved-solids loads at the studied sites show that Green River near Jensen, UT, accounts for 64 percent of the load in the river at Green River, UT, while the Duchesne River and White River contribute 10 and 13 percent, respectively. Annual streamflow and modeled dissolved-solids loads at the gaging stations were balanced between upstream and downstream sites to determine how much water and dissolved solids were transported to the Duchesne River and a section of the Green River, and how much was picked up in each drainage area. Mass-balance calculations of WY 1989–2013 mean annual dissolved-solids loads at the studied sites show that Green River near Jensen, UT, accounts for 64 percent of the load in the river at Green River, UT, while the Duchesne River and White River contribute 10 and 13 percent, respectively. The flow-normalized dissolved-solids loads estimated at Duchesne River near Randlett, UT, and White River near Watson, UT, decreased by 68,000 and 55,300 tons, or 27.8 and 20.8 percent respectively, when comparing 1989 to 2013. The drainage basins for both rivers have undergone salinity-control projects since the early 1980s to reduce the dissolved-solids load entering the Colorado River. Approximately 19 percent of the net change in flow-normalized load at Green River at Green River, UT, is from changes in load modeled at Duchesne River near Randlett, UT, and 16 percent from changes in load modeled at White River near Watson, UT. The net change in flow-normalized load estimated at Green River near Greendale, UT, for WY 1989–2013 accounts for about 45 percent of the net change estimated at Green River at Green River, UT. Mass-balance calculations of WY 1989–2013 mean annual dissolved-solids loads at the studied sites in the Duchesne River Basin show that 75,400 tons or 44 percent of the load at the Duchesne River near Randlett, UT, gaging station was not accounted for at any of the upstream gages. Most of this unmonitored load is derived from tributary inflow, groundwater discharge, unconsumed irrigation water, and irrigation tail water. A mass balance of WY 1989–2013 flow-normalized loads estimated at sites in the Duchesne River Basin indicates that the flow-normalized load of unmonitored inflow to the Duchesne River between the Myton and Randlett gaging stations decreased by 38 percent. The total net decrease in flow-normalized load calculated for unmonitored inflow in the drainage basin accounts for 94 percent of the decrease in WY 1989–2013 flow-normalized load modeled at the Duchesne River near Randlett, UT, gaging station. Irrigation improvements in the drainage basin have likely contributed to the decrease in flow-normalized load. Reductions in dissolved-solids load estimated by the Natural Resources Conservation Service (NRCS) and the Bureau of Reclamation (Reclamation) from on- and off-farm improvements in the Uinta Basin totaled about 135,000 tons in 2013 (81,900 tons from on-farm improvements and 53,300 tons from off-farm improvements). The reduction in dissolved-solids load resulting from on- and off-farm improvements facilitated by the NRCS and Reclamation in the Price River Basin from 1989 to 2013 was estimated to be 64,800 tons. The amount of sprinkler-irrigated land mapped in the drainage area or subbasin area for a gaging station was used to estimate the reduction in load resulting from the conversion from flood to sprinkler irrigation. Sprinkler-irrigated land mapped in the Uinta Basin totaled 109,630 acres in 2012. Assuming conversion to wheel-line sprinklers, a reduction in dissolved-solids load in the Uinta Basin of 95,800 tons in 2012 was calculated using the sprinkler-irrigation acreage and a pre-salinity-control project dissolved-solids yield of 1.04 tons per acre. A reduction of 72,800 tons in dissolved-solids load from irrigation improvements was determined from sprinkler-irrigated lands in the Ashley Valley and Jensen, Pelican Lake, and Pleasant Valley areas (mapped in 2012); and in the Price River Basin (mapped in 2011). This decrease in dissolved-solids load is 8,800 tons more than the decrease in unmonitored flow-normalized dissolved-solids load (-64,000 tons) determined for the Green River between the Jensen and Green River gaging stations. The net WY 1989–2013 change in flow-normalized dissolved-solids load at the Duchesne River near Randlett, UT, and the Green River between the Jensen and Green River, UT, gaging stations determined from mass-balance calculations was compared to reported reductions in dissolved-solids load from on- and off-farm improvements and estimated reductions in load determined from mapped sprinkler-irrigated areas in the Duchesne River Basin and the area draining to the Green River between the Jensen and Green River gaging stations. The combined NRCS and Reclamation estimates of reduction in dissolved-solids load from on- and off-farm improvements in the study area (200,000 tons) is more than the reduction in load estimated using the acreage with sprinkler improvements (136,000 tons) or the mass-balance of flow-normalized load (132,000 tons).

Utah

Simulated groundwater flow in the Ogallala and Arikaree aquifers, Rosebud Indian Reservation area, South Dakota – Revisions with data through water year 2008 and simulations of potential future scenarios

The Ogallala and Arikaree aquifers are important water resources in the Rosebud Indian Reservation area and are used extensively for irrigation, municipal, and domestic water supplies. Drought or increased withdrawals from the Ogallala and Arikaree aquifers in the Rosebud Indian Reservation area have the potential to affect water levels in these aquifers. This report documents revisions and recalibration of a previously published three-dimensional, numerical groundwater-flow model for this area. Data for a 30-year period (water years 1979 through 2008) were used in steady-state and transient numerical simulations of groundwater flow. In the revised model, revisions include (1) extension of the transient calibration period by 10 years, (2) the use of inverse modeling for steady-state calibration, (3) model calibration to base flow for an additional four surface-water drainage basins, (4) improved estimation of transient aquifer recharge, (5) improved delineation of vegetation types, and (6) reduced cell size near large capacity water-supply wells. In addition, potential future scenarios were simulated to assess the potential effects of drought and increased groundwater withdrawals. The model comprised two layers: the upper layer represented the Ogallala aquifer and the lower layer represented the Arikaree aquifer. The model’s grid had 168 rows and 202 columns, most of which were 1,640 feet (500 meters) wide, with narrower rows and columns near large water-supply wells. Recharge to the Ogallala and Arikaree aquifers occurs from precipitation on the outcrop areas. The average recharge rates used for the steady-state simulation were 2.91 and 1.45 inches per year for the Ogallala aquifer and Arikaree aquifer, respectively, for a total rate of 255.4 cubic feet per second (ft 3 /s). Discharge from the aquifers occurs through evapotranspiration, discharge to streams as base flow and spring flow, and well withdrawals. Discharge rates for the steady-state simulation were 171.3 ft 3 /s for evapotranspiration, 74.4 ft 3 /s for net outflow to streams and springs, and 11.6 ft 3 /s for well withdrawals. Estimated horizontal hydraulic conductivity used for the numerical model ranged from 0.2 to 84.4 feet per day (ft/d) in the Ogallala aquifer and from 0.1 to 4.3 ft/d in the Arikaree aquifer. A uniform vertical hydraulic conductivity value of 4.2x10 -4 ft/d was estimated for the Ogallala aquifer. Vertical hydraulic conductivity was estimated for five zones in the Arikaree aquifer and ranged from 8.8x10 -5 to 3.7 ft/d. Average rates of recharge, maximum evapotranspiration, and well withdrawals were included in the steady-state simulation, whereas the time-varying rates were included in the transient simulation. Inverse modeling techniques were used for steady-state model calibration. These methods were designed to estimate parameter values that are, statistically, the most likely set of values to result in the smallest differences between simulated and observed hydraulic heads and base-flow discharges. For the steady-state simulation, the root mean square error for simulated hydraulic heads for all 383 wells was 27.3 feet. Simulated hydraulic heads were within ±50 feet of observed values for 93 percent of the wells. The potentiometric surfaces of the two aquifers calculated by the steady-state simulation established initial conditions for the transient simulation. For the transient simulation, the difference between the simulated and observed means for hydrographs was within ±40 feet for 98 percent of 44 observation wells. A sensitivity analysis was used to examine the response of the calibrated steady-state model to changes in model parameters including horizontal and vertical hydraulic conductivity, evapotranspiration, recharge, and riverbed conductance. The model was most sensitive to recharge and maximum evapotranspiration and least sensitive to riverbed and spring conductances. To simulate a potential future drought scenario, a synthetic recharge record was created, the mean of which was equal to 64 percent of the average estimated recharge rate for the 30-year calibration period. This synthetic recharge record was used to simulate the last 20 years of the calibration period under drought conditions. Compared with results of the calibrated model, decreases in hydraulic-head values for the drought scenario at the end of the simulation period were as much as 39 feet for the Ogallala aquifer. To simulate the effects of potential increases in pumping, well withdrawal rates were increased by 50 percent from those estimated for the 30-year calibration period for the last 20 years of the calibration period. Compared with results of the calibrated model, decreases in hydraulic-head values for the scenario of increased pumping at the end of the simulation period were as much as 13 feet for the Ogallala aquifer. This numerical model is suitable as a tool to help understand the flow system, to help confirm that previous estimates of aquifer properties were reasonable, and to estimate aquifer properties in areas without data. The model also is useful to help assess the effects of drought and increases in pumping by simulations of these scenarios, the results of which are not precise but may be considered when making water management decisions.

South Dakota