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,369 records · Page 76Linked to original sources

Mapping the altitude of the top of the Dockum Group and paleochannel analysis using surface geophysical methods on and near Cannon Air Force Base in Curry County, New Mexico, 2020

The hydrogeology on and near Cannon Air Force Base (AFB) in eastern New Mexico was assessed to gain a better understanding of preferential groundwater flow paths through paleochannels. In and near the study area, paleochannels incised the top surface of the Dockum Group (Chinle Formation) and were subsequently filled in with electrically resistive coarse-grained sediments of the overlying Ogallala Formation, resulting in a preferential groundwater flow path in the form of a paleochannel network. A better understanding of the spatial characteristics of this preferential groundwater flow path is needed to support ongoing efforts to remediate groundwater contamination at Cannon AFB. Therefore, the U.S. Geological Survey, in cooperation with the U.S. Air Force Civil Engineer Center, used surface geophysical resistivity methods and data compiled from previous studies to better understand the spatial distribution and characteristics of the paleochannel network incised into the top of the Dockum Group. Previous studies have shown these paleochannels incised into the top of the Dockum Group with increasing resolution, but limited borehole data on and near Cannon AFB continued to make accurately mapping the top of Dockum Group challenging. For this study, surface geophysical resistivity measurements in the form of time-domain electromagnetic soundings made by the U.S. Geological Survey were used in conjunction with data previously published by Architecture, Engineering, Construction, Operations, and Management and borehole data compiled from the New Mexico Water Rights Reporting System database to prepare an updated map of the top of the Dockum Group that includes the location and characteristics of paleochannels incised into the top of the Dockum Group (Chinle Formation). A total of 149 borehole picks (determinations of the tops and bases of geologic units and their hydrogeologic-unit equivalents) were obtained from previous studies, along with 72 additional borehole picks from the New Mexico Water Rights Reporting System database and 43 picks from newly collected time-domain electromagnetic soundings. The data were gridded and contoured using Oasis Montaj v. 9.8.1. The updated map of the top of Dockum Group has many areas of uncertainty greater than 20 feet, because there are not enough data for the gridding process to reliably determine a value. However, this interpretation of the altitude of the top of the Dockum Group represents a substantial improvement in data resolution compared to previous studies. Two methodologies were used to evaluate paleochannels incised in the top of the Dockum Group across the study area: (1) trend-removal grid analysis and (2) analysis with Esri’s ArcMap Hydrology toolset. These two paleochannel analysis techniques show groundwater flow direction as well as areas having the deepest saturated thickness. Hydrologically, these techniques show where aquifer storage is highest (in the deepest parts of the paleochannel network), as well as the spatial distribution of preferential groundwater flow paths (the paleochannels). The analyses indicate a large paleochannel trending to the southeast, with smaller channels feeding in from the west. Areas where groundwater management could be more beneficial are indicated by locations where these flow lines intersect the deeper parts of the paleochannel.

New Mexico↗

Porphyry copper assessment of western Central Asia

The U.S. Geological Survey conducted an assessment of resources associated with porphyry copper deposits in the western Central Asia countries of Kyrgyzstan, Uzbekistan, Kazakhstan, and Tajikistan and the southern Urals of Kazakhstan and Russia as part of a global mineral resource assessment. The purpose of the study was to (1) delineate permissive areas (tracts) for undiscovered porphyry copper deposits; (2) compile a database of known porphyry copper deposits and significant prospects; (3) where data permit, estimate numbers of undiscovered deposits within those permissive tracts; and (4) provide probabilistic estimates the amounts of copper (Cu), molybdenum (Mo), gold (Au), and silver (Ag) that could be contained in those undiscovered deposits. Western Central Asia, a region diverse in its geologic complexity, is situated north of the Tarim and North China tectonic blocks and sandwiched between the East European and Siberian cratons. The Ural Mountains form the western margin of the region; the southern margin is formed by the high-standing ranges that make up the Tian Shan mountain range in the border regions of Kazakhstan, Kyrgyzstan, and western China, where the effects of collisional tectonics are well displayed. The tectonic collage that makes up the core of western Central Asia is perhaps the geologically least understood part of the region. There is broad agreement that the early Paleozoic is made up of tectonically juxtaposed blocks that vary from Precambrian-cored microcontinents to magmatic arc and related complexes to subduction-related accretionary complexes. The rudiments of an incipient, contiguous single Kazakhstan block were formed by the end of the Silurian. In the middle to late Paleozoic, the block was unconformably superposed by two large, nested magmatic-arc belts, one Devonian, the other Carboniferous. Both magmatic-arc complexes were folded into a horseshoe-shaped, southeast-opening orocline in response to the final collisions of the various surrounding cratonic blocks with the Kazakhstan block. Additional deformation in the upper Cenozoic derived from the collision of India and China significantly redistributed fragments of the various mosaicked blocks, particularly in the central and southern parts of the western Central Asian region. Porphyry copper deposits are associated with many of the magmatic-arc fragments and belts throughout the geologically complex region, and economically important deposits are found in arc sequences of all Paleozoic Periods. The economically most productive arcs are Carboniferous. The assessment includes a discussion of the tectonic and geologic setting of porphyry copper deposits in western Central Asia (chapter 1), an application of remote sensing data for hydrothermal alteration mapping as a tool for porphyry copper assessment in the region (chapter 2), and a probabilistic assessment of undiscovered porphyry copper resources in four areas that represent Ordovician and Late Paleozoic (Carboniferous-Permian) magmatic arcs (chapter 3). The principal litho-tectonic terrane concept used to delineate permissive tracts was that of a magmatic arc that formed in the subduction boundary zone above a subducting plate. Eight permissive tracts are delineated on the basis of mapped and inferred subsurface distributions of igneous rocks assigned to magmatic arcs of specified age ranges that define areas where the occurrence of porphyry copper deposits within 1 kilometer of the Earth’s surface is possible. These tracts range in area from about 8,000 to 200,000 square kilometers and host 18 known porphyry copper deposits that contain about 54 million metric tons of copper. Available data included geologic maps, the distribution of significant porphyry copper occurrences and potentially related deposit types, the distribution of hydrothermal alteration patterns that are consistent with porphyry copper mineralization, and information on possible subsurface extensions of permissive rocks. On the basis of analyses of these data, the assessment team estimated a mean of 25 undiscovered porphyry copper deposits for the study area. Estimates of numbers of undiscovered deposits were combined with grade and tonnage models in a Monte Carlo simulation to yield a mean estimate of 95 million metric tons of copper in undiscovered porphyry copper deposits; this represents about twice the amount of identified porphyry copper resources (54 million metric tons). Detailed descriptions of each permissive tract, including the rationales for delineation and assessment, are given in appendixes, along with a geographic information system (GIS) that includes permissive tract boundaries, point locations of known porphyry copper deposits and significant occurrences, and hydrothermal alteration data based on analysis of remote sensing data.

