Search

Groundwater-flow models require the spatial distribution of the hydraulic conductivity parameter. One approach to defining this spatial distribution in groundwater-flow model grids is to map the electrical resistivity distribution by airborne electromagnetic (AEM) survey and establish a petrophysical relation between mean resistivity calculated as a nonlinear function of the resistivity layering and thicknesses of the layers and aquifer transmissivity compiled from historical aquifer tests completed within the AEM survey area. The petrophysical relation is used to transform AEM resistivity to transmissivity and to hydraulic conductivity over areas where the saturated thickness of the aquifer is known. The US Geological Survey applied this approach to a gain better understanding of the aquifer properties of the Mississippi River Valley alluvial aquifer. Alluvial-aquifer transmissivity data, compiled from 160 historical aquifer tests in the Mississippi Alluvial Plain (MAP), were correlated to mean resistivity calculated from 16,816 line-kilometers (km) of inverted resistivity soundings produced from a frequency-domain AEM survey of 95,000 km2 of the MAP. Correlated data were used to define petrophysical relations between transmissivity and mean resistivity by omitting from the correlations the aquifer-test and AEM sounding data that were separated by distances greater than 1 km and manually calibrating the relation coefficients to slug-test data. The petrophysical relation yielding the minimum residual error between simulated and slug-test data was applied to 2,364 line-km of AEM soundings in the 1,000-km2 Shellmound (Mississippi) study area to calculate hydraulic property distributions of the alluvial aquifer for use in future groundwater-flow models.

Within the State of Arkansas, there is an increasing number of aging dams and levees that have little to no documentation concerning their construction or composition. Surface geophysical surveys offer a non-intrusive method for investigating these structures to describe their lithologic makeup, evaluate the materials constructed upon, and identify potential flow paths through them. Techniques such as electrical resistivity tomography, seismic refraction, and electromagnetic induction have been used to image dams and levees. They require additional information from geologic outcrops, geotechnical borings, or drill cores to make informed geologic interpretations of the geophysical models. These geologic models then allow the owners of these structures to make more informed decisions about their operation and maintenance. Between 2011 and 2018, the U.S. Geological Survey conducted geophysical and geotechnical investigations of three earthen structures in Arkansas. Electrical and electromagnetic geophysical data were used to develop lithologic models of these structures and characterize the underlying geology. Self-potential surveys were utilized to detect the movement of water through these structures and identify any possible seepage pathways. Geotechnical methods such as electric and hydraulic direct-push well logs and cores acted as a control on the geophysical interpretations and a confirmation of anomalies. This integrated approach detected the lack of an impermeable core within a levee, imaged a change in lithology of the bedrock forming the seal beneath a gravity dam, and identified a potential seepage feature within the core of an earthen dam. These results further support that this method of extending known lithologic features via surface and borehole geophysics is a useful approach for characterizing earthen water-control structures.

The Mississippi Alluvial Plain hosts one of the most prolific shallow aquifer systems in the United States but is experiencing chronic groundwater decline. The Reelfoot rift and New Madrid seismic zone underlie the region and represent an important and poorly understood seismic hazard. Despite its societal and economic importance, the shallow subsurface architecture has not been mapped with the spatial resolution needed for effective management. Here, we present airborne electromagnetic, magnetic, and radiometric observations, measured over more than 43,000 flight-line-kilometers, which collectively provide a system-scale snapshot of the entire region. We develop detailed maps of aquifer connectivity and shallow geologic structure, infer relationships between structure and groundwater age, and identify previously unseen paleochannels and shallow fault structures. This dataset demonstrates how regional-scale airborne geophysics can close a scale gap in Earth observation by providing observational data at suitable scales and resolutions to improve our understanding of subsurface structures. Integration of large-scale airborne geophysical surveys can enable the system-scale assessment of aquifer characteristics and improved mapping of shallow subsurface structures, according to a comprehensive investigation of the Mississippi Alluvial Plain

