USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water Resources Update”

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.

359 records · Page 20Linked to original sources

User's guide to the wetland creation/restoration data base, version 2

Wetland creation or restoration projects are frequently proposed as mitigation for unavoidable wetland losses, as components of wetland enhancement programs, and as tools to accomplish specific objectives such as waterfowl production or flood control. There is considerable controversy concerning the effectiveness of such projects as well as the most appropriate and efficient techniques to employ. The importance of the resource and the long time scales involved in fully evaluating a creation or restoration effort make it imperative to consider existing information as fully as possible in the development and evaluation of wetland creation or restoration proposals. To aid in the evaluation of wetland/creation efforts, the U.s. Fish and Wildlife Service (FWS), National Ecology Research Center, has developed the Wetland Creation/Restoration (WCR) Data Base. The data base is a highly indexed or keyworded bibliography of wetland creation or restoration articles. ("Articles" refers to any type of publication that deals specifically with wetland creation/restoration projects or studies.) The scope of the articles is international, although most of them are concerned with projects conducted in the United States. Information coded for each article includes author; citation; type of wetland and its location in terms of state, ecoregion, and FWS region; type of study undertaken; objectives in creating or restoring the wetland; actions performed to realize those objectives; length of time encompassed by the study; evaluation of results and responses to the wetland creation/restoration actions; and a listing of plant species significant to the project. A brief annotation summarizes the article and includes any significant additional information that may not be adequately reflected in the above described fields. Many of these articles describe only one or two components of a total wetland restoration effort. Planning a project that is designed to restore a wetland system (including at least some of its functions) is similar to constructing a picture from a number of puzzle pieces--missing pieces represent data gaps or information that is not available. Articles range from specific case studies, to overviews of restoration methods and techniques, to planning restoration projects and assessing programmatic and administrative backgrounds and interactions. In this data base, the term "restoration" is applied loosely to include rehabilitation of wetlands. It may refer to a number of situations or actions including, but not limited to: 1. breaching dikes or plugging drains; 2. water pollution clean-up; 3. conversion of eutrophic conditions; 4. wastewater treatment; 5. recolonization of previously disturbed or denuded areas; 6. amelioration of adverse conditions (erosion, wave, or wind action); 7. soil treatment --mulching, fertilization; 8. rerouting streams --may include construction of meander patterns; 9. monitoring natural vegetation; or 0. excluding grazers (geese, cattle) and monitoring results. This report describes the format and content of Version 2 of the WCR data base. Version 2 differs from the previous version described in SchnellerMcDonald et al. (1988): several fields have been dropped and condensed and new records have been added. Version 2 includes all records distributed with the earlier version and its updates. We recommend you replace any previous version with Version 2.

Book↗

Peak-flow frequency relations and evaluation of the peak-flow gaging network in Nebraska

Estimates of peak-flow magnitude and frequency are required for the efficient design of structures that convey flood flows or occupy floodways, such as bridges, culverts, and roads. The U.S. Geological Survey, in cooperation with the Nebraska Department of Roads, conducted a study to update peak-flow frequency analyses for selected streamflow-gaging stations, develop a new set of peak-flow frequency relations for ungaged streams, and evaluate the peak-flow gaging-station network for Nebraska. Data from stations located in or within about 50 miles of Nebraska were analyzed using guidelines of the Interagency Advisory Committee on Water Data in Bulletin 17B. New generalized skew relations were developed for use in frequency analyses of unregulated streams. Thirty-three drainage-basin characteristics related to morphology, soils, and precipitation were quantified using a geographic information system, related computer programs, and digital spatial data.For unregulated streams, eight sets of regional regression equations relating drainage-basin to peak-flow characteristics were developed for seven regions of the state using a generalized least squares procedure. Two sets of regional peak-flow frequency equations were developed for basins with average soil permeability greater than 4 inches per hour, and six sets of equations were developed for specific geographic areas, usually based on drainage-basin boundaries. Standard errors of estimate for the 100-year frequency equations (1percent probability) ranged from 12.1 to 63.8 percent. For regulated reaches of nine streams, graphs of peak flow for standard frequencies and distance upstream of the mouth were estimated.The regional networks of streamflow-gaging stations on unregulated streams were analyzed to evaluate how additional data might affect the average sampling errors of the newly developed peak-flow equations for the 100-year frequency occurrence. Results indicated that data from new stations, rather than more data from existing stations, probably would produce the greatest reduction in average sampling errors of the equations.

Nebraska↗

Analysis of projected water availability with current basin management plan, Pajaro Valley, California

