Search USGS⌕ Search

Geology topics

Neil C. Terry

Publications and source records attributed to Neil C. Terry.

At least 19 recordsLinked to original sources

Springs of the Catskill Mountains: A conduit flow system in predominantly siliciclastic strata

The Catskill Mountains of New York state are composed of thick sequences of interbedded mudstone, sandstone, and siltstone that are intercalated with thin beds of paludal limestone and calcareous mudstone. Sandstones often contain rip-up clasts of calcareous material at their bases and contain minor amounts of carbonate cement throughout their mass. This minor amount of calcareous material weathers out to produce thin zones of solutional porosity and permeability integrated within perched aquifers that feed perennial hillslope springs. Solutional features are evident within thin carbonate lenses and along fractures in sandstones, occasionally resulting in explorable cave-sized conduits. This study focuses on the behavior of bedrock springs and seeps in the headwaters of the Neversink River basin above the Neversink Reservoir. Bedrock springs are numerous throughout the Catskill Mountains and serve to provide baseflow to the headwater streams that tends to be more persistent during dry times than other baseflow sources. We observed dissolution-influenced conduits within the bedrock, and suffosional conduits within unconsolidated materials. Continuous water temperature at six springs was monitored at hourly intervals from May to October 2022 to investigate the thermal response of groundwater to daily and seasonal changes in air temperature along with event-based meteorological events. Temperature monitoring was augmented at some springs with time-lapse photography to evaluate persistence of flow. The thermal responses of springs to precipitation events were highly variable across sites, with some springs exhibiting flashy responses indicating a strong connection to the surface, while other springs showed muted or no thermal response to events. The seasonal thermal response also varied across sites having lag times of 23-48 days between peak air temperature and peak spring water temperature. The water temperatures and geochemical compositions of these bedrock springs support a conceptual model of multiple residence times, showing both event and seasonal groundwater storage within the conduit flow system that feeds the springs of the Catskill Mountains.

New York↗

Shallow hydrogeologic framework of the Tully Valley mudboil area, Onondaga County, New York

Mudboils have been documented in the Tully Valley in southern Onondaga County, New York, since the late 1890s. Sediment-laden water from the mudboils flows into Onondaga Creek, which empties into Onondaga Lake at Syracuse 15 miles to the north. Turbidity from the mudboils has degraded the water quality of Onondaga Creek despite a series of mitigation efforts that began in the early 1990s. Turbidity mitigation actions presently (2025) being considered include creek relocation and offline sediment settling. In support of these proposed actions during 2021–23, the U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation, U.S. Environmental Protection Agency, Onondaga Nation, Onondaga Environmental Institute, and Central New York Regional Planning and Development Board, collected and analyzed geologic, hydrologic, geophysical, and geotechnical data to characterize the shallow hydrogeology along four proposed creek-relocation paths and in the proposed offline settling basin area. The investigation indicated that the four proposed creek-relocation paths, two east of Onondaga Creek and two west of Onondaga Creek, are underlain by sediments including muck, alluvium, mudboil deposits, alluvial-fan sand and gravel, and lacustrine fines. The proposed excavations would penetrate partially to fully saturated conditions: generally, the water table is shallow near the creek and deep on the alluvial fans. The shallowest excavation, about 5 feet below land surface, would be near the creek and primarily in alluvium, and the deepest excavation, as much as 30 feet below land surface, would be in the alluvial-fan deposits. Brackish waters would be penetrated by proposed channel excavations on the eastern side of Onondaga Creek in an area downgradient from a potentially leaking historical salt-exploration borehole and near the main mudboil area. Excavation in these areas likely would provide a continuous source of brackish groundwater to the relocated creek. Proposed channel excavations of muck, soft to very soft lacustrine fines, and mudboil-type sediments in mudboil and suspected mudboil areas would pose an excavation and slope stability challenge and would have the greatest potential to create new mudboils. Proposed channel excavations below the water table on the Rattlesnake Gulf and Rainbow Creek alluvial fans would intercept groundwater and make the constructed streambank susceptible to seepage-induced slope instability. The substantial water-level fluctuation in the sediments of both alluvial fans would aggravate the stability condition. In addition, excavation on the Rattlesnake Gulf alluvial fan would have the potential to affect water-supply springs at the toe of the fan. The proposed offline settling basin area is in the northern part of the Rattlesnake Gulf alluvial fan. Natural and man-made diversions of Rattlesnake Gulf have resulted in saturated conditions in the general area of the proposed basin. The proposed offline settling basin would be excavated in, and berms would be constructed on, alluvial-fan deposits and lacustrine fines. In the proposed basin area, the alluvial deposits overlying the lacustrine fines are less than 10 feet thick. Excavation, berm construction, and loading of the saturated, soft to very soft lacustrine fines may be problematic and require soil strengthening.