The Mississippi Delta, a portion of the Mississippi Alluvial Plain (MAP) located in northwest Mississippi (USA), is an area dense with industrial-level agriculture sustained by groundwater-dependent irrigation supplied by the Mississippi River Valley alluvial aquifer. Observed declines in groundwater-level elevations and streamflow, contemporaneous with increases in irrigation, have raised concerns about future groundwater availability and the effects of groundwater withdrawals on streamflow. To quantify the impacts of groundwater withdrawals on streamflow and increase understanding of groundwater and surface-water interaction in the MAP, hydrograph-separation techniques were used to estimate baseflow and identify statistical streamflow trends. The analysis was conducted using the US Geological Survey Groundwater Toolbox open-source software and daily hydrologic data provided by a spatially distributed network of paired groundwater wells and streamgage sites. This study found that statistically significant reductions in stream baseflow occurred in areas with substantial groundwater-level declines. The use of hydrograph-separation and trend analyses to quantify the impacts of groundwater withdrawals and the use of streamflow as a proxy for changes in groundwater availability may be applicable in other altered environments. Characterizing and defining hydrologic relations between groundwater and surface water will help scientists and water-resource managers refine a regional groundwater-flow model that includes the Mississippi Delta, which will be used to aid water-resource managers in future decisions concerning the alluvial aquifer.

Critical groundwater resources and hidden seismic hazards underly much of the Mississippi Alluvial Plain. Spanning nearly 100,000 square kilometers across seven states, this region hosts one of the most prolific shallow aquifer systems in the United States that supports a $12 billion agricultural economy amidst chronic groundwater decline. Further, underlying fault structures of the Reelfoot Rift and New Madrid Seismic Zone represent an important and poorly understood hazard with a complex pattern of historical impacts. Despite its societal and economic importance, mapping of shallow subsurface architecture with spatial resolution needed for effective management is insufficient. Here, we report the results of 40,000 flight-line-kilometers of electromagnetic, magnetic, and radiometric data collectively providing a system-scale snapshot of an entire aquifer system, the first such effort in the United States. This survey enables new understanding of the regional hydrogeology while also revealing previously unseen large vertical displacements (exceeding 50 m) in the uppermost Tertiary units within the New Madrid Seismic Zone.

The hydrogeology below large surface water features such as rivers and estuaries is universally under-informed at the long reach to basin scales (tens of km+). This challenge inhibits the accurate modeling of fresh/saline groundwater interfaces and groundwater/surface water exchange patterns at management-relevant spatial extents. Here we introduce a towed, floating transient electromagnetic (TEM) system (i.e. FloaTEM) for rapid (up to 15 km/h) high resolution electrical mapping of the subsurface below large water bodies to depths often a factor of 10 greater than other towed instruments. The novel FloaTEM system is demonstrated at a range of diverse 4th through 6th-order riverine settings across the United States including 1) the Farmington River, near Hartford, Connecticut; 2) the Upper Delaware River near Barryville, New York; 3) the Tallahatchie River near Shellmound, Mississippi; and, 4) the Eel River estuary, on Cape Cod, near Falmouth, Massachusetts. Airborne frequency-domain electromagnetic and land-based towed TEM data are also compared at the Tallahatchie River site, and streambed geologic scenarios are explored with forward modeling. A range of geologic structures and pore water salinity interfaces were identified. Process-based interpretation of the case study data indicated FloaTEM can resolve varied sediment-water interface materials, such as the accumulation of fines at the bottom of a reservoir and permeable sand/gravel riverbed sediments that focus groundwater discharge. Bedrock layers were mapped at several sites, and aquifer confining units were defined at comparable resolution to airborne methods. Terrestrial fresh groundwater discharge with flowpaths extending hundreds of meters from shore was also imaged below the Eel River estuary, improving on previous hydrogeological characterizations of that nutrient-rich coastal exchange zone. In summary, the novel FloaTEM system fills a critical gap in our ability to characterize the hydrogeology below surface water features and will support more accurate prediction of groundwater/surface water exchange dynamics and fresh-saline groundwater interfaces.

Longitudinal conductance and transverse resistance are bulk electrical properties of a geoelectric formation and can be correlated with the bulk transmissivity of an aquifer within the formation. In aquifers that are electrically-resistive relative to adjacent layers in a horizontally layered Earth, bulk transverse resistance has been shown to correlate with aquifer transmisivity. Conversely, in aquifers that are electrically-conductive relative to adjacent layers, a correlation between bulk longitudinal conductance and aquifer transmissivity has been demonstrated. In both cases, previous investigations have relied on small datasets that have yielded coefficients of determination (R2) that are typically in the range of 0.6 to 0.7 to substantiate these relations. This paper explores electric-hydraulic relations of the surficial aquifer of Emirate Abu Dhabi, United Arab Emirates, by using a larger dataset compared to previous studies of these relations. Geophysical logs from 26 wells are used to correlate the bulk electrical properties of the aquifer, which are computed from deep electromagnetic (EM) induction resistivity logs, to the bulk transmissivity of the surficial aquifer, which is modeled from neutron porosity logs. The bulk longitudinal conductance was highly correlated (R2 = 0.8387) with the bulk transmissivity of the surficial aquifer. The bulk Dar-Zarrouk resistivity is negatively correlated with the aquifer transmissivity; an R2 value of −0.5604 is obtained.

