Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electromagnetic Surveys”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

A ModEx Framework for Watershed Subsurface Investigation With Limited Geophysical Data Using Machine Learning and Hydrologic Modeling

Abstract Subsurface heterogeneity influences watershed hydrology strongly but remains difficult to characterize at catchment scales with sparse and costly field data. Geophysical surveys such as electromagnetic induction (EMI) provide local spatial subsurface images yet scaling them to watershed scales and converting EMI‐derived resistivity into hydraulic properties remains a challenge. We present a Model–Experiment (ModEx) framework that integrates limited EMI data with machine learning (ML) and hydrologic modeling to improve process representation and guide field investigations. Sparse EMI surveys were scaled to the catchment scale using a Random Forest model, and the resulting resistivity fields were combined with nearby borehole constraints to parameterize a hydrologic model. The EMI‐informed hydrological simulations improved predictions of streamflow sustained by subsurface flow and shallow saturation patterns. By combining EMI data and ML with hydrologic modeling, the ModEx framework guides future subsurface surveys, providing a transferable and efficient strategy for data–model integration across diverse watersheds. Plain Language Summary Mapping the underground network of soil and rock that controls water is essential for predicting floods and droughts, but seeing underground is difficult and expensive. We cannot drill everywhere, so scientists use geophysical tools to scan broad areas. There are two key challenges: these geophysical scans are often sparse across the whole watershed, and the geophysical data is hard to translate into water‐related properties. We used artificial intelligence to solve these problems. We taught a computer to find patterns linking the limited geophysical data to the land surface properties. This allowed it to fill in the gaps and create a complete, useful subsurface map for the entire watershed. This new map improves hydrologic simulations, leading to more accurate predictions of water movement in the watershed. It also helps scientists build better models with less data and generates a priority map showing where to measure next, making future investigations more efficient. Key Points Limited EMI scaled with ML improves catchment‐scale subsurface parameterization for hydrologic models The framework integrates hydrologic modeling with limited geophysical data to support subsurface investigation design ModEx framework offers a transferable data–model integration strategy that quantifies and reduces uncertainty guiding watershed studies

Chen, Hang↗

Utilizing 3D Mechanical Earth Models for Calibration and Validation in a Large-Scale Carbon Capture and Storage Project in North Dakota

Conference paper presented at 17th International Conference on Greenhouse Gas Control Technologies (GHGT-17), Calgary, Alberta, Canada, October 20–24, 2024. Three-dimensional (3D) mechanical earth models (MEMs) are pivotal in assessing geological sites for carbon dioxide (CO 2 ) storage and mitigating potential risks associated with storage and injection. The Energy & Environmental Research Center is exploring innovative carbon storage monitoring techniques at a CO 2 storage site near Beulah, North Dakota. These methods aim to provide lower-impact, faster feedback for commercial carbon capture and storage (CCS) projects. The research includes monitoring CO 2 injection using various approaches: 1) an automated, integrated modular monitoring station; 2) time-lapse electromagnetic (EM) field surveys; 3) drone-based surveillance; 4) time-lapse seismic techniques; and 5) advanced wellbore monitoring.

02 PETROLEUM↗

Geothermal Fault Zone and Fluid Imaging through Joint Airborne ZTEM and Ground MT Data Inversion Analysis

