Search NASA⌕ Search

SEARCH · Search NASA

Results for “FATES”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A geospatial model for the analysis of time-dependent land subsidence induced by reservoir depletion

Land subsidence due to fluid depletion is an outcome of physical processes operating across a wide range of time and length scales. Although geomechanical models are crucial to simulate reservoir compaction and predict its long-term fate, their use across large regions often bears prohibitive computational costs. To overcome this obstacle, this paper proposes a simplified modelling framework consisting of (i) a near-field numerical solver simulating the coupling between fluid flow and rock deformation with reference to a simplified one-dimensional geometry and (ii) a far-field geospatial algorithm mapping ground settlements across a region through the superposition of poroelastic computations at multiple wells. The model computes the delay between depletion history and reservoir compaction in proximity of a producing well by assuming basal depletion of a fluid-saturated deformable disk, while the Geertsma solution of nucleus of strain is used to extrapolate the impact of such time-varying reservoir compaction around the well. This approach has been used to back-analyze the spatio-temporal progression of subsidence at the Groningen gas field. The results are compared against measurements collected over 50 year of production at 25 benchmark locations scattered over an area of 900 km 2 . It is shown that coupled simulations based on average values of rock compressibility and permeability lead to nonlinear trends of subsidence evolution in good agreement with field measurements, while uncoupled analyses overpredict settlements by more than 70%. Lastly, synthetic forecasts based on different rates of depletion were provided. The results suggest that slower depletion rates lead to lower subsidence at a given time, and that residual subsidence may continue to develop for several decades after interruption of production activities.

58 GEOSCIENCES↗
DOE OSTI2022

Sulfides mediate the migration of nanoparticle mass out of nanocomposite plastics and into aqueous environments

We show that inorganic sulfides strongly influence transfer (migration) of nanoparticle mass out of polymer nanocomposites (PNCs) and into aqueous environments. We first manufactured two families of PNCs: one incorporating silver nanoparticles (AgNPs) and one incorporating CdSe quantum dots (QDs). Then, we assessed migration out of these PNCs and into aqueous media containing Na 2 S at concentrations ranging from 0 to 10 -4 M. Results show that Na 2 S strongly suppressed migration of Ag from AgNP-based PNCs: the migration into water spiked with 10 -6 M Na 2 S was 79% less than migration into water without Na 2 S, and no migration was detected (LOD ≈ 0.01 ng/cm 2 ) in water spiked with Na 2 S at 10 -5 M or 10 -4 M. With CdSe QD-based PNCs, Na 2 S suppressed Cd migration but enhanced Se migration, resulting in only a small net effect on the total QD migration but a large shift of the leachate composition (from favoring Cd by an average of 5.8 to 1 in pure water to favoring Se 9.4 to 1 when Na 2 S was present at 10 -4 M). These results show that common inorganic substances like sulfides may play a strong role in determining the environmental fate of polymer-dispersed nanoparticles and imply that migration tests conducted in purified water may not always accurately reflect migration into real environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Efficient bi-directional coupling of 3D computational fluid-particle dynamics and 1D Multiple Path Particle Dosimetry lung models for multiscale modeling of aerosol dosimetry

The development of predictive aerosol dosimetry models has been a major focus of environmental toxicology and pharmaceutical health research for decades. One-dimensional (1D) models successfully predict overall deposition averages but fail to accurately predict local deposition. Computational fluid-particle dynamics (CFPD) models provide site-specific predictions but at a computational cost that prohibits whole lung predictions. Thus, there is a need for developing multiscale strategies to provide a realistic subject-specific picture of the fate of inhaled aerosol in the lungs. CT-based 3D/CFPD models of the large airways were bidirectionally coupled with individualized 1D Navier-Stokes airflow and particle transport based upon the widely used Multiple Path Particle Dosimetry Model (MPPD). Distribution of airflows among lobes was adjusted by measured lobar volume changes observed in CT images between FRC and FRC + 1.5 L. Additionally, as a test of the effectiveness of the coupling procedures, deposition modeling of previous 1 µm aerosol exposure studies was performed. The complete coupled model was run for 3 breaths, with the computation-intense portion being the 3D CFD Lagrangian particle tracking calculation. The average deposition per breath was 11% in the combined multiscale model with site-specific doses available in the CFPD portion of the model and airway- or region-specific deposition available for the MPPD portion. In conclusion, the key methods developed in this study enable predictions of ventilation heterogeneities and aerosol deposition across the lungs that are not captured by 3D or 1D models alone. Overall, these methods can be used as the foundation for multi-scale modeling of the full respiratory system.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