New York↗

A targeted approach for mapping groundwater discharge to surface water and fish thermal refuge in four Lake Ontario tributaries

The duration, magnitude, and frequency of heatwaves are predicted to increase in the coming decades, a combination that can reduce the survival of many fish species. Across the world, there is broad interest in identifying thermal refuge for heat-intolerant fish species and exploring opportunities to enhance or protect these areas. Because deeper groundwater maintains a relatively constant temperature, groundwater-influenced areas along streams can provide cool-water refuge for fish during periods of extreme heat. A targeted approach was developed for identifying existing cold-water zones and areas of substantial groundwater discharge in four high priority Lake Ontario tributaries. Our approach included: (1) predicting where groundwater discharge is most likely with a simple geospatial model and (2) using model predictions to select field sites for intensive high-resolution study, including ground-based mapping of groundwater features (springs, seeps, tributaries) as well as drone-based optical and thermal infrared surveys. Results from field sites were used to both verify model performance and map different types and aerial extents of thermal anomalies. Geospatial modelling successfully predicted regions of widespread groundwater upwelling, later verified and mapped by field and drone surveys. Comparison of model and field survey results further highlighted specific geospatial layers, such as soil/bedrock types and topographic wetness index, as being particularly useful for predicting groundwater influence on streams in the study area. In addition, a comparison of geospatial model results with a model of fish abundances along the studied streams showed significant positive correlations for many heat-intolerant fish species over a wide geographic area. The approach developed in this study can be applied to other watersheds to highlight areas of probable groundwater discharge and could be used by fishery and water resource managers to support cold-water fish habitat management decision-making and resource conservation.

New York↗

Multi-scale geophysical mapping of the brine and bedrock surfaces along the Dolores River, Paradox Valley, Colorado, December 2023

Total dissolved solids derived from salt dome–sourced brine in the underlying alluvial aquifer substantially increase with distance in the reach of the Dolores River that passes through Paradox Valley in southwestern Colorado. The area has been the site of salinity control operations since the 1990s to reduce salt loading to the downstream Colorado River. Previous airborne and ground/water-based electromagnetic (EM) geophysical data have successfully mapped the top of the brine surface, albeit with relatively coarse near-surface resolution and limited spatial coverage. This present December 2023 study used ground-based high-resolution EM and passive seismic (horizontal-to-vertical spectral ratio, HVSR) tools to map in detail the depth and thickness of the brine zone in the alluvial aquifer (top of the brine down to bedrock contact) in areas immediately surrounding the Dolores River where previous airborne EM (AEM) results indicated brine within 10 m of land surface. Results indicate the deepest bedrock is generally associated with the shallowest brine and local depressions in the collapse breccia (caprock to the Paradox Formation salt) may facilitate vertical migration of brine into the alluvial aquifer. Additionally, the ground-based EM mapping corroborated general patterns in depth to brine that were observed in previous AEM results while also revealing additional detail, including suspected focused brine discharge zones to the Dolores River. A river-based survey wherein EM data, channel depth and river water–specific conductance information were collected augmented these findings and indicated specific areas on both the western and eastern sides of the river where focused brine discharge may occur. This study comprises a large-scale, ground- and water-based geophysical mapping effort, including hundreds of HVSR soundings and 100s of kilometres of EM data, which were successfully translated into digital brine and bedrock surfaces that could be incorporated into groundwater modelling efforts, future well siting or other decision-making.

