USGS ScienceSearch

SEARCH · USGS Science

Results for “Water Resources Management”

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

History of irrigation and characteristics of streamflow in Nebraska, part of the North and South Platte River basins

More than 80 years of streamflow records have been collected in Nebraska. Such records are essential in planning water developments and managing water resources. Their value can be enhanced if the records are analyzed statistically. This report presents statistics on the flow of streams in the North and South Platte River basins. Major irrigation projects have been constructed in the basins since gaging of streams was begun, and the records of stream discharge at gaging stations show the effects on streamflow of those developments. Streamflow records presented in this report were computed from data that represent 1973 conditions. All reservoirs on the North Platte River above Lake McConaughy are in Wyoming. The operation of those reservoirs and irrigation projects has a direct effect on the characteristics of streamflow in the Nebraska part of the basin. A brief description of the present river system is presented.

Nebraska

Water quality and bathymetry of Sand Lake, Anchorage, Alaska

Sand Lake, a dimictic lowland lake in Anchorage, Alaska, has recently become as urban lake. Analyses indicate that the lake is oligotrophic, having low dissolved solids and nutrient concentrations. Snowmelt runoff from an adjacent residential area, however, has a dissolved-solids concentration 10 times that of the main body of Sand Lake. Lead concentrations in the runoff exceed known values from other water in the ANchorage area, including water samples taken beneath landfills. The volume of the snowmelt runoff has not been measured. The data presented can be used as a baseline for water-resource management. (Woodard-USGS)

Open-File Report

Equations for estimating synthetic unit-hydrograph parameter values for small watersheds in Lake County, Illinois

Design hydrographs computed from design storms, simple models of abstractions (interception, depression storage, and infiltration), and synthetic unit hydrographs provide vital information for stormwater, flood-plain, and water-resources management throughout the United States. Rainfall and runoff data for small watersheds in Lake County collected between 1990 and 1995 were studied to develop equations for estimation of synthetic unit-hydrograph parameters on the basis of watershed and storm characteristics. The synthetic unit-hydrograph parameters of interest were the time of concentration (TC) and watershed-storage coefficient (R) for the Clark unit-hydrograph method, the unit-graph lag (UL) for the Soil Conservation Service (now known as the Natural Resources Conservation Service) dimensionless unit hydrograph, and the hydrograph-time lag (TL) for the linear-reservoir method for unit-hydrograph estimation. Data from 66 storms with effective-precipitation depths greater than 0.4 inches on 9 small watersheds (areas between 0.06 and 37 square miles (mi2)) were utilized to develop the estimation equations, and data from 11 storms on 8 of these watersheds were utilized to verify (test) the estimation equations. The synthetic unit-hydrograph parameters were determined by calibration using the U.S. Army Corps of Engineers Flood Hydrograph Package HEC-1 (TC, R, and UL) or by manual analysis of the rainfall and run-off data (TL). The relation between synthetic unit-hydrograph parameters, and watershed and storm characteristics was determined by multiple linear regression of the logarithms of the parameters and characteristics. Separate sets of equations were developed with watershed area and main channel length as the starting parameters. Percentage of impervious cover, main channel slope, and depth of effective precipitation also were identified as important characteristics for estimation of synthetic unit-hydrograph parameters. The estimation equations utilizing area had multiple correlation coefficients of 0.873, 0.961, 0.968, and 0.963 for TC, R, UL, and TL, respectively, and the estimation equations utilizing main channel length had multiple correlation coefficients of 0.845, 0.957, 0.961, and 0.963 for TC, R, UL, and TL, respectively. Simulation of the measured hydrographs for the verification storms utilizing TC and R obtained from the estimation equations yielded good results without calibration. The peak discharge for 8 of the 11 storms was estimated within 25 percent and the time-to-peak discharge for 10 of the 11 storms was estimated within 20 percent. Thus, application of the estimation equations to determine synthetic unit-hydrograph parameters for design-storm simulation may result in reliable design hydrographs; as long as the physical characteristics of the watersheds under consideration are within the range of those for the watersheds in this study (area: 0.06-37 mi2, main channel length: 0.33-16.6 miles, main channel slope: 3.13-55.3 feet per mile, and percentage of impervious cover: 7.32-40.6 percent). The estimation equations are most reliable when applied to watersheds with areas less than 25 mi2.

Open-File Report

Water supply in the conterminous United States, Alaska, Hawaii, and Puerto Rico, water years 2010–20