The influence of cooling rate on condensation of iron, aluminum, and uranium oxide nanoparticles

Fundamental observations of particle size distributions are needed to develop models that predict the fate and transport of radioactive materials in the atmosphere following a nuclear incident. The extent of material transport is influenced by the time scales of particle formation processes (e.g., condensation, coagulation). In this study, we investigated the influence of cooling time scales on size distributions of uranium, aluminum, and iron oxide particles that are synthesized separately under identical run conditions inside the controlled environment of an argon plasma flow reactor. Two distinct temperature distributions are imposed along the flow reactor by varying the argon flow rate downstream of the plasma torch. The vaporized reactants of uranium, aluminum, and iron are cooled from about 5000K to 1000K before they are collected on silicon wafers for ex situ scanning electron microscope analysis. The microscope images show that the sizes of the largest aluminum and iron oxide particles heavily depend on the cooling time scales, whereas significant size variation with cooling rate is not observed for uranium oxide particles. In addition, the size distribution of aluminum oxide particles exhibits the broadest range among all three metal oxides studied. We performed simulations of particle size distributions using a kinetic model that couples gas phase oxidation chemistry with particle formation processes, including nucleation, condensation, and coagulation. The model results demonstrate the strong sensitivity of particle size distribution to different cooling histories (i.e., temperature vs residence time) along the flow reactor. In conclusion, the kinetic model also helps identify directions for future research to improve the predictions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Hygroscopic growth of UO 2 F 2 nanoparticles

Hygroscopicity is an important physicochemical property of aerosol that describes the ability of a particle to uptake water. The hygroscopic properties of uranyl fluoride (UO 2 F 2 ) aerosol generated from a UF 6 hydrolysis reactor was investigated for the first time using a custom-built Humidified Tandem Differential Mobility Analyzer (HTDMA). The HTDMA is capable of measuring UO 2 F 2 nanoparticle growth determined by mobility size over a wide range of atmospheric humidity from dry conditions at <10% relative humidity (RH) to 85% RH. The hygroscopic properties were determined for nanoparticles as small as 3.5 nm in this study. Although the largest size of UO 2 F 2 nanoparticles was 80 nm, monodisperse aerosol with a mobility diameter of up to approximately 500 nm can be investigated using the HTDMA. Anhydrous UO 2 F 2 nanoparticles with a mobility diameter of 3.5 nm were shown to be highly hygroscopic with a deliquescence relative humidity (DRH) of 10%. Hydrates with a larger mobility diameter from 10 to 80 nm were non-hygroscopic with no observable DRH and limited water uptake up to 85% RH. Here, these results demonstrate the hygroscopic properties of UO 2 F 2 nanoparticles are highly variable and based on both the mobility size and hydration state. Hygroscopicity affects the physicochemical properties of UO 2 F 2 nanoparticles, including the aerosol phase state and viscosity, with impacts on aerosol growth, coagulation, and deposition that is critical for understanding the fate and transport of UO 2 F 2 nanoparticles in the atmosphere.

Hygroscopicity↗
DOE OSTI2021

Stratigraphic Identification with Airborne Electromagnetic Methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell’s equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Life cycle analysis of gasification and Fischer-Tropsch conversion of municipal solid waste for transportation fuel production