Colorado↗

An evaluation of the effects of different deicing salt application rates on three watersheds in Essex County, New York

The U.S. Geological Survey, in cooperation with the New York State Department of Transportation, evaluated the effects of different deicing salt application rates on surface water, groundwater, and highway runoff quality near State highways in northern New York. Three reaches of State highways were tested with different deicing treatments between October 2019 and November 2022: a salt-sand mixture (Treatment A), a salt mixture applied at a lower rate (Treatment B), and a control mixture consistent with typical deicing salt amounts and application rates. Data on pavement conditions and the quality of surface water, highway runoff, and groundwater were collected. Surface electromagnetic data were also collected. Surface-water and groundwater quality downgradient from the State highways were compared with water quality at upgradient locations. The percentage of snow or ice coverage was used to evaluate the effectiveness of the salt applications. This report provides an overview of the transport of deicing salt. The Treatment B watershed had deicing mixture applied more frequently than other highway reaches, which caused it to have the highest annual total chloride application. Despite differences in chloride application, flow-weighted mean chloride concentrations in highway runoff were comparable across treatments. Chloride concentrations were elevated in surface water and groundwater downgradient from highways relative to chloride concentrations upgradient from highways. A chloride mass balance, calculated for one treatment watershed, indicated that groundwater affected by legacy deicing practices may be contributing additional chloride to surface water. Spatial patterns from electromagnetic surveys show a shallow saline plume alongside the highway in that area. Differences in winter severity and pavement-surface conditions drove deicing salt applications in the treatment areas. This study found that several factors affect chloride loads in the watersheds, including variable winter conditions, adaptive snow and ice management, legacy management practices, and area-specific aquifer and groundwater conditions.

New York↗

Exploring the use of non-invasive drone-based Ground Penetrating Radar (GPR) to characterize biogenic gas dynamics in subtropical peat soils

Peat soils are a critical component of the global carbon cycle as natural producers of biogenic greenhouse gases (e.g., methane and carbon dioxide) that accumulate within the soil matrix and are episodically released to the atmosphere. Previous studies have showed the unique ability of ground-based minimally-invasive geophysical methods like ground-penetrating radar (GPR) to characterize carbon dynamics in peat soils. However, the GPR methods are still limited by the scale of measurements and soil disturbance potentially altering biogenic gas releases during deployment. Here we explore the potential of drone-based GPR for identification of hot spots and hot moments of gas accumulation and release in subtropical peat soils. We collected drone-based GPR datasets across two grids (~ 17,500 m2) in the Florida Everglades during January, September, and November of 2023 to characterize peat thickness and seasonal variability of gas content. Results show that drone-based GPR is effective and efficient for: 1) capturing the temporal variation of in-situ biogenic gas content in peat soils with changes between 1-25% volumetric gas content over repeatable grids; 2) inferring peat thickness from discrete measurements ranging between 0.8-1.2 m; and 3) estimating methane flux releases of 63 and 135 mg CH4 m-2 day-1 for specific locations and time periods that are strikingly consistent with our coincident gas trap measurements. This work also indicates that (1) spatial distribution of gas content in the Everglades is strongly controlled by landscape morphology such as ridges and sloughs and (2) the temporal variation of gas content is seasonal with increased gas production during the wet season.

Florida↗

Multi-scale geophysical imaging of a hydrothermal system in Yellowstone National Park, USA

