Geology topics

Burke J. Minsley

Publications and source records attributed to Burke J. Minsley.

At least 19 recordsLinked to original sources

USGS2026

Subsurface vertical connectivity shapes solute transport to montane streams: Insights from watershed-scale geophysics-informed modeling

Headwater catchments provide essential water and nutrients to downstream ecosystems. The timing and evolution of their transport is shaped in part by the “invisible” subsurface structure but have been largely unexplored at the watershed scale. This raises a fundamental question: how and to what extent does subsurface structure influence streamflow and solute transport in montane watersheds? Here, we test controls of three-dimensionally resolved subsurface structure on solute transport, which reflects source waters and water residence time. We integrate the geophysically mapped architecture of Coal Creek, a headwater catchment in the Colorado River Basin, into numerical flow-transport models. Our results reveal greater vertical connectivity elongates flow paths and enhances deep groundwater contributions to streams, therefore substantially influencing the timing and magnitude of solute transport. In contrast, watershed-integrated outlet discharge can be predicted without high-resolution subsurface information. These findings demonstrate incorporating subsurface architecture is critical for predicting ecosystem health and contaminant transport.

Colorado↗

USGS2026

The saline groundwater legacy of a large buried coastal paleo-estuary

Elevated groundwater salinity in coastal regions threatens the beneficial use of fresh groundwater. Coastal groundwater management typically focuses on preventing intrusion from modern sources of seawater; however, past geological processes can also leave a legacy of saline groundwater now hidden in the subsurface. Here, multiple extensive airborne electromagnetic surveys provide detailed evidence of residual salinity from a paleo-estuary filling a late Pleistocene incised valley impacting more than 10,000 km 2 that is now hidden beneath coastal Louisiana’s deltaic plain. Our results show that the three-dimensional pattern of saline groundwater beneath Louisiana mimics that of near-surface aquifers surrounding the modern Delaware Bay estuary, fingerprinting the signature of the past drowning of a large, incised valley of the Mississippi River following post-glacial sea-level rise. These findings demonstrate a new framework for understanding legacy sources of saltwater critical for managing stressed water resources along global coastlines.

Louisiana↗

USGS2026

Computational electromagnetic geophysics for groundwater system studies: A review on established practices and recent advances

Identifying effective solutions for locating groundwater resources and ensuring the quality of drinking water is increasingly urgent, given the challenges posed by climate change and population growth. This review investigates electromagnetic geophysical imaging techniques, in both time- and frequency-domain, that can provide valuable insights for groundwater assessment. We explore computational electromagnetic methods used to evaluate electromagnetic data and several recent hydrogeophysical case studies. As open-source frameworks for modeling electromagnetic geophysical problems become available, a broader range of researchers can interpret their data with computationally advanced software. We provide an overview of documented open-source codes for evaluating electromagnetic data and analyze various hydrological targets in relation to their electromagnetic surveying technique and the computational method applied. Furthermore, we evaluate the potential of advanced computational techniques, including three-dimensional modeling, non-deterministic inversion and machine learning, to couple geophysical with numerical groundwater modeling and apply it in groundwater system studies. Despite obstacles such as complexity and resource demands, our findings indicate that the quantification and integration of predictive uncertainties from both electromagnetic and hydrological data and simulations would significantly improve the reliability of hydrogeophysical models. This can lead to a deeper understanding of groundwater systems and improved management practices.

Journal of Hydrology↗

USGS2025

Quantitative subsurface characterization illuminates the origin of the Quaternary Mississippi River Valley alluvial aquifer

The Mississippi River Valley alluvial aquifer (MRVA) is vital to U.S. food security and global agricultural supply. However, quantitative understanding of its Quaternary origin, architecture, and hydrologic function remains incomplete. Here we develop a three-dimensional hydrostratigraphic model to characterize the deposition of clay and silt, fine-medium sands, and graveliferous sands using lithologic data from 75,000 boreholes compiled across the Lower Mississippi Valley and a geostatistical method—interval kriging. We find that cyclic glacial entrenchments, evidenced by remnants of pre-Wisconsinan postglacial sediments, alongside geodynamic activities shaped the MRVA basal configuration. Stratal weakening from faulting and salt diapirism enhanced glacial incision and thereby produced abrupt aquifer thickening. We demarcate the top of graveliferous sands as the regional marker of the Pleistocene-Holocene transition. The MRVA hydrostratigraphy reveals hydrologic function and geologic controls on groundwater storage and quality, advancing the assessment of aquifer sustainability under a changing climate, with implications for alluvial aquifers globally.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗

USGS2025

Uncertainty quantification of geophysical and hydrologic parameters estimated from borehole nuclear magnetic resonance data

Borehole nuclear magnetic resonance (bNMR) data are typically used to infer in situ hydrologic properties. Partial water content as a function of pore size is estimated by fitting the measured NMR response to a multi-exponential T 2 distribution, and the sum of estimated T 2 amplitudes equals the total volumetric water content. From these estimated parameters, several empirical relationships are commonly used to infer hydraulic conductivity from the NMR-estimated water content and T 2 distribution. Often, parameters are estimated through deterministic inversion methods that produce a single best-fit estimate, but do not reflect uncertainties in model parameters. Here, a Bayesian Markov chain Monte Carlo (McMC) approach for analyzing bNMR data is developed that allows for comprehensive uncertainty quantification of NMR parameters and derived hydrologic properties. The underlying model that describes the T 2 distribution is defined by a set of spline interpolation points. The number of interpolation points is allowed to vary in a trans-dimensional algorithm that naturally favors simple models with fewer interpolation points, allowing the data to inform the necessary level of model complexity. Additionally, data error is estimated as an unknown parameter. Analysis of the ensemble of models output from the McMC algorithm provides useful details on the range of plausible T 2 distributions that can fit a measured bNMR decay curve, as well as uncertainty estimates of total water content. The ensemble of NMR parameters can also be propagated through commonly used relationships to produce uncertainty estimates on derived parameters such as bound/capillary/mobile water content or hydraulic conductivity.

JGR Machine Learning and Computation↗

USGS2025

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗

USGS2025

Airborne geophysical analysis to decipher salinization for coastal Louisiana

Coastal Louisiana is known for saltwater intrusion that threatens wetlands, aquifers, and rivers. However, the extent of saltwater intrusion is not well understood. This study develops an innovative framework with airborne electromagnetic (AEM) data to map chloride concentration distributions for wetlands in the Mississippi River deltaic plain and Chenier plain as well as for the Mississippi River Valley alluvial aquifer (MRVA) and Chicot aquifer. Moreover, the framework maps chloride concentrations along the Mississippi River and Atchafalaya River. Key components in the framework include the establishment of resistivity-to-chloride concentration transformation, 3D resistivity architecture building through geostatistics, and the employment of a lithologic model. The transformation functions correlate AEM resistivity data with porewater salinity measurements and groundwater and river chloride samples. The results show that AEM data reliably infers soil water chloride concentrations and correlates well with the distribution of various marsh types. AEM data reveals extensive saltwater presence at depth and near the coast, originating from salt domes and the Gulf of Mexico, respectively. The saltwater upconing pattern in the Chicot aquifer is likely due to excessive groundwater withdrawals. The AEM data also confirms a distinct tongue of saltwater intruding into the Atchafalaya Basin from the Gulf. The AEM data helps to identify faults that are obscured or eroded at the surface, which appear as leaky barriers in the subsurface where dramatic changes in chloride concentration are apparent. Finally, this study uses the AEM data to infer the presence of an extensive seawater wedge in the Mississippi River and Atchafalaya River.

Louisiana↗

USGS2024

Beyond the wedge: Impact of tidal streams on salinization of groundwater in a coastal aquifer stressed by pumping and sea-level rise