Ural Mountains↗

Assessment of Ambystomatid salamander populations and their breeding habitats in the Delaware Water Gap National Recreation Area

This report presents abundance and occurrence data for three species of ambystomad salamanders ( Ambystoma maculatum, A. jeffersonianum, and A. opacum ) collected over a 3-year period (2000, 2001, and 2002) at 200 potentional breeding sies within the Delaware Water Gap National Recreation Area (DEWA). In addition, numerous measures of inpond, near-pond, and landscape attributes were measured and used to inform statistical models to determine species-habitat relationships in the DEWA. The results of a 3-year study of ambystomatid salamander breeding habits and habitats in the (DEWA) that was conducted by the U.S. Geological Survey, in cooperation with the National Park Service, are described in the report. The objectives of the study were to document the population status and critical breeding habitats of the three species of ambystomatid salamanders known to be present in the DEWA— Ambystoma maculatum (spotted salamander), A. opacum (marbled salamander), and A. jeffersonianum (Jefferson salamander). DEWA managers are interested in ecological information on these species for several reasons. First, at the time the study began, there was little known regarding the status of pond-breeding amphibians and their habitats in the DEWA. Second, because they require undegraded habitats in both terrestrial and aquatic habitats to successfully complete their life cycles, the status of ambystomatid salamanders is widely viewed as indicative of overall ecosystem health. Third, because ambystomatid salamanders and other pond-breeding amphibians have been observed in numerous artificial impoundments with the DEWA, park managers would like to assess whether dismantling or discontinuing maintenance of artificial impoundments could affect pond-breeding amphibians and possibly other species that use pond or wetland habitats in the Park. In 2001, 2002, and 2003, the size and location of 200 wetlands, ponds, and artificial impoundments, and related landscape positions (Ridge versus Valley; Pennsylvania side versus New Jersey side of the Delaware river) were mapped, and site habitat data relating to salamander occurrence and abundance patterns were collected. The data collected during this study provide important new baseline information on ambystomatid salamanders and wetland habitats in the DEWA that will enhance long-term inventory and monitoring efforts. In addition, breeding habitat assessments indicate that ambystomatid salamanders may be sensitive to a wide variety of stresses important in the DEWA and in the region. In particular, recent trends in development (for example, roads) in and near the DEWA, regional increases in the acidity of precipitation, and predicted long-term warming trends for the region could be detrimental to pond-breeding salamander populations because of their effects on breeding site quality and quantity, and on the integrity of migration corridors. In contrast, the results of the study indicate management plans to eliminate small impoundments are not likely to adversely affect salamanders in DEWA, at least in the short-term. However, it is possible that these small impoundments may offer stable habitats that provide a rescure effect during long-term droughts.

New Jersey, Pennsylvania↗

Small basin annual yield and percentage of snowmelt runoff in North Dakota, 1931–2016

The North Dakota hydrology manual prepared by the U.S. Department of Agriculture, Soil Conservation Service, presents methodologies primarily used for developing hydrology for onfarm conservation practices, watershed projects, Resource Conservation and Development project measures, and river basin studies. The manual includes data necessary for determining hydrologic factors and developing a design discharge for a given site and intended purpose. The U.S. Geological Survey, in cooperation with the North Dakota Natural Resources Conservation Service, developed methods to reproduce and update the annual yield maps for chapter 7 of the North Dakota hydrology manual. Annual yields, in acre-feet per square mile, for the 50- and 80-percent exceedance probabilities and expected percentage of snowmelt runoff isolines were estimated using U.S. Geological Survey streamflow data from 1931 to 2016 for 71 selected streamgages with drainage areas of 505 square miles or less. An application of a modified Maintenance of Variance Extension Type III was used to estimate missing annual streamflow volumes. An alternate expected percentage of snowmelt runoff isolines was estimated using High Plains Climatic Center precipitation and snowmelt data from 1931 to 2016 for 85 selected sites. The final expected percentage of snowmelt runoff isolines was estimated using streamflow data instead of precipitation and snowfall depth data. A snowmelt runoff seasonal period of March–May produced better isoline slopes than a November–May runoff seasonal period. Slopes of the expected percentage of snowmelt runoff isolines were sensitive to amounts of missing record. Suitable isoline slopes appeared when the missing record was set to 50 percent (43 years) and 66 percent (57 years) for the 86-year period of 1931–2016.

North Dakota↗

Simulated effects of projected 2014–40 withdrawals on groundwater flow and water levels in the New Jersey Coastal Plain

Abstract Groundwater flow between 2014 through 2040 was simulated in the New Jersey Coastal Plain based on three withdrawal scenarios. Two of the scenarios were based on projected population trends and the assumption of water conservation; the nominal water-loss scenario projected a status quo in the efficiency of water loss in the delivery systems whereas the optimal water-loss scenario projected a better water-loss efficiency resulting in less withdrawals. The third scenario assumes that all wells will withdraw water at their full allocation level which is generally much more than reported withdrawals in 2013 or projected under the other two scenarios. Maps and summaries of heads and drawdowns are presented for nine confined aquifers. All the aquifers have areas with heads below sea level by 2040. Of the three scenarios, the drawdowns are most extreme in the full allocation scenarios; there are large areas of head decline greater than 20 feet in 5 of the 9 confined aquifers. The exceptions are the Vincentown aquifer, despite some areas of large drawdown in the vicinity of wells, and the three Potomac-Raritan-Magothy (PRM) aquifers where withdrawals are regulated by Critical Area restrictions. The nominal and optimal water-loss scenarios have some areas of head declines; most are less than 15 feet. The simulation of these scenarios shows some extensive areas of head recovery as well—especially in the aquifers that are regulated by the Critical Area restrictions. Budgets of inflow and outflow components were calculated for 44 hydrologic budget areas (HBAs). The budget analysis shows that the water movement is complex and varies based on the aquifer geometry and location of pumping wells. Flow components between the unconfined and confined parts of the system were summarized by HUC11 (hydrologic unit code 11) basins.

