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

Hydrologic assessment of the Edwin B. Forsythe National Wildlife Refuge

The Edwin B. Forsythe National Wildlife Refuge (hereafter Forsythe refuge or the refuge) is situated along the central New Jersey coast and provides a mixture of freshwater and saltwater habitats for numerous bird, wildlife, and plant species. Little data and information were previously available regarding the freshwater dynamics that support the refuge’s ecosystems. In cooperation with the U.S. Fish and Wildlife Service, the U.S. Geological Survey conducted an assessment of the hydrologic resources and processes in the refuge and surrounding areas to provide baseline information for evaluating restoration projects and future changes in the hydrologic system associated with climate change and other anthropogenic stressors. During spring 2015, water levels were measured at groundwater and surface-water sites in and near the Forsythe refuge. These water-level measurements, along with surface-water elevations obtained from digital elevation models, were used to construct water-table-elevation and depth-to-water maps of the refuge and surrounding areas. Water-table elevations in the refuge ranged from sea level to approximately 65 feet above sea level; in most of the refuge, the water-table elevation was within 3 feet of sea level. The water-table-elevation map indicates that the direction of shallow groundwater flow at the regional scale is generally from west to east (much of it from the northwest to the southeast), and groundwater moves downgradient from the uplands toward major groundwater discharge areas consisting of coastal streams and wetlands. The depth to water is estimated to be less than 2 feet for approximately 86 percent of the refuge, which coincides closely with the percentage of wetland area in the refuge. Depth to water in excess of 20 feet below land surface is limited to higher elevation areas of the refuge. Streamflow data collected at continuous-record streamgages and partial-record stations within the Mullica-Toms Basin were summarized. Hydrograph separation of streamflow data for eight streamgages (2004–13) reveals that base flow accounts for 68–94 percent of streamflow in basins upstream from the refuge. The high base-flow inputs underscore the importance of groundwater as a source of freshwater that supports both the streams that flow into the refuge and the hydroecology of the contributing basins. Mean annual flow typically ranged from 1.7 to 2.1 cubic feet per second per square mile at the streamgages (2004–13) and between 1.2 and 2.3 cubic feet per second per square mile at the partial-record stations (1965–2015) but was notably greater or lower than these ranges at several stations. Mean annual water budgets were estimated for multiple regions of the refuge for 2004–13 using data compiled from nearby meteorological stations and groundwater flows derived from previously calibrated groundwater-flow models. Precipitation, groundwater recharge, and evapotranspiration were estimated from available data; direct runoff was calculated as the residual component of the water balance. Groundwater recharge rates were greatest in the upland-dominated areas of the refuge with estimates of 14.4 to 18.9 inches per year, which are equivalent to 30 to 40 percent of precipitation. Groundwater recharge rates were nearly zero in the central coastal areas because these areas are major groundwater discharge zones, the water table is near land surface, the subsurface is close to saturation and cannot accept much recharge, and much of the area is underlain by thick marsh deposits likely with low permeability. Estimates of evapotranspiration varied from about 26 inches per year in the upland-dominated areas to more than 35 inches per year in the coastal wetlands, equivalent to 55–79 percent of mean annual precipitation, indicating that it is a major component of the hydrodynamics of the Forsythe refuge. On the basis of output from previously calibrated groundwater-flow models, nearly all of the groundwater exiting the surficial aquifer system in the central coastal areas of the refuge is discharged to wetlands, which highlights the importance of groundwater discharge in supporting the ecosystems of the Forsythe refuge. In the central coastal areas, horizontal flow contributes more than 90 percent of the groundwater flow to the surficial system, indicating that the upbasin areas are a substantial source of water that ultimately discharges to streams and wetlands in the refuge.

New Jersey↗

Spatial Distribution of Ground-Water Recharge Estimated with a Water-Budget Method for the Jordan Creek Watershed, Lehigh County, Pennsylvania

This report presents the results of a study by the U.S. Geological Survey, in cooperation with the Pennsylvania Geological Survey, to illustrate a water-budget method for mapping the spatial distribution of ground-water recharge for a 76-square-mile part of the Jordan Creek watershed, northwest of Allentown, in Lehigh County, Pennsylvania. Recharge was estimated by using the Hydrological Evaluation of Landfill Performance (HELP) water-budget model for 577 landscape units in Jordan Creek watershed, delineated on the basis of their soils, land use/land cover, and mean annual precipitation during 1951-2000. The water-budget model routes precipitation falling on each landscape unit to components of evapotranspiration, surface runoff, storage, and vertical percolation (recharge) for a five-layer soil column on a daily basis. The spatial distribution of mean annual recharge during 1951-2000 for each landscape unit was mapped by the use of a geographic information system. Recharge simulated by the water-budget model in Jordan Creek watershed during 1951-2000 averaged 12.3 inches per year and ranged by landscape unit from 0.11 to 17.05 inches per year. Mean annual recharge during 1951-2000 simulated by the water-budget model was most sensitive to changes to input values for precipitation and runoff-curve number. Mean annual recharge values for the crop, forest, pasture, and low-density urban land-use/land-cover classes were similar (11.2 to 12.2 inches per year) but were substantially less for high-density urban (6.8 inches per year), herbaceous wetlands (2.5 inches per year), and forested wetlands (1.3 inches per year). Recharge rates simulated for the crop, forest, pasture, and low-density urban land-cover classes were similar because those land-use/land-cover classes are represented in the model with parameter values that either did not significantly affect simulated recharge or tended to have offsetting effects on recharge. For example, for landscapes with forest land cover, values of runoff-curve number assigned to the model were smaller than for other land-use/land-cover classes (causing more recharge and less runoff), but the maximum depth of evapotranspiration was larger than for other land-use/ land-cover classes because of deeper root penetration in forests (causing more evapotranspiration and less recharge). The smaller simulated recharge for high-density urban and wetland land-use/land-cover classes was caused by the large values of runoff-curve number (greater than 90) assigned to those classes. The large runoff-curve number, however, certainly is not realistic for all wetlands; some wetlands act as areas of ground-water discharge and some as areas of recharge. Simulated mean annual recharge computed by the water-budget model for the 53-square-mile part of the watershed upstream from the streamflow-gaging station near Schnecksville was compared to estimates of recharge and base flow determined by analysis of streamflow records from 1967 to 2000. The mean annual recharge of 12.4 inches per year simulated by the water-budget method for 1967-2000 was less than estimates of mean annual recharge of 19.3 inches per year computed from the RORA computer program and base flow computed by the PART computer program (15.1 inches per year). In theory, the water-budget method provides a practical tool for estimating differences in recharge at local scales of interest, and the watershed- average recharge rate of 12.4 inches per year computed by the method is reasonable. However, the mean annual surface runoff of 4.5 inches per year simulated by the model is unrealistically small. The sum of surface runoff and recharge simulated by the water-budget model (16.9 inches per year) is 7 inches per year less than the streamflow measured at the gaging station near Schnecksville (23.9 inches per year) during 1967-2000, indicating that evapotranspiration is overestimated by the water-budget model by that amount. This discrepancy ca

Scientific Investigations Report↗

Hydrogeologic framework and occurrence, movement, and chemical characterization of groundwater in Dixie Valley, west-central Nevada

