Search USGS⌕ Search

SEARCH · Search USGS

Results for “Scientific Investigations Map”

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 1,477 records · Page 82Linked to original sources
USGS2008

Recovery of Ground-Water Levels from 1988 to 2003 and Analysis of Effects of 2003 and Full-Allocation Withdrawals in Critical Area 2, Southern New Jersey

Water levels in the Potomac-Raritan-Magothy aquifer system within Water Supply Critical Area 2 in the southern New Jersey Coastal Plain have recovered as a result of reductions in ground-water withdrawals initiated in the early 1990s. The Critical Area consists of the depleted zone and the threatened margin. The Potomac-Raritan-Magothy aquifer system consists of the Upper, Middle, and Lower aquifers. Generally, ground-water withdrawals from these aquifers declined 5 to 10 Mgal/d (million gallons per day) and water levels recovered 0 to 40 ft (foot) from 1988 to 2003. In order to reevaluate water-allocation restrictions in Critical Area 2 in response to changes in the ground-water-flow system and demands for additional water supply due to increased development, the New Jersey Department of Environmental Protection (NJDEP) needs information about the effects of changes in those allocations. Therefore, the U.S. Geological Survey (USGS), in cooperation with the NJDEP, used an existing ground-water-flow model of the New Jersey Coastal Plain to evaluate the effects of withdrawal alternatives on hydraulic heads in the Potomac-Raritan-Magothy aquifer system in Critical Area 2. The U.S. Geological Survey Regional Aquifer System Analysis model was used to simulate steady-state ground-water flow. Two withdrawal conditions were tested by using the model to evaluate hydraulic heads and differences in heads in these aquifers: 2003 withdrawals and full-allocation withdrawals (17.4 Mgal/d greater than 2003 withdrawals). Model results are presented using head maps and head-difference maps that compare 2003 to full-allocation withdrawals. Mandated hydrologic conditions for Critical Area protection are that the simulated -30-ft head contour not extend beyond the boundary of the depleted zone and (or) be at least 5 mi (miles) updip from the 250-mg/L (milligram per liter) isochlor in all three aquifers. Simulation results indicate that, for 2003 withdrawals, the simulated -30-ft head contour in all three aquifers is generally within the boundary of the depleted zone, except in the Lower aquifer in northern Camden and northwestern Burlington Counties, and is generally 1 to 10 mi downdip from the 250-mg/L isochlor. (Corresponding observed data indicate that the -30-ft water-level contour extends beyond the southwest boundary of the depleted zone in the Upper and Middle aquifers, and is generally 5 to 20 mi downdip from the 250-mg/L isochlor in all three aquifers.) The area in which heads are below -30 ft ranges from 389 mi2 (square miles) in the Middle aquifer to 427 mi2 in the Lower aquifer. For full-allocation withdrawals, the simulated -30-ft head contour extends beyond the boundary of the depleted zone in all three aquifers in northern Camden and northwestern Burlington Counties and in the Upper aquifer in Gloucester and Salem Counties, and is generally 5 to 15 mi downdip from the 250-mg/L isochlor. The area in which heads are below -30 ft ranges from 616 mi2 in the Upper aquifer to 813 mi2 in the Lower aquifer. These results and observed data indicate that any increase in withdrawals from 2003 values would likely cause heads in the three aquifers to decline below the minimum values mandated by the NJDEP for the Critical Area.

Scientific Investigations Report↗
USGS2021

A comparison of Landsat 8 Operational Land Imager and Provisional Aquatic Reflectance science product, Sentinel–2B, and WorldView–3 imagery for empirical satellite-derived bathymetry, Unalakleet, Alaska

Satellite-derived bathymetry (SDB) based upon an empirical band ratio method is a cost-effective means for mapping nearshore bathymetry in coastal areas vulnerable to natural hazards. This is particularly important for the low-lying coastal community of Unalakleet, Alaska, that has been negatively affected not only by flooding, storm surge, and historically strong storms but also by high erosion rates stemming from the Unalakleet River and Norton Sound. The purpose of this study was to assess the viability of different satellite imagery, including Landsat 8 (L8) Operational Land Imager, Sentinel–2B, WorldView–3, and L8 Provisional Aquatic Reflectance science product, for deriving SDB for Unalakleet, Alaska. Correlations were performed between satellite imagery band ratios and topobathymetric (topobathy) light detection and ranging (lidar) and in situ single-beam sound navigation and ranging (sonar). The satellite imagery correlations with topobathy lidar did not yield as high of a linear relation with water depths as the satellite imagery correlations with the single-beam sonar. An extinction depth, where light no longer attenuates through the water column, was not identified because of the shallow depths within the topobathy lidar and single-beam sonar datasets. Although some single-beam soundings measured at 7 meters deep, the correlations with the SDB band ratios did not yield a strong linear relation. Satellite imagery band ratio correlations with Electronic Navigational Chart soundings did not yield a strong linear relation because of older source data. Less than optimal linear regressions were most likely due to the geography of Unalakleet, Alaska, a low-lying coastal community subject to high erosion rates from surrounding waters. This study is one of the first attempts to compare different satellite imagery band ratio correlations with topobathy lidar and in situ sonar to assess the viability for nearshore SDB for coastal Unalakleet, Alaska.

Alaska↗
USGS2006

Proceedings of the U.S. Geological Survey Sixth Biennial Geographic Information Science Workshop, Denver, Colorado, April 24-28, 2006

The U.S. Geological Survey's (USGS) Sixth Biennial Geographic Information Science Workshop April 24 - 28, 2006, at the Denver Federal Center in Denver, Colorado, provides a unique opportunity for multi-disciplinary Geographic Information Systems (GIS) and associated scientific professionals to share, learn, present, and discuss a wide variety of geospatial-related topics. Information is exchanged through a series of plenary sessions, hands-on technical workshops, user and commercial vendor demonstrations, lecture and poster sessions, and specialty meetings. Workshop attendance is limited to Department of Interior (DOI) and USGS employees and contractors; however, plenary, hands-on, and lecture sessions have presenters from other Federal agencies, numerous commercial vendors, universities, and several consortia. Over 175 participants are expected to attend the Workshop. Several prominent speakers are featured at this Workshop. Monday evening Star Guest Speaker and National Aeronautics and Space Administration (NASA) Astronaut Captain Dominic Gorie will talk about his experiences as a veteran of three space flights and over 32 days in space, including the NASA Space Shuttle Radar Topography Mission that mapped more than 47 million miles of the Earth's land surface. Selected as an astronaut candidate by NASA in December 1994, Captain Gorie is currently Chief of the Astronaut Shuttle Branch. Monday evening also features a town hall meeting with Geographic Information Office (GIO) leaders Karen Siderelis, Kevin Gallagher, Bob Pierce, Steve Guptill, Mark DeMulder, John Mahoney, and Mark Negri, who will discuss changes and activities within the GIO in an open discussion format. Tuesday plenary sessions feature keynote speaker Dr. P. Patrick Leahy, Acting USGS Director. Dr. Leahy holds undergraduate and graduate degrees in geology (1968) and geophysics (1970) from Boston College. In 1979, he received his doctorate in geology from Rensselaer Polytechnic Institute, where he specialized in regional ground-water studies and hydraulics. Dr. Leahy has been with the USGS since 1974, and has held various technical and managerial positions, including Associate Director for Geology and Chief of the National Water-Quality Assessment (NAWQA) Program. Dr. Leahy will be discussing a broad range of topics including Bureau accomplishments, initiatives, and budgets. The purpose of this proceedings volume is to serve as an activity reference for Workshop attendees as well as an archive of technical abstracts submitted, presented, and discussed at the Workshop. Author, co-author, and presenter names, affiliations, and contact information are listed with presentation titles along with submitted abstracts. Some hands-on sessions are offered twice. In these instances, abstracts submitted for publication are presented in the proceedings on both days they are offered. All acronyms used in these proceedings are explained in the text of each abstract. The term 'ArcGIS' refers to an integrated collection of GIS software products produced by Environmental Systems Research Institute, Inc. (ESRI).