New Jersey↗

Strandlines from large floods on the Colorado River in Grand Canyon National Park, Arizona

Strandlines of peak-stage indicators (such as driftwood logs, woody debris, and trash) provide valuable data for understanding the maximum stage and extent of inundation during floods. A series of seven strandlines have been preserved along the Colorado River in Grand Canyon National Park, Arizona, USA. A survey and analysis of these strandlines was completed from the Colorado River at Lees Ferry, Ariz., gaging station to the Colorado River near Grand Canyon, Ariz., gaging station. Owing to the longitudinally discontinuous nature of the strandlines, several lines of evidence were used to determine the year of the flood associated with each strandline segment. This evidence included strandline relative vertical position, degree of peak-stage indicator weathering, datable trash drift, and map-view location. The seven distinct strandlines identified were deposited during floods with the following peak discharges (in cubic feet per second [ft 3 /s]) at the Colorado River at Lees Ferry, Ariz., gaging station (year of flood in parentheses): 210,000 ft 3 /s (1884), 170,000 ft 3 /s (1921), 125,000 ft 3 /s (1957), 108,000 ft 3 /s (1958), 97,000 ft 3 /s (1983), 52,500 ft 3 /s (1986), and 45,000 ft 3 /s (multiple events between 1996 and 2012). Stage-discharge relations were developed in areas where all, or most of the strandlines were present, and were compared to predicted stage-discharge relations from a one-dimensional flow model. River width exerted a strong control on these relations, with much greater stage change occurring for a given discharge change in narrower bedrock-dominated reaches than in wider reaches with more extensive channel-margin alluvium. This comprehensive dataset allows for the verification of model-predicted flood stage along the Colorado River in Grand Canyon National Park.

Arizona↗

Geologic studies in Alaska by the U.S. Geological Survey, 1997

The eight papers that follow continue the series of U.S. Geological Survey (USGS) reports on investigations in the geologic sciences in Alaska. The series presents new and sometimes preliminary findings that are of interest to earth scientists in academia, government, and industry; to land and resource managers; and to the general public. Reports presented in Geologic Studies in Alaska cover a broad spectrum of topics from all parts of the State (fig. 1), which serves to emphasize the diversity of USGS efforts to meet the Nation's needs for earth-science information in Alaska. The papers in this volume are organized under the topics Resources, Geologic Framework, and Environment and Climate. Such an organization is intended to reflect the scope and objectives of USGS programs currently active in Alaska. Resource papers include one that presents detailed observations from a Mississippian Zn-Pb-Cu-Ag occurrence in the Brooks Range (Werdon). Mineralogic, chemical, and isotopic data provide the basis for a proposed relationship between this vein-breccia deposit and the shale-hosted massive sulfide deposit type, which includes the active Red Dog mine in the western Brooks Range. Also included under the topic of Resources is a paper that presents geochemical and isotopic data from the Greens Creek and Woewodski Island volcanogenic massive sulfide deposits in southeastern Alaska (Newberry and Brew). The depositional environment and sedimentological setting of Tertiary coal beds in the Matanuska and Susitna Valleys is the focus of a third paper under the topic of Resources (Flores and others). Geologic Framework studies provide background information that is the scientific basis for present and future studies of the environment, mineral and energy resources, paleoclimate, and hazards in Alaska. One paper presents the results of sedimentologic and paleontologic comparisons of lower Paleozoic, deep-water-facies rock units in central Alaska (Dumoulin and others). The authors show which of these units are likely to correlate with one another, suggest likely source regions, and provide a structural restoration of units that have been fragmented by large fault motions. A second framework paper provides a map, rock descriptions, and chemical compositions of volcanic rocks in a newly recognized, geologically young volcanic center in the Aleutian volcanic arc (Hildreth and others). A third paper presents an interesting summary of gravity changes that occurred in south-central Alaska during the great earthquake of 1964 and for the following 25 years (Barnes). Gravity changes correlate with land-elevation changes in some cases, but not in others, which means that different processes are responsible for the gravity changes.

Alaska↗

Simulation of groundwater flow in the regional aquifer system on Long Island, New York, for pumping and recharge conditions in 2005–15

A three-dimensional groundwater-flow model was developed for the aquifer system of Long Island, New York, to evaluate (1) responses of the hydrologic system to changes in natural and anthropogenic hydraulic stresses, (2) the subsurface distribution of groundwater age, and (3) the regional-scale distribution of groundwater travel times and the source of water to fresh surface waters and coastal receiving waters. The model also provides the groundwater flow components used to define model boundaries for possible inset models used for local-scale analyses. The three-dimensional, groundwater flow model developed for this investigation uses the numerical code MODFLOW–NWT to represent steady-state conditions for average groundwater pumping and aquifer recharge for 2005–15. The particle-tracking algorithm MODPATH, which simulates advective transport in the aquifer, was used to estimate groundwater age, delineate the areas at the water table that contribute recharge to coastal and freshwater bodies, and estimate total travel times of water from the water table to discharge locations. A three-dimensional, 1-meter (3.3-foot) topobathymetric model was used to determine land-surface altitudes for the island and seabed altitudes for the surrounding coastal waters. The mapped extents and surface altitudes of major geologic units were compiled and used to develop a three-dimensional hydrogeologic framework of the aquifer system, including aquifers and confining units. Lithologic data from deep boreholes and previous aquifer-test results were used to estimate the three-dimensional distribution of hydraulic conductivity in principal aquifers. Natural recharge from precipitation was estimated for 2005–15 using a modified Thornthwaite-Mather methodology as implemented in a soil-water balance model. Components of anthropogenic recharge—wastewater return flow, storm water inflow, and inflow from leaky infrastructure—also were estimated for 2005–15. Groundwater withdrawals for various sources, including public water supply, industrial, remediation, and agricultural, were compiled or estimated for the same period. These data were incorporated into the model development to represent the aquifer system geometry, boundaries, and initial hydraulic properties of the regional aquifers and confining units within the Long Island aquifer system. Average hydraulic conditions—water levels and streamflows—for 2005–15 were estimated using existing data from the U.S. Geological Survey National Water Information System database. Model inputs were adjusted to best match average hydrologic conditions using inverse methods as implemented in the parameter-estimating software PEST. The calibrated model was used to simulate average hydrologic conditions in the aquifer system for 2005–15. About 656 cubic feet per second of water was withdrawn on average annually for 2005–15 for water supply and an average of about 349 cubic feet per second of water recharged the aquifer annually from return flow and leaky infrastructure. Parts of New York City have drawdowns exceeding 25 feet, mostly because of urbanization and associated large decreases in recharge rates. Large areas in the western part of the island have drawdowns exceeding 10 feet, mostly from large groundwater withdrawals and sewering, which largely eliminates wastewater return flow. Water-table altitudes in eastern parts of the island increased by more than 2 feet in some areas as a result of wastewater return flow in unsewered areas and changes in land use. Changes in streamflows show a similar pattern as water-table altitudes. Streamflows generally decrease in western parts of the island where there are large drawdowns and increase in eastern parts of the island where water-table altitudes increase.