The Mississippi Alluvial Plain (MAP) aquifer system is a vital resource that supports agriculture in one of the most productive regions of the country. The U.S. Geological Survey Water Availability and Use Science Program (WAUSP) is conducting a multi-discipline investigation of the MAP aquifer system. The investigation is utilizing borehole, surface, and airborne geophysical methods to improve the characterization and understanding of the aquifer. The combination of geophysical data collected over a range of spatial scales with varying depths of investigation and resolution is key to determining the distribution of sand and clay within the aquifer. Mobile geophysical methods that enable continuous measurements over large areas improve aquifer characterization with their increased spatial coverage. In support of the MAP investigation, a new towed Time-Domain Electromagnetic (tTEM) imaging system developed by Aarhus University was used near Shellmound, Mississippi (MS), to delineate the distribution of coarse-and fine-grained sediments underlying the site. The tTEM results compare favorably with the results of airborne EM (AEM) surveys flown in the study and improve the resolution of sand and gravel distribution within the tTEM depth of investigation.

Large sinkholes (50 to 200 m in diameter) formed over Permian bedded salt near Wink in western Texas in June 1980 and May 2002 and on the flank of a salt dome at Daisetta, Texas in May 2008. Residents, government officials, and industry representatives would like to better understand the potential for sinkhole formation and growth. We are applying several geophysical methods to better understand sinkhole precursors and assess the potential for future sinkhole development. At Wink, limited data on vertical ground movement from standard surveying has been greatly extended by satellite-based radar interferometry, which we are using to delineate areas having recent movement and determine rates of movement. Results from interferometry are guiding acquisition of high-resolution gravity data, which has identified shallow-source mass deficits that indicate potential for continued subsidence or sinkhole formation. At Daisetta, we are using interferometry to determine whether vertical movement preceded sinkhole collapse and gravimetry to identify areas where shallow mass deficits exist across the salt dome. These initial data are helping to design a more comprehensive subsurface investigation and develop a monitoring program based on radar interferometry and geodetic GPS measurements.

A reconnaissance-level geophysical study was conducted at Amistad Reservoir Dam to provide technical assistance for the horizontal and vertical placement of piezometers on the eastern embankment of the 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 directcurrent (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. 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. Because the project was a reconnaissance-level study, the report is not a comprehensive geologic or geophysical evaluation of the study area.

A geophysical characterization of a portion of American River levees in Sacramento, California was conducted in May 2007. Targets of interest included sand lenses that underlay the levees and the depth to a clay unit that underlies the sand. The concern is that the erosion of these sand lenses can lead to levee failure in highly populated areas of Sacramento. DC resistivity and electromagnetic surveys were conducted over a 6 mile length of the levee on roads and bicycle and horse trails. 2-D inversions were conducted on all the geophysical data. The OhmMapper and SuperSting surveys produced consistent inversion results that characterized the targets of interest. GEM-2 apparent resistivity data were consistent with the DC inversion results. However, the GEM-2 data could not be inverted due to large system drifts. While this would not be as large a problem in conductive terrains, it is a problem for a small induction number electromagnetic profiling system such as the GEM-2 in a resistive terrain (the sand lenses). Despite issues with the GEM-2 inversion, this geophysical investigation was successful in detection of the sand lenses and the depth to the clay zone of interest.

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 ( Stanton and others, 2007, Stanton and others, 2008).

In the North Platte River Basin, a ground-water model is being developed to evaluate the effectiveness of using leakage of water from selected irrigation canals to enhance ground-water recharge. The U.S. Geological Survey, in cooperation with the North Platte Natural Resources District, used land-based capacitively coupled and water-borne direct-current continuous resistivity profiling techniques to map the lithology of the upper 8 meters and interpret the relative canal leakage potential of 110 kilometers of the Interstate and Tri-State Canals in western Nebraska and eastern Wyoming. Lithologic descriptions from 25 test holes were used to evaluate the effectiveness of both techniques for indicating relative grain size. An interpretive color scale was developed that symbolizes contrasting resistivity features indicative of different grain-size categories. The color scale was applied to the vertically averaged resistivity and used to classify areas of the canal as having either high, moderate, or low canal leakage potential. When results were compared with the lithologic descriptions, both land-based and water-borne continuous resistivity profiling techniques were determined to be effective at differentiating coarsegrained from fine-grained sediments. Both techniques were useful for producing independent, similar interpretations of canal leakage potential.