Williston Basin Associated CO 2 Storage Field Laboratory Final Closeout Presentation
Presentation for Final U.S. Department of Energy/NETL Closeout Meeting, virtual, October 2, 2025.
SEARCH · Engineering Papers
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.
Presentation for Final U.S. Department of Energy/NETL Closeout Meeting, virtual, October 2, 2025.
Here, we discuss new quantum technologies being adapted to detect extraordinarily small magnetic fields. The signals are derived from the stimulated emission of light from single-atom defects in diamonds. The intensity of light reports the spin state of the electrons in the color center. The physics is relevant to geochemists because it points to a means of conducting NMR experiments on solutions in diamond-anvil cells. Because the method is so sensitive, it can detect signals in subnanoliter volumes, making work at elevated temperatures and pressures possible, and also on solutions that are dangerous in larger volumes.
Not provided.
Exchanges of groundwater and surface-water are fundamental to a wide range of water-supply and water-quality management issues but challenging to map beyond the reach scale. Waterborne gradient self-potential (SP) measurements are directly sensitive to water flow through riverbed sediments and can be used to infer exchange locations, direction (gain versus loss), scale, and relative changes, but to date applications to river corridor hydrology are limited. Numerical modeling and field experiments were therefore performed herein, each emphasizing waterborne gradient SP logging for identifying and locating focused vertical groundwater discharge (surface-water gain) and recharge (surface-water loss) in a river. Two and three-dimensional numerical models were constructed to simulate the polarities, appearances, and peak amplitudes of streaming-potential and electric-field anomalies on a riverbed and in the surface-water that were attributable to steady-state vertical fluxes of groundwater through high-permeability conduits in the riverbed. Effects of varied hydraulic length-scale of exchange and surface-water depth were tested through numerical modeling. Modeling results aided in data acquisition and interpretation for three repeated field experiments performed along a 1.5–2.0 km reach of the Quashnet River in Cape Cod, Massachusetts, where focused, meter-scale groundwater discharges occur at discrete locations within otherwise ubiquitous and more diffuse groundwater upwelling conditions. Strong gradient SP anomalies were repeatedly measured in the Quashnet River at previously confirmed locations of focused groundwater discharge, showing the efficacy of waterborne gradient SP logging in identifying and characterizing groundwater/surface water exchange dynamics at multiple river network scales.
The Joint Task Force, Science Monitoring And Reliable Telecommunications (JTF SMART) Subsea Cables, is working to integrate environmental sensors for ocean bottom temperature, pressure, and seismic acceleration into submarine telecommunications cables. The purpose of SMART Cables is to support climate and ocean observation, sea level monitoring, observations of Earth structure, and tsunami and earthquake early warning and disaster risk reduction, including hazard quantification. Recent advances include regional SMART pilot systems that are the first steps to trans-ocean and global implementation. Examples of pilots include: InSEA wet demonstration project off Sicily at the European Multidisciplinary Seafloor and water column Observatory Western Ionian Facility; New Caledonia and Vanuatu; French Polynesia Natitua South system connecting Tahiti to Tubaui to the south; Indonesia starting with short pilot systems working toward systems for the Sumatra-Java megathrust zone; and the CAM-2 ring system connecting Lisbon, Azores, and Madeira. This paper describes observing system simulations for these and other regions. Funding reflects a blend of government, development bank, philanthropic foundation, and commercial contributions. In addition to notable scientific and societal benefits, the telecommunications enterprise’s mission of global connectivity will benefit directly, as environmental awareness improves both the integrity of individual cable systems as well as the resilience of the overall global communications network. SMART cables support the outcomes of a predicted, safe, and transparent ocean as envisioned by the UN Decade of Ocean Science for Sustainable Development and the Blue Economy. As a continuation of the OceanObs’19 conference and community white paper (Howe et al., 2019, 10.3389/fmars.2019.00424), an overview of the SMART programme and a description of the status of ongoing projects are given.
The ocean plays an essential role in regulating Earth’s climate, influencing weather conditions, providing sustenance for large populations, moderating anthropogenic climate change, encompassing massive biodiversity, and sustaining the global economy. Human activities are changing the oceans, stressing ocean health, threatening the critical services the ocean provides to society, with significant consequences for human well-being and safety, and economic prosperity. Effective and sustainable monitoring of the physical, biogeochemical state and ecosystem structure of the ocean, to enable climate adaptation, carbon management and sustainable marine resource management is urgently needed. The Argo program, a cornerstone of the Global Ocean Observing System (GOOS), has revolutionized ocean observation by providing real-time, freely accessible global temperature and salinity data of the upper 2,000m of the ocean (Core Argo) using cost-effective simple robotics. For the past 25 years, Argo data have underpinned many ocean, climate and weather forecasting services, playing a fundamental role in safeguarding goods and lives. Argo data have enabled clearer assessments of ocean warming, sea level change and underlying driving processes, as well as scientific breakthroughs while supporting public awareness and education. Building on Argo’s success, OneArgo aims to greatly expand Argo’s capabilities by 2030, expanding to full-ocean depth, collecting biogeochemical parameters, and observing the rapidly changing polar regions. Providing a synergistic subsurface and global extension to several key space-based Earth Observation missions and GOOS components, OneArgo will enable biogeochemical and ecosystem forecasting and new long-term climate predictions for which the deep ocean is a key component. Driving forward a revolution in our understanding of marine ecosystems and the poorly-measured polar and deep oceans, OneArgo will be instrumental to assess sea level change, ocean carbon fluxes, acidification and deoxygenation. Emerging OneArgo applications include new views of ocean mixing, ocean bathymetry and sediment transport, and ecosystem resilience assessment. Implementing OneArgo requires about $100 million annually, a significant increase compared to present Argo funding. OneArgo is a strategic and cost-effective investment which will provide decision-makers, in both government and industry, with the critical knowledge needed to navigate the present and future environmental challenges, and safeguard both the ocean and human wellbeing for generations to come.
In this study, Central Appalachian coal-derived fly ashes from two power plants were beneficiated in a pilot-scale facility in order to produce a product with a relatively consistent concentration of rare earth elements (REE). The <200-mesh final fly ash product was produced by removing the carbon- and Fe-rich particles prior to screening at 200 mesh (75 µm). The Plant D fly ash had high concentrations of CaO and SO3, which were diminished through the two months when the ash was being beneficiated, representing a consequence of the heat, humidity, and excessive rainfall in the Kentucky summer. The high CaO and SO3 concentrations through the early runs likely contributed to the lower REE in the <200-mesh products of those runs. Of the non-REE minor elements, Ba, V, Mn, Zn, and As showed the greatest between-run variations within the runs for each plant. The overall REE concentrations proved to be similar, both on a between-run basis for the individual fly ash sources and on a between-plant basis. Variations in fly ash quality will occur in larger-scale operations, so on-going attention to the fly ash quality and the response of the fly ash to beneficiation is necessary. Changes in the Plant D fly ash with time imply that both the freshness of the original ash and the length and conditions of its storage at the site of beneficiation could be factors in the quality and consistency of the processed fly ash.
Recovery of rare earth elements (REEs) from coal samples collected from the Fire Clay coal seam using diluted mineral acid solutions was investigated. The initial processing step was coal recovery using conventional froth flotation which concentrated the REEs in tailing material resulting in an upgrade to values around 700 ppm on a dry whole mass basis. Leaching experiments were performed on the flotation tailings material using a 1.2 M sulfuric acid solution adjusted to a temperature of 75 °C to study the extractability of REEs from coal material. The effect of particle size, leaching time, leaching temperature, and solid concentration on REE leaching recovery were evaluated. The kinetic data obtained from leaching over a range of temperatures suggested that the leaching process follows the shrinking core model with possibly a mixed control mechanism that may be a result of several heterogenous materials leaching simultaneously. Leaching recovery increased rapidly at the beginning of the reaction then slowed as the system reached equilibrium. The apparent activation energy determined from test data obtained over a range of temperatures using 1 M sulfuric acid was 36 kJ/mol for the first 20 min of reaction time and 27 kJ/mol for the leaching period between 20 and 120 min. The leaching of light REEs during the initial stage was determined to be driven by a chemical reaction, followed by the formation of a product layer, which required lower activation energy in the later stage of leaching. In regards to the heavy REEs, the major mechanism for leaching is desorption and the product layer formation does not affect the heavy REEs significantly.
The increasing industrial demand for rare earths requires new or alternative sources to be found. Within this context, there have been studies validating the technical feasibility of coal and coal byproducts as alternative sources for rare earth elements. Nonetheless, radioactive materials, such as thorium and uranium, are frequently seen in the rare earths’ mineralization, and causes environmental and health concerns. Consequently, there exists an urgent need to remove these radionuclides in order to produce high purity rare earths to diversify the supply chain, as well as maintain an environmentally-favorable extraction process for the surroundings. In this study, an experimental design was generated to examine the effect of zeolite particle size, feed solution pH, zeolite amount, and contact time of solid and aqueous phases on the removal of thorium and uranium from the solution. The best separation performance was achieved using 2.50 g of 12-µm zeolite sample at a pH value of 3 with a contact time of 2 h. Under these conditions, the adsorption recovery of rare earths, thorium, and uranium into the solid phase was found to be 20.43 wt%, 99.20 wt%, and 89.60 wt%, respectively. The Freundlich adsorption isotherm was determined to be the best-fit model, and the adsorption mechanism of rare earths and thorium was identified as multilayer physisorption. Further, the separation efficiency was assessed using the response surface methodology based on the development of a statistically significant model.
The gold contained in copper ores is an important resource for the gold industry. In some cases, elemental gold is present and can be recovered by selective flotation. It has been reported that the gold grade and recovery can be increased, without sacrificing the copper recovery, by replacing AERO 3477 (diisobutyl dithiophosphate (DTP)) with AERO 7249 (mixture of diisobutyl monothiophosphate (MTP) and diisobutyl dithiophosphate (DTP)) as the main collector. The fundamental understanding of the improvement in selectivity with the addition of MTP in the flotation of elemental gold from pyrite is limited and is the subject of this paper. In this regard, the hydrophobicity and selectivity of DTP and MTP in the flotation of pyrite and gold are compared and discussed. Density functional theory (DFT) was used to examine the electron density, reactivity, highest occupied molecular orbital (HOMO), and lowest unoccupied molecular orbital (LUMO) of the MTP and DTP collectors. The interaction energies for the adsorption of MTP and DTP from fresh pyrite, oxidized pyrite and gold surfaces were calculated and discussed with respect to the experimental results reported in the literature. Molecular dynamics simulation (MDS) was used to examine the adsorption state of MTP and DTP on the pyrite (100) and Au (111) surfaces.
Low-temperature plasma (LTP) oxidation has been widely used to study the mineralogy of the mineral matter existing in coal sources. The current study investigated the potential of LTP oxidation as a pre-treatment method to improve rare earth element (REE) leachability from coal and its by-products. Representative density-fractionated samples of Baker and Fire Clay coarse refuse seam materials were ground to a top size of 180 µm and subjected to low-temperature plasma oxidation. Subsequently, the treated samples were leached at 1% w/v solids concentration and 75 °C for 5 h using (i) de-ionized (DI) water, (ii) 0.1 mol/L of ammonium sulfate, and (iii) 1.2 mol/L of sulfuric acid. It was determined that LTP treatment improved REE leaching characteristics, especially the leaching of heavy REEs (HREE), existing in the lighter density fractions of the Baker seam coarse refuse material. For instance, the HREE recovery for the 1.6 specific gravity (SG) float fraction increased from 8% to 33% using 0.1 mol/L of ammonium sulfate solution after 32 h of LTP treatment. This finding indicated that HREEs associated with the organic matter were released by the LTP treatment and adsorbed onto the surfaces of highly negative charged mineral matter and was exchanged with ammonium to allow their recovery. Similarly, when using 1.2 mol/L of sulfuric acid, the HREE recovery increased from 23% to 53% for the 1.6 SG float fraction. Interestingly, LTP oxidation did not provide significant improvement in REE recovery from the 2.2 sink density fractions, which was likely due to its lower organic content. No significant benefits were observed when treating the Fire Clay coarse refuse material, which was likely due to the lack of organic affinity and the difficult-to-leach REE minerals associated with the coal source such as monazite, xenotime, and zircon. Conversely, high-temperature oxidation within a temperature range of 600–750 °C significantly improved REE leaching characteristics for both coal sources. Improvement in REE recovery was due to decarbonization of the material, clay dehydroxylation and subsequent conversion of liberated REE-bearing minerals into a more leachable form. However, increasing the temperature above 800 °C decreased REE recovery due to the conversion of meta-kaolinite into mullite, which is chemically stable.
In a study funded by the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL), the University of North Dakota (UND) Energy and Environmental Research Center (EERC) performed a first-of-its-kind round-robin interlaboratory study (RRIS) to determine lab-to-lab and method-to-method variability in analyzing the rare earth element (REE) content of domestic resources. Analyses of REEs on eleven different materials were accomplished by the laboratories using four different procedures: ASTM D6357 (Procedure A), ASTM D4503 (Procedure B), an alternate in-house procedure for digestion and REE analysis (Procedure C), and neutron activation analysis (NAA). The results of the RRIS suggest that NAA is the most accurate and reliable method for many of the REEs in these types of materials; however, the method is limited in that it is able to determining only ten of the sixteen REEs. Five of the seven labs reporting data for Procedure A and three of the five labs reporting for Procedure B showed excellent performance in terms of repeatability, reproducibility, agreement with NAA, and SRM recoveries based on AOAC International guidelines on method performance. This indicates that when strictly followed these methods are suitable for REE determination in most materials, although are subject to the overall capabilities and experience of individual laboratories.
Rare earth element (REE) permanent magnets (NdFeB) are a critical element in a vast and growing number of industrial applications. In consumer electronics, a broad category encompassing computer, CD, and DVD hard drives, in addition to the ubiquitous cell phones, the nominal NdFeB magnet content may be small, but the global market share for this sector accounts for almost 30% of NdFeB demand, due to a large and continually increasing consumer base. It is estimated that wind turbines that primarily employ permanent magnets will add roughly 110 GW annually of on- and off-shore capability over the next few years. Electric vehicles (EVs) and E-bicycles (EBs) equipped with permanent magnet motors comprise the transportation contribution. Permanent magnet motors have garnered nearly 100% of the market share among EV manufacturers worldwide. Industrial, professional service, and personal robots, most using permanent magnets, are also included in the projected global need for rare earths, particularly Nd and Dy. The sector projects significant growth of approximately 10% across robotic categories. In this paper, we calculate the future demand for Nd and Dy through 2050 across these sectors using a compounded annual growth rate coupled with magnet weight and rare earth content. Uncertainties in the estimates, such as the true global production of Nd, a range of end-product scales and/or unit types in each sector, varied magnet compositions, and the variety of uses within a sector, are all considered.
The D″ region of the lower mantle, which lies just above the core–mantle boundary, is distinct from the bulk of the lower mantle in that it exhibits complex seismic heterogeneity and seismic anisotropy. Seismic anisotropy in this region is likely to be largely due to the deformation-induced texture (crystallographic preferred orientation) development of the constituent mineral phases. Thus, seismic anisotropy can provide a marker for deformation processes occurring in this dynamic region of the Earth. Post-perovskite-structured (Mg,Fe)SiO3 is believed to be the dominant mineral phase in many regions of the D”. As such, understanding deformation mechanisms and texture development in post-perovskite is important for the interpretation of observed seismic anisotropy. Here, we report on high-pressure diamond anvil cell deformation experiments on NaMgF3 neighborite (perovskite structure) and post-perovskite. During deformation, neighborite develops a 100 texture, as has been previously observed, both in NaMgF3 and MgSiO3 perovskite. Upon transformation to the post-perovskite phase, an initial texture of {130} at high angles to compression is observed, indicating that the {100} planes of perovskite become the ~{130} planes of post-perovskite. Further compression results in the development of a shoulder towards (001) in the inverse pole figure. Plasticity modeling using the elasto-viscoplastic self-consistent code shows this texture evolution to be most consistent with deformation on (001)[100] with some contribution of glide on (100)[010] and (001)<110> in NaMgF3 post-perovskite. The transformation and deformation mechanisms observed in this study in the NaMgF3 system are consistent with the behavior generally observed in other perovskite–post-perovskite systems, including the MgSiO3 system. This shows that NaMgF3 is a good analog for the mantle bridgmanite and MgSiO3 post-perovskite.
Surface wettability is determined by the attraction of a liquid phase to a solid surface. It is typically quantified by using contact angle measurements at mineral surfaces in the case of the flotation of mineral particles. Contact angle research to describe wettability has been investigated at different scales by sessile drop measurements, molecular dynamic simulation, and atomic force microscopy. In this study, the density functional theory (DFT) was employed for predicting the surface free energy and contact angles of a well-known hydrophobic phyllosilicate mineral talc and a well-known hydrophilic phyllosilicate mineral muscovite based on the calculated interfacial energy and surface charge. The results revealed that the predicted contact angle at the atomic scale was larger than the experimental value, and identified two interactions: electrostatic interaction and hydrogen bonding, between the hydrophilic muscovite surface and the water layer, while a water-exclusion zone of 3.346 Å was found between the hydrophobic talc surface and the first water layer. This investigation gives a new perspective for wettability determination at the atomic scale.
This paper investigates the feasibility of using a linear current sensing (LCS) technique integrated on an unmanned aerial system (UAS) for detecting and mapping underground infrastructure rapidly and cost-effectively. The LCS technique is based on data from a wide band of electromagnetic induction frequencies (50 kHz to 2 MHz) using a vector magnetic field gradiometer. This technique takes advantage of a slowly decaying secondary magnetic field in order to achieve greater standoff detection distance ($\frac{1}{R^2}$ vs. $\frac{1}{R^6}$ for compact metallic targets during EMI sensing, where R is the distance from a target to the sensor). These secondary magnetic fields are produced by the excite current on long conductors, allowing detection at a distance of 10 meters or more. The system operates between tens of kHz to a few MHz and uses either an active EMI source or existing EM fields to excite this linear current on a long metallic subsurface target. Once excited, these linear currents produce a secondary magnetic field that is detected with an above ground triaxial magnetic field gradiometer. By moving and tracking its geolocation, the system outputs rich datasets sufficient to support high-fidelity forward and inverse EMI models for estimating the depth and orientation of deep underground long linear metallic infrastructure. The system’s hardware and its integration to a UAS system are outlined, along with the formulation of LCS theory, and numerical and experimental data are presented. The results illustrate that the LCS technique offers large standoff detection, is adaptable to UAS, and could be used effectively for detecting deep underground infrastructure such as wires and pipes.
We analyze observations from Robo-AO’s field M dwarf survey taken on the 2.1 m Kitt Peak telescope and perform a multiplicity comparison with Gaia DR2. Through its laser-guided, automated system, the Robo-AO instrument has yielded the largest adaptive optics M dwarf multiplicity survey to date. After developing an interface to visually identify and locate stellar companions, we selected 11 low-significance Robo-AO detections for follow-up on the Keck II telescope using NIRC2. In the Robo-AO survey we find 553 candidate companions within 4″ around 534 stars out of 5566 unique targets, most of which are new discoveries. Using a position cross-match with DR2 on all targets, we assess the binary recoverability of Gaia DR2 and compare the properties of multiples resolved by both Robo-AO and Gaia. The catalog of nearby M dwarf systems and their basic properties presented here can assist other surveys which observe these stars, such as the NASA TESS mission.
Seismic coda-wave analysis is a well-developed method for detecting subtle physical changes in complex media by measuring arrival times in the late-arriving energy from multiply-scattered or reflected waves. However, a challenge arises when multiply-scattered waves are not sufficiently separated in time from the direct arrivals to provide a clear coda wave train. Additional complications for monitoring changes in fracture systems arise when the signals originate from unsynchronized internal sources, such as natural induced seismicity, from acoustic emission, or from transportable intra-fracture sources (chattering dust), that generate uncontrolled signals that vary in arrival time, amplitude and frequency content. Here, we use a twin neural network (also known as a Siamese neural network) for dimensionality reduction to analyze signals from chattering dust to classify the fluid saturation state of a synthetic fracture system. The twin neural network with shared weights generates a low-dimensional representation of the data input by minimizing contrastive loss, serving as the input to a multiclass classifier that accurately classifies whether multiple fractures in a fracture system are fully saturated or partially saturated, or whether a change in saturation has occurred in different fractures in the system. These results show that information buried in unresolved codas from uncontrolled sources can be extracted using machine learning to monitor the evolution of fracture systems caused by physical and chemical processes even when the scattered and direct wave fields overlap.