Little is known about the local plumbing systems that fuel Yellowstone’s famous hot springs, geysers and mud pots. A multi-method, multi-scale geophysical investigation was carried out in the Obsidian Pool Thermal Area (OPTA) to: (i) delineate the lateral extent of the hydrothermal area and associated surface features; (ii) estimate the dimensions of the upflow zone and identify its main controlling structures; (iii) assess fluids circulation pathways from depth to surface. Ground and airborne geophysical data were acquired to connect local and regional scales, from shallow to large depths. Maps of surface electrical resistivity show a strong correlation with hydrothermal features. At in-termediate depths, electrical resistivity permits delineating the upper limit of the upflow zone, while Poisson’s ratio highlights differences in subsurface fluid content. Combining these results with surface observations and topographic information, we speculate that differential mixing of hydrothermal and fresh water could explain the wide diversity of features observed at OPTA. Low electrical resistivity observed at large depths also suggest that a vast upflow zone, controlled by rhyolite flows and conjugate faults, underlies the OPTA. We speculate that hydrothermal fluids rise along fractures and reach the surface in topographic lows to form hydrothermal features. Our results show that synoptic, multi-scale geophysical measurements provide a roadmap for understanding where and how geologic heterogeneity, topography, fluid-gas separation, and the mixing of thermal and meteoric waters conspire to produce the wide variety of Yellowstone’s renowned hydrothermal features.

Wyoming↗

Fiber-optic distributed temperature sensing of hydrologic processes—Diverse deployments and new applications by the U.S. Geological Survey

Fiber-optic distributed temperature sensing instruments harness the temperature-dependent properties of glass to measure temperature continuously along optical fibers by using precise pulses of laser light. In the mid-2000s, this technology was refined for environmental monitoring purposes such as snowpack-air exchange, groundwater/surface-water exchange, and lake-water stratification. Fiber-optic distributed temperature sensing has revealed unprecedented details about preferential flow processes; however, the method is labor intensive and requires specific training, resulting in limited use by the broader water community. With the ongoing national implementation of the U.S. Geological Survey Next Generation Water Observing System, there has been renewed interest in harnessing the unique spatiotemporal monitoring capabilities of fiber-optic distributed temperature sensing. This fact sheet briefly describes this technology, highlights uses by the U.S. Geological Survey, and discusses current applications and future opportunities.

Fact Sheet↗

ARCHI: A new R package for automated imputation of regionally correlated hydrologic records

Missing data in hydrological records can limit resource assessment, process understanding, and predictive modeling. Here, we present ARCHI (Automated Regional Correlation Analysis for Hydrologic Record Imputation), a new, open-source software package in R designed to aggregate, impute, cluster, and visualize regionally correlated hydrologic records. ARCHI imputes missing data in “target” records by linear regression using more complete “reference” records as predictors. Automated imputation is implemented using a novel, iterative algorithm that allows each site to be considered a target or reference for regression, growing the pool of complete references with each imputed record until viable gap-filling ceases. Users can limit artifacts from spurious correlations by specifying model-acceptance criteria and applying geospatial, correlation, and group-based filters to control reference selection. ARCHI provides additional functions for visualizing results, clustering records with similar correlation structures, evaluating holdout data, and interactive parameterization with an accessible and intuitive graphical user interface (GUI). This methods brief provides an overview of the ARCHI package, modeling guidelines, and benchmarking on two regional groundwater-level datasets from the Central Valley, CA and Long Island, NY. We evaluate ARCHI alongside widely used multivariate imputation software to highlight and contextualize its computational efficiency, imputation accuracy, and model transparency when applied to large, groundwater-level datasets.

California, New York↗

Geophysical identification of potential groundwater – surface water interactions in the Bonita Peak Mining District, San Juan County, Colorado