Scientific Investigations Report↗
USGS2007

Sedimentation History of Halfway Creek Marsh, Upper Mississippi River National Wildlife and Fish Refuge, Wisconsin, 1846-2006

The history of overbank sedimentation in the vicinity of Halfway Creek Marsh near La Crosse, Wis., was examined during 2005?06 by the U.S. Geological Survey and University of Wisconsin?Madison as part of a broader study of sediment and nutrient loadings to the Upper Mississippi River bottomlands by the U.S. Environmental Protection Agency, U.S. Fish and Wildlife Service, and U.S. Geological Survey. Historical sedimentation patterns and rates were interpreted from field-scale topographic surveys and sediment cores collected from the marsh and upstream flood plains. Historical maps and aerial photographs were used to establish the timing of disturbances and to document changes in channel patterns after Euro-American settlement (post 1846). Episodic overbank sedimentation patterns and rates were linked to watershed agricultural activity, large floods, artificial levee construction, channel alterations, and dam failures over the past 160 years. These forces affected sedimentation on and between levees, the development of alluvial fans and flood-plain splays, and the general pattern of flood-plain sedimentation through the upper and lower marsh. Historical overbank deposits, episodically deposited after about 1860, are as much as 6 feet thick in the upper marsh and as much as 4 feet thick in the lower marsh, representing a total volume of approximately 1.8 million cubic yards. These stratified deposits consist of multiple layers of silt and clay, very fine to fine sand, and some medium to very coarse sand. Coarse-grained deposits are associated with flood-plain splays caused by breaches in artificial levees during large floods. Estimated sedimentation rates were highest from 1919 to 1936 [26,890 cubic yards per year (yd3/yr)] and exceeded by about 30 times the 1846?85 rate of 920 yd3/yr and exceeded by 7 times the 1994?2006 rate of 3,740 yd3/yr. The 1994?2006 sedimentation rate was the lowest since Euro-American settlement, but natural levees along the 1994?2006 channel of Halfway Creek through the lower marsh continued to form and are currently (2006) about 1 foot higher than the surrounding marsh. Natural levee building in the lower marsh from 1994?2006 was accentuated by the lack of overbank sediment storage in the upper marsh. The historical storage of sediment in the upper and lower marsh affects modern streamflow and sediment transport processes of Halfway Creek and Sand Lake Coulee through the marsh, and it also affects marsh vegetation and wildlife habitat. Results from this investigation will help improve the understanding of how past overbank sedimentation patterns continue to influence modern and future water quality, sediment transport, nutrient loads, and water-related resources in riparian habitats common to the Upper Mississippi River National Wildlife and Fish Refuge.

Scientific Investigations Report↗
USGS2012

Sandstone copper assessment of the Chu-Sarysu Basin, Central Kazakhstan: Chapter E in Global mineral resource assessment

Mineral resource assessments represent a synthesis of available information to estimate the location, quality, and quantity of undiscovered mineral resources in the upper part of the Earth’s crust. This report presents a probabilistic mineral resource assessment of undiscovered sandstone copper deposits within the late Paleozoic Chu-Sarysu Basin in central Kazakhstan by the U.S. Geological Survey as a contribution to a global assessment of mineral resources. The purposes of this study are to: (1) provide a database of known sandstone copper deposits and significant prospects in this area, (2) delineate permissive areas (tracts) for undiscovered sandstone copper deposits within 2 km of the surface at a scale of 1:1,000,000, (3) estimate numbers of undiscovered deposits within these permissive tracts at several levels of confidence, and (4) provide probabilistic estimates of amounts of copper (Cu), silver (Ag), and mineralized rock that could be contained in undiscovered deposits within each tract. The assessment uses the three-part form of mineral resource assessment based on mineral deposit models (Singer, 1993; Singer and Menzie, 2010). Delineation of permissive tracts for resources is based on the distribution of a Carboniferous oxidized nonmarine clastic (red bed) stratigraphic sequence that lies between overlying Permian and underlying Devonian evaporite-bearing sequences. Subsurface information on the extent and depth of this red bed sequence and structural features that divide the basin into sub-basins was used to define four permissive tracts. Structure contour maps, mineral occurrence databases, drill hole lithologic logs, geophysical maps, soil geochemical maps, locations of producing gas fields, and evidence for former gas accumulations were considered in conjunction with descriptive deposit models and grade and tonnage models to guide the assessment team’s estimates of undiscovered deposits in each tract. The four permissive tracts are structural sub-basins of the Chu-Sarysu Basin and range in size from 750 to 65,000 km². Probabilistic estimates of numbers of undiscovered sandstone copper deposits were made for the four tracts by a group of experts. Using these probabilistic estimates, Monte Carlo simulation was used to estimate the amount of metal contained within each tract. The results of the simulation serve as the basis for estimates of the metal endowment. The team estimates that 26 undiscovered deposits occur within the Chu-Sarysu Basin, and that these deposits contain an arithmetic mean of at least 21.5 million metric tons (Mt) of copper and 21,900 metric tons (t) of silver. The undiscovered deposits are in addition to the 7 known deposits that contain identified resources of 27.6 Mt of copper. Sixty percent of the estimated mean undiscovered copper resources are associated with the two permissive tracts that contain the identified resources; the remaining estimated resources are associated with the two tracts that have no known deposits. For the three tracts that contain 95 percent of the estimated undiscovered copper resources, the probability that each tract contains its estimated mean or more is about 40 percent. For the tract with 5 percent of the estimated undiscovered cop-per resources, the probability that it contains that amount or more is 25 percent.

Chu-Sarysu Basin↗
USGS2012

Isotropic, anisotropic, and borehole washout analyses in Gulf of Mexico Gas Hydrate Joint Industry Project Leg II, Alaminos Canyon well 21-A