Non-recyclable municipal solid waste (MSW) can be used as feedstock for liquid fuel production via gasification followed by Fischer-Tropsch (FT) processes. Given the heterogeneity of MSW material composition and variation in material properties, its convertibility to liquid hydrocarbon fuels could vary widely, affecting the sustainability of utilizing non-recyclable MSW for fuel production. This study evaluates the life cycle greenhouse gas (GHG) emissions (carbon intensities [CIs]) of FT fuels from non-recyclable MSW. Key issues that could greatly affect the CIs were examined, including fossil carbon content of the MSW, emission implications of diverting non-recyclable MSW from landfills to fuel production, and conversion efficiency. Results show that the CIs of fuels produced from various waste streams range 80–105 gCO2e/MJ, which may exceed the CI of petroleum fuels. Higher fossil carbon content in the MSW feedstock tends to incur higher GHG emissions as biogenic carbon emissions are considered carbon neutral. Meanwhile, diverting different fractions of non-recyclable MSW, such as food waste and low-quality paper, from landfills may result in GHG emissions that may include the potential avoidance of methane emissions and potential sequestration of biogenic carbon that is foregone. To reduce GHG emissions, a carbon capture and sequestration option in the fuel production stage is considered, which could reduce the CI by 53–64 gCO2e/MJ. Carbon fates of different non-recyclable MSW in landfills are further evaluated to determine how they vary and impact the CIs of MSW-derived fuels.

02 PETROLEUM↗
DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

The adsorption of As(V) on poorly crystalline Fe oxyhydroxides, revisited: Effect of the reaction media and the drying treatment

Arsenic (As) adsorbed on Fe oxyhydroxides (adsorbent) is widely occurring in many environmental settings such as in acid mine drainage systems or in the hydrometallurgical operations to form Fe-As coprecipitates. However, the influence of the reaction media and the drying treatment on the microstructure of the directly formed adsorbents at various pHs was still not fully understood. In this work, As adsorption behaviors on various forms of Fe oxyhydroxides were systematically investigated by using XRD, FTIR, Raman, XANES, and HRTEM. Overall, the results revealed that at weak acidic pH, more As could adsorbed on the suspension adsorbent formed in sulfate and chloride media than that in nitrate media, possibly due to the microstructure alteration of the adsorbent in the presence of sulfate and chloride. Besides, the increasing crystallinity of the Fe oxyhydroxides and the aggregation effect after drying were the major reasons why less As could be hold by the dried adsorbents than that of the corresponding suspension adsorbents. These findings could shed more light on the nature of the Fe oxyhydroxides which may be helpful for more precisely predicting the fate of some toxic metal(loid)s in the environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Spatial and temporal baseline of perfluorooctanesulfonic acid retained in sediment core samples from Puget Sound, Washington, USA

Per- and poly-fluorinated alkyl substances (PFAS) are a group of highly persistent synthetic chemicals utilized in many industrial and consumer products, and — most significantly toward introduction to the marine environment — in fire-fighting foams. Recently, PFAS have been linked to adverse health effects, prompting the need to understand transport and fate. Puget Sound’s main basin, hosting densely populated urban and industrial areas, is known to be contaminated with various persistent organic pollutants and provides a central focus for ecosystem monitoring. Perfluorooctanesulfonic acid, a manufactured PFAS and biodegradation product, partitions to marine sediments and thus can be used as a primary indicator toward regulatory efforts. The current study offers a spatial and temporal analysis of the Puget Sound from cores collected adjacent, Tacoma and Seattle, WA, as well as cores from central Hood Canal and Carr Inlet.

Per- and poly-fluorinated alkyl substances (PFAS),↗
DOE OSTI2021

Monte Carlo simulation of the passage of $\gamma$-rays and $\alpha$-particles in CsI

In this work, theoretical and computational methods for simulating the creation of ionization tracks by fast ions in solids were applied to the passage of a-particles in CsI, an inorganic scintillator commonly used for radiation detection. The methods were implemented in a Monte Carlo program to simulate the interaction of $\alpha$-particles, with incident energies of up to 1 MeV, with CsI. The simulations followed the fate of individual electron-hole pairs and thus allowed for a detailed description of the microscopic structure of ionization tracks created by incident radiation. Simulations were also performed with $\gamma$-rays of the same energy to compare and contrast the ionization tracks obtained with both types of particle. Intrinsic properties such as the mean energy per electron-hole pair, Fano factor, maximum theoretical light yield, and spatial distributions of electron-hole pairs were computed for both $\alpha$-particles and $\gamma$-rays. $\alpha$-particles created cylindrical tracks that were initially aligned with the incident direction and with initial radii of a few nanometers, whereas $\gamma$-rays showed significant scattering, resulting in probability distributions with lower intensities and much greater radial extents.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Sequentially coupled flow and geomechanical simulation with a discrete fracture model for analyzing fracturing fluid recovery and distribution in fractured ultra-low permeability gas reservoirs