We present an assessment of water supply across the conterminous United States (CONUS), Alaska, Hawaii, and Puerto Rico covering water years 2010–20. Our analysis drew on two national hydrologic models, the National Hydrologic Model Precipitation-Runoff Modeling System and the Weather Research and Forecasting model hydrologic modeling system. Both models produced estimates of streamflow, evapotranspiration, soil moisture, snow water equivalent, and other hydrologic states and fluxes. The models were driven by the bias-adjusted 4-kilometer-resolution, long-term regional hydroclimate simulation over the conterminous United States dataset (CONUS404). We assessed spatial and temporal error distributions by comparing monthly simulations at the 12-digit hydrologic unit code and regional scale from both models against external benchmarking datasets. Results showed that average annual rainfall across the CONUS was 857 millimeters per year for the period of analysis, with water year 2012 the driest year (729 millimeters) and water year 2019 the wettest year (995 millimeters). Key interannual variability results included the following: (1) the California–Nevada hydrologic region had the highest variability in precipitation and snow accumulation, and (2) the Texas hydrologic region was among hydrologic regions with the highest variability in precipitation. We related interannual variability in precipitation to storage volumes in soil moisture, snow water equivalent, and lakes and reservoirs to highlight areas with little storage and large year-to-year variability in precipitation. These areas included the Southern High Plains, Central High Plains, Texas, Souris–Red–Rainy, Mississippi Embayment, and Midwest regions. Our analysis of groundwater-level data showed that several of these areas overlap aquifers where groundwater levels were considerably lower than historical averages, including the Colorado Plateaus aquifers, the Rio Grande aquifer system, and the Central and Southern regions of the High Plains aquifer. Many of these lowered groundwater levels are continuations of decades-long declines from overpumping that started well before the assessment period. The resulting water budgets and their analyses provide a high-resolution foundational assessment of the mean state and variability of the terrestrial hydrologic cycle across the CONUS and Alaska, Hawaii, and Puerto Rico to support a wide range of water resource management applications.

Professional Paper

Hydrogeologic and geospatial data for the assessment of focused recharge to the carbonate-rock Aquifer in Genesee County, New York

Existing hydrogeologic and geospatial data useful for the assessment of focused recharge to the carbonate-rock aquifer in the central part of Genesee County, NY, were compiled from numerous local, State, and Federal agency sources. Data sources utilized in this pilot study include available geospatial datasets from Federal and State agencies, interviews with local highway departments and the Genesee County Soil and Water Conservation District, and an initial assessment of karst features through the analysis of ortho-photographs, with minimal field verification. The compiled information is presented in a series of county-wide and quadrangle maps. The county-wide maps present generalized hydrogeologic conditions including distribution of geologic units, major faults, and karst features, and bedrock-surface and water-table configurations. Ten sets of quadrangle maps of the area that overlies the carbonate-rock aquifer present more detailed and additional information including distribution of bedrock outcrops, thin and (or) permeable soils, and karst features such as sinkholes and swallets. Water-resource managers can utilize the information summarized in this report as a guide to their assessment of focused recharge to, and the potential for surface contaminants to reach the carbonate-rock aquifer.

New York

Bathymetry and capacity of Shawnee Reservoir, Oklahoma, 2016

Shawnee Reservoir (locally known as Shawnee Twin Lakes) is a man-made reservoir on South Deer Creek with a drainage area of 32.7 square miles in Pottawatomie County, Oklahoma. The reservoir consists of two lakes connected by an equilibrium channel. The southern lake (Shawnee City Lake Number 1) was impounded in 1935, and the northern lake (Shawnee City Lake Number 2) was impounded in 1960. Shawnee Reservoir serves as a municipal water supply, and water is transferred about 9 miles by gravity to a water treatment plant in Shawnee, Oklahoma. Secondary uses of the reservoir are for recreation, fish and wildlife habitat, and flood control. Shawnee Reservoir has a normal-pool elevation of 1,069.0 feet (ft) above North American Vertical Datum of 1988 (NAVD 88). The auxiliary spillway, which defines the flood-pool elevation, is at an elevation of 1,075.0 ft. The U.S. Geological Survey (USGS), in cooperation with the City of Shawnee, has operated a real-time stage (water-surface elevation) gage (USGS station 07241600) at Shawnee Reservoir since 2006. For the period of record ending in 2016, this gage recorded a maximum stage of 1,078.1 ft on May 24, 2015, and a minimum stage of 1,059.1 ft on April 10–11, 2007. This gage did not report reservoir storage prior to this report (2016) because a sufficiently detailed and thoroughly documented bathymetric (reservoir-bottom elevation) survey and corresponding stage-storage relation had not been published. A 2011 bathymetric survey with contours delineated at 5-foot intervals was published in Oklahoma Water Resources Board (2016), but that publication did not include a stage-storage relation table. The USGS, in cooperation with the City of Shawnee, performed a bathymetric survey of Shawnee Reservoir in 2016 and released the bathymetric-survey data in 2017. The purposes of the bathymetric survey were to (1) develop a detailed bathymetric map of the reservoir and (2) determine the relations between stage and reservoir storage capacity and between stage and reservoir surface area. The bathymetric map may serve as a baseline to which temporal changes in storage capacity, due to sedimentation and other factors, can be compared. The stage-storage relation may be used in the reporting of real-time Shawnee Reservoir storage capacity at USGS station 07241600 to support water-resource management decisions by the City of Shawnee.