Through the use of three-dimensional seismic amplitude mapping, several gas hydrate prospects were identified in the Alaminos Canyon area of the Gulf of Mexico. Two of the prospects were drilled as part of the Gulf of Mexico Gas Hydrate Joint Industry Program Leg II in May 2009, and a suite of logging-while-drilling logs was acquired at each well site. Logging-while-drilling logs at the Alaminos Canyon 21–A site indicate that resistivities of approximately 2 ohm-meter and P-wave velocities of approximately 1.9 kilometers per second were measured in a possible gas-hydrate-bearing target sand interval between 540 and 632 feet below the sea floor. These values are slightly elevated relative to those measured in the hydrate-free sediment surrounding the sands. The initial well log analysis is inconclusive in determining the presence of gas hydrate in the logged sand interval, mainly because large washouts in the target interval degraded well log measurements. To assess gas-hydrate saturations, a method of compensating for the effect of washouts on the resistivity and acoustic velocities is required. To meet this need, a method is presented that models the washed-out portion of the borehole as a vertical layer filled with seawater (drilling fluid). Owing to the anisotropic nature of this geometry, the apparent anisotropic resistivities and velocities caused by the vertical layer are used to correct measured log values. By incorporating the conventional marine seismic data into the well log analysis of the washout-corrected well logs, the gas-hydrate saturation at well site AC21–A was estimated to be in the range of 13 percent. Because gas hydrates in the vertical fractures were observed, anisotropic rock physics models were also applied to estimate gas-hydrate saturations.

Scientific Investigations Report↗
USGS2007

Relation of Lake-Floor Characteristics to the Distribution of Variable Leaf Water-Milfoil in Moultonborough Bay, Lake Winnipesaukee, New Hampshire, 2005

Geophysical, water, and sediment surveys were done to characterize the effects of surficial geology, water and sediment chemistry, and surficial-sediment composition on the distribution of variable leaf water-milfoil in Moultonborough Bay, Lake Winnipesaukee, New Hampshire. Geophysical surveys were conducted in a 180-square-kilometer area, and water-quality and sediment samples were collected from 24 sites in the survey area during July 2005. Swath-bathymetric data revealed that Moultonborough Bay ranged in depth from less than 1 meter (m) to about 15 m and contained three embayments. Seismic-reflection profiles revealed erosion of the underlying bedrock and subsequent deposition of glaciolacustrine and Holocene lacustrine sediments within the survey area. Sediment thickness ranged from 5 m along the shoreward margins to more than 15 m in the embayments. Data from sidescan sonar, surficial-sediment samples, bottom photographs, and video revealed three distinct lake-floor environments: rocky nearshore, mixed nearshore, and muddy basin. Rocky nearshore environments were found in shallow water (less than 5 m deep) and contained sediments ranging from coarse silt to very coarse sand. Mixed nearshore environments also were found in shallow water and contained sediments ranging from silt to coarse sand with different densities of aquatic vegetation. Muddy basin environments contained the finest-grained sediments, ranging from fine to medium silt, and were in the deepest waters of the bay. Acoustic Ground Discrimination Systems (AGDS) survey data revealed that 86 percent of the littoral zone (the area along the margins of the bay and islands that extends from 0 to 4.3 m in water depth) contained submerged aquatic vegetation (SAV) in varying densities: approximately 36 percent contained SAV bottom cover of 25 percent or less, 43 percent contained SAV bottom cover of more than 25 and less than 75 percent, and approximately 7 percent contained SAV bottom cover of more than 75 percent. SAV included variable leaf water-milfoil, native milfoil, bassweed, pipewort, and other species, which were predominantly found near shoreward margins and at depths ranging from less than 1 to 4 m. AGDS data were used in a Geographic Information System to generate an interpolated map that distinguished variable leaf water-milfoil from other SAV. Furthermore, these data were used to isolate areas susceptible to variable leaf water-milfoil growth. Approximately 21 percent of the littoral zone contained dense beds (more than 59 percent bottom cover) of variable leaf water-milfoil, and an additional 44 percent was determined to be susceptible to variable leaf water-milfoil infestation. Depths differed significantly between sites with variable leaf water-milfoil and sites with other SAV (p = 0.04). Variable leaf water-milfoil was found at depths that ranged from 1 to 4 m, and other SAV had a depth range of 1 to 2 m. Although variable leaf water-milfoil was observed at greater depths than other SAV, it was not observed below the photic zone. Analysis of constituent concentrations from the water column, interstitial pore water, and sediment showed little correlation with the presence of variable leaf water-milfoil, with two exceptions. Iron concentrations were significantly lower at variable leaf water-milfoil sites than at other sampling sites (p = 0.04). Similarly, the percentage of total organic carbon also was significantly lower at the variable leaf water-milfoil sites than at other sampling sites (p = 0.04). Surficial-sediment-grain size had the greatest correlation to the presence of variable leaf water-milfoil. Variable leaf water-milfoil was predominantly growing in areas of coarse sand (median grain-size 0.62 millimeters). Surficial-sediment-grain size was also correlated with total ammonia plus organic nitrogen (Rho = 0.47; p = 0.02) and with total phosphorus (Rho = 0.44; p = 0.05) concentrations in interstitial pore-water samples.

Scientific Investigations Report↗
USGS2016

Groundwater-flow model of the northern High Plains aquifer in Colorado, Kansas, Nebraska, South Dakota, and Wyoming

The High Plains aquifer is a nationally important water resource underlying about 175,000 square miles in parts of eight states: Colorado, Kansas, Oklahoma, Nebraska, New Mexico, South Dakota, Texas, and Wyoming. Droughts across much of the Northern High Plains from 2001 to 2007 have combined with recent (2004) legislative mandates to elevate concerns regarding future availability of groundwater and the need for additional information to support science-based water-resource management. To address these needs, the U.S. Geological Survey began the High Plains Groundwater Availability Study to provide a tool for water-resource managers and other stakeholders to assess the status and availability of groundwater resources. A transient groundwater-flow model was constructed using the U.S. Geological Survey modular three-dimensional finite-difference groundwater-flow model with Newton-Rhapson solver (MODFLOW–NWT). The model uses an orthogonal grid of 565 rows and 795 columns, and each grid cell measures 3,281 feet per side, with one variably thick vertical layer, simulated as unconfined. Groundwater flow was simulated for two distinct periods: (1) the period before substantial groundwater withdrawals, or before about 1940, and (2) the period of increasing groundwater withdrawals from May 1940 through April 2009. A soil-water-balance model was used to estimate recharge from precipitation and groundwater withdrawals for irrigation. The soil-water-balance model uses spatially distributed soil and landscape properties with daily weather data and estimated historical land-cover maps to calculate spatial and temporal variations in potential recharge. Mean annual recharge estimated for 1940–49, early in the history of groundwater development, and 2000–2009, late in the history of groundwater development, was 3.3 and 3.5 inches per year, respectively. Primary model calibration was completed using statistical techniques through parameter estimation using the parameter estimation suite of software with Tikhonov regularization. Calibration targets for the groundwater model included 343,067 groundwater levels measured in wells and 10,820 estimated monthly stream base flows at streamgages. A total of 1,312 parameters were adjusted during calibration to improve the match between calibration targets and simulated equivalents. Comparison of calibration targets to simulated equivalents indicated that, at the regional scale, the model correctly reproduced groundwater levels and stream base flows for 1940–2009. This comparison indicates that the model can be used to examine the likely response of the aquifer system to potential future stresses. Mean calibrated recharge for 1940–49 and 2000–2009 was smaller than that estimated with the soil-water-balance model. This indicated that although the general spatial patterns of recharge estimated with the soil-water-balance model were approximately correct at the regional scale of the Northern High Plains aquifer, the soil-water-balance model had overestimated recharge, and adjustments were needed to decrease recharge to improve the match of the groundwater model to calibration targets. The largest components of the simulated groundwater budgets were recharge from precipitation, recharge from canal seepage, outflows to evapotranspiration, and outflows to stream base flow. Simulated outflows to irrigation wells increased from 7 percent of total outflows in 1940–49 to 38 percent of 1970–79 total outflows and 49 percent of 2000–2009 total outflows.