Dixie Valley, a primarily undeveloped basin in west-central Nevada, is being considered for groundwater exportation. Proposed pumping would occur from the basin-fill aquifer. In response to proposed exportation, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation and Churchill County, conducted a study to improve the understanding of groundwater resources in Dixie Valley. The objective of this report is to characterize the hydrogeologic framework, the occurrence and movement of groundwater, the general water quality of the basin-fill aquifer, and the potential mixing between basin-fill and geothermal aquifers in Dixie Valley. Various types of geologic, hydrologic, and geochemical data were compiled from previous studies and collected in support of this study. Hydrogeologic units in Dixie Valley were defined to characterize rocks and sediments with similar lithologies and hydraulic properties influencing groundwater flow. Hydraulic properties of the basin-fill deposits were characterized by transmissivity estimated from aquifer tests and specific-capacity tests. Groundwater-level measurements and hydrogeologic-unit data were combined to create a potentiometric surface map and to characterize groundwater occurrence and movement. Subsurface inflow from adjacent valleys into Dixie Valley through the basin-fill aquifer was evaluated using hydraulic gradients and Darcy flux computations. The chemical signature and groundwater quality of the Dixie Valley basin-fill aquifer, and potential mixing between basin-fill and geothermal aquifers, were evaluated using chemical data collected from wells and springs during the current study and from previous investigations. Dixie Valley is the terminus of the Dixie Valley flow system, which includes Pleasant, Jersey, Fairview, Stingaree, Cowkick, and Eastgate Valleys. The freshwater aquifer in the study area is composed of unconsolidated basin-fill deposits of Quaternary age. The basin-fill hydrogeologic unit can be several orders of magnitude more transmissive than surrounding and underlying consolidated rocks and Dixie Valley playa deposits. Transmissivity estimates in the basin fill throughout Dixie Valley ranged from 30 to 45,500 feet squared per day; however, a single transmissivity value of 0.1 foot squared per day was estimated for playa deposits. Groundwater generally flows from the mountain range uplands toward the central valley lowlands and eventually discharges near the playa edge. Potentiometric contours east and west of the playa indicate that groundwater is moving eastward from the Stillwater Range and westward from the Clan Alpine Mountains toward the playa. Similarly, groundwater flows from the southern and northern basin boundaries toward the basin center. Subsurface groundwater flow likely enters Dixie Valley from Fairview and Stingaree Valleys in the south and from Jersey and Pleasant Valleys in the north, but groundwater connections through basin-fill deposits were present only across the Fairview and Jersey Valley divides. Annual subsurface inflow from Fairview and Jersey Valleys ranges from 700 to 1,300 acre-feet per year and from 1,800 to 2,300 acre-feet per year, respectively. Groundwater flow between Dixie, Stingaree, and Pleasant Valleys could occur through less transmissive consolidated rocks, but only flow through basin fill was estimated in this study. Groundwater in the playa is distinct from the freshwater, basin-fill aquifer. Groundwater mixing between basin-fill and playa groundwater systems is physically limited by transmissivity contrasts of about four orders of magnitude. Total dissolved solids in playa deposit groundwater are nearly 440 times greater than total dissolved solids in the basin-fill groundwater. These distinctive physical and chemical flow restrictions indicate that groundwater interaction between the basin fill and playa sediments was minimal during this study period (water years 2009–11). Groundwater in Dixie Valley generally can be characterized as a sodium bicarbonate type, with greater proportions of chloride north of the Dixie Valley playa, and greater proportions of sulfate south of the playa. Analysis of major ion water chemistry data sampled during the study period indicates that groundwater north and south of Township 22N differ chemically. Dixie Valley groundwater quality is marginal when compared with national primary and secondary drinking-water standards. Arsenic and fluoride concentrations exceed primary drinking water standards, and total dissolved solids and manganese concentrations exceed secondary drinking water standards in samples collected during this study. High concentrations of boron and tungsten also were observed. Chemical comparisons between basin-fill and geothermal aquifer water indicate that most basin-fill groundwater sampled could contain 10–20 percent geothermal water. Geothermal indicators such as high temperature, lithium, boron, chloride, and silica suggest that mixing occurs in many wells that tap the basin-fill aquifer, particularly on the north, south, and west sides of the basin. Magnesium-lithium geothermometers indicate that some basin-fill aquifer water sampled for the current study likely originates from water that was heated above background mountain-block recharge temperatures (between 3 and 15 degrees Celsius), highlighting the influence of mixing with warm water that was possibly derived from geothermal sources.

Nevada↗

Hydrogeology of the Mogollon Highlands, central Arizona

The Mogollon Highlands, 4,855 square miles of rugged, mountainous terrain at the southern edge of the Colorado Plateau in central Arizona, is characterized by a bedrock-dominated hydrologic system that results in an incompletely integrated regional ground-water system, flashy streamflow, and various local water-bearing zones that are sensitive to drought. Increased demand on the water resources of the area as a result of recreational activities and population growth have made necessary an increased understanding of the hydrogeology of the region. The U.S. Geological Survey conducted a study of the geology and hydrology of the region in cooperation with the Arizona Department of Water Resources under the auspices of the Arizona Rural Watershed Initiative, a program launched in 1998 to assist rural areas in dealing with water-resources issues. The study involved the analysis of geologic maps, surface-water and ground-water flow, and water and rock chemical data and spatial relationships to characterize the hydrogeologic framework. The study area includes the southwestern corner of the Colorado Plateau and the Mogollon Rim, which is the eroded edge of the plateau. A 3,000- to 4,000-foot sequence of early to late Paleozoic sedimentary rocks forms the generally south-facing scarp of the Mogollon Rim. The area adjacent to the edge of the Mogollon Rim is an erosional landscape of rolling, step-like terrain exposing Proterozoic metamorphic and granitic rocks. Farther south, the Sierra Ancha and Mazatzal Mountain ranges, which are composed of various Proterozoic rocks, flank an alluvial basin filled with late Cenozoic sediments and volcanic flows. Eight streams with perennial to intermittent to ephemeral flow drain upland regions of the Mogollon Rim and flow into the Salt River on the southern boundary or the Verde River on the western boundary. Ground-water flow paths generally are controlled by large-scale fracture systems or by karst features in carbonate rocks. Stream channels are also largely controlled by structural features, such as regional joint or fault systems. Precipitation, which shows considerable variability in amount and intensity, recharges the ground-water system along the crest of the Mogollon Rim and to a lesser extent along the crests and flanks of the rim and the Mazatzal Mountains and Sierra Ancha. Flashy runoff in the mainly bedrock stream channels is typical. Springs are distributed throughout the region, typically discharging at or above the contact of variably permeable formations along the face of the Mogollon Rim with a scattering of low-discharge springs in the Proterozoic rocks below the rim. The surface of the Colorado Plateau is the primary recharge area for the C aquifer in which ground-water flows north toward the Little Colorado River and south toward the Mogollon Highlands. Within the study area, flow from the C aquifer primarily discharges from large, stable springs in the upper East Verde River, Tonto Creek, and Canyon Creek Basins along the top of the Mogollon Rim and to the west as base flow in West Clear Creek. On the basis of chemical evidence and the distribution and flow characteristics of springs and perennial streams, the C aquifer is also the source of water for the limestone aquifer that discharges from carbonate rocks near the base of the Mogollon Rim. Vertical flow from the C aquifer, the base of which is in the Schnebly Hill Formation, recharges the limestone aquifer that discharges mainly at Fossil Springs in the western part of the study area and as base flow in Cibecue Creek on the eastern edge of the study area. Local, generally shallow aquifers of variable productivity occur in plateau and mesa-capping basalts in the sedimentary rocks of the Schnebly Hill and Supai Formations, in fractured zones of the Proterozoic Payson granite, and in the alluvium of the lower Tonto Creek Basin. Where time series data exist, such water-bearing zones are shown to be sensitive to short-

Scientific Investigations Report↗

Occurrence of Pesticides in Ground Water of Wyoming, 1995-2006