Oklahoma

Bathymetry of Clear Creek Reservoir, Chaffee County, Colorado, 2016

To better characterize the water supply capacity of Clear Creek Reservoir, Chaffee County, Colorado, the U.S. Geological Survey, in cooperation with the Pueblo Board of Water Works and Colorado Mountain College, carried out a bathymetry survey of Clear Creek Reservoir. A bathymetry map of the reservoir is presented here with the elevation-surface area and the elevation-volume relations. The bathymetry survey was carried out June 6–9, 2016, using a man-operated boat-mounted, multibeam echo sounder integrated with a Global Positioning System and a terrestrial survey using real-time kinematic Global Navigation Satellite Systems. The two collected datasets were merged and imported into geographic information system software. The equipment and methods used in this study allowed water-resource managers to maintain typical reservoir operations, eliminating the need to empty the reservoir to carry out the survey.

Colorado

Periphyton and Macroinvertebrate Communities at Five Sites in the San Joaquin River Basin, California, during June and September, 2001

The effects of agriculture, particularly from the use of pesticides, on aquatic ecosystems in the San Joaquin River Basin concern many aquatic resource managers, water quality managers, and water users. A total of five sites were sampled once in June 2001 and once in September 2001 to document the periphyton (attached algae) community, the benthic macroinvertebrate (insects and non-insects) community, and stream habitat conditions. The purposes of the study were to document existing conditions and, to the extent possible, relate the periphyton and macroinvertebrate community condition to environmental conditions. A total of 161 taxa of algae were collected during the study. Samples from the richest targeted habitat, woody debris, included 109 taxa. In both the June and September samples, greater than 95 percent of the taxa collected were diatoms. Cluster analysis and detrended correspondence analysis of sample data showed that Orestimba Creek had a very different periphyton community than the Merced and Tuolumne Rivers. Salt Slough and the San Joaquin River had community compositions that were intermediate between the two extremes. A total of 126 taxa of macroinvertebrates were collected during the study. Samples from woody debris included 59 taxa. The samples included a variety of both insect and non-insect taxa. Cluster analysis and detrended correspondence analysis of sample data showed that Orestimba Creek was very different from the Merced River and Tuolumne River, similar to the results for periphyton. Orestimba Creek was dominated by non-insects, while the Merced and Tuolumne Rivers were dominated by insects. Salt Slough was more similar to Orestimba Creek because of the abundance of non-insects. The San Joaquin River was more similar to the Merced and Tuolumne Rivers. There was no evidence of major differences between June and September samples for either the periphyton or macroinvertebrate communities. Specific conductance (a surrogate for salinity) and several habitat measures were associated with differences in the periphyton and macroinvertebrate communities at the five sites. Additional sampling at more sites over a longer period of time will likely be necessary before the effects of water quality and habitat conditions on aquatic communities are fully understood in the San Joaquin River Basin.

Scientific Investigations Report

Conceptual model and numerical simulation of the ground-water-flow system in the unconsolidated deposits of the Colville River Watershed, Stevens County, Washington

Increased use of ground- and surface-water supplies in watersheds of Washington State in recent years has created concern that insufficient instream flows remain for fish and other uses. Issuance of new ground-water rights in the Colville River Watershed was halted by the Washington Department of Ecology due to possible hydraulic continuity of the ground and surface waters. A ground-water-flow model was developed to aid in the understanding of the ground-water system and the regional effects of ground-water development alternatives on the water resources of the Colville River Watershed. The Colville River Watershed is underlain by unconsolidated deposits of glacial and non-glacial origin. The surficial geologic units and the deposits at depth were differentiated into aquifers and confining units on the basis of areal extent and general water-bearing characteristics. Five principal hydrogeologic units are recognized in the study area and form the basis of the ground-water-flow model. A steady-state ground-water-flow model of the Colville River Watershed was developed to simulate September 2001 conditions. The simulation period represented a period of below-average precipitation. The model was calibrated using nonlinear regression to minimize the weighted differences or residuals between simulated and measured hydraulic head and stream discharge. Simulated inflow to the model area was 53,000 acre-feet per year (acre-ft/yr) from precipitation and secondary recharge, and 36,000 acre-ft/yr from stream and lake leakage. Simulated outflow from the model was primarily through discharge to streams and lakes (71,000 acre-ft/yr), ground-water outflow (9,000 acre-ft/yr), and ground-water withdrawals (9,000 acre-ft/yr). Because the period of simulation, September 2001, was extremely dry, all components of the ground-water budget are presumably less than average flow conditions. The calibrated model was used to simulate the possible effects of increased ground-water pumping. Although the steady-state model cannot be used to predict how long it would take for effects to occur, it does simulate the ultimate response to such changes relative to September 2001 (relatively dry) conditions. Steady-state simulations indicated that increased pumping would result in decreased discharge to streams and lakes and decreased ground-water outflow. The location of the simulated increased ground-water pumping determined the primary source of the water withdrawn. Simulated pumping wells in the northern end of the main Colville River valley diverted a large percentage of the pumpage from ground-water outflow. Simulated pumping wells in the southern end of the main Colville River valley diverted a large percentage of the pumpage from flow to rivers and streams. The calibrated steady-state model also was used to simulate predevelopment conditions, during which no ground-water pumping, secondary recharge, or irrigation application occurred. Cumulative streamflow in the Colville River Watershed increased by 1.1 cubic feet per second, or about 36 percent of net ground-water pumping in 2001. The model is intended to simulate the regional ground-water-flow system of the Colville River Watershed and can be used as a tool for water-resource managers to assess the ultimate regional effects of changes in stresses. The regional scale of the model, coupled with relatively sparse data, must be considered when applying the model in areas of poorly understood hydrology, or examining hydrologic conditions at a larger scale than what is appropriate.