This project has aimed to achieve detailed electrical resistivity resolution at geothermal reservoir scales by combining airborne natural electromagnetic (EM) field surveying (ZTEM) with ground magnetotelluric (MT) measurements to approximate an airborne MT geophysical method. MT alone is relatively expensive and may have permitting challenges in sensitive areas. Airborne ZTEM field data contains only the magnetic field, requires a background assumption, and has been limited to relatively high frequencies, thus suffering uniqueness problems. Based on proto-type 2D simulations, ZTEM ambiguities may be reduced through formal incorporation with possibly sparse ground MT soundings, which we pursued in full 3D for this project. The methodology was tested at the high-temperature Roosevelt Hot Springs geothermal system, Utah, which was considered advantageous given the near total exposure of crystalline reservoir rocks across the project area. ZTEM and ground MT survey data were acquired in 2017, subcontracted to outside parties with which we have worked in the past. These included 80 remote-referenced tensor MT soundings over the Mineral Mountains and adjacent Roosevelt Hot Spring producing geothermal system. These MT stations abut later coverage of a similar number of MT stations taken for the Utah FORGE project providing excellent total data aperture to re-solve structure beneath both project areas better than either set alone. The airborne ZTEM survey covered 704 line kilometers in E-W flight lines with a 250 m line spacing. Although this survey was timed during a maintenance-related shutdown of power production at the Roosevelt Hot Springs, other noise sources difficult to identify but including two high-voltage state-scale transmission lines compromised the ZTEM survey badly leading to unusable responses. Thus, with DOE management concurrence, the project proceeded to emphasize inversion and interpretation of the joint SubTER-FORGE MT data sets with regard to the Roosevelt Hot Springs reservoir recharge and to deep heat sources for both it and the Utah FORGE EGS project area. We also investigated the joint ZTEM-MT sampling concept with data sets from the Eleven Mile Canyon prospect area donated by the U.S. Navy (A. Sabin, PoC). Inversion of the SubTER-FORGE MT data using the HexMT 3D finite element algorithm reveals a large, low-resistivity anomaly extending sub-vertically through the depth range of the crust beneath the western Mineral Mountains. The steep conductive zone connects in the lower crust to a more tabular conductor characteristic of much of the Great Basin that generally is ascribed to current mafic magmatic underplating, hybridization and fluid release. The location of the resolved anomaly relative to the recent (0.5-0.8 Ma) eruptive centers of the Mineral Mountains implicates it as remnants of the magma body which fed these centers. This structure appears to be currently feeding heat and fluids upward into the Roosevelt Hot Springs hydrothermal system, as well as heat laterally to the FORGE project area. Separate and joint inversion models were carried out for the donated Eleven Mile Canyon MT-ZTEM data set to demonstrate concept. ZTEM only inversion showed two main alteration zones in the western portion of the project area known from geological mapping. Joint inversion including an E-W profile of MT soundings sharpened these features considerably. It also resolved in much greater detail the graben related normal faulting structure of the central project area which lies at depths exceeding the sensitivity of ZTEM alone. The sparse number of MT da-ta relative to the ZTEM required upweighting the former by a factor of several, but an exact procedure awaits future research. Our final impression is that sparse MT data can improve resolution of the subsurface over that of ZTEM alone. However, well sampled MT data are to be preferred and offer the simplicity of interpreting just one data type, and possess the superior resolution capability coming with the electric field everywhere, and from their high bandwidth.

15 GEOTHERMAL ENERGY↗

Findings on subtask 1.6 – basin electric carbon storage research project: novel monitoring techniques

The Energy & Environmental Research Center (EERC) led a study to validate novel and emerging technologies as commercial monitoring techniques for application in carbon dioxide (CO 2 ) injection operations. This applied research was conducted at Basin Electric Power Cooperative’s (Basin Electric’s) active CO 2 -injection operations in Beulah, North Dakota, in two phases. The EERC previously completed a set of baseline (preinjection) activities in Phase 1, which included 1) design of an automated, integrated, modular (AIM) monitoring station; 2) time-lapse electromagnetic (EM) field surveys; 3) drone-based surveillance studies; 4) time-lapse monitoring with seismic methods; and 5) advanced wellbore-monitoring methods.

54 ENVIRONMENTAL SCIENCES↗

SUBTASK 1.6 – BASIN ELECTRIC CARBON STORAGE RESEARCH PROJECT: NOVEL MONITORING TECHNIQUES