The projection and analysis of the Pajaro Valley Hydrologic Model (PVHM) 34 years into the future using MODFLOW with the Farm Process (MF-FMP) facilitates assessment of potential future water availability. The projection is facilitated by the integrated hydrologic model, MF-FMP that fully couples the simulation of the use and movement of water from precipitation, streamflow, runoff, groundwater flow, and consumption by natural and agricultural vegetation throughout the hydrologic system at all times. MF-FMP allows for more complete analysis of conjunctive-use water-resource systems than previously possible with MODFLOW by combining relevant aspects of the landscape with the groundwater and surface-water components. This analysis is accomplished using distributed cell-by-cell supply-constrained and demand-driven components across the landscape within “water-balance subregions” (WBS) comprised of one or more model cells that can represent a single farm, a group of farms, watersheds, or other hydrologic or geopolitical entities. Analysis of conjunctive use would be difficult without embedding the fully coupled supply-and-demand into a fully coupled simulation, and are difficult to estimate a priori. The analysis of projected supply and demand for the Pajaro Valley indicate that the current water supply facilities constructed to provide alternative local sources of supplemental water to replace coastal groundwater pumpage, but may not completely eliminate additional overdraft. The simulation of the coastal distribution system (CDS) replicates: 20 miles of conveyance pipeline, managed aquifer recharge and recovery (MARR) system that captures local runoff, and recycled-water treatment facility (RWF) from urban wastewater, along with the use of other blend water supplies, provide partial relief and substitution for coastal pumpage (aka in-lieu recharge). The effects of these Basin Management Plan (BMP) projects were analyzed subject to historical climate variations and assumptions of 2009 urban water demand and land use. Water supplied directly from precipitation, and indirectly from reuse, captured local runoff, and groundwater is necessary but inadequate to satisfy agricultural demand without coastal and regional storage depletion that facilitates seawater intrusion. These facilities reduce potential seawater intrusion by about 45% with groundwater levels in the four regions served by the CDS projected to recover to levels a few feet above sea level. The projected recoveries are not high enough to prevent additional seawater intrusion during dry-year periods or in the deeper aquifers where pumpage is greater. While these facilities could reduce coastal pumpage by about 55% of the historical 2000–2009 pumpage for these regions, and some of the water is delivered in excess of demand, other coastal regions continue to create demands on coastal pumpage that will need to be replaced to reduce seawater intrusion. In addition, inland urban and agricultural demands continue to sustain water levels below sea level causing regional landward gradients that also drive seawater intrusion. Seawater intrusion is reduced by about 45% but it supplies about 55% of the recovery of groundwater levels in the coastal regions served by the CDS. If economically feasible, water from summer agricultural runoff and tile-drain returnflows could be another potential local source of water that, if captured and reused, could offset the imbalance between supply and demand as well as reducing discharge of agricultural runoff into the National Marine Sanctuary of Monterey Bay. A BMP update (2012) identifies projects and programs that will fund a conservation program and will provide additional, alternative water sources to reduce or replace coastal and inland pumpage, and to replenish the aquifers with managed aquifer recharge in an inland portion of the Pajaro Valley.

California↗

Geologic map of Wupatki National Monument and vicinity, Coconino County, Northern Arizona

Introduction The geologic map of Wupatki National Monument is a cooperative effort between the U.S. Geological Survey, the National Park Service, and the Navajo Nation to provide geologic information for resource management officials of the National Park Service, U.S. Forest Service, Navajo Indian Reservation (herein the Navajo Nation), and visitor information services at Wupatki National Monument, Arizona. Funding for the map was provided in part by the Water Rights Branch of the Water Resources Division of the National Park Service. Field work on the Navajo Nation was conducted under a permit from the Navajo Nation Minerals Department. Any persons wishing to conduct geologic investigations on the Navajo Nation must first apply for, and receive, a permit from the Navajo Nation Minerals Department, P.O. Box 1910, Window Rock, Arizona 86515, telephone (928)-871-6587. Wupatki National Monument lies within the USGS 1:24,000-scale Wupatki NE, Wupatki SE, Wupatki SW, Gray Mountain, East of SP Mountain, and Campbell Francis Wash quadrangles in northern Arizona. The map is bounded approximately by longitudes 111? 16' to 111? 32' 30' W. and latitudes 35? 30' to 35? 37' 40' N. The map area is in Coconino County on the southern part of the Colorado Plateaus geologic province (herein Colorado Plateau). The map area is locally subdivided into three physiographic parts, the Coconino Plateau, the Little Colorado River Valley, and the San Francisco Volcanic Field as defined by Billingsley and others (1997) [fig. 1]. Elevations range from 4,220 ft (1,286 m) at the Little Colorado River near the northeast corner of the map area to about 6,100 ft (1,859 m) at the southwest corner of the map area. The small community of Gray Mountain is about 16 mi (26 km) northwest of Wupatki National Monument Visitor Center, and Flagstaff, Arizona, the nearest metropolitan area, is about 24 mi (38 km) southwest of the Visitor Center (fig. 1). U.S. Highway 89 provides access to the west entrance of Wupatki National Monument. A paved Coconino County road provides a loop from Wupatki National Monument south to Sunset Crater National Monument and back to U.S. Highway 89 about 10 mi (16 km) north of Flagstaff, Arizona. Access to Coconino National Forest is via dirt roads maintained by the National Forest Service. Several unimproved dirt roads on Babbitt Ranch lands provide limited access to remote areas north of Wupatki National Monument. Travel is mostly restricted to paved roads within Wupatki National Monument, and a dirt road that crosses the Little Colorado River provides access to the Navajo Nation area east and northeast of the Little Colorado River. The Little Colorado River crossing is not bridged and can be impassable when the river is flowing. Four-wheel-drive vehicles are recommended but not necessary for travel in remote parts of the Navajo Nation. Extra food and water are highly recommended for travel in this sandy area. Land ownership north of Wupatki National Monument forms a checkerboard pattern between private and State land. Coconino National Forest manages lands south of Wupatki National Monument and the National Park Service manages Wupatki National Monument. The Leupp and Tolani Lake Chapters of the Navajo Nation manage the area northeast and east of the Little Colorado River (see land management boundaries on map). The geologic map of Wupatki National Monument provides updated geologic framework information for this part of the Colorado Plateau. The geologic information supports Federal, State, and private land managers when conducting geologic, biologic, and hydrologic investigations and will support future and ongoing geologic and associated scientific investigations of all disciplines within the Wupatki National Monument area.

Arizona↗

Simulation of reservoir storage and firm yields of three surface-water supplies, Ipswich River Basin, Massachusetts