More accurate characterization and prediction of the in-situ distribution of fracturing fluid in fractured reservoirs are needed for enhancing well productivity. In this study, an implicit-sequentially coupled flow/geomechanics simulator incorporating an efficient discrete fracture model is developed to model fluid distribution and recovery performance of ultra-low permeability gas reservoirs. The finite-volume and finite-element methods are used for space discretization of the flow and geomechanics equations, respectively, while the backward Euler method is employed for time discretization. The flow and geomechanics equations are solved sequentially based on fixed-stress splitting. An efficient discrete-fracture model is used to explicitly model the fractured system. Flexible unstructured gridding is employed to model arbitrarily-oriented fractures. The interrelations among pore volume, permeability and geomechanical conditions are considered dynamically using two-way coupled flow and geomechanics computations. The geometry of fracture (networks) due to hydraulic fracturing has significant impacts on the fracturing fluid recovery efficiency and ensuing fluid distribution. Under the same injection volume, the fracturing fluid recovery is higher when the fracture geometry is planar. Fluid recovery is relatively lower whenever natural fractures are activated during fracturing treatments; flowback time is also shortened when complex fracture network with enlarged fracture interface is present. Fracturing fluid in hydraulic fractures may leak off into the natural fractures and subsequently imbibes into the surrounding matrix due to capillarity effects. The fracturing fluid recovery and in-situ fluid distribution are sensitive to the shut-in duration and fracture closure behavior. This study analyzes the coupled flow-geomechanical responses of fractured gas reservoirs during the post-fracturing periods. Understanding the fate of the fracturing fluid can provide insights on, to some extent, the stimulated fracture volume, size of the water invasion zone, and efficiency of the fracturing design. The simulation predictions can also provide more accurate initial reservoir conditions (e.g. distributions of different phases and pressure) for long-term well performance estimation.

42 ENGINEERING↗
DOE OSTI2024

Probing O 2 -dependence of cyclopentyl reactions via isomer-resolved speciation

Modeling chemical kinetics relevant to low-temperature combustion requires complete description of reactions involving critical species such as hydroperoxyalkyl radicals, Q̇OOH, which undergo competing unimolecular reactions and bimolecular reactions with O 2 . The balance of flux across the two pathways affects rates of chain-branching and depends on temperature, pressure, and oxygen concentration. Accordingly, the influence of [O 2 ] on product formation from alkyl + O 2 reactions and the subsequent fate of Q̇OOH and related products is central to the development of an accurate chemical kinetics mechanism. Furthermore, chemical reactions consuming Q̇OOH-mediated species are often simplified to such a degree that mechanism truncation error (uncertainty derived from incomplete reaction networks) becomes significant.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Chemistry and P-V-T equation of state of FeO 2 H x at the base of Earth’s lower mantle and their geophysical implications