Colorado, Kansas, Nebraska, South Dakota, Wyoming↗
USGS2007

Investigation of wind and water level for the Giacomini Wetland Restoration Project, Point Reyes National Seashore

Point Reyes National Seashore (PRNS), comprising unique elements of geological, biological, and historical interest, is located on the central California coast approximately 60 km northwest of San Francisco. The National Seashore contains nearly 130 km of exposed and protected shorelines, spectacular coastal cliffs and headlands, lagoons, open grasslands, bushy hillsides, and forested ridges. Approximately 30 km of the shoreline are coastal-dune habitat that supports 11 federally listed species, including the threatened western snowy plover and the endangered plants Tidestrom's lupine ( Lupinus tidestromii ) and beach layia ( Layia carnosa ). The San Andreas Fault, a right-lateral strike-slip fault, trends northwest along the northeastern side of the park. Tomales Bay, which is straight, long, narrow, and shallow, runs along the northeastern boundary of PRNS. The Bay, which fills the northwestern end of a rift valley at the intersection of the San Andreas Fault with the coastline, is approximately 20 km long, 2 km wide, and 6 m deep with mountainous terrain to the southwest and rolling hills to the northeast. Tomales Bay is one of the cleanest estuaries on the West Coast. In winter, approximately 17,000 to 20,000 shorebirds inhabit Tomales Bay and Bodega Bay, which lies directly to the north. At the head of Tomales Bay, the Giacomini Ranch comprises 563 acres of pastureland currently being used for grazing dairy cattle. After more than 50 years of operation as a dairy, the National Park Service acquired the Giacomini property with the intention to restore most of it and the nearby Olema Marsh to tidal wetland. Restoration will add approximately 4% to the existing coastal wetlands in California. The project will return the headwaters of Tomales Bay and two major stream intersections to an intertidal marsh environment, enhancing habitat for both wildlife and fish populations and contributing to the long-term health of Tomales Bay. Prior to the establishment of the ranch, the area was primarily salt marsh that formed as the delta of Lagunitas Creek expanded into Tomales Bay. In converting the salt marsh to dairy land, levees and tide gates were constructed to prevent tidal incursion and stream flooding. Those levees have significantly altered the patterns of estuarine circulation and sediment deposition. To restore natural hydrologic processes within the area and to promote the return of ecological functions and processes, the levees will have to be breached or removed. Developing a successful restoration strategy requires knowledge of elevations within the pastureland and the range of water depths that can be expected from tidal, river, and wind action. In support of the restoration program, the USGS provides technical assistance to PRNS in the form of a scientific study focusing on understanding the physical processes that could affect the Giacomini wetland restoration. The study will yield scientific products that NPS resource managers can use in designing and implementing the restoration project. Research elements include: - Develop a Geodetic Control Network (GCN) throughout PRNS that meets the standards specified National Geodetic Survey data base (the NGS "Bluebook"). The grid will allow this and future studies to be conducted to a precision commensurate with the expressed goals of PRNS. The survey will consist of three steps: (1) verify existing GPS control monuments in the area; (2) tie control monuments in the study areas to the GPS control monuments; and (3) establish NAVD88 elevations using a digital electronic level. - Conduct a detailed survey of the Giacomini site to produce an accurate topographic map of the property. The site survey can be coupled with on-site water-level measurements to produce an empirical flooding model. - Measure water level and wind regime at the Giacomini site. The water-level range is critical to determining the wetland types based on the elevation of the dairy land. Water level at Sacramento Landing, in central Tomales Bay, will also be measured for comparison. As of November 2005, we have created a GCN, produced a detailed topographic map of the Giacomini site, and collected approximately three years of water-level and wind data at the Giacomini site and over one year of usable water-level data at the Sacramento Landing pier.

California↗
USGS2023

Interaction of a legacy groundwater contaminant plume with the Little Wind River from 2015 through 2017, Riverton Processing site, Wyoming

The Riverton Processing site was a uranium mill 4 kilometers southwest of Riverton, Wyoming, that prepared uranium ore for nuclear reactors and weapons from 1958 to 1963. The U.S. Department of Energy completed surface remediation of the uranium tailings in 1989; however, groundwater below and downgradient from the tailings site and nearby Little Wind River was not remediated. Beginning in 2010, a series of floods along the Little Wind River began to mobilize contaminants in the unsaturated zone, resulting in substantial increases of uranium and other contaminants of concern in monitoring wells completed inside the contaminant plume. In 2011, the U.S. Department of Energy started a series of university and Government agency retrospective and field investigations to understand the processes controlling contaminant increases in the groundwater plume. The goals of the field investigations described in this report were to (1) identify and quantify the contaminant flux and potential associated biological effects from groundwater associated with the legacy plume as it enters a perennial stream reach, and (2) assess chemical exposure and potential effects to biological receptors from the interaction of the contaminant plume and the river. Field investigations along the Little Wind River were completed by the U.S. Geological Survey during 2015–17 in cooperation with the U.S. Department of Energy Office of Legacy Management to characterize: (1) seepage areas and seepage rates; (2) pore-water and bed sediment chemistry and hyporheic exchange and reactive loss; and (3) exposure pathways and biological receptors. All data collected during the study are contained in two U.S. Geological Survey data releases, available at https://doi.org/10.5066/F7BR8QX4 and https://doi.org/10.5066/P9J9VJBR . A variety of tools and methods were used during the field characterizations. Streambed temperature mapping, electrical resistivity tomography, electromagnetic induction, fiber-optic distributed temperature sensing, tube seepage meters, vertical thermal sensor arrays, and an environmental tracer (radon) were used to identify areas of groundwater seepage and associated seepage rates along specific sections of the study reach of the river. Drive points, minipiezometers, diffusive equilibrium in thin-film/diffusive gradients in thin-film probes, bed-sediment samples, and equal discharge increment sampling methods were used to characterize pore-water chemistry, estimate hyporheic exchange and reactive loss of selected chemical constituents, and quantify contaminant loadings entering the study reach. Sampling and analysis of surface sediments, filamentous algae, periphytic algae, and macroinvertebrates were used to characterize biological exposure pathways, metal uptake, and receptors. Areas of focused groundwater discharge identified by the fiber-optic distributed temperature sensing surveys corresponded closely with areas of elevated electrical conductivity identified by the electromagnetic induction survey results in the top 5 meters of sediment. During three monitoring periods in 2016, the mean vertical seepage rate measured with tube seepage meters was 0.45 meter per day, ranging from −0.02 to 1.55 meters per day. Five of the 11 locations where vertical thermal profile data were collected along the study reach during August 2017 indicated mean upwelling values ranging from 0.11 to 0.23 meter per day. Radon data collected from the Little Wind River during June, July, and August 2016 indicated a consistent inflow of groundwater to the central part of the study reach, in the area congruous with the center of the previously mapped groundwater plume discharge zone. During August 2017, the greatest attenuation of uranium from reactive loss in pore-water samples was observed at three locations along the study reach, at depths between 6 and 15 centimeters, and similar trends in molybdenum attenuation were also observed. Bed-sediment concentration profiles collected during 2017 also indicated attenuation of uranium and molybdenum from groundwater during hyporheic mixing of surface water with the legacy plume during groundwater upwelling into the river. Streamflow measurements combined with equal discharge increment water sampling along the study reach indicated an increase in dissolved uranium concentrations in the downstream direction during 2016 and 2017. Net uranium load entering the Little Wind River study reach was about 290 and 435 grams per day during 2016 and 2017, respectively. Biological samples indicated that low levels of uranium and molybdenum exposure were confined to the benthos in the Little Wind River within and immediately downstream from the perimeter of the groundwater plume. Concentrations of molybdenum and uranium in filamentous algae were consistently low at all sites in the study reach with no indication of increased exposure of dissolved bioavailable molybdenum or uranium at sites next to or downstream from the groundwater plume. Comparison of the August 2017 results from electromagnetic induction, tube seepage meters, vertical thermal profiling, and pore-water chemistry surveys were in general agreement in identifying areas with upwelling groundwater conditions along the study reach. However, the electroconductivity values measured with electromagnetic induction in the top 100 centimeters of sediment did not agree with sodium concentrations measured in pore-water samples collected at similar streambed depths. Differences and similarities between multiple methods can result in additional insights into hydrologic and biogeochemical processes that may be occurring along a reach of a river system interacting with shallow groundwater inputs. It may be advantageous to apply a variety of geophysical, geochemical, hydrologic, and biological tools at other Uranium Mill Tailings Remedial Action ( https://www.energy.gov/sites/prod/files/2014/10/f19/UMTRCA.pdf ) sites during the investigation of legacy contaminant plume interactions with surface-water systems.