Scientific Investigations Report

Regionalized equations for bankfull-discharge and channel characteristics of streams in New York State—Hydrologic Region 5 in central New York

Equations that relate drainage area to bankfull discharge and channel dimensions (width, depth, and cross-sectional area) at gaged sites are needed to define bankfull discharge and channel dimensions at ungaged sites and to provide information for the design of stream-restoration projects. Such equations are most accurate if derived from streams within an area of uniform hydrologic, climatic, and physiographic conditions and applied only within that region. A study to develop equations to predict bankfull data for ungaged streams in New York established eight regions that coincided with previously defined hydrologic regions. This report presents drainage areas and bankfull characteristics (discharge and channel dimensions) for streams in central New York (Region 5) selected for this pilot study. Stream-survey data and discharge records from seven active (currently gaged) sites and nine inactive (discontinued gaged) sites were used in regression analyses to relate size of drainage area to bankfull discharge and bankfull channel width, depth, and cross-sectional area. The resulting equations are: (1) bankfull discharge, in cubic feet per second = 45.3*(drainage area, in square miles) 0.856 ; (2) bankfull channel width, in feet = 13.5*(drainage area, in square miles) 0.449 ; (3) bankfull channel depth, in feet = 0.801*(drainage area, in square miles) 0.373 ; and (4) bankfull channel cross-sectional area, in square feet = 10.8*(drainage area, in square miles) 0.823 . The high correlation coefficients (R 2 ) for these four equations (0.96, 0.92, 0.91, 0.98, respectively) indicate that much of the variation in the variables is explained by the size of the drainage area. Recurrence intervals for the estimated bankfull discharge of each stream ranged from 1.11 to 3.40 years; the mean recurrence interval was 1.51 years. The 16 surveyed streams were classified by Rosgen stream type; most were mainly C-type reaches, with occasional B- and F-type reaches. The Region 5 equation was compared with equations developed for six other large areas in the Northeast. The major differences among results indicate a need to refine equations so they can be applied by water-resources managers to local planning and design efforts.

Scientific Investigations Report

Regionalized equations for bankfull-discharge and channel characteristics of streams in New York State—Hydrologic Region 6 in the Southern Tier of New York

Equations that relate bankfull discharge and channel characteristics (width, depth, and cross-sectional area) to drainage-area size at gaged sites are needed to define bankfull discharge and channel dimensions at ungaged sites and to provide information for watershed assessments, stream-channel classification, and the design of stream-restoration projects. Such equations are most accurate if derived from streams within an area of uniform hydrologic, climatic, and physiographic conditions and applied only within that region. In New York State, eight hydrologic regions were previously defined on the basis of similar high-flow (flood) characteristics. This report presents drainage areas and associated bankfull characteristics (discharge and channel dimensions) for surveyed streams in southwestern New York (Region 6). Stream-survey data and discharge records from 11 active (currently gaged) sites and 3 inactive (discontinued) sites were used in regression analyses to relate bankfull discharge and bankfull channel width, depth, and cross-sectional area to the size of the drainage area. The resulting equations are: (1) bankfull discharge, in cubic feet per second = 48.0*(drainage area, in square miles) 0.842 ; (2) bankfull channel width, in feet = 16.9*(drainage area, in square miles) 0.419 ; (3) bankfull channel depth, in feet = 1.04*(drainage area, in square miles) 0.244 ; and (4) bankfull channel cross-sectional area, in square feet = 17.6*(drainage area, in square miles) 0.662 . The coefficient of determination (R 2 ) for these four equations were 0.90, 0.79, 0.64, and 0.89, respectively. The high correlation coefficients for bankfull discharge and cross-sectional area indicate that much of the variation in these variables is explained by the size of the drainage area. The smaller correlation coefficients for bankfull channel width and depth indicate that other factors also affect these relations. Recurrence intervals for the estimated bankfull discharge of each stream ranged from 1.01 to 2.35 years; the mean recurrence interval was 1.54 years. The 14 surveyed streams were classified by Rosgen stream type; most were C-type reaches, with occasional B-type reaches. The Region 6 equation (curve) for bankfull discharge was compared with equations previously developed for four other large areas in New York State and southeastern Pennsylvania. The differences among results indicate that, although the equations need to be refined by region before being applied by water-resources managers to local planning and design efforts, similar regions have similar relations between bankfull discharge and channel characteristics.

