Geochemical Modeling: EQ3/6 and PHREEQC
Explore the source record for details and available documents.
Domski, Paul Stanley↗
SEARCH · Search NASA
Results for “PHREEQC”
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.
Data for Zheng et al. (2025), "AquaMEND: Reconciling multiple impacts of salinization on soil carbon biogeochemistry"
Soil salinization, exacerbated by climate change, poses a global threat to coastal ecosystems and soil function. Salinity affects soil carbon cycling by directly impacting microbial activity and indirectly altering soil physicochemical properties, but current models inadequately represent these complexities. This dataset contains the observational and modeling data from Zheng et al. (2025), which described a process-based modeling framework that couples soil solution chemistry with microbial carbon cycling reactions to study the impacts of soil salinization. This conceptual model is implemented numerically into the open-source geochemical program PHREEQC 3.0 (Parkhurst and Appelo, 2013). This dataset consists of: - Figure2_AquaMEND_salinity_buffer: Contains model simulation outputs to assess the impact of three different cation exchange and surface complexation processes on salinity buffering (Fig. 2 from Zheng et al. 2025). - Figure3_Salinity_function: Contains salinity function fitting for literature data (Fig. 3 from Zheng et al. 2025). - Figure4_AquaMEND_microbial_mechanisms: Contains model simulation outputs for testing various microbial process-based hypotheses related to soil salinization, including microbial mortality, carbon use efficiency (CUE), extracellular enzyme activity, and other microbial mechanisms (Fig. 4 from Zheng et al. 2025). - Figure5_AquaMEND_Redox: Contains on model simulation outputs to evaluate shifts among key redox processes, such as aerobic respiration, sulfate reduction, and methanogenesis (Fig.5 from Zheng et al. 2025). - Figure6_AquaMEND_sorption: Contains on model simulation outputs for investigating the effects of salinity on dissolved organic matter (DOM) sorption and desorption processes (Fig. 6 from Zheng et al. 2025). - Figure7_AquaMEND_process_couple: Contains on model simulation outputs for exploring coupled biotic-abiotic processes and their interactions (Fig. 7 from Zheng et al. 2025). - data: Includes datasets used to develop salinity response functions and evaluate salinity buffering capacity. Datasets for MEND model calibration. - database: Contains the `.dat` file required by PHREEQC for model execution. - README.md: A Markdown plain text file describing the computational tools and directories. Files are a mixture of plain text CSV (comma-separated value) and plain text *.dat files written by the model; no special software is required to read them.
EARTH SCIENCE > AGRICULTURE > SOILS > SOIL SALINIT↗
Modeling Anaerobic Soil Organic Carbon Decomposition in Arctic Polygon Tundra: Insights into Soil Geochemical Influences on Carbon Mineralization: Modeling Archive
This anaerobic carbon decomposition model framework was developed with explicit representation of fermentation, methanogenesis and iron reduction by combining three well-known modeling approaches developed in different disciplines. A pool-based model to represent upstream carbon transformations and replenishment of dissolved organic carbon (DOC) pool, a thermodynamically-based model to calculate rate kinetics and biomass growth for methanogenesis and Fe(III) reduction, and a humic ion-binding model for aqueous phase speciation and pH calculations. The combined models are implemented in the open source geochemical model PHREEQC (V3.0). Installation of PHREEQC is required to run this model. This dataset contains one zipped file and one pdf.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).
54 ENVIRONMENTAL SCIENCES↗
An Optimization-Based Law of Mass Action Precipitation/Dissolution Model
Rare earth elements (REE) and many other critical minerals are necessary for the manufacturing of modern everyday technologies, including microchips, batteries and electric motors. Recovery of these materials typically involves aqueous systems which can be modeled as chemical equilibrium problems. One common method for solving these problems is the law of mass action approach (LMA), where a system of non-linear equations involving the equilibrium constants is solved. However, despite being theoretically simple, these problems are in practice very difficult to solve. Currently, the use of iterative heuristics based on saturation indices to decide on which species and reactions to include in the calculations is the state of the art to arrive at a solution. Here, we present an optimization-based alternative to solve chemical equilibria problems involving precipitation/dissolution reactions without the need for such heuristics. Our approach is first validated against the LMA software MINTEQ and PHREEQC for a number of case studies, and then applied to a novel REE recovery process reported in the literature. Overall, our approach was found to have close agreement with MINTEQ and PHREEQC, and we were able to successfully replicate the reported yield and purity for the published REE process.
42 ENGINEERING↗
Plutonium Speciation and Oxidation State Distributions in the Presence of Citrate
We explored the speciation and kinetics of the Pu(VI)-citrate and Pu(III)-citrate systems (pH m = 2.5–11.0, I = 0.1 M NaCl, T = 23 °C, O 2 (g) < 2 ppm) using ultraviolet–visible–near-infrared (UV–vis–NIR) spectrophotometry, solvent extraction, and PHREEQC modeling. Formation constants were determined for PuO 2 (HcitH)(aq) (log K° 1,1 = 1.09 ± 0.05) and PuO 2 (HcitH)(citH) 3– (log K° 1,2 = -0.20 ± 0.07), and evidence for (PuO 2 ) m (citH –k ) n (OH) x 2m(3+k)n–x was identified under alkaline conditions. Pu(VI) species were found to be less stable in the presence of citrate than in the absence of citrate (t ≤ 168 days); the rate of reduction increased with increasing pH. Further, the direct reduction of Pu(VI) to Pu(IV) was required to fit experimental data in the presence of citrate but did not improve the fit for Pu in the absence of citrate. We also observed increased Pu(III) stability in the presence of citrate (t ≤ 293 days), with higher concentrations of Pu(III) favored at lower pH. Finally, we provide evidence of a radiolysis-driven mechanism for the citrate-mediated reduction of plutonium that involves electron transfer from the oxidative breakdown of citrate. Our work highlights the need to investigate the redox effect of organic ligands on plutonium oxidation states under repository-relevant conditions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The potential for New Mexico basalts to sequester CO2 with a focus on critical element mobility
This study investigates the potential for New Mexico basalts to serve as geologic CO₂ sequestration media through mineral carbonation. Batch-type experiments were conducted using crushed basalts from the Carrizozo and Taos regions at low PCO₂ (0.7–2 bar) and temperatures up to 40°C over 180 days. The experiments evaluated the kinetics of CO₂-water-rock interaction, cation release, and critical element mobility (Cu, Ni, Co, Zn). Solid and fluid samples were analyzed using SEM, EMPA, ICP-OES, and ICP-MS, and geochemical modeling was performed with GEM-Selektor and PHREEQC. Results indicated that up to 80% of the injected CO₂ was mineralized within 100 days, forming secondary carbonates such as siderite and ankerite. Additionally, significant mobilization of critical metals was observed, which has implications for resource recovery and mine-waste valorization in future CO₂ sequestration projects.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Mercury Complexation with Dissolved Organic Matter Released from Thirty-Six Types of Biochar
Previous studies show mercury (Hg) can be effectively removed from solution by biochar, but limited attention was paid on the complexation between Hg and components released from biochars, e.g. dissolved organic matter (DOM). Here, aqueous data from batch-style experiments were modeled using PHREEQC, incorporating thermodynamic constants between Hg and DOM, which was assumed to be composed of thiol, carboxylic, and phenolic functional groups. Modelling results suggest that > 99% Hg complexed with thiol groups in DOM. The modelled concentrations of Hg–DOM complexes from low-T (low-temperature, 300°C) biochars were greater than from high-T (600°C) biochars. The concentrations of Hg–DOM complexes were lower in wood-based than in agricultural residue- and manure-based biochars. Hg–DOM complexes may affect Hg speciation, bioavailability, transport, and methylation processes. This research describes a method to evaluate Hg–DOM interactions, and the results indicate extra caution regarding Hg–DOM complex formation is required in the selection of biochar for Hg remediation.
54 ENVIRONMENTAL SCIENCES↗
Speciation in the Fe(III)-Cl(I)-H 2 O System at 298.15 K, 313.15 K, and 333.15 K (25 °C, 40 °C, and 60 °C)
This article presents the results on speciation of ferric iron generated by the dissolution of chemical reagent hydromolysite (ferric chloride hexahydrate, FeCl 3 :6H 2 O) in water at 298.15 K, 313.15 K, and 333.15 K (25 °C, 40 °C, and 60 °C). Experiments were performed with a thermoregulated system up to the equilibrium point, as manifested by solution pH. Solution samples were analyzed in terms of concentration, pH, and electrical conductivity. Measurements of density and refractive index were obtained at different temperatures and iron concentrations. A decrease of pH was observed with the increase in the amount of dissolved iron, indicating that ferric chloride is a strong electrolyte that reacts readily with water. Experimental results were modeled using the hydrogeochemical code PHREEQC in order to obtain solution speciation. Cations and neutral and anion complexes were simultaneously present in the system at the studied conditions according to model simulations, where dominant species included Cl - , FeCl 2+ , FeCl 2 + , FeOHCl 2 0 , and H + . A decrease in the concentration of Cl - and Fe 3+ ions took place with increasing temperature due to the association of Fe-Cl species. Standard equilibrium constants for the formation of FeOHCl 2 0 obtained in this study were log K f 0 = -0.8 ± 0.01 at 298.15 K (25 °C), -0.94 ± 0.02 at 313.15 K (40 °C), and -1.03 ± 0.01 at 333.15 K (60 °C).
36 MATERIALS SCIENCE↗
Laboratory column experiments and transport modeling to evaluate retardation of uranium in an aquifer downgradient of a uranium in-situ recovery site
The purpose of this study was to determine the attenuation potential and retardation of uranium in sediments taken from boreholes at the Smith-Ranch Highland in-situ recovery (ISR) site. Five column experiments with four different sediments were conducted to study the effects of variable mineralogy and alkalinity on uranium breakthrough. Uranium transport was modeled with PHREEQC using a generalized composite surface complexation model (GC SCM) with one, two, and, three generic surfaces, respectively. Reactive surface areas were approximated with PEST using BET derived surface areas to constrain fitting parameters. Uranium breakthrough was delayed by a factor of 1.68, 1.69 and 1.47 relative to the non-reactive tracer for three of the 5 experiments at an alkalinity of 540 mg/l. A sediment containing smectite and kaolinite retained uranium by a factor of 2.80 despite a lower measured BET surface area. Decreasing alkalinity to 360 mg/l from 540 mg/l increased retardation by a factor of 4.26. Model fits correlated well to overall BET surface area in the three columns where clay content was less than 1%. For the sediment with clay, models consistently understated uranium retardation when reactive surface sites were restricted by BET results. Calcite saturation was shown to be a controlling factor for uranium desorption as the pH of the system changes. A pH of 6 during a secondary background water flush remobilized previously sorbed uranium resulting in a secondary uranium peak at twice the influent concentrations. Furthermore, this study demonstrates the potential of GC SCM models to predict uranium transport in sediments with homogenous mineral composition, but highlights the need for further research to understand the role of sediment clay composition and calcite saturation in uranium transport.
58 GEOSCIENCES↗
Chemical and environmental stability of monazite-cheralite solid solutions Ln 1-2 x Ca x Th x PO 4 (Ln = Pr, Nd; x = 0–0.15): A thermodynamic study
Monazite-cheralite ceramics are a promising waste form for actinides. To elucidate the long-term behavior of this matrix in aqueous solutions, this study measured thermodynamic data for Th-rhabdophanes Ln 1-2 x Ca x Th x PO 4 ·nH 2 O (with Ln = Pr, Nd; x = 0–0.15) and the associated anhydrous monazite-cheralites Ln 1-2 x Ca x Th x PO 4 . Further, solubility experiments at 298K and high temperature oxide melt solution calorimetry were combined for calculation of ΔG$^°_f$, ΔH$^°_f$ and S$^°_m$ of Th-rhabdophanes and associated monazite-cheralites. Standard solubility constants were employed in a geochemical simulation using the PHREEQC software, the results of which confirmed the high chemical stability of the monazite-cheralite phases and supported their use as a specific conditioning matrix for the long-term immobilization of actinides.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
A chemistry-informed hybrid machine learning approach to predict metal adsorption onto mineral surfaces
Historically, surface complexation model (SCM) constants and distribution coefficients (K d ) have been employed to quantify mineral-based retardation effects controlling the fate of metals in subsurface geologic systems. Our recent SCM development workflow, based on the Lawrence Livermore National Laboratory Surface Complexation/Ion Exchange (L-SCIE) database, illustrated a community FAIR data approach to SCM development by predicting uranium(VI)-quartz adsorption for a large number of literature-mined data. Here, we present an alternative hybrid machine learning (ML) approach that shows promise in achieving equivalent high-quality predictions compared to traditional surface complexation models. At its core, the hybrid random forest (RF) ML approach is motivated by the proliferation of incongruent SCMs in the literature that limit their applicability in reactive transport models. Our hybrid ML approach implements PHREEQC-based aqueous speciation calculations; values from these simulations are automatically used as input features for a random forest (RF) algorithm to quantify adsorption and avoid SCM modeling constraints entirely. Named the LLNL Speciation Updated Random Forest (L-SURF) model, this hybrid approach is shown to have applicability to U(VI) sorption cases driven by both ion-exchange and surface complexation, as is shown for quartz and montmorillonite cases. The approach can be applied to reactive transport modeling and may provide an alternative to the costly development of self-consistent SCM reaction databases.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
Chemical degradation of fly ash blended concrete with the seasonal variation of rainwater in a radioactive waste repository: A thermodynamic modeling approach
Highlights: • The thermodynamic model for the cementitious system was successfully developed. • The chemical degradation of concrete was greatly affected by rainwater types. • Inorganic carbon species in rainwater reduced the chemical resistance of concrete. • The durability of concrete to rainwater decreased with increasing fly ash content. This study presents a long-term degradation behavior of fly ash blended concrete in a vault type low- and intermediate-level radioactive waste (LILW) repository by thermodynamic equilibrium calculations using PHREEQC combined with CEMDATA18 database. Since rainwater plays a predominant role as leachate for the concrete waste container in the repository, the effect of chemical properties of seasonal rainwater on the chemical degradation of concrete was analyzed. In addition, the impact of the blending with fly ash on the chemical degradation was evaluated through a gradual replacement of ordinary Portland cement (OPC) by fly ash. Regardless of the replacement level, autumn rainwater showed the greatest influence on the concrete degradation owing to the presence of inorganic carbon species, which induce the carbonation. The resistance of concrete to chemical degradation by rainwater was persistently reduced due to the decreasing formation of portlandite and calcium silicate hydrate (C-S-H) with increasing replacement level.
36 MATERIALS SCIENCE↗
Oil/water/rock wettability: Influencing factors and implications for low salinity water flooding in carbonate reservoirs
Wettability of the oil/brine/rock system is an essential petro-physical parameter which governs subsurface multiphase flow behaviour and the distribution of fluids, thus directly affecting oil recovery. Recent studies [1–3] show that manipulation of injected brine composition can enhance oil recovery by shifting wettability from oil-wet to water-wet. However, what factor(s) control system wettability has not been completely elucidated due to incomplete understanding of the geochemical system. To isolate and identify the key factors at play we used in this paper SO 4 2 —free solutions to examine the effect of salinity (formation brine/FB, 10 times diluted formation brine/10 dFB, and 100 times diluted formation brine/100 dFB) on the contact angle of oil droplets at the surface of calcite. We then compared contact angle results with predictions of surface complexation by low salinity water using PHREEQC software. We demonstrate that the conventional dilution approach likely triggers an oil-wet system at low pH, which may explain why the low salinity water EOR-effect is not always observed by injecting low salinity water in carbonated reservoirs. pH plays a fundamental role in the surface chemistry of oil/brine interfaces, and wettability. Our contact angle results show that formation brine triggered a strong water-wet system (35°) at pH 2.55, yet 100 times diluted formation brine led to a strongly oil-wet system (contact angle = 175°) at pH 5.68. Surface complexation modelling correctly predicted the wettability trend with salinity; the bond product sum ([>CaOH 2 + ][–COO - ] + [>CO 3 - ][–NH + ] + [>CO 3 - ][–COOCa + ]) increased with decreasing salinity. Finally, at pH < 6 dilution likely makes the calcite surface oil-wet, particularly for crude oils with high base number. Yet, dilution probably causes water wetness at pH > 7 for crude oils with high acid number.
02 PETROLEUM↗
Periodic changes in effluent chemistry at cold-water geyser: Crystal geyser in Utah
Crystal geyser is a CO 2 -driven cold-water geyser which was originally drilled in the late 1930’s in Green River, Utah. By utilizing a suite of temporal groundwater sample datasets, in situ monitoring of temperature, pressure, pH and electrical conductivity from multiple field trips to Crystal geyser from 2007 to 2014, periodic trends in groundwater chemistry from the geyser effluent were identified. Based on chemical characteristics, the primary sourcing aquifers are characterized to be both the Entrada and Navajo Sandstones with a minor contribution from Paradox Formation brine. The single eruption cycle at Crystal geyser lasted over four days and was composed of four parts: Minor Eruption (mEP), Major Eruption (MEP), Aftershock Eruption (Ae) and Recharge (R). During the single eruption cycle, dissolved ionic species vary 0–44% even though the degree of changes for individual ions are different. Generally, Na + , K + , Cl - and SO 4 2- regularly decrease at the onset and throughout the MEP. These species then increase in concentration during the mEP. In contrast, Ca 2+ , Mg 2+ , Fe 2+ and Sr 2+ increase and decrease in concentration during the MEP and mEP, respectively. The geochemical inverse modeling with PHREEQC was conducted to characterize the contribution from three end-members (Entrada Sandstone, Navajo Sandstone and Paradox Formation brine) to the resulting Crystal geyser effluent. Results of the inverse modeling showed that, during the mEP, the Navajo, Entrada and brine supplied 62–65%, 36–33% and 1–2%, respectively. During the MEP, the contribution shifted to 53–56%, 45–42% and 1–2% for the Navajo, Entrada and Paradox Formation brine, respectively. Finally, these changes in effluent characteristics further support the hypothesis by Watson et al. (2014) that the mEP and MEP are driven by different sources and mechanisms.
54 ENVIRONMENTAL SCIENCES↗
Community Data Mining Approach for Surface Complexation Database Development
This paper presents a comprehensive data-to-model workflow, including a findable, accessible, interoperable, reusable (FAIR) community sorption database (newly developed LLNL Surface Complexation/Ion Exchange (L-SCIE) database) along with a data fitting workflow to efficiently optimize surface complexation reaction constants with multiple surface complexation model (SCM) constructs. This workflow serves as a universal framework to mine, compile, and analyze large numbers of published sorption data as well as to estimate reaction constants for parameterizing reactive transport models. Here the framework includes (1) data digitization from published papers, (2) data unification including unit conversions, and (3) data-model integration and reaction constant estimation using geochemical software PHREEQC coupled with the universal parameter estimation code PEST. We demonstrate our approach using an analysis of U(VI) sorption to quartz based on a first L-SCIE implementation, concluding that a multisite SCM construct with carbonate surface species yielded the best fit to community data. Surface complexation reaction constants extracted from this approach captured all available sorption data available in the literature and provided insight into previously published reaction constants and surface complexation model constructs. The L-SCIE sorption database presented herein allows for automating this approach across a wide range of metals and minerals and implementing novel machine learning approaches to reactive transport in the future.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
The Uranium from Seawater Program at the Pacific Northwest National Laboratory: Overview of Marine Testing, Adsorbent Characterization, Adsorbent Durability, Adsorbent Toxicity, and Deployment Studies
The Pacific Northwest National Laboratory’s (PNNL) Marine Science Laboratory (MSL) located along the coast of Washington State is evaluating the performance of uranium adsorption materials being developed for seawater extraction under realistic marine conditions with natural seawater. Two types of exposure systems were employed in this program: flow-through columns for testing of fixed beds of individual fibers and pellets and a recirculating water flume for testing of braided adsorbent material. Testing consists of measurements of the adsorption of uranium and other elements from seawater as a function of time, typically 42 to 56 day exposures, to determine the adsorbent capacity and adsorption rate (kinetics). Analysis of uranium and other trace elements collected by the adsorbents was conducted following strong acid digestion of the adsorbent with 50% aqua regia using either Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES) or Inductively Coupled Plasma Mass Spectrometer (ICP-MS). The ORNL 38H adsorbent had a 56 day adsorption capacity of 3.30 ± 0.68 g U/ kg adsorbent (normalized to a salinity of 35 psu), a saturation adsorption capacity of 4.89 ± 0.83 g U/kg of adsorbent material (normalized to a salinity of 35 psu) and a half-saturation time of 28 ± 10 days. The AF1 adsorbent material had a 56 day adsorption capacity of 3.9 ± 0.2 g U/kg adsorbent material (normalized to a salinity of 35 psu), a saturation capacity of 5.4 ± 0.2 g U/kg adsorbent material (normalized to a salinity of 35 psu) and a half saturation time of 23 ± 2 days. The ORNL amidoxime-based adsorbent materials are not specific for uranium, but also adsorb other elements from seawater. The major doubly charged cations in seawater (Ca and Mg) account for a majority of the cations adsorbed (61% by mass and 74% by molar percent). For the ORNL AF1 adsorbent material, U is the 4th most abundant element adsorbed by mass and 7th most abundant by molar percentage
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Coupled Lattice Boltzmann Modeling Framework for Pore-Scale Fluid Flow and Reactive Transport
In this paper, we propose a modeling framework for pore-scale fluid flow and reactive transport based on a coupled lattice Boltzmann model (LBM). We develop a modeling interface to integrate the LBM modeling code parallel lattice Boltzmann solver and the PHREEQC reaction solver using multiple flow and reaction cell mapping schemes. The major advantage of the proposed workflow is the high modeling flexibility obtained by coupling the geochemical model with the LBM fluid flow model. Consequently, the model is capable of executing one or more complex reactions within desired cells while preserving the high data communication efficiency between the two codes. Meanwhile, the developed mapping mechanism enables the flow, diffusion, and reactions in complex pore-scale geometries. We validate the coupled code in a series of benchmark numerical experiments, including 2D single-phase Poiseuille flow and diffusion, 2D reactive transport with calcite dissolution, as well as surface complexation reactions. The simulation results show good agreement with analytical solutions, experimental data, and multiple other simulation codes. In addition, we design an AI-based optimization workflow and implement it on the surface complexation model to enable increased capacity of the coupled modeling framework. Compared to the manual tuning results proposed in the literature, our workflow demonstrates fast and reliable model optimization results without incorporating pre-existing domain knowledge.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Uncertainty propagation in pore water chemical composition calculation using surrogate models
Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