A Hydrologic Simulation Program FORTRAN (HSPF) model previously developed for the Ipswich River Basin was modified to simulate the hydrologic response and firm yields of the water-supply systems of Lynn, Peabody, and Salem-Beverly. The updated model, expanded to include a portion of the Saugus River Basin that supplies water to Lynn, simulated reservoir system storage over a 35-year period (1961-95) under permitted withdrawals and hypothetical restrictions designed to maintain seasonally varied streamflow for aquatic habitat. A firm yield was calculated for each system and each withdrawal restriction by altering demands until the system failed. This is considered the maximum withdrawal rate that satisfies demands, but depletes reservoir storage. Simulations indicate that, under the permitted withdrawals, Lynn and Salem-Beverly were able to meet demands and generally have their reservoir system recover to full capacity during most years; reservoir storage averaged 83 and 82 percent of capacity, respectively. The firm yields for the Lynn and Salem-Beverly systems were 11.4 and 12.2 million gallons per day (Mgal/d), respectively, or 8 and 21 percent more than average 1998-2000 demands, respectively. Under permitted withdrawals and average 1998-2000 demands, the Peabody system failed in all years; thus Peabody purchased water to meet demands. The firm yield for the Peabody system is 3.70 Mgal/d, or 37 percent less than the average 1998-2000 demand. Simulations that limit withdrawals to levels recommended by the Ipswich River Fisheries Restoration Task Group (IRFRTG) indicate that under average 1998-2000 demands, reservoir storage was depleted in each of the three systems. Reservoir storage under average 1998-2000 demands and IRFRTG-recommended streamflow requirements averaged 15, 22, and 71 percent of capacity for the Lynn, Peabody, Salem-Beverly systems, respectively. The firm-yield estimates under the IRFRTG-recommended streamflow requirements were 6.02, 1.94, and 7.69 Mgal/d or 43, 64, and 34 percent less than the average 1998-2000 demands for the Lynn, Peabody, and Salem-Beverly systems, respectively. Simulations that limit withdrawals from the Saugus River to a less stringent set of restrictions (based on an Instream Flow Incremental Methodology study) than those previously simulated indicate that the firm yield of the Lynn system is about 31 percent less than the average 1998-2000 withdrawals (7.31 Mgal/d).

Massachusetts↗

Analysis of the magnitude and frequency of floods in Colorado

Regionalized flood-frequency relations need to be updated on a regular basis (about every 10 years). The latest study on regionalized flood-frequency equations for Colorado used data collected through water year 1981. A study was begun in 1994 by the U.S. Geological Survey, in cooperation with the Colorado Department of Transportation and the Bureau of Land Management, to include streamflow data collected since water year 1981 in the regionalized flood-frequency relations for Colorado. Longer periods of streamflow data and improved statistical analysis methods were used to define regression relations for estimating peak discharges having recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years for unregulated streams in Colorado. The regression relations can be applied to sites of interest on gaged and ungaged streams. Ordinary least-squares regression was used to determine the best explanatory basin or climatic characteristic variables for each peak-discharge characteristic, and generalized least-squares regression was used to determine the best regression relation. Drainage-basin area, mean annual precipitation, and mean basin slope were determined to be statistically significant explanatory variables in the regression relations. Separate regression relations were developed for each of five distinct hydrologic regions in the State. The mean standard errors of estimate and average standard error of prediction associated with the regression relations generally ranged from 40 to 80 percent, except for one hydrologic region where the errors ranged from about 200 to 300 percent. Methods are presented for determining the magnitude of peak discharges for sites located at gaging stations, for sites located near gaging stations on the same stream when the ratio of drainage-basin areas is between about 0.5 and 1.5, and for sites where the drainage basin crosses a flood-region boundary or a State boundary. Methods are presented for determining the magnitude of peak discharges for sites located at gaging stations, for sites located near gaging stations on the same stream when the ratio of drainage-basin areas is between about 0.5 and 1.5, and for sites where the drainage basin crosses a flood-region boundary or a State boundary.

Colorado↗

A deposit model for Mississippi Valley-Type lead-zinc ores

This report is a descriptive model of Mississippi Valley-Type (MVT) lead-zinc deposits that presents their geological, mineralogical and geochemical attributes and is part of an effort by the U.S. Geological Survey Mineral Resources Program to update existing models and develop new models that will be used for an upcoming national mineral resource assessment. This deposit modeling effort by the USGS is intended to supplement previously published models for use in mineral-resource and mineral-environmental assessments. Included in this report are geological, geophysical and geochemical assessment guides to assist in mineral resource estimation. The deposit attributes, including grade and tonnage of the deposits described in this report are based on a new mineral deposits data set of all known MVT deposits in the world. Mississippi Valley-Type (MVT) lead-zinc deposits are found throughout the world but the largest, and more intensely researched deposits occur in North America. The ores consist mainly of sphalerite, galena, and generally lesser amounts of iron sulfides. Silver is commonly an important commodity, whereas Cu is generally low, but is economically important in some deposits. Gangue minerals may include carbonates (dolomite, siderite, ankerite, calcite), and typically minor barite. Silicification of the host rocks (or quartz gangue) is generally minor, but may be abundant in a few deposits. The deposits have a broad range of relationships with their host rocks that includes stratabound, and discordant ores; in some deposits, stratiform and vein ore are important. The most important characteristics of MVT ore deposits are that they are hosted mainly by dolostone and limestone in platform carbonate sequences and usually located at flanks of basins, orogenic forelands, or foreland thrust belts inboard of the clastic rock-dominated passive margin sequences. They have no spatial or temporal relation to igneous rocks, which distinguishes them from skarn or other intrusive rock-related Pb-Zn ores. Abundant evidence has shown that the ore fluids were derived mainly from evaporated seawater and were driven within platform carbonates by large-scale tectonic events. MVT deposits formed mainly during the Phanerozoic with more than 80 percent of the deposits hosted in Phanerozoic rocks and less than 20 percent in Precambrian rocks. Phanerozoic-hosted MVT deposits also account for 94 percent of total MVT ore, and 93 percent of total MVT lead and zinc metal. Many MVT deposits formed during Devonian to Permian time, corresponding to a series of intense tectonic events during assimilation of Pangea. The second most important period for MVT deposit genesis was Cretaceous to Tertiary time when microplate assimilation affected the western margin of North America and Africa-Eurasia. Many subtypes or alternative classifications have been applied to MVT deposits. These alternative classifications reflect geographic and or specific geological features that some workers believe set them apart as unique (for example, Appalachian-, Alpine-, Reocin-, Irish-, Viburnum trend-types). However, we do not consider these alternative classifications or sub-types to be sufficiently different to warrant using them. This report also describes the geoenvironmental characteristic of MVT deposits. The response of MVT ores in the supergene environment is buffered by their placement in carbonate host rocks which commonly results in near-neutral associated drainage water. The geoenvironmental features and anthropogenic mining effects presented in this report illustrates this important environmental aspect of MVT deposits which separates them from other deposit types (especially coal, VHMS, Cu-porphyry, SEDEX, acid-sulfate polymetallic vein).