Scientific Investigations Report

Suspended-sediment yields and stream-channel processes on Judy's Branch watershed in the St. Louis Metro East region in Illinois

Judy's Branch watershed, a small basin (8.64 square miles) in the St. Louis Metro East region in Illinois, was selected as a pilot site to determine suspended-sediment yields and stream-channel processes in the bluffs and American Bottoms (expansive low-lying valley floor in the region). Suspended-sediment and stream-chan-nel data collected and analyzed for Judy's Branch watershed are presented in this report to establish a baseline of data for water-resource managers to evaluate future stream rehabilitation and manage-ment alternatives. The sediment yield analysis determines the amount of sediment being delivered from the watershed and two subwatersheds: an urban tributary and an undeveloped headwater (pri-marily agricultural). The analysis of the subwater-sheds is used to compare the effects of urbanization on sediment yield to the river. The stream-channel contribution to sediment yield was determined by evaluation of the stream-channel processes operat-ing on the streambed and banks of Judy's Branch watershed. Bank stability was related to hydrologic events, bank stratigraphy, and channel geometry through model development and simulation. The average suspended-sediment yield from two upland subwatersheds (drainage areas of 0.23 and 0.40 sq.mi. was 1,163 tons per square mile per year (tons/sq.mi.-year) between July 2000 and June 2004. The suspended-sediment yield at the Route 157 station was 2,523 tons/sq.mi.-year, near the outlet of Judy's Branch watershed (drainage area = 8.33 sq.mi.). This is approximately 1,360 tons/sq.mi.-year greater than the average at the upland stations for the same time period. This result is unexpected in that, generally, the suspended-sediment yield decreases as the watershed area increases because of sediment stored in the channel and flood plain. The difference indicates a possible increase in yield from a source, such as bank retreat, and supports the concept that land-use changes increase stream-flows that may in turn result in higher rates of bank retreat. Utilizing both bank-rod data and resurveyed cross-section data, it was determined that approxi-mately half of the suspended- sediment yield at Route 157 during July 2000-June 2004 came from bank retreat. Given that bank retreat can be a substantial portion of the sediment yield, understanding bank stability processes is important. Bank stability can be assessed mathematically by computing the factor of safety, which is defined by the ratio of the shear strength (resisting force) along the failure surface and the shear stress (driving gravitational force). Once the factor of safety falls below one, the bank theoretically becomes unstable. Bank-stability conditions were related to hydrologic events, bank type, and channel geometry through model develop-ment and simulation. The most common type of bank in the watershed consists of cohesive alluvial soil deposits overlying a stiff glacial till. A stabil-ity chart for different bank types was developed using a bank-stability analysis. Banks steeper than 70 degrees and higher than from 10 to 11.5 feet (depending on bank type) become at risk for mass failure in the watershed under conditions that pro-mote saturation of the bank and a sudden drop in the river level.

Scientific Investigations Report

Groundwater conditions in the Brunswick-Glynn County area, Georgia, 2009