The Energy & Environmental Research Center (EERC) conducted baseline activities associated with an applied research project at Basin Electric Power Cooperative’s (Basin’s) carbon capture and storage (CCS) site in Beulah, North Dakota, to establish novel carbon storage-monitoring techniques as commercial methods under Cooperative Agreement No. DE-FE0024233, Subtask 1.6. The following report summarizes the baseline activities performed and briefly describes the subsequent (operational monitoring) activities that have been proposed to the U.S. Department of Energy (DOE) as part of the overall project to develop and demonstrate novel monitoring techniques at North America’s largest permitted CCS operation. Dakota Gasification Company (DGC), a wholly owned subsidiary of Basin, owns and operates the Great Plains Synfuels Plant (GPSP) approximately 5 miles northwest of the town of Beulah, North Dakota (Figure 1). In 2023, DGC received approval from the North Dakota Industrial Commission (NDIC) to develop a storage facility on-site for injecting a stream of carbon dioxide (CO2) captured from GPSP. DGC will transport the captured CO2 stream with approximately 6.8 miles of transmission lines that extend north of GPSP and inject >1 million tonnes (MMt) of CO2 annually (>1 MMt/yr) over a 12-year period with up to six underground injection control (UIC) Class VI-compliant injection wells completed in the Broom Creek Formation, a predominantly sandstone reservoir and saline aquifer underlying GPSP. The Broom Creek Formation lies approximately 5900 feet (ft) below ground surface (bgs) at GPSP. The commercial scale (i.e., >1 MMt/yr) of DGC’s permitted carbon storage project is ideal for developing and testing the novel monitoring techniques included within Subtask 1.6. The goals of this project are to demonstrate 1) the cost-effectiveness of novel monitoring technologies included as part of this research, 2) technology capability for tracking the CO2 plume and/or associated pressure response in the subsurface and monitoring out-of-zone migration, and 3) compliance with UIC Class VI program requirements. The research activities proposed for the overall project include 1) design of an automated, integrated, modular (AIM) monitoring station; 2) time-lapse electromagnetic (EM) field surveys; 3) drone-based surveillance studies; 4) time-lapse monitoring with seismic methods; 5) advanced wellbore-monitoring methods; 6) deployment of an AIM monitoring network; 7) EM monitoring of CO2 with real-time data processing; 8) continued seasonal drone-based surveillance studies; 9) seismic monitoring with passive and active surveys; and 10) wellbore monitoring with nuclear magnetic resonance (NMR) for near-surface characterization. Completion of Activities 1.0–5.0 (baseline activities) are described in this report. Upon authorization of funding by DOE, the EERC will initiate Activities 6.0– 10.0 (operational monitoring activities). Current state-of-the-art (SOA) carbon storage-monitoring techniques require countless labor hours dedicated to the acquisition of data. Once data are gathered, these SOA techniques often rely on commercial facilities to process raw data from the field. However, it is anticipated that next-generation monitoring techniques, such as those being demonstrated, will lower acquisition footprints, be less operationally intensive, and improve data acquisition efficiencies. These new techniques are more conducive to the application of machine learning, artificial intelligence, and automation, thus providing a pathway for integration into active control systems, informing site operability, and improving the integration of data for future CCS projects across the United States. Additionally, reclaimed and active mining lands are present within the project site, creating a unique opportunity to demonstrate the effectiveness of remote sensing and surface-based geophysics monitoring techniques at similar project sites that may include disturbed, unconsolidated, or actively excavated near-surface environments. The efforts included in the overall project will produce necessary designs, learnings, and data acquired during the baseline and operational monitoring periods that are necessary for time-lapse demonstration and validation of the described monitoring techniques. In addition, it is anticipated that the monitoring technologies included in this study will be compliant with UIC Class VI requirements to enable the potential for implementation at other CCS sites across the United States.

42 ENGINEERING↗

IceCube-Gen2: the window to the extreme Universe

The observation of electromagnetic radiation from radio to γ-ray wavelengths has provided a wealth of information about the Universe. However, at PeV (10 15 eV) energies and above, most of the Universe is impenetrable to photons. New messengers, namely cosmic neutrinos, are needed to explore the most extreme environments of the Universe where black holes, neutron stars, and stellar explosions transform gravitational energy into non-thermal cosmic rays. These energetic particles havemillions of times higher energies than those produced in the most powerful particle accelerators on Earth. As neutrinos can escape from regions otherwise opaque to radiation, they allow an unique view deep into exploding stars and the vicinity of the event horizons of black holes. The discovery of cosmic neutrinos with IceCube has opened this new window on the Universe. IceCube has been successful in finding first evidence for cosmic particle acceleration in the jet of an active galactic nucleus. Yet, ultimately, its sensitivity is too limited to detect even the brightest neutrino sources with high significance, or to detect populations of less luminous sources. In thiswhite paper, we present an overview of a next-generation instrument, IceCube-Gen2, which will sharpen our understanding of the processes and environments that govern the Universe at the highest energies. IceCube-Gen2 is designed to: (a) Resolve the high-energy neutrino sky from TeV to EeV energies (b) Investigate cosmic particle acceleration through multi-messenger observations (c) Reveal the sources and propagation of the highest energy particles in the Universe (d) Probe fundamental physics with high-energy neutrinos IceCube-Gen2 will enhance the existing IceCube detector at the South Pole. It will increase the annual rate of observed cosmic neutrinos by a factor of ten compared to IceCube, and will be able to detect sources five times fainter than its predecessor. Furthermore, through the addition of a radio array, IceCube- Gen2 will extend the energy range by several orders of magnitude compared to IceCube. Construction will take 8 years and cost about $350M. The goal is to have IceCube-Gen2 fully operational by 2033. IceCube-Gen2 will play an essential role in shaping the new era of multimessenger astronomy, fundamentally advancing our knowledge of the highenergy Universe. Overall, this challenging mission can be fully addressed only through the combination of the information from the neutrino, electromagnetic, and gravitational wave emission of high-energy sources, in concert with the new survey instruments across the electromagnetic spectrum and gravitational wave detectors which will be available in the coming years.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