New York↗

Hydrogeologic Framework and Ground-Water Budget of the Spokane Valley-Rathdrum Prairie Aquifer, Spokane County, Washington, and Bonner and Kootenai Counties, Idaho

The U.S. Geological Survey, in cooperation with the Idaho Department of Water Resources and Washington State Department of Ecology, investigated the hydrogeologic framework and ground-water budget of the Spokane Valley-Rathdrum Prairie (SVRP) aquifer located in northern Idaho and northeastern Washington. Descriptions of the hydrogeologic framework, water-budget components, and further data needs are provided. The SVRP aquifer, which covers about 370 square miles including the Rathdrum Prairie, Idaho, and the Spokane Valley and Hillyard Trough, Washington, is the sole source of drinking water for more than 500,000 residents. Continued growth, water-management issues, and potential effects on water availability and water quality in the aquifer and in the Spokane and Little Spokane Rivers have illustrated the need to better understand and manage the region's water resources. The SVRP aquifer consists mostly of gravels, cobbles, and boulders - deposited during a series of outburst floods resulting from repeated collapse of the ice dam that impounded ancient Glacial Lake Missoula. In most places, the SVRP aquifer is bounded by bedrock of pre-Tertiary granite or metasedimentary rocks, or Miocene basalt and associated sedimentary deposits. Discontinuous fine-grained layers are scattered throughout the SVRP aquifer at considerably different altitudes and with considerably different thicknesses. In the Hillyard Trough and the Little Spokane River Arm of the aquifer, a massive fine-grained layer with a top altitude ranging from about 1,500 to 1,700 feet and thickness ranging from about 100 to 200 feet separates the aquifer into upper and lower units. Most of the Spokane Valley part of the aquifer is devoid of fine-grained layers except near the margins of the valley and near the mouths of lakes. In the Rathdrum Prairie, multiple fine-grained layers are scattered throughout the aquifer with top altitudes ranging from about 1,700 to 2,400 feet with thicknesses ranging from 1 to more than 135 feet. The altitude of the base of the aquifer ranges from less than 1,800 feet near Lake Pend Oreille to less than 1,200 feet near the aquifer's outlet near Long Lake. The thickness of the aquifer is more than 800 feet in the northwestern part of the northern Rathdrum Prairie, through the West Channel area, and through the west-central part of the Rathdrum Prairie. In Washington, the areas of greatest thickness, more than 600 feet, are mapped in the central parts of the Spokane Valley, Spokane, and the Hillyard Trough. Recharge or inflow to the SVRP aquifer occurs from six main sources: the Spokane River, lakes, infiltration from precipitation over the aquifer, tributaries, infiltration from landscape irrigation and septic systems, and subsurface inflow. Discharge or outflow from the SVRP aquifer occurs from five main sources: the Spokane River, the Little Spokane River, pumpage, subsurface discharge to Long Lake, and infiltration of ground water to sewers. Total estimated mean annual inflow to and outflow from the SVRP aquifer is about 1,470 cubic feet per second. Several data needs were identified during this investigation that would improve the definition of the hydrogeologic framework and ground-water budget components for the SVRP aquifer study area. Deep drilling along the axis of the aquifer could determine the depth to the bottom of the aquifer where data are currently unavailable as well as identify the presence of fine-grained layers and their thickness. A more detailed analysis of the geologic and hydrologic setting near the southern ends of Spirit and Hoodoo Valleys could help determine the location of the ground-water divide between the two valleys and the Rathdrum Prairie. Better estimates of seepage into the aquifer from Coeur d'Alene Lake and Lake Pend Oreille and underflow from the aquifer to Long Lake would strengthen the recharge and discharge estimates of the aquifer. A hydrochemical study incorporating analyses of envi

Scientific Investigations Report↗

Hydrogeology and steady-state numerical simulation of groundwater flow in the Lost Creek Designated Ground Water Basin, Weld, Adams, and Arapahoe Counties, Colorado