Scientific Investigations Report↗

Determination of upstream boundary points on southeastern Washington streams and rivers under the requirements of the Shoreline Management Act of 1971

Regulation of the shorelines of the State of Washington, as mandated by the Shoreline Management Act of 1971, requires knowledge of the locations on streams and river reaches where specific regulatory criteria are satisfied. The U.S. Geological Survey conducted a study in 1971 to determine the upstream boundary points of these reaches for many of the State's streams and rivers. Updated upstream boundary points were determined in the current study for all the streams and rivers in southeastern Washington that fall under the jurisdiction of the Shoreline Management Act of 1971. Upstream boundary point locations where the mean annual discharge equals 20 cubic feet per second were determined for 149 streams. In addition, upstream boundary point locations where the mean annual discharge equals 200 cubic feet per second or the drainage area equals 300 square miles were determined for 22 rivers. Boundary point locations were determined by application of multiple-linear-regression equations that relate mean annual discharge to drainage area and mean annual precipitation. Southeastern Washington was divided into five hydrologically distinct regions, and a separate regression equation was developed for each region. The regression equations are based on data for gaging stations with at least 10 years of record. The number of stations in the regression analysis for each of the five regions ranged from 5 to 33. The coefficient of determination, R2, of the regression equations ranged from 0.953 to 0.997. The equation for the Upper Yakima region had the lowest standard error, ranging from -7 to +9 percent for a regression estimate of 20 cubic feet per second. The equation for the Columbia Basin to Palouse region had the highest standard error, ranging from -36 to +55 percent for a regression estimate of 20 cubic feet per second. The approximate error in the location of an upstream boundary point can be calculated using the variables mean annual precipitation of the basin upstream from a boundary point and average basin width in the vicinity of the boundary point. The calculation gives only a rough estimate of the error of the boundary point location, because of the uncertainty in estimating average basin width.

Water-Resources Investigations Report↗

Revised techniques for estimating peak discharges from channel width in Montana