Wyoming↗
USGS2018

Evaluating the potential for near-shore bathymetry on the Majuro Atoll, Republic of the Marshall Islands, using Landsat 8 and WorldView-3 imagery

Satellite-derived near-shore bathymetry (SDB) is becoming an increasingly important method for assessing vulnerability to climate change and natural hazards in low-lying atolls of the northern tropical Pacific Ocean. Satellite imagery has become a cost-effective means for mapping near-shore bathymetry because ships cannot collect soundings safely while operating close to the shore. Also, green laser light detection and ranging (lidar) acquisitions are expensive in remote locations. Previous research has demonstrated that spectral band ratio-based techniques, commonly called the natural logarithm approach, may lead to more precise measurements and modeling of bathymetry because of the phenomenon that different substrates at the same depth have approximately equal ratio values. The goal of this research was to apply the band ratio technique to Landsat 8 at-sensor radiance imagery and WorldView-3 atmospherically corrected imagery in the coastal waters surrounding the Majuro Atoll, Republic of the Marshall Islands, to derive near-shore bathymetry that could be incorporated into a seamless topobathymetric digital elevation model of Majuro. Attenuation of light within the water column was characterized by measuring at-sensor radiance and reflectance at different depths and calculating an attenuation coefficient. Bathymetric lidar data, collected by the U.S. Naval Oceanographic Office in 2006, were used to calibrate the SDB results. The bathymetric lidar yielded a strong linear relation with water depths. The Landsat 8-derived SDB estimates derived from the blue/green band ratio exhibited a water attenuation extinction depth of 6 meters with a coefficient of determination R 2 =0.9324. Estimates derived from the coastal/red band ratio had an R 2 =0.9597. At the same extinction depth, SDB estimates derived from WorldView-3 imagery exhibited an R 2 =0.9574. Because highly dynamic coastal shorelines can be affected by erosion, wetland loss, hurricanes, sea-level rise, urban development, and population growth, consistent bathymetric data are needed to better understand sensitive coastal land/water interfaces in areas subject to coastal disasters.

Majuro Atoll↗
USGS2010

StreamStats in Oklahoma– Drainage-basin characteristics and peak-flow frequency statistics for ungaged streams

The USGS Streamflow Statistics (StreamStats) Program was created to make geographic information systems-based estimation of streamflow statistics easier, faster, and more consistent than previously used manual techniques. The StreamStats user interface is a map-based internet application that allows users to easily obtain streamflow statistics, basin characteristics, and other information for user-selected U.S. Geological Survey data-collection stations and ungaged sites of interest. The application relies on the data collected at U.S. Geological Survey streamflow-gaging stations, computer aided computations of drainage-basin characteristics, and published regression equations for several geographic regions comprising the United States. The StreamStats application interface allows the user to (1) obtain information on features in selected map layers, (2) delineate drainage basins for ungaged sites, (3) download drainage-basin polygons to a shapefile, (4) compute selected basin characteristics for delineated drainage basins, (5) estimate selected streamflow statistics for ungaged points on a stream, (6) print map views, (7) retrieve information for U.S. Geological Survey streamflow-gaging stations, and (8) get help on using StreamStats. StreamStats was designed for national application, with each state, territory, or group of states responsible for creating unique geospatial datasets and regression equations to compute selected streamflow statistics. With the cooperation of the Oklahoma Department of Transportation, StreamStats has been implemented for Oklahoma and is available at http://water.usgs.gov/osw/streamstats/. The Oklahoma StreamStats application covers 69 processed hydrologic units and most of the state of Oklahoma. Basin characteristics available for computation include contributing drainage area, contributing drainage area that is unregulated by Natural Resources Conservation Service floodwater retarding structures, mean-annual precipitation at the drainage-basin outlet for the period 1961-1990, 10-85 channel slope (slope between points located at 10 percent and 85 percent of the longest flow-path length upstream from the outlet), and percent impervious area. The Oklahoma StreamStats application interacts with the National Streamflow Statistics database, which contains the peak-flow regression equations in a previously published report. Fourteen peak-flow (flood) frequency statistics are available for computation in the Oklahoma StreamStats application. These statistics include the peak flow at 2-, 5-, 10-, 25-, 50-, 100-, and 500-year recurrence intervals for rural, unregulated streams; and the peak flow at 2-, 5-, 10-, 25-, 50-, 100-, and 500-year recurrence intervals for rural streams that are regulated by Natural Resources Conservation Service floodwater retarding structures. Basin characteristics and streamflow statistics cannot be computed for locations in playa basins (mostly in the Oklahoma Panhandle) and along main stems of the largest river systems in the state, namely the Arkansas, Canadian, Cimarron, Neosho, Red, and Verdigris Rivers, because parts of the drainage areas extend outside of the processed hydrologic units.