The Lost Creek Designated Ground Water Basin (Lost Creek basin) is an important alluvial aquifer for irrigation, public supply, and domestic water uses in northeastern Colorado. Beginning in 2005, the U.S. Geological Survey, in cooperation with the Lost Creek Ground Water Management District and the Colorado Water Conservation Board, collected hydrologic data and constructed a steady-state numerical groundwater flow model of the Lost Creek basin. The model builds upon the work of previous investigators to provide an updated tool for simulating the potential effects of various hydrologic stresses on groundwater flow and evaluating possible aquifer-management strategies. As part of model development, the thickness and extent of regolith sediments in the basin were mapped, and data were collected concerning aquifer recharge beneath native grassland, nonirrigated agricultural fields, irrigated agricultural fields, and ephemeral stream channels. The thickness and extent of regolith in the Lost Creek basin indicate the presence of a 2- to 7-mile-wide buried paleovalley that extends along the Lost Creek basin from south to north, where it joins the alluvial valley of the South Platte River valley. Regolith that fills the paleovalley is as much as about 190 ft thick. Average annual recharge from infiltration of precipitation on native grassland and nonirrigated agricultural fields was estimated by using the chloride mass-balance method to range from 0.1 to 0.6 inch, which represents about 1-4 percent of long-term average precipitation. Average annual recharge from infiltration of ephemeral streamflow was estimated by using apparent downward velocities of chloride peaks to range from 5.7 to 8.2 inches. Average annual recharge beneath irrigated agricultural fields was estimated by using passive-wick lysimeters and a water-balance approach to range from 0 to 11.3 inches, depending on irrigation method, soil type, crop type, and the net quantity of irrigation water applied. Estimated average annual recharge beneath irrigated agricultural fields represents about 0-43 percent of net irrigation. The U.S. Geological Survey modular groundwater modeling program, MODFLOW-2000, was used to develop a steady-state groundwater flow model of the Lost Creek basin. Groundwater in the basin is simulated generally to flow from the basin margins toward the center of the basin and northward along the paleovalley. The largest source of inflow to the model occurs from recharge beneath flood- and sprinkler-irrigated agricultural fields (14,510 acre-feet per year [acre-ft/yr]), which represents 39.7 percent of total simulated inflow. Other substantial sources of inflow to the model are recharge from precipitation and stream-channel infiltration in nonirrigated areas (13,810 acre-ft/yr) seepage from Olds Reservoir (4,280 acre-ft/yr), and subsurface inflow from ditches and irrigated fields outside the model domain (2,490 acre-ft/yr), which contribute 37.7, 11.7, and 6.8 percent, respectively, of total inflow. The largest outflow from the model occurs from irrigation well withdrawals (26,760 acre-ft/yr), which represent 73.2 percent of total outflow. Groundwater discharge (6,640 acre-ft/yr) at the downgradient end of the Lost Creek basin represents 18.2 percent of total outflow, and evapotranspiration (3,140 acre-ft/yr) represents about 8.6 percent of total outflow.

Scientific Investigations Report↗

Flood-frequency characteristics of Wisconsin streams

Flood-frequency characteristics for 360 gaged sites on unregulated rural streams in Wisconsin are presented for percent annual exceedance probabilities ranging from 0.2 to 50 using a statewide skewness map developed for this report. Equations of the relations between flood-frequency and drainage-basin characteristics were developed by multiple-regression analyses. Flood-frequency characteristics for ungaged sites on unregulated, rural streams can be estimated by use of the equations presented in this report. The State was divided into eight areas of similar physiographic characteristics. The most significant basin characteristics are drainage area, soil saturated hydraulic conductivity, main-channel slope, and several land-use variables. The standard error of prediction for the equation for the 1-percent annual exceedance probability flood ranges from 56 to 70 percent for Wisconsin Streams; these values are larger than results presented in previous reports. The increase in the standard error of prediction is likely due to increased variability of the annual-peak discharges, resulting in increased variability in the magnitude of flood peaks at higher frequencies. For regulated streams, a graphical method for estimating flood-frequency characteristics was developed from the relations of discharge and drainage area for selected annual exceedance probabilities. Graphs for the major regulated streams in Wisconsin are presented in the report.

Wisconsin↗

Sequence stratigraphy, seismic stratigraphy, and seismic structures of the lower intermediate confining unit and most of the Floridan aquifer system, Broward County, Florida

Deep well injection and disposal of treated wastewater into the highly transmissive saline Boulder Zone in the lower part of the Floridan aquifer system began in 1971. The zone of injection is a highly transmissive hydrogeologic unit, the Boulder Zone, in the lower part of the Floridan aquifer system. Since the 1990s, however, treated wastewater injection into the Boulder Zone in southeastern Florida has been detected at three treated wastewater injection utilities in the brackish upper part of the Floridan aquifer system designated for potential use as drinking water. At a time when usage of the Boulder Zone for treated wastewater disposal is increasing and the utilization of the upper part of the Floridan aquifer system for drinking water is intensifying, there is an urgency to understand the nature of cross-formational fluid flow and identify possible fluid pathways from the lower to upper zones of the Floridan aquifer system. To better understand the hydrogeologic controls on groundwater movement through the Floridan aquifer system in southeastern Florida, the U.S. Geological Survey and the Broward County Environmental Planning and Community Resilience Division conducted a 3.5-year cooperative study from July 2012 to December 2015. The study characterizes the sequence stratigraphy, seismic stratigraphy, and seismic structures of the lower part of the intermediate confining unit aquifer and most of the Floridan aquifer system. Data obtained to meet the study objective include 80 miles of high-resolution, two-dimensional (2D), seismic-reflection profiles acquired from canals in eastern Broward County. These profiles have been used to characterize the sequence stratigraphy, seismic stratigraphy, and seismic structures in a 425-square-mile study area. Horizon mapping of the seismic-reflection profiles and additional data collection from well logs and cores or cuttings from 44 wells were focused on construction of three-dimensional (3D) visualizations of eight sequence stratigraphic cycles that compose the Eocene to Miocene Oldsmar, Avon Park, and Arcadia Formations. The mapping of these seismic-reflection and well data has produced a refined Cenozoic sequence stratigraphic, seismic stratigraphic, and hydrogeologic framework of southeastern Florida. The upward transition from the Oldsmar Formation to the Avon Park Formation and the Arcadia Formation embodies the evolution from (1) a tropical to subtropical, shallow-marine, carbonate platform, represented by the Oldsmar and Avon Park Formations, to (2) a broad, temperate, mixed carbonate-siliciclastic shallow marine shelf, represented by the lower part of the Arcadia Formation, and to (3) a temperate, distally steepened carbonate ramp represented by the upper part of the Arcadia Formation. In the study area, the depositional sequences and seismic sequences have a direct correlation with hydrogeologic units. The approximate upper boundary of four principal permeable units of the Floridan aquifer system (Upper Floridan aquifer, Avon Park permeable zone, uppermost major permeable zone of the Lower Floridan aquifer, and Boulder Zone) have sequence stratigraphic and seismic-reflection signatures that were identified on cross sections, mapped, or both, and therefore the sequence stratigraphy and seismic stratigraphy were used to guide the development of a refined spatial representation of these hydrogeologic units. In all cases, the permeability of the four permeable units is related to stratiform megaporosity generated by ancient dissolution of carbonate rock associated with subaerial exposure and unconformities at the upper surfaces of carbonate depositional cycles of several hierarchical scales ranging from high-frequency cycles to depositional sequences. Additionally, interparticle porosity also contributes substantially to the stratiform permeability in much of the Upper Floridan aquifer. Information from seismic stratigraphy allowed 3D geomodeling of hydrogeologic units—an approach never before applied to this area. Notably, the 3D geomodeling provided 3D visualizations and geocellular models of the depositional sequences, hydrostratigraphy, and structural features. The geocellular data could be used to update the hydrogeologic structure inherent to groundwater flow simulations that are designed to address the sustainability of the water resources of the Floridan aquifer system. Two kinds of pathways that could enable upward cross-formational flow of injected treated wastewater from the Boulder Zone have been identified in the 80 miles of high-resolution seismic data collected for this study: a near-vertical reverse fault and karst collapse structures. The single reverse fault, inferred to be of tectonic origin, is in extreme northeastern Broward County and has an offset of about 19 feet at the level of the Arcadia Formation. Most of the 17 karst collapse structures identified manifest as columniform, vertically stacked sagging seismic reflections that span early Eocene to Miocene age rocks equivalent to much of the Floridan aquifer system and the lower part of the overlying intermediate confining unit. In some cases, the seismic-sag structures extend upward into strata of Pliocene age. The seismic-sag structures are interpreted to have a semicircular shape in plan view on the basis of comparison to (1) other seismic-sag structures in southeastern Florida mapped with two 2D seismic cross lines or 3D data, (2) comparison to these structures located in other carbonate provinces, and (3) plausible extensional ring faults detected with multi-attribute analysis. The seismic-sag structures in the study area have heights as great as 2,500 vertical feet, though importantly, one spans about 7,800 feet. Both multi-attribute analysis and visual detection of offset of seismic reflections within the seismic-sag structures indicate faults and fractures are associated with many of the structures. Multi-attribute analysis highlighting chimney fluid pathways also indicates that the seismic-sag structures have a high probability for potential vertical cross-formational fluid flow along the faulted and fractured structures. A collapse of the seismic-sag structures within a deep burial setting evokes an origin related to hypogenic karst processes by ascending flow of subsurface fluids. In addition, paleo-epigenic karst related to major regional subaerial unconformities within the Florida Platform generated collapse structures (paleo-sinkholes) that are much smaller in scale than the cross-formational seismic-sag structures.