Little existing information was available describing pesticide occurrence in ground water of Wyoming, so the U.S. Geological Survey, in cooperation with the Wyoming Department of Agriculture and the Wyoming Department of Environmental Quality on behalf of the Wyoming Ground-water and Pesticides Strategy Committee, collected ground-water samples twice (during late summer/early fall and spring) from 296 wells during 1995-2006 to characterize pesticide occurrence. Sampling focused on the State's ground water that was mapped as the most vulnerable to pesticide contamination because of either inherent hydrogeologic sensitivity (for example, shallow water table or highly permeable aquifer materials) or a combination of sensitivity and associated land use. Because of variations in reporting limits among different compounds and for the same compound during this study, pesticide detections were recensored to two different assessment levels to facilitate qualitative and quantitative examination of pesticide detection frequencies - a common assessment level (CAL) of 0.07 microgram per liter and an assessment level that differed by compound, referred to herein as a compound-specific assessment level (CSAL). Because of severe data censoring (fewer than 50 percent of the data are greater than laboratory reporting limits), categorical statistical methods were used exclusively for quantitative comparisons of pesticide detection frequencies between seasons and among various natural and anthropogenic (human-related) characteristics. One or more pesticides were detected at concentrations greater than the CAL in water from about 23 percent of wells sampled in the fall and from about 22 percent of wells sampled in the spring. Mixtures of two or more pesticides occurred at concentrations greater than the CAL in about 9 percent of wells sampled in the fall and in about 10 percent of wells sampled in the spring. At least 74 percent of pesticides detected were classified as herbicides. Considering only detections using the CAL, triazine pesticides were detected much more frequently than all other pesticide classes, and the number of different pesticides classified as triazines was the largest of all classes. More pesticides were detected at concentrations greater than the CSALs in water from wells sampled in the fall (28 different pesticides) than in the spring (21 different pesticides). Many pesticides were detected infrequently as nearly one-half of pesticides detected in the fall and spring at concentrations greater than the CSALs were detected only in one well. Using the CSALs for pesticides analyzed for in 11 or more wells, only five pesticides (atrazine, prometon, tebuthiuron, picloram, and 3,4-dichloroaniline, listed in order of decreasing detection frequency) were each detected in water from more than 5 percent of sampled wells. Atrazine was the pesticide detected most frequently at concentrations greater than the CSAL. Concentrations of detected pesticides generally were small (less than 1 microgram per liter), although many infrequent detections at larger concentrations were noted. All detected pesticide concentrations were smaller than U.S. Environmental Protection Agency (USEPA) drinking-water standards or applicable health advisories. Most concentrations were at least an order of magnitude smaller; however, many pesticides did not have standards or advisories. The largest percentage of pesticide detections and the largest number of different pesticides detected were in samples from wells located in the Bighorn Basin and High Plains/ Casper Arch geographic areas of north-central and southeastern Wyoming. Prometon was the only pesticide detected in all eight geographic areas of the State. Pesticides were detected much more frequently in samples from wells located in predominantly urban areas than in samples from wells located in predominantly agricultural or mixed areas. Pesticides were detected distinctly less often in sa

Scientific Investigations Report↗

Aquifer geometry, lithology, and water levels in the Anza–Terwilliger area—2013, Riverside and San Diego Counties, California

The population of the Anza–Terwilliger area relies solely on groundwater pumped from the alluvial deposits and surrounding bedrock formations for water supply. The size, characteristics, and current conditions of the aquifer system in the Anza–Terwilliger area are poorly understood, however. In response to these concerns, the U.S. Geological Survey, in cooperation with the High Country Conservancy and Rancho California Water District, undertook a study to (1) improve mapping of groundwater basin geometry and lithology and (2) to resume groundwater-level monitoring last done during 2004–07 in the Anza–Terwilliger area. Inversion of gravity data, including new data collected for this study, was done to estimate the thickness of the alluvial deposits that form the Cahuilla and Terwilliger groundwater basins and to understand the geometry of the underlying basement complex. After processing of the gravity data, the thickness of the alluvial aquifer materials was modeled by using all available lithology, density, and geophysical data. The thickest alluvial deposits (greater than 500 feet) are in the northern part of the study area along the south side of the San Jacinto fault zone, in the southern part of the Cahuilla groundwater basin, and in the western part of the Terwilliger groundwater basin. Through most of the area of alluvial materials, the thickness of the alluvium estimated from gravity data is less than 400 feet. Analysis of more than 900 drillers’ logs indicated that in areas having relatively thick alluvium, particularly along the San Jacinto fault zone and in the Terwilliger Valley, the alluvium is predominantly composed of sands and gravels. Fine-textured sediments appeared to be discontinuous rather than forming laterally extensive, low-permeability layers. More than 500 drillers’ logs indicated only bedrock is present, indicating that the fractured bedrock is an important source of groundwater, primarily for domestic use, in the study area. The depths of the holes drilled into the bedrock indicated that fractures potentially supplying water to wells persist in the upper few hundred feet and that the permeable zone of the fractured bedrock extends to depths greater than weathered zones in the upper part of the basement complex. Water-level data were collected from 59 wells during fall 2013. These data indicated that hydraulic head did not vary substantially with well depth and that the measured water levels in bedrock and alluvium were similar. Large offsets in groundwater altitude across the San Jacinto fault zone indicated that the fault zone is a barrier to groundwater flow in the northeastern part of the Anza Valley. On the basis of data from 33 wells, water levels mostly declined between the fall of 2006 and the fall of 2013; the median decline was 5.1 feet during this period, for a median rate of decline of about 0.7 feet/year. Based on data from 40 wells, water-level changes between fall 2004 and fall 2013 were variable in magnitude and trend, but had a median decline of 2.4 feet and a median rate of decline of about 0.3 feet/ year. These differences in apparent rates of groundwater-level change highlight the value of ongoing water-level measurements to distinguish decadal, or longer term, trends in groundwater storage often associated with climatic variability and trends. Fifty-four long-term hydrographs indicated the sensitivity of groundwater levels to climatic conditions; they also showed a general decline in water levels across the study area since 1986 and, in some cases, dating back to the 1950s.

California↗

Multi-resource analysis: A proof of concept study of natural resource tradeoffs in the Piceance Basin, Colorado, using the net resources assessment (NetRA) decision support tool

Executive Summary The U.S. Geological Survey (USGS) is developing a multi-resource analysis (MRA) line of products to inform land-use decision makers. Specifically, MRA products will integrate scientific information, include considerations for natural resource interrelations, and quantify the effects of resource management decisions in biophysical, economic, and societal terms. As part of the establishment of the MRA, the USGS, in collaboration with the University of New Mexico, has developed the Net Resources Assessment (NetRA) decision support tool. As a proof of concept analysis, the NetRA was applied to the Piceance basin in Colorado in a hypothetical example to illustrate how resource managers could use the NetRA to consider tradeoffs of natural resources among alternative development plans and land cover patterns within a geographic region. The NetRA is a policy-relevant approach to assess the availability of multiple natural resources. It is an analytical toolset that may be used to examine the spatiotemporal relations between development of energy and mineral resources and delivery of biological natural resources. The NetRA operates at multiple map scales and contains a set of integrated, compatible submodels with specific data requirements for natural resource stocks, engineering economics, biophysical, and ecological data for ecosystem services stocks, market prices, regulations, and nonmarket values. The NetRA includes an explicit process to consider the interdependence between development and conservation, which is a crucial consideration in land-management and land-use decisions. The NetRA is used to estimate an expected net resource value (NRV). The NRV is the expected, present value, economic benefit from the extraction of a resource (for example, natural gas) minus the total cost of production, which is the aggregation of the development, production, and social costs. Social costs include private costs plus any external costs. There can be external social benefits associated with natural gas production, such as increased demand for locally produced goods and increased employment in the local area through backward and forward linkages of natural gas production. The NRV is used to compare development outcomes (scenarios) from a range of exploration and development plans for cumulative energy production. The Piceance basin application of the NetRA uses the NRV to assess the tradeoff between continuous natural gas extraction and the effects to the local populations of Odocoileus hemionus (mule deer) and aquatic species and to consumptive water uses for an area the size and resolution of a USGS energy resource assessment unit. In the proof of concept simulation, the 2.9-square-mile-area of USGS oil and gas assessment unit 50200263 (Piceance basin continuous gas unit of the Mesaverde Total Petroleum System) was gridded into 588 cells. From this area, seven clusters with potential for development and three that cannot be developed were identified; the three clusters that cannot be developed were identified as wilderness study areas, areas of critical environmental concern, and national forests. On the basis of these criteria, there are 118 cells unsuitable for development in the oil and gas assessment unit: 84 are in national forests, 23 are areas of critical environmental concern, and 11 are wilderness study areas. The remaining cells in the oil and gas assessment unit can be developed on both private and public lands. Two scenarios were considered that are distinguished as plan 1 and plan 2. Plan 1 keeps the amount of land disturbance unchanged and limits the number of development locations to 140 grid cells for the production period, which constrains the amount of the energy resources available for development; the plan requires the usage of the Bureau of Land Management (BLM) unsuitability criteria. Plan 2 also limits the number of development locations to 140 grid cells for the production period but provides a constant volume of energy production by increasing the density of well pads within the cells. The effects of plan 2 to the NRV when there are five wells per pad and five pads per square mile happen mostly in the first 5 years of development, even though the effects on the population of mule deer continue in later years. This outcome is the result of the upfront development and investment costs and the initial effect to the ecosystem services.