Oklahoma↗
USGS2020

Hydrogeologic and geochemical characterization of groundwater resources in Pine and Wah Wah Valleys, Iron, Beaver, and Millard Counties, Utah

Pine and Wah Wah Valleys are neighboring structural basins that encompass about 1,330 square miles in Beaver, Iron, and Millard Counties in Utah, approximately 50 miles northwest of Cedar City, Utah, and 50 miles southeast of Baker, Nevada. Perennial streamflow is limited and only exists in higher-altitude reaches of small mountain streams in both basins. Groundwater is in unconsolidated basin-fill aquifers and bedrock mountain aquifers. Groundwater in Pine and Wah Wah Valleys is being targeted for large-scale groundwater extraction and export to provide municipal supply to the growing population in Iron County, Utah. Concern about declining groundwater levels and spring flows from proposed groundwater withdrawals has increased interest in an improved understanding of the groundwater system. Previous studies have indicated that an average of 28,000 acre-feet per year of recharge occurs mostly as infiltration of precipitation in high-altitude regions in the two basins. Groundwater discharge in the mountain hydrologic systems was estimated to average 8,500 acre-feet per year and is assumed to be consumed before subsequently recharging the valley basin-fill aquifers. Subsurface groundwater outflow moves from basin-fill aquifers in Pine and Wah Wah Valleys northward to adjacent regional basins and was estimated to average 19,500 acre-feet per year. An updated water-level map for the basin-fill aquifers in Pine and Wah Wah Valleys indicates that groundwater moves northward along the lengths of both valleys toward adjacent basins. Measured depths to water range from about 210 to 750 feet below land surface in Wah Wah Valley, and from about 300 to 620 feet below land surface in Pine Valley. Long-term water levels at seven wells completed in the basin-fill aquifers of Pine and Wah Wah Valleys with records spanning more than 40 years are generally stable with observed fluctuations of less than 5 feet. Observed discharge from two springs monitored between 2013 and 2016 also is generally stable. Groundwater leaving Pine and Wah Wah Valleys through the subsurface moves northward, converges with regional groundwater flow, and discharges by evapotranspiration at regional groundwater discharge areas, likely Tule Valley, Utah. In this study, basin-scale groundwater discharge was estimated by (1) mapping the groundwater discharge areas in each valley; (2) evaluating the 2005–11 summer multispectral satellite images against the Basin and Range carbonate-rock aquifer system study evapotranspiration measurements to select scenes broadly representative of average conditions in the study area and partitioning the groundwater discharge areas into evapotranspiration units using the selected satellite images and field reconnaissance; and (3) scaling evapotranspiration to the evapotranspiration units using evapotranspiration-rate estimates from several studies in the Great Basin. The resulting updated estimates of average annual groundwater evapotranspiration in the Tule Valley and Sevier Lake groundwater discharge areas were 35,000 and 10,500 acre-feet per year, respectively, with a likely uncertainty of plus or minus 35 percent. Groundwater samples from 13 sites in Pine Valley and 11 sites in Wah Wah Valley were analyzed for major ions and nutrients, to characterize geochemistry and water quality. Groundwater samples also were analyzed for the stable isotopes of oxygen, hydrogen, and carbon, the radioactive isotopes of carbon and hydrogen, and dissolved noble gases including helium-3, helium-4, neon, argon, krypton and xenon. Groundwater sampling sites included 12 wells and 12 springs. Carbon-14 and tritium/helium groundwater age dating indicate that groundwater in the basin-fill aquifers is typically thousands to tens of thousands of years older than groundwater in the shallow mountain aquifers. Dissolved-solids concentrations are lower and noble-gas temperatures are warmer in the valley wells compared to almost all groundwater sampled from wells and springs in the surrounding mountains. These results indicate a hydraulic discontinuity between the mountain and valley aquifers throughout much of the study area, and that much of the valley recharge is not derived from direct infiltration of precipitation in the mountains.

Utah↗
USGS2007

Hydrology of Polk County, Florida

Local water managers usually rely on information produced at the State and regional scale to make water-resource management decisions. Current assessments of hydrologic and water-quality conditions in Polk County, Florida, commonly end at the boundaries of two water management districts (South Florida Water Management District and the Southwest Florida Water Management District), which makes it difficult for managers to determine conditions throughout the county. The last comprehensive water-resources assessment of Polk County was published almost 40 years ago. To address the need for current countywide information, the U.S. Geological Survey began a 3?-year study in 2002 to update information about hydrologic and water-quality conditions in Polk County and identify changes that have occurred. Ground-water use in Polk County has decreased substantially since 1965. In 1965, total ground-water withdrawals in the county were about 350 million gallons per day. In 2002, withdrawals totaled about 285 million gallons per day, of which nearly 95 percent was from the Floridan aquifer system. Water-conservation practices mainly related to the phosphate-mining industry as well as the decrease in the number of mines in operation in Polk County have reduced total water use by about 65 million gallons per day since 1965. Polk County is underlain by three principal hydrogeologic units. The uppermost water-bearing unit is the surficial aquifer system, which is unconfined and composed primarily of clastic deposits. The surficial aquifer system is underlain by the intermediate confining unit, which grades into the intermediate aquifer system and consists of up to two water-bearing zones composed of interbedded clastic and carbonate rocks. The lowermost hydrogeologic unit is the Floridan aquifer system. The Floridan aquifer system, a thick sequence of permeable limestone and dolostone, consists of the Upper Floridan aquifer, a middle semiconfining unit, a middle confining unit, and the Lower Floridan aquifer. The Upper Floridan aquifer provides most of the water required to meet demand in Polk County. Data from about 300 geophysical and geologic logs were used to construct hydrogeologic maps showing the tops and thicknesses of the aquifers and confining units within Polk County. Thickness of the surficial aquifer system ranges from several feet thick or less in the extreme northwestern part of the county and along parts of the Peace River south of Bartow to more than 200 feet along the southern part of the Lake Wales Ridge in eastern Polk County. Thickness of the intermediate aquifer system/intermediate confining unit is highly variable throughout the county because of past erosional processes and sinkhole formation. Thickness of the unit ranges from less than 25 feet in the extreme northwestern part of the county to more than 300 feet in southwestern Polk County. The altitude of the top of the Upper Floridan aquifer in the county ranges from about 50 feet above National Geodetic Vertical Datum of 1929 (NGVD 29) in the northwestern part to more than 250 feet below NGVD 29 in the southern part. Water levels in the Upper Floridan aquifer fluctuate seasonally, increasing during the wet season (June through September) and decreasing during the rest of the year. Water levels in the Upper Floridan aquifer also can change from year to year, depending on such factors as pumpage and climatic variations. In the southwestern part of the county, fluctuations in water use related to phosphate mining have had a major impact on ground-water levels. Hydrographs of selected wells in southwestern Polk County show a general decline in water levels that ended in the mid-1970s. This water-level decline coincides with an increase in water use associated with phosphate mining. A substantial increase in water levels that began in the mid-1970s coincides with a period of decreasing water use in the county. Despite reductions in water use since 1970, howev