Saltwater intrusion (SWI) is a well-studied phenomenon that threatens the freshwater supplies of coastal communities around the world. The development and advancement of numerical models has led to improved assessment of the risk of salinization. However, these studies often fail to include the impact of surface waters as potential sources of aquifer salinity and how they may impact SWI. Based on field-collected data, we developed a regional, variable-density groundwater model using SEAWAT for east Dover, Delaware. In this location, major users of groundwater from the surficial aquifer are the City of Dover and irrigation for agriculture. Our model includes salinized marshland and tidal streams, along with irrigation and municipal pumping wells. Model scenarios were run for 100 years and included changes in pumping rates and sea-level rise (SLR). We examined how these drivers of SWI affect the extent and location of salinization in the surficial aquifer by evaluating differences in chloride concentration near surface waters and the subsurface freshwater-saltwater interface. We found the presence of the marsh inverts the typical freshwater-saltwater wedge interface and that the edge of the interface did not migrate farther inland. Additionally, we found that tidal streams are the dominant pathways of SWI at our site with salinization from streams being exacerbated by SLR. Our results also show that spatial distribution of pumping affects both the magnitude and extent of salinization, with an increase in concentrated pumping leading to more intensive salinization than a more widely distributed increase of the same total pumping volume.

Delaware↗

USGS2024

Upland Yedoma taliks are an unpredicted source of atmospheric methane

Landscape drying associated with permafrost thaw is expected to enhance microbial methane oxidation in arctic soils. Here we show that ice-rich, Yedoma permafrost deposits, comprising a disproportionately large fraction of pan-arctic soil carbon, present an alternate trajectory. Field and laboratory observations indicate that talik (perennially thawed soils in permafrost) development in unsaturated Yedoma uplands leads to unexpectedly large methane emissions (35–78 mg m −2 d −1 summer, 150–180 mg m −2 d −1 winter). Upland Yedoma talik emissions were nearly three times higher annually than northern-wetland emissions on an areal basis. Approximately 70% emissions occurred in winter, when surface-soil freezing abated methanotrophy, enhancing methane escape from the talik. Remote sensing and numerical modeling indicate the potential for widespread upland talik formation across the pan-arctic Yedoma domain during the 21 st and 22 nd centuries. Contrary to current climate model predictions, these findings imply a positive and much larger permafrost-methane-climate feedback for upland Yedoma.

Alaska, Siberia↗

USGS2023

Quantifying salinity in the layered coastal aquifers underlying and adjacent to Delaware Bay USA using AEM-derived resistivity

Airborne electromagnetic (AEM) methods are particularly well suited to coastal aquifer salinity studies, yet the quantitative translation from bulk resistivity to fluid salinity carries uncertainty that can impact mapped salinity distributions and interpretations of the freshwater-saline interface and hydrostratigraphic layers. A recent AEM survey of the region near the Delaware Bay, USA highlights several challenges common to coastal hydrogeologic settings that may influence both qualitative and quantitative interpretation. We use a Bayesian inversion to estimate geophysical parameter uncertainty, and results are integrated with hydrogeologic measurements to develop quantitative interpretations of salinity across the freshwater-saline interface in stacked aquifers.

Delaware, New Jersey↗

USGS2023

Estimating streambed hydraulic conductivity for selected streams in the Mississippi Alluvial Plain using continuous resistivity profiling methods—Delta region

Introduction The Mississippi Alluvial Plain is one of the most important agricultural regions in the United States, and crop productivity relies on groundwater irrigation from an aquifer system whose full capacity is unknown. Groundwater withdrawals from the Mississippi River Valley alluvial aquifer have resulted in substantial groundwater-level declines and reductions in base flow in streams within the Mississippi Alluvial Plain. These effects are limiting well production and threatening future water availability in the region. A comprehensive assessment of water availability in the Mississippi Alluvial Plain is critically important for making well-informed management decisions about sustainability, establishing best practices for water use, and predicting changes to water levels in the Mississippi Alluvial Plain over the next 50–100 years. The first step in the new regional modeling effort was to run the existing Mississippi Embayment Regional Aquifer Study (MERAS) model and perform data-worth and uncertainty analyses to prioritize data collection efforts to improve model forecasts. Parameter estimation indicated that streambed conductance was one of the variables that the model was most sensitive to, but little data were available to constrain those general estimates. From this characterization of the existing data, a map of the streams that the MERAS model was most sensitive to was created by the U.S. Geological Survey to guide the collection of 862 kilometers of waterborne resistivity surveys within the Delta region of Mississippi to characterize streambed lithology. This technique characterizes the streambed itself and the 15–30 meters below the streambed that control the exchange of water between the stream and the alluvial aquifer. These data can be used to map changes in the lithology of the streambed and identify areas of potential groundwater/surface-water exchange. Additionally, electrical and nuclear well logs from the study area were compared to facilitate the development of a petrophysical relation between the waterborne resistivity data and hydraulic conductivity. Resistivity values may then be used as a cost-effective way to approximate aquifer hydraulic conductivity distributions for use in regional groundwater models.