Colorado↗

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri and Mississippi Rivers near St. Louis, Missouri, August 3–10, 2020

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near 15 bridges at 10 highway crossings of the Missouri and Mississippi Rivers near Washington, Louisiana, and St. Louis, Missouri, on August 3–10, 2020. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches about 1,640 to 1,970 feet longitudinally and generally extending laterally across the active channel from bank to bank during moderate flood-flow conditions. These surveys provided channel geometry and hydraulic conditions at the time of the surveys and provided characteristics of scour holes that may be useful in developing predictive guidelines or equations for computing potential scour depth. These data also may be useful to the Missouri Department of Transportation as a low to moderate flood-flow assessment of the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every in-channel pier. Scour holes were present at most piers for which bathymetry could be obtained, except those on banks or surrounded by riprap. All the bridge sites in this study were previously surveyed and documented in previous studies, including the two new bridge structures at Louisiana and Washington (structures A8141 and A8504, sites 22 and 32, respectively). Comparisons between bathymetric surfaces from the previous surveys and those of the current (2020) study do not indicate any consistent correlation in channel-bed elevations with streamflow conditions. The comparisons of the 2020 surveys to two previous surveys at the new bridge structure A8141 at Washington (site 22) resulted in net erosion of the channel bed in both comparisons, despite the 2020 streamflow being less than either previous survey. Alternatively, there was a net gain of sediment at new bridge structure A8504 at Louisiana (site 32) between 2014 and 2020, which was the most substantial increase in the surveys detailed in this report; substantially less flow in 2020 than in 2014 or changes to the channel and spur dikes near the bridge may have contributed to the observed sediment gain. Pier size, nose shape, and skew to approach flow had a substantial effect on the size of the scour hole observed at a given pier. Larger and deeper scour holes were present at piers with wide or blunt noses caused by exposed footings, seal courses, or caissons. When a pier was skewed to primary approach flow, the scour hole was generally deeper and larger than at a similar pier without skew; however, the shape of the scour hole near skewed piers in this study generally was longer and deeper on the leeward side, contrary to the general shape of scour holes for skewed piers. However, this phenomenon has been observed historically at these sites, and likely is exacerbated by debris rafts or other turbulence-inducing features near the atypical scour holes. A substantial scour hole was observed near pier 11 of structure A6500 (site 33), which was deeper than in the 2016 survey. The scour holes observed at pier 17 of structure L0561 (site 25) and piers 3 and 4 of structure A1500 (site 34) also were slightly deeper and wider in 2020 than in 2016. At new bridge structures A8141 at Washington (site 22) and A8504 at Louisiana (site 32), the smaller cross-sectional area and configuration of the piers of the new bridges resulted in substantially less scour than with the wider old piers.

Illinois, Missouri↗

Potentiometric surfaces in the Cockfield and Wilcox aquifers of southern and northeastern Arkansas, 2003

This report presents the results of water-level measurements made at wells in the Cockfield Formation and Wilcox Group of southern and northeastern Arkansas during 2003, and the water levels are displayed in potentiometric-surface maps and hydrographs. During March and April 2003, the water level was measured at 55 wells completed in the Cockfield aquifer, 13 wells completed in the Wilcox aquifer of southern Arkansas, and 43 wells completed in the Wilcox aquifer of northeastern Arkansas. The Cockfield Formation generally consists of discontinuous sand units interbedded with silt, clay, and lignite in southeastern Arkansas. Sand beds near the base of the Cockfield Formation constitute most of the Cockfield aquifer. Withdrawals from the Cockfield aquifer in the study area during 2000 totaled about 9 million gallons per day. The potentiometric surface of the Cockfield aquifer constructed from the 2003 water levels shows that regional direction of ground-water flow generally is towards the east and southeast, away from the outcrop, except in areas of intense ground-water withdrawals. Some local ground-water flow in the outcrop area is toward rivers that have eroded into the Cockfield Formation and deposited alluvium in south Bradley and Calhoun Counties (Ouachita River), and in north Dallas County (Saline River). An evaluation of 20 wells with water-level data from 1983 to 2003 shows that water levels in 15 wells have declined at a rate of -0.04 to -0.97 feet per year, and water levels in 5 wells have risen at a rate of 0.07 to 0.32 feet per year. An evaluation of the same 20 wells from 2000 to 2003 shows that water levels have declined in only 8 wells, and water levels have risen in 12 wells. The Wilcox Group is distributed throughout most of southern and eastern Arkansas. There are two study areas in southern and northeastern Arkansas. The Wilcox Group of the southern study area consists of interbedded clay, sandy clay, sand, and lignite. Thin discontinuous sand units constitute the Wilcox aquifer in the southern study area. Withdrawals from the aquifer in the southern study area were about 1 million gallons per day during 2000. The potentiometric surface of the Wilcox aquifer in the southern study area shows that regional ground-water flow generally is south and east, except in Clark County where flow is towards the Ouachita River. The Wilcox Group in the northeastern study area consists of thin interbedded lignitic sand and clays. A sand bed of about 200 feet thick in the middle to lower part of the Wilcox Group constitutes the major producing unit of the Wilcox aquifer in the northeastern study area. Withdrawals from the aquifer in the northeastern study area were about 23 million gallons per day during 2000. The potentiometric surface of the Wilcox aquifer in the northeastern study area shows that ground-water flow generally is south and east, except where ground-water withdrawals may have altered the natural direction of flow near the centers of pumping at Paragould and West Memphis. An evaluation of 27 wells with water-level data from 1983 to 2003 in the northeastern study area shows that water levels in all 27 wells have been declining at a rate of -0.17 to -1.73 feet per year. An evaluation of the same 27 wells from 2000 to 2003 shows that water levels in 18 wells have risen and in 9 wells have declined.

Scientific Investigations Report↗

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri and Mississippi Rivers on the periphery of Missouri, June 13–22, 2022

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near seven bridges at six highway crossings of the Missouri and Mississippi Rivers on the periphery of Missouri from June 13–22, 2022. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches about 1,640 feet longitudinally and generally extending laterally across the active channel from bank to bank during minor flood-flow conditions. These surveys provided channel geometry and hydraulic conditions at the time of the surveys and provided characteristics of scour holes that may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. These data also may be useful to the Missouri Department of Transportation as a minor flood-flow assessment of the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every in-channel pier. Scour holes were present at most piers for which bathymetry could be obtained, except those on banks or surrounded by riprap. Occasionally, scour holes were minor and difficult to discern from nearby dunes and ripples. All bridge sites in this study were surveyed and documented in previous studies. Although partial exposure of substructural support elements was observed at several piers, at most sites the exposure most likely is minimal compared to the overall substructure that remains buried in bed material at these piers. The notable exceptions are piers 12 and 13 at structure L0135 on State Highway 51 at Chester, Illinois, where the bedrock material was fully exposed around the piers. The average difference between the bathymetric surfaces between 2022 and 2018 varied from 0.41 foot higher to 1.86 feet lower. Between 2022 and 2014, the average difference between the bathymetric surfaces varied from 1.02 feet higher to 4.69 feet lower. Only the two sites on the Missouri River and the Caruthersville site were surveyed in 2011; for those sites, the average difference between the bathymetric surfaces varied from 5.83 feet higher to 1.34 feet lower. The most substantial overall net gain of sediment in a reach was between 2011 and 2022 at structure A1700 near Caruthersville, Mo. (site 38). This result was expected because structure A1700 is downstream from the confluences of the Missouri and Ohio Rivers, and therefore subject to the largest streamflows, the largest streamflow fluctuations, and the most substantial sediment flux, as has historically been observed at this site. The presence of riprap blankets, pier size and nose shape, and alignment to flow had a substantial effect on the size of the scour hole observed for a given pier. Piers that were surrounded by riprap blankets had scour holes that were substantially smaller (to nonexistent) compared to piers at which no rock or riprap were present. New riprap blankets were surveyed at pier 3 of structure L0098 at Brownville, Nebraska, and at piers 15–18 of structure A1700 near Caruthersville, Mo., that effectively mitigated the scour holes historically observed at these piers. Narrow piers having round or sharp noses that were aligned with flow often had scour holes that were difficult to discern from nearby bed features, whereas piers having wide or blunt noses resulted in larger, deeper scour holes. Several of the structures had piers that were skewed to primary approach flow. Scour holes near these piers consistently displayed greater depth on the side of the pier with impinging flow and deposition on the leeward side of the pier.