Scientific Investigations Report↗
USGS2012

Hydrogeology of the stratified-drift aquifers in the Cayuta Creek and Catatonk Creek valleys in parts of Tompkins, Schuyler, Chemung, and Tioga Counties, New York

The surficial deposits, areal extent of aquifers, and the water-table configurations of the stratified-drift aquifer systems in the Cayuta Creek and Catatonk Creek valleys and their large tributary valleys in Tompkins, Schuyler, Chemung, and Tioga Counties, New York were mapped in 2009, in cooperation with the New York State Department of Environmental Conservation. Well and test-boring records, surficial deposit maps, Light Detection and Ranging (LIDAR) data, soils maps, and horizontal-to-vertical ambient-noise seismic surveys were used to map the extent of the aquifers, construct geologic sections, and determine the depth to bedrock (thickness of valley-fill deposits) at selected locations. Geologic materials in the study area include sedimentary bedrock, unstratified drift (till), stratified drift (glaciolacustrine and glaciofluvial deposits), and recent alluvium. Stratified drift consisting of glaciofluvial sand and gravel is the major component of the valley fill in this study area. The deposits are present in sufficient amounts in most places to form extensive unconfined aquifers throughout the study area and, in some places, confined aquifers. Stratified drift consisting of glaciolacustrine fine sand, silt, and clay are present locally in valleys underlying the surficial sand and gravel deposits in the southern part of the Catatonk Creek valley. These unconfined and confined aquifers are the source of water for most residents, farms, and businesses in the valleys. A generalized depiction of the water table in the unconfined aquifer was constructed using water-level measurements made from the 1950s through 2010, as well as LIDAR data that were used to determine the altitudes of perennial streams at 10-foot contour intervals and water surfaces of ponds and wetlands that are hydraulically connected to the unconfined aquifer. The configuration of the water-table contours indicate that the general direction of groundwater flow within Cayuta Creek and Catatonk Creek stratified-drift aquifers is predominantly from the valley walls toward the main streams in the valleys. The groundwater discharges from the aquifer system to the main-stem streams in the valleys. Locally, the direction of groundwater flow is radially away from groundwater mounds that have formed beneath upland tributaries that typically lose water where they flow on alluvial fans in the valleys. In some places, groundwater that would normally flow toward streams is intercepted by pumping wells.

New York↗
USGS2022

Assessment of habitat availability for juvenile Chinook salmon (Oncorhynchus tshawytscha) and steelhead (O. mykiss) in the Willamette River, Oregon

The Willamette River, Oregon, is home to two salmonid species listed as threatened under the Endangered Species Act, Upper WIllamette River spring Chinook salmon ( Oncorhynchus tshawytscha ) and Upper Willamette River winter steelhead ( O. mykiss ). Streamflow in the Willamette River is regulated by upstream dams, 13 of which are operated by the U.S. Army Corps of Engineers (USACE) as part of the Willamette Valley Project. In 2008, these dams were determined to have a deleterious effect on Endangered Species Act-listed salmonids, resulting in USACE taking actions to mitigate those effects. Mitigation actions included setting seasonal streamflow targets at various locations along the river to improve survival and migration of juvenile salmonids. Although these targets were established with the best available information at the time, recent data and models have advanced understanding of Willamette River bathymetric, hydraulic, and thermal conditions, allowing for a more robust analysis of the effect of streamflow on downstream habitat. This study integrates those recent advances to build high-resolution models of usable habitat for juvenile Chinook salmon and steelhead to assess variation in spatial and seasonal patterns of habitat availability. Specifically, this study develops detailed maps of habitat availability for juvenile Chinook salmon and steelhead for two size classes (fry and pre-smolt). Habitat availability is modeled in a three-step process whereby (1) two-dimensional hydraulic models are paired with literature-supplied data on habitat preferences to create spatially explicit maps of rearing habitats for a wide range of streamflows; (2) reach-specific relations between streamflow and habitat area are developed and paired with streamgage records to create habitat time series for 2011, 2015, and 2016, which reflect “cool and wet,” “hot and dry,” and “warm but average precipitation” conditions, respectively; (3) temperature models are coupled with literature-based thermal thresholds to determine time periods and locations along the river corridor when rearing habitat has optimal, harmful, or lethal temperature conditions; (4) finally, habitat availability is summarized at several spatial scales to characterize longitudinal and seasonal patterns. Findings show that modeled area of rearing habitat for Chinook salmon and steelhead responds non-uniformly to streamflow, where habitat in some reaches of the Willamette River consistently increase with additional streamflow, while in other reaches, habitat area decreases when streamflows increase from low to moderate flows. Modeled differences in flow-habitat relations are primarily explained by local geomorphology in each reach and resulting hydraulic conditions that arise with different streamflows. These are most pronounced when comparing laterally active, multi-channel reaches upstream from Corvallis with downstream reaches that are laterally stable with single-channel planforms. The reaches upstream from Corvallis generally have more habitat available per unit stream distance than downstream reaches, but all reaches display greatest amounts of habitat at the highest streamflows. Finally, results show that warm water temperature in summer greatly decreases the utility of habitat available to the focal species, particularly downstream from Corvallis. Together, these findings serve to inform flow management by characterizing spatial and seasonal patterns of habitat availability for juvenile spring Chinook salmon and winter steelhead and provide a quantitative assessment of the effects of streamflow on rearing habitat.

Oregon↗
USGS2011

Geomorphic Classification and Evaluation of Channel Width and Emergent Sandbar Habitat Relations on the Lower Platte River, Nebraska

This report presents a summary of geomorphic characteristics extracted from aerial imagery for three broad segments of the Lower Platte River. This report includes a summary of the longitudinal multivariate classification in Elliott and others (2009) and presents a new analysis of total channel width and habitat variables. Three segments on the lower 102.8 miles of the Lower Platte River are addressed in this report: the Loup River to the Elkhorn River (70 miles long), the Elkhorn River to Salt Creek (6.9 miles long), and Salt Creek to the Missouri River (25.9 miles long). The locations of these segments were determined by the locations of tributaries potentially significant to the hydrology or sediment supply of the Lower Platte River. This report summarizes channel characteristics as mapped from July 2006 aerial imagery including river width, valley width, channel curvature, and in-channel habitat features. In-channel habitat measurements were not made under consistent hydrologic conditions and must be considered general estimates of channel condition in late July 2006. Longitudinal patterns in these features are explored and are summarized in the context of the longitudinal multivariate classification in Elliott and others (2009) for the three Lower Platte River segments. Detailed descriptions of data collection and classification methods are described in Elliott and others (2009). Nesting data for the endangered interior least tern (Sternula antillarum) and threatened piping plover (Charadrius melodus) from 2006 through 2009 are examined within the context of the multivariate classification and Lower Platte River segments. The widest reaches of the Lower Platte River are located in the segment downstream from the Loup River to the Elkhorn River. This segment also has the widest valley and highest degree of braiding of the three segments and many large vegetated islands. The short segment of river between the Elkhorn River and Salt Creek has a fairly low valley width and high channel sinuosities at larger scales. The segment from Salt Creek to the Missouri River has narrow valleys and generally low channel sinuosity. Tern and plover nest sites from 2006 through 2009 in the multi-scale multivariate classification indicated relative nesting selection of cluster 2 reaches among the four-cluster classification and reaches containing clusters 2, 3, and 6 from the seven-cluster classification. These classes, with the exception of cluster 6 are common downstream from the Elkhorn River. Trends in total channel width indicated that reaches dominated by dark vegetation (islands) are the widest on the Lower Platte River. Reaches with high percentages of dry sand and dry sand plus light vegetation were the narrowest reaches. This suggests that narrow channel reaches have sufficient transport capacity to maintain sandbars under recent (2006) flow regimes and are likely to be most amenable to maintaining tern and plover habitat in the Lower Platte River. Further investigations into the dynamics of emergent sandbar habitat and the effects of bank stabilization on in-channel habitats will require the collection and analysis of new data, particularly detailed elevation information and an assessment of existing bank stabilization structures.