Arkansas, Louisiana, Mississippi↗

USGS2023

A model of transmissivity and hydraulic conductivity from electrical resistivity distribution derived from airborne electromagnetic surveys of the Mississippi River Valley Alluvial Aquifer, Midwest USA

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 km 2 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-km 2 Shellmound (Mississippi) study area to calculate hydraulic property distributions of the alluvial aquifer for use in future groundwater-flow models.

Mississippi River Valley Alluvial Aquifer↗

USGS2023

System-scale airborne electromagnetic surveys in the lower Mississippi River Valley support multidisciplinary applications

The lower Mississippi River Valley spans over 200,000 square kilometres in parts of seven states, encompassing areas of critical groundwater supplies, natural hazards, infrastructure, and low-lying coastal regions. From 2018 - 2022, the U.S. Geological Survey acquired over 82,000 line-kilometres of airborne electromagnetic, radiometric, and magnetic data over this region to provide comprehensive and systematic information about subsurface geologic and hydrologic properties that support multiple scientific and societal interests. Most of the data were acquired on a regional grid of west-east flight lines separated by 3 - 6 kilometres; however, several high-resolution inset grids with line spacing as close as 200 m were acquired in targeted areas of interest. Approximately 8,000 line-kilometres were acquired along streams and rivers to characterise the potential for surface water-groundwater connection, and another 6,000 line-kilometres were acquired along the Mississippi and Arkansas River levees to characterise this critical infrastructure. Here, we present a summary of the data along with several examples of how they are being used to inform regional groundwater model development, inferences of groundwater salinity, identification of faults in the New Madrid seismic zone, and levee infrastructure.

lower Mississippi River Valley↗

USGS2023

Using airborne electromagnetics to improve depth to bedrock estimates in Wisconsin

Depth to bedrock is an important factor in hydrologic systems. The thickness of glacial sediments overlying bedrock in Wisconsin controls the routing of groundwater in surficial aquifers and its Depth to bedrock is often an important factor in hydrologic systems because hydraulic properties of bedrock and overlying sediments are typically appreciably different. For example, the thickness of glacial sediments overlying bedrock in Wisconsin controls the routing of groundwater in surficial aquifers and its connection with surface water bodies such as lakes and wetlands. In fractured bedrock environments, shallow bedrock can be vulnerable to degraded water quality when contaminants at the surface infiltrate quickly through permeable formations. Here, airborne electromagnetic surveys were acquired in three different parts of Wisconsin, totalling more than 5,700 flight line-kilometers, to improve understanding of depth to bedrock, the lithologic composition of overlying sediments, and as input structure for groundwater model development.

Wisconsin↗

USGS2022

Rapid and gradual permafrost thaw: A tale of two sites

Warming temperatures and increasing disturbance by wildfire and extreme weather events is driving permafrost change across northern latitudes. The state of permafrost varies widely in space and time, depending on landscape, climate, hydrologic, and ecological factors. Despite its importance, few approaches commonly measure and monitor the changes in deep (>1 m) permafrost conditions with high spatial resolution. Here, we use electrical resistivity tomography surveys along two transects in interior Alaska previously disturbed by wildfire and more recently by warming temperatures and extreme precipitation. Long-term point observations of permafrost depth, temperature, and water content inform geophysical measurements which, in turn, are used to extrapolate interpretations over larger areas and with high spatial fidelity. We contrast gradual loss of recently formed permafrost driven by warmer temperatures and increased snowfall, with rapid permafrost loss driven by changes in air temperature, snow depth, and extreme summer precipitation in 2014.

Alaska↗

USGS2022

GSPy: A new toolbox and data standard for Geophysical Datasets