Missouri↗

Hydrogeologic characterization of the Cahuilla Valley and Terwilliger Valley Groundwater Basins, Riverside County, California

The relation between the groundwater and the amount of natural recharge to the Cahuilla Valley and Terwilliger Valley groundwater basins is not well understood. During the 20th century, the reliance on groundwater near Anza, California, used for agricultural, domestic, and municipal reasons has increased, and there is the potential for changes in groundwater availability related to climate change. Several types of existing data were evaluated, and new data were collected for this study, with the goal of characterizing the region’s hydrogeology. The study’s scope included constructing a geologic framework model to show where the groundwater-bearing units are present and their relation to each other, estimating the major components of the groundwater budget, and understanding local short-term and regional long-term groundwater flow and how that has changed since the early 1900s. Two electrical resistivity tomography surveys were done in the Durasno Valley about 2,150 feet apart to identify the thickness of the alluvium, its horizontal extent, and the depth-to-basement along two profiles perpendicular to Cahuilla Creek. The subsurface sediments were mostly horizontally layered and the transitional boundary between the alluvium and basement was thinner and shallower along the upgradient profile where the depth-to-basement was about 70 feet below land surface; the depth-to-basement at the downgradient profile was more than about 140 feet below land surface. The results from the surveys were used to place four monitoring wells at two sites along the survey profiles. Artesian flow from the deepest well at the downgradient site indicated that the decomposed and competent basement likely contributed some groundwater to the overlying alluvium, laterally, from below, or both. A digital three-dimensional geologic framework model was constructed using EarthVision software to represent the subsurface geometry of the alluvium, decomposed basement, and competent basement. Maps and cross sections of the modeled thicknesses of the alluvium and decomposed basement, and the modeled elevation of the top of the competent basement, were made to show the subsurface geometry of vertical faults, selected wells, and the groundwater-bearing units. Because natural recharge is related to the variable cycles of precipitation, estimates are difficult to quantify. Recharge and runoff have extreme interannual variability in the study area; recharge and runoff can be sporadic, and a substantive amount may not occur in some years. Estimates of recharge from a previous study and the regional-scale Basin Characterization Model for California for four different periods ranged from 3,800 acre-feet/year for 1897–1947 to 5,900 acre-feet/year for 1971–2000. Potential recharge from the disposal of domestic septic systems may have been as much as 500 acre-feet in 2020. It was estimated that between about 400 and 2,400 acre-feet/year of groundwater is lost through evapotranspiration by vegetation and evaporation from open water bodies, but the main source of discharge is through pumpage, mainly used for agriculture from the alluvium in the Cahuilla Valley and Terwilliger Valley groundwater basins. The estimated total pumpage for 1991–2021 ranged from about 1,140 acre-feet in 2019 to about 3,450 acre-feet in 1994. When summed, the cumulative amount of estimated pumpage between 1991 and 2021 was about 81,400 acre-feet. The general direction of groundwater flow is from the northeast along the San Jacinto fault zone at the headwaters of Cahuilla and Hamilton Creeks, to the surface-water outlets at the west and southeast parts of the study area. Groundwater-level data from the 1950s and earlier indicate that there was a natural groundwater divide between the Cahuilla Valley and Terwilliger Valley groundwater basins, but the changing magnitude and extent of the groundwater depressions caused by pumping since about 1950 indicate that the location of the natural groundwater boundary between the Cahuilla Valley and Terwilliger Valley groundwater basins has migrated over time. Flow from the upper to the lower parts of the Cahuilla Valley groundwater basin roughly follows the course of Cahuilla Creek through the narrow Durasno Valley where an estimated volume of flow in April 2019 was about 10–150 acre-feet/year. Short-term trends in groundwater levels, particularly in wells where groundwater is shallow and in the basement unit, show how some areas respond quickly to recharge and discharge. Wells located further to the east within the Cahuilla Valley groundwater basin in the alluvium show much less of a response to recharge events; areas of sustained pumpage from the alluvium, primarily for agriculture, show long-term declines in groundwater levels and generally do not show the effects of storm events or recent runoff. Groundwater levels in wells that are farthest from where most of the recharge occurs and where pumping has been the greatest, had some of the largest long-term groundwater-level declines at a rate of about 0.8 foot/year between 1971 and 2021.

California↗

Origins and delineation of saltwater intrusion in the Biscayne aquifer and changes in the distribution of saltwater in Miami-Dade County, Florida

Intrusion of saltwater into parts of the shallow karst Biscayne aquifer is a major concern for the 2.5 million residents of Miami-Dade County that rely on this aquifer as their primary drinking water supply. Saltwater intrusion of this aquifer began when the Everglades were drained to provide dry land for urban development and agriculture. The reduction in water levels caused by this drainage, combined with periodic droughts, allowed saltwater to flow inland along the base of the aquifer and to seep directly into the aquifer from the canals. The approximate inland extent of saltwater was last mapped in 1995. An examination of the inland extent of saltwater and the sources of saltwater in the aquifer was completed during 2008–2011 by using (1) all available salinity information, (2) time-series electromagnetic induction log datasets from 35 wells, (3) time-domain electromagnetic soundings collected at 79 locations, (4) a helicopter electromagnetic survey done during 2001 that was processed, calibrated, and published during the study, (5) cores and geophysical logs collected from 8 sites for stratigraphic analysis, (6) 8 new water-quality monitoring wells, and (7) analyses of 69 geochemical samples. The results of the study indicate that as of 2011 approximately 1,200 square kilometers (km 2 ) of the mainland part of the Biscayne aquifer were intruded by saltwater. The saltwater front was mapped farther inland than it was in 1995 in eight areas totaling about 24.1 km 2 . In many of these areas, analyses indicated that saltwater had encroached along the base of the aquifer. The saltwater front was mapped closer to the coast than it was in 1995 in four areas totaling approximately 6.2 km 2 . The changes in the mapped extent of saltwater resulted from improved spatial information, actual movement of the saltwater front, or a combination of both. Salinity monitoring in some of the canals in Miami-Dade County between 1988 and 2010 indicated influxes of saltwater, with maximum salinities ranging from 1.4 to 32 practical salinity units (PSU) upstream of the salinity control structures. Time-series electromagnetic induction log data from monitoring wells G–3601, G–3608, and G–3701, located adjacent to the Biscayne, Snapper Creek, and Black Creek Canals, respectively, and upstream of the salinity control structures, indicated shallow influxes of conductive water in the aquifer that likely resulted from leakage of brackish water or saltwater from these canals. The determination that saltwater influxes were recent is supported by the similarity in the oxygen and hydrogen stable isotope composition in samples from the Snapper Creek Canal, 1.6 kilometers (km) inland of a salinity control structure, and in samples from well G–3608, which is adjacent to the canal, as well as by the relative ages of the water sampled from well G–3608 and other wells open to the aquifer below the saltwater interface. Historical and recent salinity information from the Card Sound Road Canal, monitoring well FKS8 located adjacent to the canal, and the 2001 helicopter electromagnetic survey indicated that saltwater may occasionally leak from this canal as far inland as 15 km. This leakage may be prevented or reduced by a salinity control structure that was installed in May 2010. Saltwater also may have leaked from the Princeton Canal. Results of geochemical sampling and analysis indicate a close correspondence between droughts and saltwater intrusion. Tritium/helium-3 apparent (piston-flow) ages determined from samples of saltwater with chloride concentrations of about 1,000 milligrams per liter (mg/L) or greater generally corresponded to a period during which droughts were frequent. Comparison of average daily air temperatures in Miami, Florida, with estimates of recharge temperatures determined from the dissolved gas composition in water samples indicated that saltwater likely entered the aquifer in April or early May when water levels are typically at their lowest during the year. Conversely, most of the samples of freshwater with chloride concentrations less than about 1,000 mg/L indicate recharge temperatures corresponding to air temperatures in mid to late May when rainfall and water levels in the aquifer increase, and the piston-flow ages of these samples correspond to wet years. The piston-flow ages of freshwater samples generally were younger than ages of samples of saltwater. Saltwater samples that were depleted in boron, magnesium, potassium, sodium, and sulfate, and enriched in calcium relative to the concentrations theoretically produced by freshwater/seawater mixing, generally were found to be associated with areas where saltwater had recently intruded. The calcium to (bicarbonate + sulfate) molar ratios (Ca/(HCO 3 +SO 4 )) of these samples generally were greater than 1. Saltwater samples from some of the monitoring wells, however, indicated little or no enrichment or depletion of these ions relative to the theoretical freshwater/seawater mixing line, and the Ca/(HCO 3 +SO 4 ) molar ratios of these samples generally were less than 1. Results indicated that aquifer materials are approaching equilibrium with seawater at these well locations.