This report includes the results and interpretation of geophysical investigations to aid in the location, identification, and conceptual site model (CSM) development of the near surface groundwater and surface water hydrology of the Upper Cement Creek and California Gulch of the Animas Rivers within the Bonita Peak Mining Districts (BPMD). Mine influenced water (MIW) from numerous mine adits within the BPMD has altered the geochemistry and ecology of these ground and surface waters. Many studies of varying scales have investigated the geology, hydrogeology, geochemistry, and aqueous chemistry of these adits and impacted areas. Field investigations using non-invasive geophysical methods were conducted in the summer of 2019 to aid in identifying potential groundwater inputs to these streams and further the CSM of the near surface within and along these streams. Two field deployments were completed in August and September of 2019. In August 2019, Phase 1 was completed using multifrequency electromagnetic induction (EMI) and thermal infrared imaging (TIR) to measure the bulk electrical conductivity, magnetic susceptibility, and the temperature response at known and suspected groundwater seeps into the Upper Cement Creek and California Gulch stretches of these alpine streams. Results and interpretation from Phase 1 were used to co-locate stretches of these streams for further detailed fiber optic distributed temperature system (FODTS) monitoring of the temperature dynamics for submeter identification of potential groundwater seeps during Phase 2. Phase 2 included additional EMI measurements, FODTS measurements, and repeated some Phase 1 measurements for improved data quality. All these data were used for the combined qualitative interpretation identifying seep locations. The data from both Phase 1 and Phase 2 are available in the USGS data release at: https://www.sciencebase.gov/catalog/item/5f18a30082cef313ed84879b. These data include detailed near-surface electrical conductivity and magnetic susceptibility mapping, discrete TIR imaging, and several kilometers of FODTS data. These data were georeferenced and used to identify potential groundwater seeps into these streams. The results indicate a complex near surface geophysical and hydrogeologic regime. Seeps are classified and spatially correlated to the EMI, TIR, and geochemical sampling results. This report is a thorough compilation of the geophysical field efforts in the summer of 2019, the data processing, and interpretations of these data. These data alone should not be used for final interpretations of the complex hydrogeology within these alpine systems. Rather, these data should be combined with other known data from various investigations to further refine the CSM and guide decisions for future investigations, watershed protective measures, and other site decisions. Geographical information system (GIS) layers of these results are available to project managers for such integration and combined interpretations.

Colorado↗

Geophysical characterization of mine influenced groundwater and surface water in the Mayflower section of the Animas River, Bonita Peak Mining District, Silverton Colorado

This report details findings from geophysical investigations to identify possible groundwater - surface water interactions near the Mayflower Section of the Animas River in Silverton, Colorado. The Mayflower Section is a mine influenced Superfund Site in Colorado. This investigation utilized electromagnetic induction (EM or EMI), magnetic, and fiber optic distributed temperature system (FODTS) geophysical methods to measure the bulk earth electrical conductivity, magnetic susceptibility, and temperature of specific surveyed volumes of the earth near the Mayflower Section of the Animas River. These physical parameters are used to understand the groundwater – surface water interactions, which can guide decision makers in their assessment of mine-impacted surface water. This report details the results from characterization and monitoring technologies to provide high data density and continuous monitoring of bulk earth electrical conductivity and temperature in the Mayflower section of the Animas River to identify zones of groundwater – surface water interactions and potential metal loading from mine influenced water. The investigation separated right and left bank characterization for each method and indicates more groundwater is entering from the right bank than the left bank and these predominantly right bank discharges potentially contain metal-rich water compared to the left bank. Results also indicate mineral veins facilitate preferential groundwater discharge to the river due to possible jointing, fractures, and permeability differences sometimes occurring along veins relative to host rock. For example, Boulder Gulch is likely groundwater dominated and may be receive water impacted by Mayflower Impoundments #1 and #2. Additionally, the beaver ponds near Blair Gulch may influence groundwater discharge and Mayflower Impoundment #4 is possibly impacting groundwater and surface water and may be connected to the wetlands to the west of the impoundment. These data could be further analyzed for smaller spatial scale analysis within areas of interest. The identification of these locations along sections of the river likely impacted by mine influenced groundwater potentially entering the Animas River and can be used by site investigators, decision makers, and stakeholders in mitigation decisions and strategies.

Colorado↗

Paired Air and Stream Temperature Analysis (PASTA) to evaluate groundwater influence on streams