Florida↗

Yield of bedrock wells in the Nashoba terrane, central and eastern Massachusetts

The yield of bedrock wells in the fractured-bedrock aquifers of the Nashoba terrane and surrounding area, central and eastern Massachusetts, was investigated with analyses of existing data. Reported well yield was compiled for 7,287 wells from Massachusetts Department of Environmental Protection and U.S. Geological Survey databases. Yield of these wells ranged from 0.04 to 625 gallons per minute. In a comparison with data from 103 supply wells, yield and specific capacity from aquifer tests were well correlated, indicating that reported well yield was a reasonable measure of aquifer characteristics in the study area. Statistically significant relations were determined between well yield and a number of cultural and hydrogeologic factors. Cultural variables included intended water use, well depth, year of construction, and method of yield measurement. Bedrock geology, topography, surficial geology, and proximity to surface waters were statistically significant hydrogeologic factors. Yield of wells was higher in areas of granites, mafic intrusive rocks, and amphibolites than in areas of schists and gneisses or pelitic rocks; higher in valleys and low-slope areas than on hills, ridges, or high slopes; higher in areas overlain by stratified glacial deposits than in areas overlain by till; and higher in close proximity to streams, ponds, and wetlands than at greater distances from these surface-water features. Proximity to mapped faults and to lineaments from aerial photographs also were related to well yield by some measures in three quadrangles in the study area. Although the statistical significance of these relations was high, their predictive power was low, and these relations explained little of the variability in the well-yield data. Similar results were determined from a multivariate regression analysis. Multivariate regression models for the Nashoba terrane and for a three-quadrangle subarea included, as significant variables, many of the cultural and hydrogeologic factors that were individually related to well yield, in ways that are consistent with conceptual understanding of their effects, but the models explained only 21 percent (regional model for the entire terrane) and 30 percent (quadrangle model) of the overall variance in yield. Moreover, most of the explained variance was due to well characteristics rather than hydrogeologic factors. Hydrogeologic factors such as topography and geology are likely important. However, the overall high variability in the well-yield data, which results from the high variability in aquifer hydraulic properties as well as from limitations of the dataset, would make it difficult to use hydrogeologic factors to predict well yield in the study area. Geostatistical analysis (variograms), on the other hand, indicated that, although highly variable, the well-yield data are spatially correlated. The spatial continuity appears greater in the northeast-southwest direction and less in the southeast-northwest direction, directions that are parallel and perpendicular, respectively, to the regional geologic structural trends. Geostatistical analysis (kriging), used to estimate yield values throughout the study area, identified regional-scale areas of higher and lower yield that may be related to regional structural features—in particular, to a northeast-southwest trending regional fault zone within the Nashoba terrane. It also would be difficult to use kriging to predict yield at specific locations, however, because of the spatial variability in yield, particularly at small scales. The regional-scale analyses in this study, both with hydrogeologic variables and geostatistics, provide a context for understanding the variability in well yield, rather a basis for precise predictions, and site-specific information would be needed to understand local conditions.

Massachusetts↗

Groundwater recharge estimates for Maine using a Soil-Water-Balance model—25-year average, range, and uncertainty, 1991 to 2015