Florida↗

Simulation and particle-tracking analysis of ground-water flow near the Savannah River site, Georgia and South Carolina, 2002, and for selected ground-water management scenarios, 2002 and 2020

Ground-water flow under 2002 hydrologic conditions was evaluated in an eight-county area in Georgia and South Carolina near the Savannah River Site (SRS), by updating boundary conditions and pumping rates in an existing U.S. Geological Survey (USGS) ground-water model. The original ground-water model, developed to simulate hydrologic conditions during 1987-92, used the quasi-three-dimensional approach by dividing the Floridan, Dublin, and Midville aquifer systems into seven aquifers. The hydrogeologic system was modeled using six active layers (A2-A7) that were separated by confining units with an overlying source-sink layer to simulate the unconfined Upper Three Runs aquifer (layer A1). Potentiometric- surface maps depicting September 2002 for major aquifers were used to update, evaluate, and modify boundary conditions used by the earlier ground-water flow model. The model was updated using the USGS finite-difference code MODFLOW-2000 for mean-annual conditions during 1987-92 and 2002. The specified heads in the source-sink layer A1 were lowered to reflect observed water-level declines during the 1998-2002 drought. These declines resulted in a decrease of 12.1 million gallons per day (Mgal/d) in simulated recharge or vertical inflow to the uppermost confined aquifer (Gordon, layer A2). Although ground-water pumpage in the study area has increased by 32 Mgal/d since 1995, most of this increase (17.5 Mgal/d) was from the unconfined Upper Three Runs aquifer (source-sink layer A1) with the remaining 14.5 Mgal/d assigned to the active layers within the model (A2-A7). The simulated water budget for 2002 shows a decrease from the 1987-92 model from 1,040 Mgal/d to 1,035 Mgal/d. The decreased ground-water inflows and increased ground-water withdrawal rates reduced the simulated ground-water outflow to river cells in the active layers of the model by 43 Mgal/d. The calibration statistics for all layers of the 2002 simulation resulted in a decrease in the root mean square (RMS) of the residuals from 10.6 to 8.0 feet (ft). The residuals indicate 83.3 percent of the values for the 2002 simulation met the calibration error criteria established in the original model, whereas 88.8 percent was within the specified range for the 1987-92 simulation. Simulated ground-water outflow to the Savannah River and its tributaries during water year 2002 was 560 cubic feet per second (ft3/s), or 86 percent of the observed gain in mean-annual streamflow between streamflow gaging stations at the Millhaven, Ga., and Augusta, Ga. At Upper Three Runs Creek, simulated ground-water discharge during 2002 was 110 ft3/s, or 83 percent of the observed streamflow at two streamflow gaging stations near the SRS. These results indicate that the constructed model calibrated to 1987-92 conditions and modified for 2002 dry conditions is still representative of the hydrologic system. The USGS particle-tracking code MODPATH was used to generate advective water-particle pathlines and their associated time-of-travel based on MODFLOW simulations for 1987-92, 2002, and each of four hypothetical ground-water management scenarios. The four hypothetical ground-water management scenarios represent hydrologic conditions for (1) reported pumping for 2002 and boundary conditions for an average year; (2) reported pumping for 2002 with SRS pumping discontinued and boundary conditions for an average year; (3) projected 2020 pumping and boundary conditions for an average year; and (4) projected 2020 pumping and boundary conditions for a dry year. The MODPATH code was used in forward-tracking mode to evaluate flowpaths from areas on the SRS and in backtracking mode to evaluate further areas of previously documented trans-river flow on the Georgia side of the Savannah River. Trans-river flow is a condition in which the local head gradients might allow migration of contaminants from the SRS into the underlying aquifers and beneath the Savannah River into Georgia. More...

Georgia, South Carolina↗

Potentiometric surface maps and groundwater-level hydrographs for confined aquifers of the New Jersey Coastal Plain, 2018

The U.S. Geological Survey, in cooperation with the New Jersey Department of Environmental Protection (NJDEP), prepared potentiometric surface maps for 10 confined aquifers of the New Jersey Coastal Plain physiographic province based on water-level measurements collected during late 2018 and early 2019 from 951 wells in New Jersey and parts of Pennsylvania and Delaware. Maps were prepared for the confined Cohansey aquifer, Rio Grande water-bearing zone, Atlantic City 800-foot sand, Piney Point aquifer, Vincentown aquifer, Wenonah-Mount Laurel aquifer, Englishtown aquifer system, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system. Potentiometric surface maps indicate regional cones of depression in the following aquifers and the counties in which they are centered: Atlantic City 800-foot sand in Atlantic County, the Piney Point aquifer in Cumberland County, the Wenonah-Mount Laurel aquifer and Englishtown aquifer system in Monmouth and Ocean Counties, the Wenonah-Mount Laurel aquifer in Camden and Gloucester County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Ocean County, and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy aquifer system in Camden County. Cones of depression with smaller areal extents were in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Atlantic City 800-foot sand centered in Cape May County, the Piney Point aquifer centered in Ocean County, the Wenonah-Mount Laurel aquifer in Salem and Burlington Counties, the Englishtown aquifer system in Camden County, the upper aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth County, and the middle aquifer of the Potomac-Raritan-Magothy aquifer system in Monmouth, Ocean, and Salem Counties. No cone of depression was interpreted in the Vincentown aquifer. Long-term hydrographs are presented for 75 wells spanning each of the 10 confined aquifers, and contain a mix of discrete water-level measurements and daily mean water levels based on continuously recorded 15-minute data. Changes of water levels during 2014–19, as indicated by the hydrographs, were compared with those of previous periods to assess any departures from historical data. During 2014–19, water levels were stable and fluctuated within similar ranges as previous periods in the following aquifers and locations: all wells in the confined Cohansey aquifer, the Rio Grande water-bearing zone, the Vincentown aquifer, the Englishtown aquifer system, the Piney Point aquifer wells in Burlington and Ocean Counties, six of eight wells in the Wenonah-Mount Laurel aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system except those within NJDEP Critical Area II. Increasing water levels in 2014–19, ongoing since historical periods, were indicated in the following aquifers and locations: Atlantic County wells in the Piney Point aquifer, all wells in the upper and lower aquifers of the Potomac-Raritan-Magothy aquifer system outside NJDEP Critical Areas, and all wells in the middle aquifer of the Potomac-Raritan-Magothy aquifer system within NJDEP Critical Area II. Water levels in the Atlantic City 800-foot sand also increased during 2014–19 in wells in Atlantic County and northern Cape May County closer to the center of the cone of depression in that aquifer, which is a response unique to this period and absent from previous periods. During 2013–19, continued decreasing water levels, ongoing since previous periods, were indicated by hydrographs of Atlantic City 800-foot sand wells in southern Cape May County, Piney Point aquifer wells in Cumberland County where the regional cone of depression is located, and two wells in the Wenonah-Mount Laurel aquifer—070478, which in 2014–19 departed from previous periods, and 330020, which continued a gradual decrease throughout its period of record.

New Jersey↗

Potential interaction of groundwater and surface water including autonomous underwater vehicle reconnaissance at Nolin River Lake, Kentucky, 2016