The Upper Floridan aquifer is contaminated with saltwater in a 2-square-mile area of downtown Brunswick, Georgia. The presence of this saltwater has limited the development of the groundwater supply in the Glynn County area. Hydrologic, geologic, and water-quality data are needed to effectively manage water resources. Since 1959, the U.S. Geological Survey (USGS) has conducted a cooperative water program with the City of Brunswick and Glynn County to monitor and assess the effect of groundwater development on saltwater intrusion within the Floridan aquifer system. The potential development of alternative sources of water in the Brunswick and surficial aquifer systems also is an important consideration in coastal areas. During calendar year 2009, the cooperative water program included continuous water-level recording of 13 wells completed in the Floridan, Brunswick, and surficial aquifer systems; collecting water levels from 46 wells to map the potentiometric surface of the Upper Floridan aquifer in Glynn County during August 2009; and collecting and analyzing water samples from 55 wells completed in the Floridan aquifer system, of which 27 wells were used to map chloride concentrations in the upper water-bearing zone of the Upper Floridan aquifer in the Brunswick area during August 2009. Periodic water-level measurements also were collected from two wells completed in the Upper Floridan aquifer and four wells completed in the Brunswick aquifer system on Jekyll Island. Equipment was installed on one well to enable real-time specific conductance monitoring in the area surrounding the chloride plume. During 2008-2009, water levels in 30 of the 32 wells monitored in the Brunswick-Glynn County area rose at a rate of 0.24 to 7.58 feet per year (ft/yr). The largest rise of 7.58 ft/yr was in the Upper Floridan aquifer. These rises corresponded to a period of above normal precipitation and decreased pumping. Declines during 2008-2009 were recorded in wells completed in the Brunswick aquifer system (0.37 ft/yr) and Lower Floridan aquifer (0.83 ft/yr). Chloride data collected by two local industrial groundwater users at their well fields since 1958 were compiled and compared with data collected by the USGS during the same period. The results indicate that chloride concentrations at the two well fields have continued to rise despite modification of production wells to eliminate deep saline zones and decreases in pumpage at both facilities. One of the industrial users, Pinova Inc., plugged the lower portions of nine production wells in the mid to late 1960s, which generally decreased chloride concentrations to less than 100 milligrams per liter (mg/L) for a period of 10 to 20 years. However, chloride concentrations eventually returned to previous levels despite decreases in pumpage. During 1990-2009, chloride concentrations at the other industrial user's well field (Georgia-Pacific Cellulose LLC) generally increased despite a 16 million gallon per day decrease in pumpage during this period. Data from the Georgia-Pacific Cellulose well field and additional chloride data from USGS observation wells located to the east indicate continued movement of chloride from the source area located southeast of the site toward the well field.

Georgia

Land subsidence along the Delta-Mendota Canal in the northern part of the San Joaquin Valley, California, 2003-10

Extensive groundwater withdrawal from the unconsolidated deposits in the San Joaquin Valley caused widespread aquifer-system compaction and resultant land subsidence from 1926 to 1970—locally exceeding 8.5 meters. The importation of surface water beginning in the early 1950s through the Delta-Mendota Canal and in the early 1970s through the California Aqueduct resulted in decreased pumping, initiation of water-level recovery, and a reduced rate of compaction in some areas of the San Joaquin Valley. However, drought conditions during 1976–77 and 1987–92, and drought conditions and regulatory reductions in surface-water deliveries during 2007–10, decreased surface-water availability, causing pumping to increase, water levels to decline, and renewed compaction. Land subsidence from this compaction has reduced freeboard and flow capacity of the Delta-Mendota Canal, the California Aqueduct, and other canals that deliver irrigation water and transport floodwater. The U.S. Geological Survey, in cooperation with the U.S. Bureau of Reclamation and the San Luis and Delta-Mendota Water Authority, assessed land subsidence in the vicinity of the Delta-Mendota Canal as part of an effort to minimize future subsidence-related damages to the canal. The location, magnitude, and stress regime of land-surface deformation during 2003–10 were determined by using extensometer, Global Positioning System (GPS), Interferometric Synthetic Aperture Radar (InSAR), spirit leveling, and groundwater-level data. Comparison of continuous GPS, shallow extensometer, and groundwater-level data, combined with results from a one-dimensional model, indicated the vast majority of the compaction took place beneath the Corcoran Clay, the primary regional confining unit. Land-surface deformation measurements indicated that much of the northern portion of the Delta-Mendota Canal (Clifton Court Forebay to Check 14) was fairly stable or minimally subsiding on an annual basis; some areas showed seasonal periods of subsidence and of uplift that resulted in little or no longer-term elevation loss. Many groundwater levels in this northern area did not reach historical lows during 2003–10, indicating that deformation in this region was primarily elastic. Although the northern portion of the Delta-Mendota Canal was relatively stable, land-surface deformation measurements indicated the southern portion of the Delta-Mendota Canal (Checks 15–21) subsided as part of a large subsidence feature centered about 15 kilometers northeast of the Delta-Mendota Canal, south of the town of El Nido. Results of InSAR analysis indicated at least 540 millimeters of subsidence near the San Joaquin River and the Eastside Bypass during 2008–10, which is part of a 3,200 square-kilometer area—including the southern part of the Delta-Mendota Canal—affected by 20 millimeters or more of subsidence during the same period. Calculations indicated that the subsidence rate doubled in 2008 in some areas. The GPS surveys done in 2008 and 2010 confirmed the high subsidence rate measured by using InSAR for the same period. Water levels in many shallow and deep wells in this area declined during 2007–10; water levels in many deep wells reached historical lows, indicating that subsidence measured during this period was largely inelastic. InSAR-derived subsidence maps for various periods during 2003–10 showed that the area of maximum active subsidence (that is, the largest rates of subsidence) shifted from its historical (1926–70) location southwest of Mendota to south of El Nido. Continued groundwater-level and land-subsidence monitoring in the San Joaquin Valley is important because (1) regulatory- and drought-related reductions in surface-water deliveries since 1976 have resulted in increased groundwater pumping and associated land subsidence, and (2) land use and associated groundwater pumping continue to change throughout the valley. The availability of surface water remains uncertain; even during record-setting precipitation years, such as 2010–11, water deliveries have fallen short of requests and groundwater pumping was required to meet the irrigation demand. Due to the expected continued demand for irrigation supply water and the limitations and uncertainty of surface-water supplies, groundwater pumping and associated land subsidence is likely to continue in the future. Spatially detailed information on land subsidence is needed to facilitate minimization of future subsidence-related damages to the Delta-Mendota Canal and other infrastructure in the San Joaquin Valley. The integration of subsidence, deformation, and water-level measurements—particularly continuous measurements—enables the analysis of aquifer-system response to increased groundwater pumping, which in turn, enables identification of the preconsolidation head and calculation of aquifer-system storage properties. This information can be used to improve numerical model simulations of groundwater flow and aquifer-system compaction and allow for consideration of land subsidence in the evaluation of water-resource management alternatives.