To address the lack of information on the spatial and temporal variability of recharge to groundwater systems in Maine, a study was initiated in cooperation with the Maine Geological Survey to use the U.S. Geological Survey Soil-Water-Balance model to evaluate annual average potential recharge across the State over a 25-year period from 1991 to 2015. The Maine Soil-Water-Balance model was calibrated using annual observations of recharge, runoff, and evapotranspiration for 32 calibration watersheds in the State during 2001–12 (902 total observations). Observations of recharge, runoff, and evapotranspiration were developed for each watershed to reduce the possibility of nonunique combinations of model parameters during the calibration. The Maine Soil-Water-Balance model was run using an optional evapotranspiration calculation method that provides more control for calibration than the standard method. The model was calibrated using the Parameter ESTimation software suite. The overall mean model error (average of all annual residuals for recharge, runoff, and precipitation) was 0.39 inch. The mean of the absolute value of the residuals, or the mean absolute error, was 2.32 inches. The root mean squared error for the calibrated model overall was 3.14 inches. Statistical tests indicated that the model residuals are normally distributed. To determine the potential uncertainty in the median annual potential recharge that results from uncertainty in the parameters as they relate to information contained in the observations, 300 alternate model realizations were run, and the standard deviation of the median potential recharge value at every pixel was calculated. Simulated 25-year median potential recharge across the State is widely variable; this variability closely follows patterns of precipitation, with additional variability contributed by the patchwork nature of the combinations of land-use class and hydrologic soil group inputs, and distribution of available water capacity in the soil across the State. Overall, the 25-year median annual potential recharge across the State is 7.5 inches, ranging from a low of about 5 inches to over 30 inches. The statewide range in the 25-year minimum values is from just over 2 inches to just over 20 inches. The statewide range in the 25-year maximum potential recharge is between 15 and 48 inches per year. The model areas with the highest simulated median potential recharge include areas underlain by type A soils (sandy and well drained), particularly those that also have land uses with low or little vegetation (blueberry barrens, developed, open space, scrub/shrub, and cropland, for example). The potential recharge values for these areas are similar to previously published values for comparable soil types. The 25-year average potential recharge grids were compared to recharge evaluated through groundwater-flow models or other methods in four hydrogeologic settings at six study areas in the State. A key factor in the ability of the Soil-Water-Balance model to reproduce the earlier study results was whether the available water-capacity data were an appropriate match for the hydrologic soil groups. The Maine Soil-Water-Balance model does a good job in representing an accurate potential recharge under circumstances where the surficial mapped soils extend below the surface to the water-table aquifer and where the available water-capacity data are in an appropriate range for the hydrologic soil group. One hydrogeologic setting that was challenging for the model was where a silt and clay layer was below a shallow soil unit that did not have available water-capacity data that were appropriate for the hydrologic soil group. In these cases, typically the available water-capacity data were very low, not accounting for the impedance of water flow provided by the underlying soil. The model also does not simulate well areas where bedrock surfaces are above the water table but below the plant rooting zone. The data products accompanying this report are intended to be used to provide first-cut estimates of recharge for geographic areas no smaller than the smallest watersheds used in the calibration of the model—or about 1.5 square miles. It is recommended that the grids are used to calculate an area-wide average potential recharge for any given area of study, and an uncertainty around the mean should be calculated from the standard deviation grid at the same time.

Maine↗

Geophysical analysis of the Salmon Peak Formation near Amistad Reservoir Dam, Val Verde County, Texas, and Coahuila, Mexico, March 2006, to aid in piezometer placement

Since 1992, numerous sinkholes have developed northwest of the Amistad Reservoir dam on the Rio Grande. Increases in the discharge of springs south of the dam, on the western side of the Rio Grande, in Coahuila, Mexico, have been documented. In 1995 the Mexico Section of the International Boundary and Water Commission (IBWC) completed a study of the western embankment (Coahuila, Mexico) of the dam that included surface geophysics, borehole geophysics, and installation of piezometers to learn more about subsurface conditions. As part of a 5-year safety inspection in 2005, technical advisors recommended that one line of similarly constructed piezometers be installed on the eastern embankment (Val Verde County, Texas) of the dam for comparison of water levels (potentiometric head) on both the western and eastern embankments of Amistad Reservoir dam. To provide technical assistance for the horizontal and vertical placement of piezometers on the eastern embankment of Amistad Reservoir dam, the U.S. Geological Survey, in cooperation with the U.S. Section of the IBWC, conducted a study along both the western and eastern embankments of Amistad Reservoir dam. The study involved an integrated approach using surface and borehole geophysical methods. In the western embankment investigation, geological and geophysical characteristics that indicate relatively large water-yielding properties of the Salmon Peak Formation were identified. The direct-current (DC) resistivity method was selected as the surface geophysical reconnaissance technique to correlate relatively large water-yielding properties of the Salmon Peak Formation, identified from analysis of borehole geophysical logs, with variations in subsurface resistivity. The dipole-dipole array and the reciprocal Schlumberger array were selected as the most applicable DC resistivity arrays. Two resistivity units were identified in both the dipole-dipole array data and the reciprocal Schlumberger array data along DC resistivity profiles on both embankments. Resistivity unit 1 generally is of relatively low resistivity, ranging from 45 to 150 ohm-meters compared with resistivity unit 2, which ranges from 120 to 345 ohm-meters (depending on the DC array type). The presence of mapped sinkholes in the reservoir north of the western embankment study area and the zone of increased water content (as indicated by zones of low neutron log count rates in nearby piezometers) leads to the conclusion that resistivity unit 1 is a preferential flow path where surface water from Amistad Reservoir is forced into the ground-water system (because of increased head from the reservoir). In the eastern embankment investigation, trends in the spatial distribution of sinkholes and the occurrence of weathered zones were identified from geologic descriptions of cores. The correlation of surface geophysical DC resistivity, historical lithologic data, and general trend of documented sinkholes along the eastern end of the eastern embankment profile were used to justify further exploration (drilling of piezometers) in the eastern expression of resistivity unit 1. The spatial location of the piezometers and the screened intervals were selected to best match the locations of the screened intervals of the western embankment piezometers. Six piezometers were installed on the eastern embankment and logged using borehole geophysical techniques. Surface DC resistivity sections superimposed on the resistivity logs for two piezometers indicate three discernible resistivity units that correlate with resistivity units 2, 1, and 2, respectively, identified in the western embankment study area. Resistivity units 1 and 2 in the DC resistivity profiles generally correspond with low and high resistivity zones, respectively, on the normal and lateral resistivity logs collected in the nearby piezometers at the time of installation.

Scientific Investigations Report↗