The U.S. Geological Survey collaborated with the U.S. Army Corps of Engineers, Louisville District, on a synoptic study of water quality at Nolin River Lake during August 2016. The purpose of the study was to develop a better understanding of the potential for interaction between groundwater and surface water at Nolin River Lake, Kentucky. Groundwater can have properties that are measurably different from those in adjacent surface water, and inflows and outflows can be an important component of water quality and quantity. An improved understanding of potential interaction of groundwater and surface water at Nolin River Lake may be used to refine lake-management strategies. This study (1) compiled and interpreted existing information to characterize the hydrogeological setting and implications for potential interaction of groundwater and surface water in the Nolin River Lake watershed; (2) collected transects of onsite water-quality parameters using an autonomous underwater vehicle (AUV) in areas with potential for interaction of groundwater and surface water, including five sites on Nolin River Lake and one site on the Nolin River; and (3) collected discrete water-quality and phytoplankton community data at the same six sites. A review of existing hydrogeologic information did not indicate the presence of karst features adjacent to or beneath Nolin River Lake that would facilitate groundwater interaction with the reservoir. Observations leading to this conclusion include (1) limestone that is adjacent to the shoreline and perhaps beneath the lake, is overlain with siliciclastic rocks and fine-grained sediment that inhibits infiltration and development of karst features that encourage rapid groundwater flow; (2) the geologic deposits surrounding the reservoir are described as having limited or no potential for development of karst features, some exceptions may exist in tributary valleys; (3) very few karst features were mapped within 1 mile of the reservoir or in the area currently occupied by the reservoir; and (4) faults that intersect the reservoir but may not possess hydraulic properties that cause the faults to be conduits for groundwater flow. Groundwater interaction with reservoir tributaries is likely more common in areas of the watershed upstream from Nolin River Lake where karst hydrogeology is prevalent. Results of water-quality surveys using an AUV from August 15 to 19, 2016, did not identify areas of anomalous values that might indicate groundwater inflows through preferential flow zones. Spatial distributions of water-quality parameters were generally uniform within each constant-depth layer. The constant-depth layers were selected to be above, within, and below the thermocline and ranged from the water surface to 25 feet. Surveys near the bottom of the reservoir that might have been more sensitive to groundwater inflows were not done because presurvey data were not available to indicate locations of obstacles that could ensnare the AUV. Water-quality data collected with the AUV did identify water-quality anomalies where stream tributaries were discharging to the reservoir. The discrete water-quality samples indicated uniformity among the five reservoir sites. The riverine site that is immediately upstream from Nolin River Lake, however, had some unique water-quality characteristics relative to sites on the reservoir. The highest concentrations of nitrate plus nitrite as nitrogen (0.145 milligrams per liter [mg/L]), total phosphorous (0.07 mg/L), chlorophyll a (36.1 micrograms per liter), and pheophytin a (10.2 micrograms per liter) were measured at the Nolin River Lake riverine site (site 2NRR20034). The concentrations of nutrients and chlorophyll a at the riverine site did exceed the 25th percentile of median concentrations measured by the U.S. Environmental Protection Agency (EPA) at other lakes and reservoirs in EPA level IV ecoregion 71a. Concentrations of most nutrients and chlorophyll a at the five reservoir sites also exceeded the 25th percentile of median concentrations in EPA level IV ecoregion 72h. The exception was the concentrations of total phosphorus as phosphorus at the reservoir sites that were at or below the 25th percentile of median concentrations measured by EPA (0.03 mg/L). Concentrations of orthophosphate as phosphorus were less than the method detection limit of 0.004 mg/L at all sites. The phytoplankton community in Nolin River Lake was almost exclusively (greater than 90 percent of total phytoplankton abundance) cyanobacteria, also known as blue-green algae. A species of Cylindrospermopsis dominated the cyanobacterial community at the five reservoir sites, while Chroococcus microscopicus was most abundant at the riverine site. Cyanobacterial cell densities ranged from 10,000 to 198,067,460 cells per liter in five areas in the reservoir and from 4,800 to 73,751,253 cells per liter at the riverine site. Multiple potential sources of water to Nolin River Lake include direct precipitation, overland flow, interflow, groundwater, and surface water. Understanding the exact contribution of each of these components to the water budget at Nolin River Lake may help the U.S. Army Corps of Engineers manage the water quality, water quantity, and biological communities in the reservoir. Additional hydrogeologic and water-quality data that builds on the results of this study may refine the inferences of this study; for example, deeper AUV surveys that target the largest fault zones might further the understanding of the potential for groundwater flow through those features. A complete understanding of the reservoir hydrology, however, may require the use of scientific methods intended for water bodies as large as Nolin River Lake, such as aerial infrared photography and imagery; water mass, chemical, and isotopic balance studies; geophysical measurements; and numerical simulations.

Kentucky↗

Bathymetric and velocimetric surveys at highway bridges crossing the Missouri River between Kansas City and St. Louis, Missouri, May 19–26, 2021

Bathymetric and velocimetric data were collected by the U.S. Geological Survey, in cooperation with the Missouri Department of Transportation, near nine bridges at eight highway crossings of the Missouri River between Kansas City and St. Louis, Missouri, from May 19 to 26, 2021. A multibeam echosounder mapping system was used to obtain channel-bed elevations for river reaches about 1,640 to 1,840 feet (ft) longitudinally and generally extending laterally across the active channel from bank to bank during low to moderate flood-flow conditions. These surveys provided channel geometry and hydraulic conditions at the time of the surveys and provided characteristics of scour holes that may be useful in developing or verifying predictive guidelines or equations for computing potential scour depth. These data also may be useful to the Missouri Department of Transportation as a low to moderate flood-flow assessment of the bridges for stability and integrity issues with respect to bridge scour during floods. Bathymetric data were collected around every in-channel pier. Scour holes were present at most piers for which bathymetry could be obtained, except those on banks or surrounded by riprap. Occasionally, scour holes were minor and difficult to discern from nearby dunes and ripples. All the bridge sites in this study were previously surveyed and documented in previous studies. Comparisons between bathymetric surfaces from the previous surveys and those of the current (2021) study do not indicate any consistent correlation between channel-bed elevations and streamflow conditions. The average difference between the bathymetric surfaces varied from 1.59 ft higher to 0.95 ft lower in 2021 than 2017, which corresponds to a gain of 100,200 cubic yards and a loss of 55,800 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 2.74 ft higher to 3.05 ft lower in 2021 than 2013, which corresponds to a gain of 111,500 cubic yards and a loss of 169,200 cubic yards, respectively. The average difference between the bathymetric surfaces varied from 4.52 ft higher to 1.38 ft lower in 2021 than 2011, which corresponds to a gain of 221,100 cubic yards and a loss of 90,300 cubic yards, respectively. The most substantial overall net gain was 221,100 cubic yards between 2011 and 2021 at structures L0550 and A4497 at Jefferson City (site 20). The large net gain likely results from a combination of the mitigation of the scour holes near pier 4 of both bridges and the substantially lower flow in 2021 than in 2011. Alternatively, the most substantial overall net loss was 169,200 cubic yards between 2013 and 2021 at structure A6288 at Hermann (site 21), despite comparable streamflows. Pier size, nose shape, and skew to approach flow had a substantial effect on the size of the scour hole observed at a given pier. Larger and deeper scour holes were present at piers with wide or blunt noses caused by exposed footings or caissons. When a pier was skewed to primary approach flow, the scour hole was generally deeper and larger than at a similar pier without skew; furthermore, the shape of the scour hole near skewed piers in this study generally was longer and deeper on the side with impinging flow. At structure A6288 at Hermann (site 21), the scour hole near pier 5 was difficult to discern from nearby dunes and ripples, whereas the upstream edge of the footing was visible at pier 4, which likely contributes to the larger scour hole near that pier; the top of the footing may blunt the horseshoe vortex at pier 5, but the exposed front of the footing may exacerbate the vortex at pier 4.

Missouri↗

Lithology, hydrologic characteristics, and water quality of the Arkansas River Valley alluvial aquifer in the vicinity of Van Buren, Arkansas