The hydrogen-absorbing ability of a mantle mineral in its structure determines the capacity of the water reservoirs hosted by the mineral. Water reservoirs at the base of Earth’s mantle directly influence the fate of water brought down by slab subduction and the seismic heterogeneity such as ultralow-velocity zones (ULVZs) at the core-mantle boundary. Pyrite-FeO 2 H x (0 ≤ x ≤ 1) presents a possibility of such reservoirs in the deep mantle. Ever since the discovery of this mineral phase, however, its chemistry at the lower mantle conditions has been debated. We conducted kinetics experiments of pyrite-FeO 2 H x dehydrogenation at 110 GPa/2100 K, 110 GPa/2300 K, and 120 GPa/2300 K and P-V-T equation of state analysis using in situ synchrotron X-ray diffraction. We found that x approaches 0.80, 0.75, and 0.79, respectively, at the above conditions. The collective P-V-T data yield K 0 = 241(13) GPa, K' = 4.2(4), dK/dT = –0.028(1) GPa/K, α 0 = 4.32(13) × 10 –5 K –1 , and α 1 = 0.31(10) × 10 –8 K –2 for the composition of x = 0.75 ± 0.04. Our first-principles calculations indicate that FeO 2 H 0.75 with a slightly distorted pyrite structure is stable at 100 GPa. These results indicate that this mineral is likely present in the deep mantle with rather a partially dehydrogenated composition than FeO 2 or FeOOH. Furthermore, the results also clarify the difference between the ULVZs originated from pyrite-FeO 2 H x and those from partial melting in terms of shear and compressional wave seismic velocity reduction ratio δlnV S /δlnV P .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Organic compounds alter the preference and rates of heavy metal adsorption on ferrihydrite

The availability of heavy metals in terrestrial environments is largely controlled by their interactions with minerals and organic matter, with iron minerals having a particularly strong role in heavy metal fate. Because soil organic matter contains a variety of compounds that differ in their chemical properties, the underlying impact organic matter-soil mineral associations bestow on heavy metal binding is still unresolved. Here, we systematically examine the binding of Cd, Zn and Ni by a suite of organic-ferrihydrite assemblages, chosen to account for various compound chemistries within soil organic matter. We posited that organic compound functionality would dictate the extent of association with the organic-ferrihydrite assemblages. Increased heavy metal binding to the assemblages was observed and attributed to the introduction of additional binding sites by the organic functional groups with differing metal affinities. Furthermore, the relative increase depended on the metal's Lewis acidity and followed the order Cd > Zn > Ni, whereas the reverse order was obtained for metal binding by pristine ferrihydrite (Ni > Zn > Cd). Citric acid-, aspartic acid- and cysteine-ferrihydrite assemblages also enhanced the metal binding rate. X-ray absorption spectroscopy revealed that the organic coating contributed significantly to Zn binding by the assemblages, despite relatively low organic surface coverage. Our findings provide valuable information on the nature of heavy metal-organic-mineral interactions and metal adsorption processes regulating their bioavailability and transport.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Natural organic matter composition and nanomaterial surface coating determine the nature of platinum nanomaterial-natural organic matter corona

Natural organic matter corona (NOM corona) is an interfacial area between nanomaterials (NMs) and the sur- rounding environment, which gives rise to NMs' unique surface identity. While the importance of the formation of natural organic matter (NOM) corona on engineered nanomaterials (NMs) to NM behavior, fate, and toxicity has been well-established, the understanding of how NOM molecular properties affect NOM corona composition remains elusive due to the complexity and heterogeneity of NOM. This is further complicated by the variation of NOMs from different origins. Here we use eight NOM isolates of different molecular composition and ultrahigh resolution Fourier-transform ion cyclotron resonance-mass spectrometry (ESI-FT-ICR-MS) to determine the mol- ecular composition of platinum NM-NOM corona as a function of NOM composition and NM surface coating. We observed that the composition of PtNM-NOM corona varied with the composition of the original NOM. The per- centage of NOM formulas that formed PVP-PtNM-NOM corona was higher than those formed citrate-PtNM-NOM corona, due to increased sorption of NOM formulas, in particular condensed hydrocarbons, to the PVP coating. The relative abundance of heteroatom formulas (CHON, CHOS, and CHOP) was higher in the PVP-PtNM-NOM corona than in citrate-PtNM-corona which was in turn higher than those in the original NOM isolate, indicating preferential partitioning of heteroatom-rich molecules to NM surfaces. The relative abundance of CHO, CHON, CHOS, CHOP and condensed hydrocarbons in PtNM-NOM corona increased with their increase in NOM isolates. Furthermore, PtNM-NOM corona is rich with compounds with high molecular weight. This study demonstrates that the composition and properties of PtNM-NOM corona depend on NOM molecular properties and PtNM sur- face coating. The results here provide evidence of molecular interactions between NOM and NMs, which are crit- ical to understanding NM colloidal properties (e.g., surface charge and stability), interaction forces (e.g., van der Waals and hydrphobic), environmental behaviors (e.g., aggregation, disso.