California

Simulation of hydrologic conditions and suspended-sediment loads in the San Antonio River Basin downstream from San Antonio, Texas, 2000-12

Suspended sediment in rivers and streams can play an important role in ecological health of rivers and estuaries and consequently is an important issue for water-resource managers. To better understand suspended-sediment loads and transport in a watershed, the U.S. Geological Survey (USGS), in cooperation with the San Antonio River Authority, developed a Hydrological Simulation Program—FORTRAN model to simulate hydrologic conditions and suspended-sediment loads during 2000–12 for four watersheds, which comprise the overall study area in the San Antonio River Basin (hereinafter referred to as the “USGS–2014 model”). The study area consists of approximately 2,150 square miles encompassing parts of Bexar, Guadalupe, Wilson, Karnes, DeWitt, Goliad, Victoria, and Refugio Counties. The USGS–2014 model was calibrated for hydrology and suspended sediment for 2006–12. Overall, model-fit statistics and graphic evaluations from the calibration and testing periods provided multiple lines of evidence indicating that the USGS–2014 model simulations of hydrologic and suspended-sediment conditions were mostly “good” to “very good.” Model simulation results indicated that approximately 1,230 tons per day of suspended sediment exited the study area and were delivered to the Guadalupe River during 2006–12, of which approximately 62 percent originated upstream from the study area. Sample data and simulated model results indicate that most of the suspended-sediment load in the study area consisted of silt- and clay-sized particles (less than 0.0625 millimeters). The Cibolo Creek watershed was the largest contributor of suspended sediment from the study area. For the entire study area, open/developed land and cropland exhibited the highest simulated soil erosion rates; however, the largest contributions of sediment (by land-cover type) were pasture and forest/rangeland/shrubland, which together composed approximately 80 percent of the land cover of the study area and generated about 70 percent of the suspended-sediment load from the study area.

Texas

Comparison of regression relations of bankfull discharge and channel geometry for the glaciated and nonglaciated settings of Pennsylvania and southern New York

Streambank erosion in areas of past glacial deposition has been shown to be a dominant source of sediment to streams. Water resource managers are faced with the challenge of developing long and short term (emergency) stream restoration efforts that rely on the most suitable channel geometry for project design. A geomorphic dataset of new (2016, n=5) and previous (1999–2006, n=96) estimates of bankfull discharge and channel dimensions at U.S. Geological Survey streamflow-gaging stations was compiled to present and contrast the glaciated and unglaciated noncarbonate settings of southern New York and Pennsylvania that included selected areas of Maryland. Empirical models were developed by using simple linear regressions that relate bankfull discharge and channel geometry to drainage area (regional curves). Significant relations (p<0.05) were able to explain variability with coefficient of determination (R2 ) values of 0.89 for bankfull discharge, 0.94 for cross-sectional area, 0.87 for bankfull width, and 0.83 for bankfull depth. These regression relations for the glaciated noncarbonate settings of northern Pennsylvania and southern New York were able to provide a slightly better fit than regional curve models developed previously for the entire noncarbonate region of Pennsylvania. Although, the analysis of covariance (ANCOVA) results for comparison between regression equations for the glaciated and unglaciated settings showed that except for the significant intercept of bankfull discharge versus drainage area (F=8.26, p-value<0.005), the regression equations are not significantly different between the glaciated and unglaciated setting of Pennsylvania and southern New York. Therefore, data stratification by glaciation does not improve regional curves relations developed previously for the noncarbonate (glaciated and unglaciated) and carbonate settings of Pennsylvania and Maryland. Further analysis that incorporates data stratification or multivariate approaches based on mean annual runoff, precipitation, slope, stream classification, or other relevant parameters may optimize the accuracy and utility of statewide models. The new estimates of bankfull discharge and channel dimensions at streamflowgaging sites and updated drainage areas from StreamStats were incorporated into previously developed regional curves to produce an updated set of regression relations of bankfull discharge and channel geometry for the noncarbonate and carbonate settings of Pennsylvania and Maryland.