This study was conducted to develop new estimating equations based on channel width and the updated flood frequency curves of previous investigations. Simple regression equations for estimating peak discharges with recurrence intervals of 2, 5, 10 , 25, 50, and 100 years were developed for seven regions in Montana. The standard errors of estimates for the equations that use active channel width as the independent variables ranged from 30% to 87%. The standard errors of estimate for the equations that use bankfull width as the independent variable ranged from 34% to 92%. The smallest standard errors generally occurred in the prediction equations for the 2-yr flood, 5-yr flood, and 10-yr flood, and the largest standard errors occurred in the prediction equations for the 100-yr flood. The equations that use active channel width and the equations that use bankfull width were determined to be about equally reliable in five regions. In the West Region, the equations that use bankfull width were slightly more reliable than those based on active channel width, whereas in the East-Central Region the equations that use active channel width were slightly more reliable than those based on bankfull width. Compared with similar equations previously developed, the standard errors of estimate for the new equations are substantially smaller in three regions and substantially larger in two regions. Limitations on the use of the estimating equations include: (1) The equations are based on stable conditions of channel geometry and prevailing water and sediment discharge; (2) The measurement of channel width requires a site visit, preferably by a person with experience in the method, and involves appreciable measurement errors; (3) Reliability of results from the equations for channel widths beyond the range of definition is unknown. In spite of the limitations, the estimating equations derived in this study are considered to be as reliable as estimating equations based on basin and climatic variables. Because the two types of estimating equations are independent, results from each can be weighted inversely proportional to their variances, and averaged. The weighted average estimate has a variance less than either individual estimate. (Author 's abstract)

Water-Resources Investigations Report↗

User's guide to revised method-of-characteristics solute-transport model (MOC--version 31)

The U.S. Geological Survey computer model to simulate two-dimensional solute transport and dispersion in ground water (Konikow and Bredehoeft, 1978; Goode and Konikow, 1989) has been modified to improve management of input and output data and to provide progressive run-time information. All opening and closing of files are now done automatically by the program. Names of input data files are entered either interactively or using a batch-mode script file. Names of output files, created automatically by the program, are based on the name of the input file. In the interactive mode, messages are written to the screen during execution to allow the user to monitor the status and progress of the simulation and to anticipate total running time. Information reported and updated during a simulation include the current pumping period and time step, number of particle moves, and percentage completion of the current time step. The batch mode enables a user to run a series of simulations consecutively, without additional control. A report of the model's activity in the batch mode is written to a separate output file, allowing later review. The user has several options for creating separate output files for different types of data. The formats are compatible with many commercially available applications, which facilitates graphical postprocessing of model results. Geohydrology and Evaluation of Stream-Aquifer Relations in the Apalachicola-Chattahoochee-Flint River Basin, Southeastern Alabama, Northwestern Florida, and Southwestern Georgia By Lynn J. Torak, Gary S. Davis, George A. Strain, and Jennifer G. Herndon Abstract The lower Apalachieola-Chattahoochec-Flint River Basin is underlain by Coastal Plain sediments of pre-Cretaceous to Quaternary age consisting of alternating units of sand, clay, sandstone, dolomite, and limestone that gradually thicken and dip gently to the southeast. The stream-aquifer system consism of carbonate (limestone and dolomite) and elastic sediments, which define the Upper Floridan aquifer and Intermediate system, in hydraulic connection with the principal rivers of the basin and other surface-water features, natural and man made. Separate digital models of the Upper Flori-dan aquifer and Intermediate system were constructed by using the U.S. Geological Survey's MODular Finite-Element model of two dimensional ground-water flow, based on concep- tualizations of the stream-aquifer system, and calibrated to drought conditions of October 1986. Sensitivity analyses performed on the models indicated that aquifer hydraulic conductivity, lateral and vertical boundary flows, and pumpage have a strong influence on groundwater levels. Simulated pumpage increases in the Upper Floridan aquifer, primarily in the Dougherty Plain physiographic district of Georgia,. caused significant reductions in aquifer discharge to streams that eventually flow to Lake Seminole and the Apalachicola River and Bay. Simulated pumpage increases greater than 3 times the October 1986 rates caused drying ofsome stream reaches and parts of the Upper Floridan aquifer in Georgia. Water budgets prepared from simulation results indicate that ground- water discharge to streams and recharge by horizontal and vertical flow are the principal mechanisms for moving water through the flow system. The potential for changes in ground-water quality is high in areas where chemical constituents can be mobilized by these mechanisms. Less than 2 percent of ground-water discharge to streams comes from the Intermediate system; thus, it plays a minor role in the hydrodynamics of the stream- aquifer system.

Water-Resources Investigations Report↗

Fate and behavior tools related to inland spill response—Workshop on the U.S. Geological Survey’s role in Federal science support

Executive Summary There is a growing body of tools available for science support for determining the fate and behavior of industrial and agricultural chemicals that are rapidly injected (“spilled”) into aquatic environments. A 2-day roundtable-style workshop was held by the U.S. Geological Survey (USGS) in Middleton, Wisconsin, in December 2017 to describe and explore existing Federal science support for spill fate and behavior tools used for inland spills, ongoing and new fate and behavior studies, and science gaps in planning and response tools as part of the USGS Midcontinent Region’s efforts to include spill response as part of its strategic plans. A total of 28 attendees representing a variety of Federal, State, and regional entities presented on programs and tools used in various aspects of spill response. Most programs and tools discussed were for spills in riverine environments but tools and applications for spills in lakes, on land surfaces, in urban storm sewer networks, and groundwater also were discussed. A primary workshop focus was to facilitate communication and increase potential for future collaboration among agencies for inland spill science support. The role and need for more USGS science support within the inland spill community was discussed. Enhanced communication is needed within the USGS and the U.S. Department of the Interior science programs, as well as within and among other agencies that do emergency planning and response. A main conclusion of the workshop was that there are untapped resources of the USGS outlined in the agency’s science strategy that could strengthen science support for fate and behavior tools in inland areas, especially in the Upper Mississippi River, Ohio River, and Great Lakes Basins where large freshwater resources overlap with dense corridors of oil and hazardous substances, with transportation networks, and with large populations centers. Fate and behavior tools are being developed quickly for inland spill response by multiple Federal agencies in partnership with local and regional entities. Applicability of these tools ranges from planning and preparedness, to the early stages of spill response for protection of human life and property, and to the application of monitoring and models to assess the long-term consequences of spills. Key findings from the workshop, with an emphasis on potential further development of USGS science support, include the following: •The national and regional response to spills occurs within an established system that must be respected by all parties involved in spill response. The USGS’s role is to support spill responders who are physically working at a spill scene, deploying booms and using other efforts to contain and recover spilled materials. •The USGS has tools that have been used throughout spill response operations, from early response to recovery and restoration. Developing a more formal role for the USGS to participate in science support for inland spills on a consistent basis is a desired outcome. This will require the USGS to improve internal and external communication and would be best accomplished by assigning one or more coordinator positions within the agency to plan and oversee USGS spill-response efforts. More involvement of the USGS on National and Regional Response Teams, especially in the realm of the Science and Technology Subcommittees, will gofar in increasing external communication and integration of fate and behavior tools. •Rapid response to spills requires modeling and mapping of plumes and associated time-of-travel estimation for a range of stream sizes across the United States. Many existing models use USGS streamgage data and the USGS National Hydrography Dataset. Nearly all existing models would benefit from updated linkages to USGS StreamStats and its soon-to-be released time-of-travel estimates,real-time velocity, stream morphology, and slope data. Integrating USGS tools with those from other agencies could be done to better serve the larger spill response community. • A problem is that existing models to rapidly predict plume extent, as well as more followup/longer-term fate and transport models, can be unknown or unavailable to spill responders. Thus, creating and strengthening linkages among USGS scientists skilled at using these tools is needed to support spill response with the on-scene responders. • Research for inland spill fate and behavior done outside of an immediate spill response can assist with spill planning and preparedness by (1) revealing sites likely to experience spills in the future (high-risk sites) and (2) understanding how a spilled substance might behave under a range of environmental conditions. However, USGS research on this topic has been scarce and subject to funding availability. Examples include the 2010 Line 6B Spill release into the Kalamazoo River in Michigan, where the USGS provided science support for a variety of fate and behavior tools for stream and impoundment environments. A long-term research site in Bemidji, Minnesota, provides important insights into transformations and longevity of spilled oil in groundwater and groundwater-surface water interactions. • Linking stream models to other components of this inland environment, including groundwater, overland flow, and karst, is needed. Stream network data can be linked to underground conduits such as storm sewers and karst groundwater systems. Stream models can also be linked with geospatial data such as that contained in U.S. Environmental Protection Agency’s interactive mapping tools. • The USGS is uniquely qualified to collect water-quality data during spills in the United States because of its many geographically dispersed water science centers, its knowledge and preparedness for flood measurement and documentation, and its cadre of skilled water-quality employees. Rapid-deployment gages, used for floods, could also be used for spills if they included spill-specific sensors. Coordinated expertise at USGS water and environmental science centers can be used for monitoring spill effects and for assessing risk to water quality and ecological communities. • Scientists at the USGS have proven capable of providing science coordination and technical assistance within the Incident Command Structure at the request of the lead on-scene coordinator. This external coordination, as well as internal communication within USGS Water, Hazards, and Ecosystems Mission Areas, could be improved by establishing and naming a USGS spills coordinator. Scott Morlock, Jo Ellen Hinck, and Faith Fitzpatrick are currently (2017) serving in informal coordination roles in addition to their traditional duties.

Open-File Report↗

Hydrogeology of the alluvial aquifers at the Pueblo Depot Activity near Pueblo, Colorado

In 1992, the U.S. Army Corps of Engineers and the U.S. Army Pueblo Depot Activity requested that the U.S. Geological Survey study the hydrogeology of the Pueblo Depot Activity, a 36-square-mile facility that has been operated by the U.S. Army since 1942. The purpose of the study was to provide an updated hydrogcological framework to facilitate the investigation of specific sites on the facility that might require remediation. This report describes the hydrogeology of the alluvial aquifers beneath the facility and the distribution of specific conductance of ground water in those aquifers. The Pueblo Depot Activity is underlain by two alluvial aquifers: (1) The terrace alluvial aquifer, which is a southernmost, downgradicnt part of an erosional remnant of an extensive terrace deposit; and (2) the Chico Creek alluvial aquifer, a smaller alluvial system along Chico Creek. These aquifers primarily consist of sand separated by clay layers and are underlain by the almost impermeable Pierre Shale of Upper Cretaceous age. The bedrock surface, which has an average slope of 28 feet per mile to the south-southeast, is relatively regular beneath the northern two-thirds of the terrace deposits at the Pueblo Depot Activity, but forms an irregular surface of troughs, hills, and ridges in the southwestern part of the terrace alluvium. Saturated thickness of the terrace aquifer ranges from 0 to about 45 feet. The bedrock surface beneath the Chico Creek aquifer slopes about 31 feet per mile to the south. Saturated thickness of the Chico Creek alluvium ranges from 0 to about 30 ft, but generally is less than 15 ft. Total thickness of the Chico Creek alluvium in the saturated area ranges from 16 to 41 ft. Water in the terrace alluvial aquifer generally flows southward, except in the southwestern part where directions of flow are complex. Measured hydraulic conductivity ranges from 0.4 to 400 feet per day (median 26 feet per day). Estimates for vertically averaged ground-water-flow velocity range from 0.02 to 3 feet per day (median 0.9 foot per day). Water in the Chico Creek alluvial aquifer generally flows southward to the Arkansas River alluvium. Measured hydraulic conductivity ranges from 14 to 310 feet per day (median 42 feet per day). Estimates for vertically averaged ground water-flow velocity range from 0.5 to 4 feet per day (median 0.7 foot per day). Specific conductance of ground water in the terrace alluvial aquifer generally is less than 800 microsiemens per centimeter; the smallest values were observed in the north-central part of the Pueblo Depot Activity. In the southwestern part of the terrace alluvial aquifer, values varied in an irregular pattern, and values as large as 3,300 microsiemens per centimeter were measured locally. Water in the terrace alluvial aquifer was dominated by the sodium cation and usually by the bicarbonate anion, and sulfate usually was present in substantial (and locally predominant) concentrations. Measured specific conductance of water in the Chico Creek alluvial aquifer ranged from 683 to 1,460 microsiemens per centimeter. This water was dominated by the sodium cation and by the bicarbonate and sulfate anions; sulfate was more predominant to the south.

Colorado↗

Hydrology of the Helena area bedrock, west-central Montana, 1993-98

The Generalized Bedrock Geologic Map of the Helena Area, West-Central Montana (plate 1 in the report) provides an intermediate-scale overview of bedrock in the Helena area. The geologic map has been compiled at a scale of 1:100,000 from the most widely available sources of geologic map information (see index to geologic mapping on pl. 1). That information has been updated by M.W. Reynolds for this report with more recent geologic mapping and field revision of published maps. All well locations and all bedrock units penetrated during drilling have been confirmed on geologic maps at the largest scale available. Source geologic maps are all at scales larger than 1:100,000 scale. Care has been taken to ensure accurate representation of the original geology at the compilation scale. However, positional accuracy of some features might be somewhat diminished at the smaller scale of the base map when compared with the original data source. Also, line thicknesses for contacts and faults necessarily assume a greater width, relative to the real geologic feature, at the scale of the generalized map than on any original map. The map is not intended for large-scale, site-specific detailed planning. Bedrock units throughout the Helena area are generally covered by young surficial deposits such as alluvium, colluvium, glacial debris, or windblown sediment. Thickness of such deposits varies from veneers through which the underlying bedrock is clearly discernible to major thicknesses that conceal all underlying bedrock and structure. Boundaries of major accumulations of surficial deposits are attributed separately from bedrock contacts. These boundaries should not be considered precise at the map scale or at larger scales. Boundaries shown may be less accurate positionally than bedrock contacts and faults because (1) surficial deposits commonly thin to a knife edge; (2) different mappers will interpret the edge differently when drawing a boundary; or (3) the original geologic map maker was concerned principally with bedrock units and structure and thus overlooked, or did not originally map as consistently, some surficial deposits. Veneers of surficial sediment, when saturated, can be local sources of recharge to underlying bedrock. Use of the generalized map to define their distribution does not substitute for site specific mapping of such deposits. Specific knowledge is needed to determine the water-bearing properties of the geologic units at and surrounding a site because the units, including the igneous and metamorphic rocks, have internal differences in stratigraphy, composition, mineralogy and grain size or crystallinity. These differences, together with structural imprints such as faults, folds, and the spacing, orientation, degree of openness of fractures, and extent and type of mineral filling in fractures and faults, all affect the ability of rocks to store and transmit water.

Montana↗

Hydraulic properties and ground-water flow in the St Peter-Prairie du Chien-Jordan aquifer, Rochester area, southeastern Minnesota

The hydraulic properties were updated and their effects on ground-water flow in the St. Peter-Prairie du Chien-Jordan aquifer in the Rochester area in southeastern Minnesota were evaluated, using new information compiled since a study by Delin (1990). Since 1988, new information on the hydrogeology of the ground-water system in the Rochester area has become available from well-drilling and construction activity associated with Rochester's rapid growth. The St. Peter-Prairie du Chien-Jordan aquifer consists of the St. Peter Sandstone, the Prairie du Chien Group (limestones and dolomites), and the Jordan Sandstone. Horizontal hydraulic conductivity and transmissivity were determined from 15 aquifer tests and specific-capacity information compiled for 310 wells. A 140-square-mile area of the aquifer bounded on the west, south, and east by a ground-water divide contributes water to the Rochester, Minnesota, municipal wells. Transmissivities for the St. Peter-Prairie du Chien-Jordan aquifer in the study area range from less than 5,000 square feet per day (ft 2 /d) to greater than 20,000 ft 2 /d. Transmissivities greater than 20,000 ft 2 /d occur in the west-central, northwestern, and east-central parts of the study area. Transmissivities of less than 5,000 ft 2 /d occur in the northern, northeastern, central, and southern parts of the study area. The areas of greatest potential well yield coincide with areas of greatest transmissivity. Delin (1990) developed a ground-water-flow model to simulate flow of ground water in the St. Peter-Prairie du Chien-Jordan aquifer in the Rochester area. The 1988 Rochester model was rerun using revised horizontal hydraulic conductivity arrays in the model, based on the transmissivity distribution determined for this study. The results of the simulations using horizontal hydraulic conductivities based on the transmissivity distribution determined for this study may indicate that transmissivity values derived from specific-capacity information generally are too high. The transmissivity distribution determined for this study, however, is valid as an indicator of the spatial variability of the relative magnitude of transmissivity and potential well yield for the St. Peter-Prairie du Chien-Jordan aquifer in the study area. Water-level changes in wells from January through February 1988 to February through March 1995 ranged from -6.8 to +15.3 feet. Water-level changes in 12 Rochester municipal wells for the same period ranged from -7.4 to +8.0 feet. Water levels in wells generally rose in the northern and eastern parts of the study area and generally declined in the southwestern and western parts. Near Rochester, water levels in wells generally declined near the city boundaries and showed little change or rose in the central part of the city. Water-level changes from 1988 to 1995 near the ground-water divide generally were less than 2 feet, resulting in no appreciable changes in the location of the divide.

Minnesota↗

Techniques for estimating the magnitude and frequency of floods in rural basins of South Carolina, 1999

Data from 167 streamflow-gaging stations in or near South Carolina with 10 or more years of record through September 30, 1999, were used to develop two methods for estimating the magnitude and frequency of floods in South Carolina for rural ungaged basins that are not significantly affected by regulation. Floodfrequency estimates for 54 gaged sites in South Carolina were computed by fitting the water-year peak flows for each site to a log-Pearson Type III distribution. As part of the computation of flood-frequency estimates for gaged sites, new values for generalized skew coefficients were developed. Flood-frequency analyses also were made for gaging stations that drain basins from more than one physiographic province. The U.S. Geological Survey, in cooperation with the South Carolina Department of Transportation, updated these data from previous flood-frequency reports to aid officials who are active in floodplain management as well as those who design bridges, culverts, and levees, or other structures near streams where flooding is likely to occur. Regional regression analysis, using generalized least squares regression, was used to develop a set of predictive equations that can be used to estimate the 2-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence- interval flows for rural ungaged basins in the Blue Ridge, Piedmont, upper Coastal Plain, and lower Coastal Plain physiographic provinces of South Carolina. The predictive equations are all functions of drainage area. Average errors of prediction for these regression equations ranged from -16 to 19 percent for the 2-year recurrence-interval flow in the upper Coastal Plain to -34 to 52 percent for the 500-year recurrenceinterval flow in the lower Coastal Plain. A region-of-influence method also was developed that interactively estimates recurrence- interval flows for rural ungaged basins in the Blue Ridge of South Carolina. The region-of-influence method uses regression techniques to develop a unique relation between flow and basin characteristics for an individual watershed. This, then, can be used to estimate flows at ungaged sites. Because the computations required for this method are somewhat complex, a computer application was developed that performs the computations and compares the predictive errors for this method. The computer application includes the option of using the region-of-influence method, or the generalized least squares regression equations from this report to compute estimated flows and errors of prediction specific to each ungaged site. From a comparison of predictive errors using the region-of-influence method with those computed using the regional regression method, the region-of-influence method performed systematically better only in the Blue Ridge and is, therefore, not recommended for use in the other physiographic provinces. Peak-flow data for the South Carolina stations used in the regionalization study are provided in appendix A, which contains gaging station information, log- Pearson Type III statistics, information on stage-flow relations, and water-year peak stages and flows. For informational purposes, water-year peak-flow data for stations on regulated streams in South Carolina also are provided in appendix D. Other information pertaining to the regulated streams is provided in the text of the report.

South Carolina↗

Simulation of the interaction of karstic lakes Magnolia and Brooklyn with the upper Floridan Aquifer, southwestern Clay County, Florida

The stage of Lake Brooklyn, in southwestern Clay County, Florida, has varied over a range of 27 feet since measurements by the U.S. Geological Survey began in July 1957. The large stage changes have been attributed to the relation between highly transient surface-water inflow to the lake and subsurface conduits of karstic origin that permit a high rate of leakage from the lake to the Upper Floridan aquifer. After the most recent and severe stage decline (1990-1994), the U.S. Geological Survey began a study that entailed the use of numerical ground-water flow models to simulate the interaction of the lake with the Upper Floridan aquifer and the large fluctuations of stage that were a part of that process. A package (set of computer programs) designed to represent lake/aquifer interaction in the U.S. Geological Survey Modular Finite-Difference Ground-Water Flow Model (MODFLOW-96) and the Three-Dimensional Method-of-Characteristics Solute-Transport Model (MOC3D) simulators was prepared as part of this study, and a demonstration of its capability was a primary objective of the study. (Although the official names are Brooklyn Lake and Magnolia Lake (Florida Geographic Names), in this report the local names, Lake Brooklyn and Lake Magnolia, are used.) In the simulator of lake/aquifer interaction used in this investigation, the stage of each lake in a simulation is updated in successive time steps by a budget process that takes into account ground-water seepage, precipitation upon and evaporation from the lake surface, stream inflows and outflows, overland runoff inflows, and augmentation or depletion by artificial means. The simulator was given the capability to simulate both the division of a lake into separate pools as lake stage falls and the coalescence of several pools into a single lake as the stage rises. This representational capability was required to simulate Lake Brooklyn, which can divide into as many as 10 separate pools at sufficiently low stage. In the first of two calibrated models, recharge to the water table, specified as a monthly rate, was set equal to 40 percent of the monthly rainfall rate. The specified rate of inflow to the uppermost stream segment was set equal to outflows from Lake Lowry estimated from lake stage and the 1994-97 rating table. Leakage to the intermediate and Upper Floridan aquifers was assumed to occur from the surficial aquifer system through the confining layers directly beneath deeper parts of the lake bottom. A leakance coefficient value of 0.001 feet per day per foot of thickness was used beneath Lake Magnolia, and a value of 0.005 feet per day per foot of thickness was used beneath most of Lake Brooklyn. With these values, the conductance through the confining layers beneath Lake Brooklyn was about 19 times that beneath Lake Magnolia. The simulated stages of Lake Brooklyn matched the measured stages reasonably well in the early (1957-72) and later (1990-98) parts of the simulation time period, but the match was unsatisfactory in an intermediate time period (1973-89). To resolve this discrepancy, the hypothesis was proposed that undocumented losses of water from Alligator Creek upstream from Lake Brooklyn or from the lake itself occurred between 1973 and 1989 when there was sufficient streamflow. The resulting simulation of lake stages matched the measured lake stages accurately during the entire simulation time period. The model was then revised to incorporate the assumption that only 20 percent of precipitation recharged the water table (the second calibrated model). Recalibration of the model required that leakance values for the confining units under deeper parts of the lakes also be reduced by nearly 50 percent. The stages simulated with the new parameter assumptions, but retaining the assumption of surface-water losses, were an excellent match of the measured values. The stage of Lake Magnolia was also simulated accurately. The results of sensitivity analyses show that simulated streamflow between Lakes Magnolia and Brooklyn tends to be water-budget controlled, and is not appreciably affected by the specified outflow altitude or channel characteristics of the receiving stream. To match heads measured in observation wells of the surficial aquifer network, the assigned hydraulic conductivity values were zoned, and ranged from a minimum of 4 feet per day to a maximum of 400 feet per day in the first calibrated model. These values were reduced by about 50 percent in the second calibrated model. Differences between observation wells were noted in the abruptness of changes of measured head values, and in the relation of the timing of peak measured heads and simulated peak heads. These differences seemed to be correlated with the depth of the water table below land surface. Spatially uniform values of transmissivity were specified for the intermediate (10,000 feet squared per day) and Upper Floridan (100,000 feet squared per day) aquifers. Simulated heads in the Upper Floridan aquifer layer follow the trend of the heads measured in a long-term observation well with data beginning in 1960. This result suggests that the observed head decline could be explained entirely in terms of the stage decline in Lake Brooklyn and may not indicate a regional trend.

Florida↗

Methods for estimating flood frequency in Montana based on data through water year 1998

Annual peak discharges having recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years (T-year floods) were determined for 660 gaged sites in Montana and in adjacent areas of Idaho, Wyoming, and Canada, based on data through water year 1998. The updated flood-frequency information was subsequently used in regression analyses, either ordinary or generalized least squares, to develop equations relating T-year floods to various basin and climatic characteristics, equations relating T-year floods to active-channel width, and equations relating T-year floods to bankfull width. The equations can be used to estimate flood frequency at ungaged sites. Montana was divided into eight regions, within which flood characteristics were considered to be reasonably homogeneous, and the three sets of regression equations were developed for each region. A measure of the overall reliability of the regression equations is the average standard error of prediction. The average standard errors of prediction for the equations based on basin and climatic characteristics ranged from 37.4 percent to 134.1 percent. Average standard errors of prediction for the equations based on active-channel width ranged from 57.2 percent to 141.3 percent. Average standard errors of prediction for the equations based on bankfull width ranged from 63.1 percent to 155.5 percent. In most regions, the equations based on basin and climatic characteristics generally had smaller average standard errors of prediction than equations based on active-channel or bankfull width. An exception was the Southeast Plains Region, where all equations based on active-channel width had smaller average standard errors of prediction than equations based on basin and climatic characteristics or bankfull width. Methods for weighting estimates derived from the basin- and climatic-characteristic equations and the channel-width equations also were developed. The weights were based on the cross correlation of residuals from the different methods and the average standard errors of prediction. When all three methods were combined, the average standard errors of prediction ranged from 37.4 percent to 120.2 percent. Weighting of estimates reduced the standard errors of prediction for all T-year flood estimates in four regions, reduced the standard errors of prediction for some T-year flood estimates in two regions, and provided no reduction in average standard error of prediction in two regions. A computer program for solving the regression equations, weighting estimates, and determining reliability of individual estimates was developed and placed on the USGS Montana District World Wide Web page. A new regression method, termed Region of Influence regression, also was tested. Test results indicated that the Region of Influence method was not as reliable as the regional equations based on generalized least squares regression. Two additional methods for estimating flood frequency at ungaged sites located on the same streams as gaged sites also are described. The first method, based on a drainage-area-ratio adjustment, is intended for use on streams where the ungaged site of interest is located near a gaged site. The second method, based on interpolation between gaged sites, is intended for use on streams that have two or more streamflow-gaging stations.

Montana↗