2D Seismic Surveys as an Alternative for CO 2 Monitoring in the North Dakota CarbonSAFE Project

Conference presentation for International Meeting for Applied Geoscience & Energy (IMAGE), Houston, TX, August 28 – September 1, 2023. As part of the North Dakota CarbonSAFE project, seismic, controlled-source electromagnetic, and microgravity baseline monitoring surveys were collected to assess the safe, permanent, commercial-scale geologic storage of CO 2 generated by the Milton R. Young coal-fired power plant. Data from five 2D seismic lines were assessed as an alternative to 3D seismic for CO 2 monitoring. 2D and 3D surveys were compared using wave equation-based (WEB)-AVO inversion. The encouraging WEB-AVO results demonstrated that 2D surveys could be included in CO 2 monitoring programs as a cost-effective method.

20 FOSSIL-FUELED POWER PLANTS↗

Phase II Field Demonstration at Lansing Smith Generating Plant, Southport, Florida (Final Report)

The Final Technical Report: Field Test Design and Pressure Management Strategies for Phase II Field Demonstration of Optimal Pressure Control, Plume Management, and Produced Water Strategies presents the culmination of multi-year efforts under the U.S. Department of Energy’s Brine Extraction and Storage Test (BEST) program, led by EPRI in partnership with Gulf Power at the Plant Smith site near Panama City, Florida. The project was designed to evaluate and demonstrate the technical feasibility of managing subsurface pressures and fluid movement associated with large-scale CO₂ injection, using low-salinity water as a proxy fluid. Through a combination of field injection testing, reservoir modeling, and optimization studies, the research team developed and refined active and passive brine extraction strategies aimed at controlling injection-induced pressure buildup, mitigating risks of fault activation, and managing plume migration. The field demonstration incorporated a new injection well (TIW-2), a new monitoring/extraction well (TEMW-A), and an existing well (TIW-1) repurposed for passive pressure relief. Complementary geophysical monitoring was designed to track plume development and evaluate the effectiveness of pressure management strategies. The report details the integrated workflow encompassing wellfield development, regulatory permitting, model calibration, and survey design. It includes a comprehensive electromagnetic (EM) modeling and inversion study used to develop a cost-effective, time-lapse geophysical monitoring plan capable of imaging the evolving low-salinity plume within the highly saline Lower Tuscaloosa formation. Reservoir simulation results guided the design of the 17-month injection program and the timing of active extraction to maintain formation pressures below a hypothetical fault reactivation threshold. Supporting analyses evaluated potential injectivity risks related to clay fines migration and geochemical incompatibility, identifying practical mitigation measures such as salinity and pH control. The final design recommends a series of crosswell EM surveys—conducted before, during, and after injection—as the most effective and economical approach for plume imaging, supplemented by continuous downhole pressure and flow monitoring. Collectively, the study provides a field-ready framework for cost-effective pressure management and monitoring in support of future CO₂ storage projects.

01 COAL, LIGNITE, AND PEAT↗

Electromagnetic inversion for a baseline conductivity model at the Kemper CarbonSAFE site

Geophysics is one component of a successful monitoring strategy for carbon storage sites. Seismic surveys are a common element of these strategies, with repeat surveys occurring about every 5 years. However, seismic is expensive, the sensitivity to CO2 saturation is limited, and it requires a large footprint on the surface (Yurikov et al., 2022). Electromagnetic (EM) methods offer repeat surveys at a significantly lower cost, the ability to utilize airborne methods to increase coverage, and potentially provide enough information to increase the time between seismic surveys. Specifically, we focus on semi-airborne EM, with a grounded dipole transmitter, and airborne magnetic field receivers as described and modeled in Kohnke et al. (2023).

Kohnke, Colton↗

Fiber Optic Distributed Acoustic Sensing in DFPS East Well

Fiber optic files in *.hdf format, acquired in the Devine Fracture Pilot Site East well, using Fotech interrogator. There was no change in channel spacing or acquisition settings during the project. Channel spacing is 0.68 m.