New York, Pennsylvania

Estimated groundwater recharge from a water-budget model incorporating selected climate projections, Island of Maui, Hawai‘i

Demand for freshwater on the Island of Maui is expected to increase by 45 percent between 2015 and 2035. Groundwater availability on Maui is affected by changes in climate and agricultural irrigation. To evaluate the availability of fresh groundwater under projected future climate conditions and changing agricultural irrigation practices, estimates of groundwater recharge are needed. A water-budget model with a daily computation interval was used to estimate the spatial distribution of recharge on Maui for one present-day and two future-climate scenarios. All three scenarios used 2017 land cover. The two future-climate scenarios, including one wetter than the present-day scenario and one drier than the present-day scenario, were developed using available high-resolution downscaled climate projections. The drier future scenario was developed using projections for a Representative Concentration Pathway warming scenario during 2071–99 with total radiative forcing of 8.5 Watts per square meter by the year 2100 (RCP8.5 2071–99 scenario), whereas the wetter future scenario was developed using projections for a “Special Report on Emissions Scenarios” A1B emission scenario during 2080–99 (A1B 2080–99 scenario). For the RCP8.5 2071–99 scenario, projected mean annual recharge decrease for Maui is about 172 million gallons per day, or about 14 percent less than present-day recharge, which is estimated to be 1,232 million gallons per day. Recharge for the RCP8.5 2071–99 scenario is projected to decrease in 22 of Maui’s 25 aquifer systems, which are defined by the Hawaiʻi Commission on Water Resource Management. For the A1B 2080–99 future scenario, projected mean annual recharge increase for Maui is about 144 million gallons per day, or about 12 percent more than present-day recharge. Recharge for the A1B 2080–99 scenario is projected to increase in 17 of Maui’s 25 aquifer systems. Between the two future scenarios, a total of 11 aquifer systems show similar direction in drying (Kahului, Kama‘ole, Lualaʻilua, Makawao, Olowalu, Pāʻia, Ukumehame, Waikapū) or wetting (Honopou, Kawaipapa, and Waikamoi) changes for recharge. Selectively modifying the climate inputs for the A1B 2080–99 scenario indicates that the projected changes in rainfall account for most of the projected changes in recharge for Maui’s 25 aquifer systems. However, projected changes in reference evapotranspiration and forest-canopy evaporation also can account for a substantial part of the projected changes in recharge where changes in reference evapotranspiration are relatively large and where changes in forest-canopy evaporation extend across large forested areas. Projected changes in daily rainfall frequency have a relatively small but non-negligible impact on recharge estimates.

Hawaii

Simulation of groundwater-level changes from projected groundwater withdrawals in the Truxton basin, northwestern Arizona

A three-dimensional, numerical groundwater flow model of the Hualapai Plateau and Truxton basin was developed to assist water-resource managers in understanding the potential effects of projected groundwater withdrawals on groundwater levels and storage in the basin. The Truxton Basin Hydrologic Model (TBHM) is a transient model that simulates the hydrologic system for the years 1976 through 2139, including hypothetical low-, medium-, and high-groundwater withdrawal scenarios beginning in 2020. The simulated effects of these withdrawal scenarios are presented as groundwater-level changes from the year 2020 to 2070, and from 2020 to 2140. Hydrologic properties in the TBHM are derived from calibration of a steady-state model of the predevelopment (before 1976) groundwater system. The future pumping scenarios are each simulated with three different interpretations of basin depth supported by geophysical data. For each of the resulting nine transient models, a Monte Carlo approach is used to produce a range of possible and probable groundwater-level changes at points throughout the basin given probabilistic ranges of hydrologically reasonable aquifer property values supported by the model calibration results. The ensemble of models that simulate the future pumping scenarios include pumping from the existing well field (three wells) plus additional pumping from a proposed new well. Simulated high future pumping increases progressively to 1,840 acre-feet per year in 2120 and produces a range of drawdowns between 20 and 39 feet (ft) near the pumping center, with a median drawdown of 28 ft. The low future pumping scenario, which increases progressively to 650 acre-ft per year in 2120, produces a range of drawdowns between 5 and 15 ft, with a median drawdown of 10 ft at the same location over the same period of time.

Arizona