The diversity of geophysical methods and datatypes, as well as the isolated nature of various specialties (e.g., electromagnetic, seismic, potential fields) leads to a profusion of separate data file formats and documentation conventions. This can hinder cooperation and reduce the impact of datasets researchers have invested in heavily to collect and prepare. An open, portable, and well-supported community data standard could greatly improve the interoperability, transferability, and long-term archival of geophysical data. Airborne geophysical methods particularly need an open and accessible data standard, and they exemplify the complexity that is common in geophysical datasets where critical auxiliary information on the survey and system parameters are required to fully utilize and understand the data. Here, we propose a new Geophysical Standard, termed the GS convention, that leverages the well-established and widely used NetCDF file format and builds on the Climate and Forecasts (CF) metadata convention. We also present an accompanying open-source Python package, GSPy, to provide methods and workflows for building the GS-standardized NetCDF files, importing and exporting between common data formats, preparing input files for geophysical inversion software, and visualizing data and inverted models. By using the NetCDF format, handled through the Xarray Python package, and following the CF conventions, we standardize how metadata is recorded and directly stored with the data, from general survey and system information down to specific variable attributes. Utilizing the hierarchical nature of NetCDF, GS-formatted files are organized with a root Survey group that contains global metadata about the geophysical survey. Data are then organized into subgroups beneath Survey and are categorized as Tabular or Raster depending on the geometry and point of origin for the data. Lastly, the standard ensures consistency in constructing and tracking coordinate reference systems, which is vital for accurate portability and analysis. Development and adoption of a NetCDF-based data standard for geophysical surveys can greatly improve how these complex datasets are shared and utilized, making the data more accessible to a broader science community. The architecture of GSPy can be easily transferred to additional geophysical datatypes and methods in future releases.

Frontiers in Earth Science - Environmental Informa↗

USGS2022

Mapped predictions of manganese and arsenic in an alluvial aquifer using boosted regression trees

Manganese (Mn) concentrations and the probability of arsenic (As) exceeding the drinking-water standard of 10 μg/L were predicted in the Mississippi River Valley alluvial aquifer (MRVA) using boosted regression trees (BRT). BRT, a type of ensemble-tree machine-learning model, were created using predictor variables that affect Mn and As distribution in groundwater. These variables included iron (Fe) concentrations and specific conductance predicted from previously developed BRT models, groundwater flux and age estimates from MODFLOW, and hydrologic characteristics. The models also included results from the first airborne geophysical survey conducted in the United States to target an entire aquifer system. Predictions of high Mn and As occurred where Fe was high. Predicted high Mn concentrations were correlated with fraction of young groundwater (less than 65 years) computed from MODFLOW results. High probabilities of As exceedance were predicted where groundwater was relatively old and airborne electromagnetic resistivity was high, typically proximal to streams. Two-variable partial-dependence plots and sensitivity analysis were used to provide insight into the factors controlling Mn and As distribution in groundwater. The maps of predicted Mn concentrations and As exceedance probabilities can be used to identify areas where these constituents may be high, and that could be targeted for further study. This paper shows that incorporation of a selected set of process-informed data, such as MODFLOW results and airborne geophysics, into a machine-learning model improves model interpretability. Incorporation of process-rich information into machine-learning models will likely be useful for addressing a wide range of problems of interest to groundwater hydrologists.

Arkansas, Louisiana, Mississippi, Tennessee↗

USGS2022

Surface parameters and bedrock properties covary across a mountainous watershed: Insights from machine learning and geophysics

Bedrock property quantification is critical for predicting the hydrological response of watersheds to climate disturbances. Estimating bedrock hydraulic properties over watershed scales is inherently difficult, particularly in fracture-dominated regions. Our analysis tests the covariability of above- and belowground features on a watershed scale, by linking borehole geophysical data, near-surface geophysics, and remote sensing data. We use machine learning to quantify the relationships between bedrock geophysical/hydrological properties and geomorphological/vegetation indices and show that machine learning relationships can estimate most of their covariability. Although we can predict the electrical resistivity variation across the watershed, regions of lower variability in the input parameters are shown to provide better estimates, indicating a limitation of commonly applied geomorphological models. Our results emphasize that such an integrated approach can be used to derive detailed bedrock characteristics, allowing for identification of small-scale variations across an entire watershed that may be critical to assess the impact of disturbances on hydrological systems.

Colorado↗