Baalousha, Mohammed↗
DOE OSTI2022

The role of dissolved pyrogenic carbon from biochar in the sorption of As(V) in biogenic iron (oxyhydr)oxides

Water contamination by arsenic (As) affects millions of people around the world, making techniques to immobilize or remove this contaminant a pressing societal need. Biochar and iron (oxyhydr)oxides [in particular, biogenic iron (oxyhydr)oxides (BIOS)] offer the possibility of stabilizing As in remediation systems. However, little is known about the potential antagonism in As sorption generated by the dissolved organic carbon (DOC) from biochar, or whether DOC affects how As(V) interacts with BIOS. For this reason, our objectives were to evaluate the i) As(V) sorption potential in BIOS when there is presence of DOC from pyrolyzed biochars at different temperatures; and ii) identify whether the presence of DOC alters the surface complexes formed by As(V) sorbed in the BIOS. We conducted As(V) sorption experiments with BIOS at circumneutral pH conditions and in the presence of DOC from sugarcane (Saccharum officinarum) straw biochar at pyrolyzed 350 (BC350) and 750 °C (BC750). The As(V) content was quantified by inductively coupled plasma mass spectrometry, and the BIOS structure and As(V) sorption mechanisms were investigated by X-ray absorption spectroscopy. In addition, the organic moieties comprising the DOC from biochars were investigated by attenuated total reflectance Fourier transform infrared spectroscopy. The addition of DOC did not change the biomineral structure or As(V) oxidation state. The presence of DOC, however, reduced by 25 % the sorption of As(V), with BC350 being responsible for the greatest reduction in As(V) sorption capacity. Structural modeling revealed As(V) predominantly formed binuclear bidentate surface complexes on BIOS. The presence of DOC did not change the binding mechanism of As(V) in BIOS, suggesting that the reduction of As(V) sorption to BIOS was due to site blocking. In conclusion, our results bring insights into the fate of As(V) in surface waters and provide a basis for understanding the competitive sorption of As(V) in environments with biochar application.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Prediction of alkaline earth metal ion adsorption on goethite for various background electrolytes with the CD-MUSIC model

As water scarcity drives the use of more saline water sources, contaminant fate and transport models must capture the impact of high concentrations of alkaline earth metal ions (AEMs) and background electrolytes in these more complex waters. By utilizing macroscopic adsorption data from various electrolyte systems, a Charge Distribution – Multisite Complexation (CD-MUSIC) model, capable of incorporating electrolyte adsorption, was able to accurately simulate the adsorption behavior of alkaline earth metal ions onto goethite. The modeling effort was guided by previous spectroscopic and surface complexation modeling of alkaline earth metal adsorption and built on previous CD-MUSIC modeling that accounted for changes in crystal face contributions to the surface site density as a function of specific surface area. Further, the model was constrained to consider only two dominant surface complex species for each metal ion adsorption reaction. These two species were selected from 44 possible species through objective curve fitting of single-solute macroscopic adsorption data. While most of the alkaline earth metal surface complexes formed outer-sphere complexes at the goethite surface, an inner-sphere species was utilized for Mg 2+ . With the surface complex species and equilibrium constants obtained from this study, the calibrated model successfully predicted alkaline earth metal ion adsorption over a wide range of solution and surface conditions; the model predictions encompassed a wide range of pH (5–11), solute/solid ratio (1.37 × 10 -5 – 8.33 × 10 -4 mol -solute /g -solid ), ionic strengths (0.01 M – 0.7 M), and background electrolytes (Na + , Cs + , Rb + , Cl - , and NO 3 - ) using the same crystal face contribution methodology for site density, capacitance values, and surface acidity constants adopted for proton and cadmium adsorption in previous work (Han and Katz, 2019). Model simulations for a range of background water chemistries demonstrated the potential for Mg 2+ to reduce Cd 2+ adsorption to goethite in model seawater and oil- and gas-produced waters.

42 ENGINEERING↗