Scientific Investigations Report↗
USGS2011

Hydrogeologic setting and simulation of groundwater flow near the Canterbury and Leadville Mine Drainage Tunnels, Leadville, Colorado

The Leadville mining district is historically one of the most heavily mined regions in the world producing large quantities of gold, silver, lead, zinc, copper, and manganese since the 1860s. A multidisciplinary investigation was conducted by the U.S. Geological Survey, in cooperation with the Colorado Department of Public Health and Environment, to characterize large-scale groundwater flow in a 13 square-kilometer region encompassing the Canterbury Tunnel and the Leadville Mine Drainage Tunnel near Leadville, Colorado. The primary objective of the investigation was to evaluate whether a substantial hydraulic connection is present between the Canterbury Tunnel and Leadville Mine Drainage Tunnel for current (2008) hydrologic conditions. Altitude in the Leadville area ranges from about 3,018 m (9,900 ft) along the Arkansas River valley to about 4,270 m (14,000 ft) along the Continental Divide east of Leadville, and the high altitude of the area results in a moderate subpolar climate. Winter precipitation as snow was about three times greater than summer precipitation as rain, and in general, both winter and summer precipitation were greatest at higher altitudes. Winter and summer precipitation have increased since 2002 coinciding with the observed water-level rise near the Leadville Mine Drainage Tunnel that began in 2003. The weather patterns and hydrology exhibit strong seasonality with an annual cycle of cold winters with large snowfall, followed by spring snowmelt, runoff, and recharge (high-flow) conditions, and then base-flow (low-flow) conditions in the fall prior to the next winter. Groundwater occurs in the Paleozoic and Precambrian fractured-rock aquifers and in a Quaternary alluvial aquifer along the East Fork Arkansas River, and groundwater levels also exhibit seasonal, although delayed, patterns in response to the annual hydrologic cycle. A three-dimensional digital representation of the extensively faulted bedrock was developed and a geophysical direct-current resistivity field survey was performed to evaluate the geologic structure of the study area. The results show that the Canterbury Tunnel is located in a downthrown structural block that is not in direct physical connection with the Leadville Mine Drainage Tunnel. The presence of this structural discontinuity implies there is no direct groundwater pathway between the tunnels along a laterally continuous bedrock unit. Water-quality results for pH and major-ion concentrations near the Canterbury Tunnel showed that acid mine drainage has not affected groundwater quality. Stable-isotope ratios of hydrogen and oxygen in water indicate that snowmelt is the primary source of groundwater recharge. On the basis of chlorofluorocarbon and tritium concentrations and mixing ratios for groundwater samples, young groundwater (groundwater recharged after 1953) was indicated at well locations upgradient from and in a fault block separate from the Canterbury Tunnel. Samples from sites downgradient from the Canterbury Tunnel were mixtures of young and old (pre-1953) groundwater and likely represent snowmelt recharge mixed with older regional groundwater that discharges from the bedrock units to the Arkansas River valley. Discharge from the Canterbury Tunnel contained the greatest percentage of old (pre-1953) groundwater with a mixture of about 25 percent young water and about 75 percent old water. A calibrated three-dimensional groundwater model representing high-flow conditions was used to evaluate large-scale flow characteristics of the groundwater and to assess whether a substantial hydraulic connection was present between the Canterbury Tunnel and Leadville Mine Drainage Tunnel. As simulated, the faults restrict local flow in many areas, but the fracture-damage zones adjacent to the faults allow groundwater to move along faults. Water-budget results indicate that groundwater flow across the lateral edges of the model controlled the majority of flow in and out of the aquifer (79 percent and 63 percent of the total water budget, respectively). The largest contributions to the water budget were groundwater entering from the upper reaches of the watershed and the hydrologic interaction of the groundwater with the East Fork Arkansas River. Potentiometric surface maps of the simulated model results were generated for depths of 50, 100, and 250 m. The surfaces revealed a positive trend in hydraulic head with land-surface altitude and evidence of increased control on fluid movement by the fault network structure at progressively greater depths in the aquifer. Results of advective particle-tracking simulations indicate that the sets of simulated flow paths for the Canterbury Tunnel and the Leadville Mine Drainage Tunnel were mutually exclusive of one another, which also suggested that no major hydraulic connection was present between the tunnels. Particle-tracking simulations also revealed that although the fault network generally restricted groundwater movement locally, hydrologic conditions were such that groundwater did cross the fault network at many locations. This cross-fault movement indicates that the fault network controls regional groundwater flow to some degree but is not a complete barrier to flow. The cumulative distributions of adjusted age results for the watershed indicate that approximately 30 percent of the flow pathways transmit groundwater that was younger than 68 years old (post-1941) and that about 70 percent of the flow pathways transmit old groundwater. The particle-tracking results are consistent with the apparent ages and mixing ratios developed from the chlorofluorocarbon and tritium results. The model simulations also indicate that approximately 50 percent of the groundwater flowing through the study area was less than 200 years old and about 50 percent of the groundwater flowing through the study area is old water stored in low-permeability geologic units and fault blocks. As a final examination of model response, the conductance parameters of the Canterbury Tunnel and Leadville Mine Drainage Tunnel were manually adjusted from the calibrated values to determine if altering the flow discharge in one tunnel affects the hydraulic behavior in the other tunnel. The examination showed no substantial hydraulic connection. The multidisciplinary investigation yielded an improved understanding of groundwater characteristics near the Canterbury Tunnel and the Leadville Mine Drainage Tunnel. Movement of groundwater between the Canterbury Tunnel and Leadville Mine Drainage Tunnel that was central to this investigation could not be evaluated with strong certainty owing to the structural complexity of the region, study simplifications, and the absence of observation data within the upper sections of the Canterbury Tunnel and between the Canterbury Tunnel and the Leadville Mine Drainage Tunnel. There was, however, collaborative agreement between all of the analyses performed during this investigation that a substantial hydraulic connection did not exist between the Canterbury Tunnel and the Leadville Mine Drainage Tunnel under natural flow conditions near the time of this investigation.

Colorado↗