Groundwater is critical for maintaining stream baseflow and thermal stability; however, the influence of groundwater on streamflow has been difficult to evaluate at broad spatial scales. Techniques such as baseflow separation necessitate streamflow records and do not directly indicate whether groundwater inflow may be sourced from more dynamic shallow flowpaths. We present a web tool application PASTA (Paired Air and Stream Temperature Analysis; https://cuahsi.shinyapps.io/pasta/ ) that capitalizes on increased public stream temperature data availability and large-scale, gridded climate observations to provide new and efficient insights regarding relative groundwater influence on streams. PASTA analyzes paired air and stream water temperature signals to evaluate spatiotemporal patterns in stream thermal sensitivity and relative groundwater influence, including inference regarding the dominant source groundwater depth (shallow or deep (i.e., approximately >6 m depth)). The tool is linked to publicly available stream temperature datasets and accepts user-uploaded datasets. As local air temperature is not often monitored, PASTA pulls daily air temperature data from the comprehensive Daymet products when directly measured data are unavailable, allowing the repurposing of existing stream temperature data. After data are selected or uploaded, the tool (a) fits sinusoidal curves of daily stream and air temperatures by year (water or calendar) to indicate groundwater influence characteristics and (b) performs linear regressions for stream versus air temperatures to indicate stream thermal sensitivity. Results are exported in ASCII file format, creating an efficient and approachable analysis tool for the adoption of newly developed heat tracing analysis from stream reach to landscape scales.

Water Resources Research↗

Field evaluation of semi-automated moisture estimation from geophysics using machine learning

Geophysical methods can provide three-dimensional (3D), spatially continuous estimates of soil moisture. However, point-to-point comparisons of geophysical properties to measure soil moisture data are frequently unsatisfactory, resulting in geophysics being used for qualitative purposes only. This is because (1) geophysics requires models that relate geophysical signals to soil moisture, (2) geophysical methods have potential uncertainties resulting from smoothing and artifacts introduced from processing and inversion, and (3) results from multiple geophysical methods are not easily combined within a single soil moisture estimation framework. To investigate these potential limitations, an irrigation experiment was performed wherein soil moisture was monitored through time, and several surface geophysical datasets indirectly sensitive to soil moisture were collected before and after irrigation: ground penetrating radar, electrical resistivity tomography (ERT), and frequency domain electromagnetics (FDEM). Data were exported in both raw and processed form, and then snapped to a common 3D grid to facilitate moisture prediction by standard calibration techniques, multivariate regression, and machine learning. A combination of inverted ERT data, raw FDEM, and inverted FDEM data was most informative for predicting soil moisture using a random regression forest model (one-thousand 60/40 training/test cross-validation folds produced root mean squared errors ranging from 0.025–0.046 cm 3 /cm 3 ). This cross-validated model was further supported by a separate evaluation using a test set from a physically separate portion of the study area. Machine learning was conducive to a semi-automated model-selection process that could be used for other sites and datasets to locally improve accuracy.

Connecticut↗

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

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

Wyoming↗

Advancing subsurface investigations beyond the borehole with passive seismic horizontal-to-vertical spectral ratio and electromagnetic geophysical methods at transportation infrastructure sites in New Hampshire

The U.S. Geological Survey (USGS), in cooperation with the New Hampshire Department of Transportation (NHDOT), surveyed transportation infrastructure sites using rapidly deployable geophysical methods to assess benefits added to a comprehensive site characterization with traditional geotechnical techniques. Horizontal-to-vertical spectral-ratio (HVSR) passive-seismic and electromagnetic-induction (EMI) methods were applied at 4 sites including a roadway-stream crossing, roadway-bridge rail-trail crossing, commuter-parking expansion, and a railroad-adjacent river-cutbank slope-failure site. Additionally, ground-penetrating-radar (GPR) was used at the slope-failure site. Typically, subsurface geotechnical properties are determined from boring data; however, borings are often spaced hundreds of feet apart, potentially missing important spatial variability between boreholes. Geotechnical site characterization including geophysical surveys helped provide a more accurate characterization by using continuous or near continuous profiling. Three-component ambient noise measured with HVSR methods were used to determine resonance frequency and estimate sediment thickness. The method works when there is a strong shear-wave acoustic impedance contrast (> 2:1) between sediment and bedrock. Sediment thickness estimates from HVSR measurements were combined with boring data to make detailed maps of the bedrock surface altitude. The bulk electrical conductivity of the subsurface was indirectly measured with EMI methods and was used to identify lithologic variations, shallow bedrock, and conductive groundwater. Ground penetrating radar, which transmits pulses of electromagnetic energy into the subsurface and records the amplitude and timing of reflected signals, was used to identify bedding and changes in lithology or water content. By combining geophysical and boring data analyses, transportation projects produced more spatially comprehensive representations of geotechnical subsurface conditions than would be determined using conventional borings alone.

New Hampshire↗

Geostatistical mapping of salinity conditioned on borehole logs, Montebello Oil Field, California

We present a geostatistics-based stochastic salinity estimation framework for the Montebello Oil Field that capitalizes on available total dissolved solids (TDS) data from groundwater samples as well as electrical resistivity (ER) data from borehole logging. Data from TDS samples ( n = 4924) was coded into an indicator framework based on falling below four selected thresholds (500, 1000, 3000, and 10,000 mg/L). Collocated TDS-ER data from the surrounding groundwater basin were then employed to produce a kernel density estimator to establish conditional probabilities for ER data ( n = 8 boreholes) falling below the selected TDS thresholds within the Montebello Oil Field area. Directional variograms were estimated from these indicator coded data, and 500 TDS realizations from conditional indicator simulation were generated for the subsurface region above the Montebello Oil Field reservoir. Simulations were summarized as 3D maps of median TDS, most likely salinity class, and probability for exceeding each of the specified TDS thresholds. Results suggested TDS was below 500 mg/L in most of the study area, with a trend toward higher values (500 to 1000 mg/L) to the southwest; consistent with the average regional groundwater flow direction. Discrete localized zones of TDS greater than 1000 mg/L were observed, with one of these zones in the greater than 10,000 mg/L range; however, these areas were not prevalent. The probabilistic approach used here is adaptable and is readily modified to include additional data and types and can be employed in time-lapse salinity modeling through Bayesian updating.

California↗

Beware of spatial autocorrelation when applying machine learning algorithms to borehole geophysical logs

Although many of the algorithms now considered to be machine learning algorithms (MLAs) have existed for nearly a century (e.g., Rosenblatt 1958 ), interest in MLAs has recently increased exponentially for solving data-driven problems across a variety of fields due to the expanded availability of large, complex datasets that may be difficult to interrogate using other methods, increases in computing power, and a growing library of easily implemented machine learning tools. While MLAs are often similar to statistical methods, there are key differences in the approach to problem solving. Namely, statistical methods are more concerned with generating informative models from “long” data (i.e., many more observations than explanatory variables), whereas MLAs are typically concerned with generating accurate predictions from “wide” data (i.e., a large number of variables with relatively fewer observations, Bzdok et al. 2018 ). In hydrogeologic studies, such wide datasets may be available from boreholes, where various types of geophysical, geochemical, and lithological information may exist. Borehole datasets are therefore a tempting target for MLAs to reveal hidden relations among gathered data and parameters of interest (e.g., contaminant concentration), and as a method of parameter reduction (e.g., reduce costs by collecting fewer datasets).

Groundwater↗

Exploring the potential of ground-penetrating radar (GPR) to measure the extent of chronic disturbance in peatlands: Examples from acid mine drainage and peat fire

Peatlands are accumulations of partially decayed organic soil that cover approximately 3% of Earth’s surface and have been shown to serve essential environmental and ecological functions such as sequestering carbon, purifying water, and providing habitat for organisms. However, peatlands are threatened by pressures from agriculture, urban development, mining, and climate change. Geophysical methods have been used in peatlands to determine peat volume and carbon stocks (e.g., Comas et al., 2017), observe differences in humification and water content (e.g., Ulriksen, 1982), guide engineering projects (e.g., Jol and Smith, 1995), learn about subsurface greenhouse gas dynamics (Wright and Comas, 2016), observe seasonal variations in pore water salinity (Walter et al., 2018), and assess hydrological processes (Hare et al., 2017). Among various geophysical methods, ground penetrating radar (GPR) is arguably the most popular for studying peat properties given the method’s sensitivity to variations in water content and ability to resolve major structural properties within the peat at high spatial resolution. Though less widely applied, frequency-domain analysis of GPR may also yield useful information.

Conference Paper↗