A study to assess the potential of the Arkansas River Valley alluvial aquifer in the vicinity of Van Buren, Arkansas, as a viable source of public-supply water was conducted by the U.S. Geological Survey in cooperation with the Little Rock, District, U.S. Army Corps of Engineers. An important study component was to identify possible changes in hydrologic conditions following installation of James W. Trimble Lock and Dam 13 (December 1969) on the Arkansas River near the study area. Data were gathered for the study in regard to the lithology, hydrologic characteristics, and water quality of the aquifer. Lithologic information was obtained from drillers’ logs of wells drilled from 1957 through 1959. Water-quality samples were collected from 10 irrigation wells and analyzed for inorganic constituents and pesticides. To evaluate the potential viability of the alluvial aquifer in the Van Buren area, these data were compared to similar stratigraphic, lithologic, and groundwater-quality data from the Arkansas River Valley alluvial aquifer at Dardanelle, Ark., where the aquifer provides a proven, productive, sole-source of public-supply water. Drillers’ logs for 59 wells in the Van Buren study area revealed well depths ranging from 25 to 52 feet (ft), with a mean depth of 42 ft. The thickness of the lower sand/gravel interval serving as the water-producing zone ranged from 5 to 47 ft, with a mean thickness of 29 ft. The presence of gravel was noted in only 4 of 59 well logs available for review from the study area. Percent sand was calculated from well logs in the study area, and these sand percentages were overlain onto an orthophotograph map to examine the areal distribution of sand percentage in relation to geomorphologic features of the flood plain in the study area. The logs denoting the greatest percent sand tend to occur in areas near to the river and on the concave (point bar) side of abandoned channels, while the lower percent sand tends to occur on the convex (channel fill and backswamp deposits) side of the abandoned channels. Comparison of hydrographs from water levels collected between 1957 and 1972 to cumulative departure from mean monthly and mean annual precipitation showed overall good fit and explained the long-term decreasing water levels from the earliest period of record through October 1967, followed by a sharp rise in water levels concurrent with rises in cumulative departure from mean monthly and mean annual precipitation. Hydrographs for four wells ranging from 0.8 to 4.5 miles upstream from the dam and potentially affected by rising river stage were compared to graphs of river stage and cumulative departure from mean monthly precipitation. Water levels for these wells showed minimal discernible effect by rising river stage following dam completion. Periods of increased precipitation compared closely to increases in water level for all hydrographs, regardless of river stage, and periods of no precipitation resulted in declining water levels, although river stage continued to slowly rise during these same periods. The Arkansas River has greater salinity than local groundwater, providing a quantitative tracer for any groundwater recharge originating from the river. Comparison of predam and postdam groundwater-chloride concentrations showed no increase in chloride concentrations after dam installation, which is consistent with hydrologic data. These data suggest that the dominant source of groundwater recharge in the Arkansas River Valley alluvial aquifer is infiltration of precipitation through proximal, coarse channel deposits, with minimal influx of river water. Groundwater-quality data collected from 10 wells in the study area indicated a calcium-bicarbonate water type. No primary drinking-water standards were exceeded for any constituents, and iron and manganese were the only constituents exceeding secondary drinking-water regulations. Six of the 10 well-water samples were analyzed for the presence of pesticides, as row-crop agriculture is the dominant land use in the study area. Six herbicide compounds and one herbicide metabolite were detected at concentrations substantially below those of the Federal primary drinking-water standards and health advisories. The hydrologic and geochemical data gathered for this study provide a qualitative assessment of the potential of the Arkansas River Valley alluvial aquifer as a source of public water supply in the Van Buren area. Results indicate minimal influx of water from the Arkansas River, and recharge to the aquifer appears to be dominantly by infiltration of precipitation through overlying alluvium. If vertical wells are used as a source of public water supply, then several wells will have to be used in combination at relatively low pumping rates and placed in areas with a greater percent sand. Use of a horizontal well configuration near the river to increase production may depend on infiltration of river water to supplement water removed from storage, especially where areas of lower permeability sediments might be encountered within the surrounding alluvium. If a poor hydraulic connection exists between the river and the alluvium, as indicated by this study, then production will depend on ample precipitation and recharge throughout the year and groundwater storage sufficient to prevent declining water levels where pumping rates exceed recharge.

Arkansas↗

Stratigraphy and tectonic history of the Tucson Basin, Pima County, Arizona, based on the Exxon state (32)-1 well

The Tucson Basin is a relatively large late Cenozoic extensional basin developed in the upper plate of the Catalina detachment fault in the southern Basin and Range Province, southeastern Arizona. In 1972, Exxon Company, U.S.A., drilled an exploration well (Exxon State (32)-1) near the center of the Tucson Basin that penetrated 3,658 m (12,001 ft) of sedimentary and volcanic rocks above granitoid basement. Detailed study of cuttings and geophysical logs of the Exxon State well has led to revision of the previously reported subsurface stratigraphy for the basin and provided new insight into its depositional and tectonic history. There is evidence that detachment faulting and uplift of the adjacent Catalina core complex on the north have affected the subsurface geometry of the basin. The gravity anomaly map of the Tucson Basin indicates that the locations of subbasins along the north-trending axis of the main basin coincide with the intersection of this axis with west-southwest projections of synforms in the adjacent core complex. In other words, the subbasins overlie synforms and the ridges between subbasins overlie antiforms. The Exxon State well was drilled near the center of one of the subbasins. The Exxon well was drilled to a total depth of 3,827 m (12,556 ft), and penetrated the following stratigraphic section: Pleistocene(?) to middle(?) Miocene upper basin-fill sedimentary rocks (0-908 m [0-2,980 ft]) lower basin-fill sedimentary rocks (908-1,880 m [2,980-6,170 ft]) lower Miocene and upper Oligocene Pantano Formation (1,880-2,516 m [6,170-8,256 ft]) upper Oligocene to Paleocene(?) volcanic and sedimentary rocks (2,516-3,056 m [8,256-10,026 ft]) Lower Cretaceous to Upper Jurassic Bisbee Group (3,056-3,658 m [10,026-12,001 ft]) pre-Late Jurassic granitoid plutonic rock (3,658-3,827 m [12,001- 12,556 ft]). Stratigraphy and Tectonic History of the Tucson Basin, Pima County, Arizona, Based on the Exxon State (32)-1 Well The 1,880 m (6,170 ft) of basin-fill sedimentary rocks consist of alluvial-fan, alluvial-plain, and playa facies. The uppermost unit, a 341-m-thick (1,120-ft) lower Pleistocene and upper Pliocene alluvial-fan deposit (named the Cienega Creek fan in this study), is an important aquifer in the Tucson basin. The facies change at the base of the alluvial fan may prove to be recognizable in well data throughout much of the basin. The well data show that a sharp boundary at 908 m (2,980 ft) separates relatively unconsolidated and undeformed upper basin fill from denser, significantly faulted lower basin fill, indicating that there were two stages of basin filling in the Tucson basin as in other basins of the region. The two stages apparently occurred during times of differing tectonic style in the region. In the Tucson area the Pantano Formation, which contains an andesite flow dated at about 25 Ma, fills a syntectonic basin in the hanging wall of the Catalina detachment fault, reflecting middle Tertiary extension on the fault. The formation in the well is 636 m thick (2,086 ft) and consists of alluvial-fan, playa, and lacustrine sedimentary facies, a lava flow, and rock- avalanche deposits. Analysis of the geophysical logs indicates that a K-Ar date of 23.4 Ma reported previously for the Pantano interval of the well was obtained on selected cuttings collected from a rock-avalanche deposit near the base of the unit and, thus, does not date the Pantano Formation. The middle Tertiary volcanic and sedimentary rocks have an aggregate thickness of 540 m (1,770 ft). We obtained a new 40Ar/ 39Ar age of 26.91+0.18 Ma on biotite sampled at a depth of 2,584-2,609 m (8,478-8,560 ft) from a 169-m-thick (554-ft) silicic tuff in this interval. The volcanic rocks probably correlate with other middle Tertiary volcanic rocks of the area, and the sedimentary rocks may correlate with the Cloudburst and Mineta Formations exposed on the flanks of the San Pedro Basin to the northeast. The Bisbee Group in the Exxon well is 602 m (1,975 f

Scientific Investigations Report↗