Potentiometric surface of the Ozark aquifer in northern Arkansas, 2010

The Ozark aquifer in northern Arkansas is composed of dolomite, limestone, sandstone, and shale of Late Cambrian to Middle Devonian age and ranges in thickness from approximately 1,100 feet to more than 4,000 feet. Hydrologically, the aquifer is complex, characterized by discrete and discontinuous flow components with large variations in permeability. The potentiometric-surface map, based on 56 well and 5 spring water-level measurements made in 2010 in Arkansas and Missouri, has a maximum water-level altitude measurement of 1,174 feet in Carroll County and a minimum water-level altitude measurement of 120 feet in Randolph County. Regionally, the flow within the aquifer is to the south and southeast in the eastern and central part of the study area and to the west, northwest, and north in the western part of the study area. Water-level altitudes changed 0.5 feet or less in 31 out of 56 wells measured between 2007 and 2010. Despite rapidly increasing population within the study area, the increase appears to have minimal effect on groundwater levels, although the effect may have been minimized by the development and use of surface-water distribution infrastructure, suggesting that most of the incoming populations are fulfilling their water needs from surface-water sources. The conversion of some users from groundwater to surface water may be allowing water levels in some wells to recover (rise) or decline at a slower rate in some areas such as in Benton, Carroll, and Washington Counties.

Arkansas↗

Physical subdivision and description of the water-bearing sediments of the Santa Clara Valley, California

A thick Quaternary alluvial section fills a sedimentary basin beneath the Santa Clara Valley, California, located within the San Andreas Fault system at the south end of San Francisco Bay. This section consists of an upper sequence about 1,000 feet thick containing eight sedimentary cycles and a lower fine-grained unit as thick as several hundred feet. Together these constitute the Quaternary Santa Clara Basin. The section overlies an irregular unconformity with more than 1,200 feet of relief cut into the underlying bedrock. This stratigraphy is determined through study of new wells and seismic reflection profiles, together with a sample of the many thousands of water wells in the valley. It represents a major change and improvement in understanding of the basin, particularly with regard to the upper cyclic sequence, which forms a large groundwater system that is an important resource in the San Francisco Bay region. Each of the eight sedimentary cycles consists of a coarse-grained bottom interval overlain by a fine-grained top, with the coarse bottom forming a permeable sheet that is more or less continuous around the basin and the fine top forming a similarly extensive, relatively impermeable confining layer. This stratigraphic organization contrasts with most previous views, which have considered the coarse sediment in the basin to occur as scattered, discrete lenses and (or) sinuous channel sands, all embedded in a predominantly fine-grained section. Temperature logs in several wells demonstrate that the fine cycle tops do limit vertical movement of groundwater, although this may not be the case where those tops are thin to perhaps locally absent around parts of the basin margin. Age control has been obtained from previous work, in which the sedimentary cycles were correlated with the marine oxygen isotope record and the ages of two deeper Quaternary unconformities were estimated, and from detailed paleomagnetic study of cores from the new wells by E.A. Mankinen. Despite careful search of the cores, very few fossils were found, and none that are helpful in subdividing the section. No tephra (volcanic ash) was recovered, and the few carbon samples found and dated radiometrically are limited to the upper 120 feet of the section. The upper cyclic section ranges in age from 0 to somewhat older than 718 thousand years (ka), and the lower fine-grained section lies between unconformities with estimated ages of 950 and 1500 ka. Reflections in the seismic profiles indicate that layering in the basin is subparallel to the ground surface, and this fact, together with the continuous stratigraphic detail provided by geophysical logs of the new wells, allows the confident interwell correlation required to delineate the sedimentary cycles. The sequence of layers within any one cycle tends to persist laterally between the wells in the dataset, which are spaced 1 to 3 km apart, with most changes occurring gradually. The eight cycles, in contrast, tend to differ from each other in the details of their internal organization. Maps and cross sections show the elevations of cycle boundaries and the underlying bedrock surface, the varying thicknesses of the cycles and of their fine tops and coarse bottoms, and the aggregate thickness of coarse layers in those bottom intervals. Coarse sediment is more abundant toward some parts of the basin margin and in the southern part of the basin. Cycle boundary surfaces are relatively smooth, and their shapes are consistent with having been intercycle topographic surfaces. The underlying bedrock surface has a relief of more than 1,200 feet and deepens toward the center of the basin and the west edge of the fault-bounded Evergreen Basin, which is concealed beneath the east side of the Quaternary basin. The absence of consistent abrupt changes in thicknesses or boundary elevations across the basin or in cross section indicates that the interior of the basin is largely unfaulted, with the Silver Creek strand of the San Andreas system at the west edge of the Evergreen Basin being the sole exception. The east and west margins of the Santa Clara Basin, in contrast, are marked by reverse and thrust fault systems.

California↗

Comparison of Surface Flow Features from Lidar-Derived Digital Elevation Models with Historical Elevation and Hydrography Data for Minnehaha County, South Dakota

The U.S. Geological Survey (USGS) has taken the lead in the creation of a valuable remote sensing product by incorporating digital elevation models (DEMs) derived from Light Detection and Ranging (lidar) into the National Elevation Dataset (NED), the elevation layer of 'The National Map'. High-resolution lidar-derived DEMs provide the accuracy needed to systematically quantify and fully integrate surface flow including flow direction, flow accumulation, sinks, slope, and a dense drainage network. In 2008, 1-meter resolution lidar data were acquired in Minnehaha County, South Dakota. The acquisition was a collaborative effort between Minnehaha County, the city of Sioux Falls, and the USGS Earth Resources Observation and Science (EROS) Center. With the newly acquired lidar data, USGS scientists generated high-resolution DEMs and surface flow features. This report compares lidar-derived surface flow features in Minnehaha County to 30- and 10-meter elevation data previously incorporated in the NED and ancillary hydrography datasets. Surface flow features generated from lidar-derived DEMs are consistently integrated with elevation and are important in understanding surface-water movement to better detect surface-water runoff, flood inundation, and erosion. Many topographic and hydrologic applications will benefit from the increased availability of accurate, high-quality, and high-resolution surface-water data. The remotely sensed data provide topographic information and data integration capabilities needed for meeting current and future human and environmental needs.

Scientific Investigations Report↗