Electrically active proppant,Electromagnetic,Fiber↗

Enhanced Multi-Dimensional Inversion Through Target-Specific Inversion Parameter Bounds With an Application to Crosswell Electromagnetic for Sequestration Monitoring

In geophysical inversions, lower and upper model parameter bounds are a means of solution stabilization. Further, constraints that intend to let only geologically plausible inverse solutions pass are amenable to lower and upper bounds. Reliable prior information is paramount to construct such bound constraints. It is common practice to narrow and widen bound intervals for regions of, respectively, more and less certain prior information. Contrary to this practice, we experiment with widened bound intervals in zones that are poorly resolved by a given survey configuration but where prior information would suggest structural anomalies of interest. The purpose of enlarged parameter bounds that correlate spatially with predefined targets is to let the inversion explore a larger solution space, thus increasing the potential to resolve otherwise hidden anomalies. Application of the method is based on a carbon-sequestration baseline (pre-injection) crosswell electromagnetic (EM) field survey at the Containment and Monitoring Institute Field Research Station (Alberta, Canada), where impeded measurements led to generally reduced sensitivities for the interwell region. Synthetic-data proofs of concept use augmented bounds designed to boost the resolution of artificial plume targets, indicating an enhanced illumination compared to constant bounds. Comparative field data inversions with spatially variable bounds constructed from prior resistivity and velocity information highlight non-horizontal baseline structures.

3-D electrical resistivity imaging↗

The Kimberlina synthetic multiphysics dataset for CO 2 monitoring investigations

Abstract We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO 2 reservoir model based on a potential CO 2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO 2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO 2 , and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO 2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO 2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.

58 GEOSCIENCES↗

El-CID: a filter for gravitational-wave electromagnetic counterpart identification

ABSTRACT As gravitational-wave (GW) interferometers become more sensitive and probe ever more distant reaches, the number of detected binary neutron star mergers will increase. However, detecting more events farther away with GWs does not guarantee corresponding increase in the number of electromagnetic counterparts of these events. Current and upcoming wide-field surveys that participate in GW follow-up operations will have to contend with distinguishing the kilonova (KN) from the ever increasing number of transients they detect, many of which will be consistent with the GW sky-localization. We have developed a novel tool based on a temporal convolutional neural network architecture, trained on sparse early-time photometry and contextual information for Electromagnetic Counterpart Identification (El-CID). The overarching goal for El-CID is to slice through list of new transient candidates that are consistent with the GW sky localization, and determine which sources are consistent with KNe, allowing limited target-of-opportunity resources to be used judiciously. In addition to verifying the performance of our algorithm on an extensive testing sample, we validate it on AT2017gfo – the only EM counterpart of a binary neutron star merger discovered to date – and AT2019npv – a supernova that was initially suspected as a counterpart of the GW event, GW190814, but was later ruled out after further analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Magnetotelluric Investigations of the Kilauea Volcano, Hawaii

In 2002 and 2003 a collaborative effort was undertaken between Lawrence Berkeley National Laboratory, Sandia National Laboratories, the U.S. Geological Survey (USGS) Menlo Park, the USGS Hawaiian Volcano Observatory, and Electromagnetic Instruments Inc. to study the Kilauea volcano in Hawaii using the magnetotelluric (MT) technique. The work was motivated by a desire to improve understanding of the magma reservoirs and conduits within Kilauea and the East and Southwest Rift zones, which has implications for understanding Kilauea's plumbing system. An improved understanding of the rift zones has implications in understanding large-scale landslides that are generated in the Hilina Slump, which produce significant impacts on coastal communities. Up to eight stations operated simultaneously, with multiple remote reference sites, and data were processed using multi-station robust processing techniques. Data was acquired at 70 sites over the Southwest and East rift zones. Good to excellent quality data were obtained even in the harshest conditions, such as those encountered on the fresh lava flows of the East Rift Zone, where electrical contact resistances are on the order of 100 kΩ. A three-dimensional (3D) MT model study was done to guide interpretation of the observed MT measurements. Synthetic modeling demonstrates that conductive bodies in the upper 3 km can be spatially resolved where MT station sampling is good. Resistivity anomalies in the 3D inversions have a high degree of spatial correlation with previously published seismic velocity anomalies beneath Kilauea. Melt fractions between 0.096 and 0.117 are calculated for the Kilauea and Pu'u'ō'ō low resistivity anomalies, respectively.

58 GEOSCIENCES↗

Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”

This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.

ATS↗