Uncertainty analysis of multi-rate kinetics of uranium desorption from sediments.

Multi-rate surface complexation models have been proposed to describe the kinetics of uranyl (U(VI)) surface complexation reactions (SCR) rate-limited by diffusive mass transfer to and from intragranular sorption sites in subsurface sediments. In this study, a Bayesian-based, Differential Evolution Markov Chain method was used to assess the uncertainty and to identify factors controlling the uncertainties of the multi-rate SCR model. The rate constants in the multi-rate SCR were estimated with and without assumption of a specified lognormal distribution to test the lognormal assumption typically used to minimize the number of the rate constants in the multi-rate model. U(VI) desorption under variable chemical conditions from a contaminated sediment at US Hanford 300 Area, Washington was used as an example. The results indicated that the estimated rate constants without a specified lognormal assumption approximately followed a lognormal distribution, indicating that the lognormal is an effective assumption for the rate constants in the multi-rate SCR model. However, those rate constants with their corresponding half-lives longer than the experimental durations for model characterization had larger uncertainties and could not be reliably estimated. The uncertainty analysis revealed that the time-scale of the experiments for calibrating the multi-rate SCR model, the assumption for the rate constant distribution, the geochemical conditions involved in predicting U(VI) desorption, and equilibrium U(VI) speciation reaction constants were the major factors contributing to the extrapolation uncertainties of the multi-rate SCR model. Overall, the results from this study demonstrated that the multi-rate SCR model with a lognormal distribution of its rate constants is an effective approach for describing rate-limited U(VI) desorption; however, the model contains uncertainties, especially for those smaller rate constants, that require careful consideration for predicting U(VI) sorption and desorption.

[1]  Simulating adsorption of U(VI) under transient groundwater flow and hydrochemistry: Physical versus chemical nonequilibrium model , 2011 .

[2]  David B. Dunson,et al.  Bayesian Data Analysis , 2010 .

[3]  Scott C Brooks,et al.  U(VI) adsorption to heterogeneous subsurface media: application of a surface complexation model. , 2002, Environmental science & technology.

[4]  G. C. Tiao,et al.  Bayesian inference in statistical analysis , 1973 .

[5]  Sylvia Richardson,et al.  Markov Chain Monte Carlo in Practice , 1997 .

[6]  Charles F. Harvey,et al.  Aquifer remediation: A method for estimating mass transfer rate coefficients and an evaluation of pulsed pumping , 1994 .

[7]  James A. Davis,et al.  Newly recognized hosts for uranium in the Hanford Site vadose zone , 2009 .

[8]  K. Williams,et al.  Rate‐limited U(VI) desorption during a small‐scale tracer test in a heterogeneous uranium‐contaminated aquifer , 2012 .

[9]  J. S. Hunter,et al.  Statistics for experimenters : an introduction to design, data analysis, and model building , 1979 .

[10]  Tapabrata Maiti,et al.  Bayesian Data Analysis (2nd ed.) (Book) , 2004 .

[11]  D. Kent,et al.  Characterizing particle‐scale equilibrium adsorption and kinetics of uranium(VI) desorption from U‐contaminated sediments , 2013 .

[12]  Teresa B. Culver,et al.  Modeling the Desorption of Organic Contaminants from Long-Term Contaminated Soil Using Distributed Mass Transfer Rates , 1997 .

[13]  John M. Zachara,et al.  Scale‐dependent desorption of uranium from contaminated subsurface sediments , 2008 .

[14]  D. Higdon,et al.  Accelerating Markov Chain Monte Carlo Simulation by Differential Evolution with Self-Adaptive Randomized Subspace Sampling , 2009 .

[15]  Steven C. Smith,et al.  Sorption of Cs+ to micaceous subsurface sediments from the Hanford site, USA , 2002 .

[16]  Wassana Yantasee,et al.  Microscopic reactive diffusion of uranium in the contaminated sediments at Hanford, United States , 2006 .

[17]  E. Ziegel.  Statistics and Data Analysis in Geology (3rd ed.) , 2005 .

[18]  S. Kerisit,et al.  Diffusion and adsorption of uranyl carbonate species in nanosized mineral fractures. , 2012, Environmental science & technology.

[19]  James A. Davis,et al.  Approaches to surface complexation modeling of Uranium(VI) adsorption on aquifer sediments , 2004 .

[20]  W. D. Ray,et al.  Statistics for Experiments. An Introduction to Design, Data Analysis and Model Building , 1979 .

[21]  S. Sorooshian,et al.  A Shuffled Complex Evolution Metropolis algorithm for optimization and uncertainty assessment of hydrologic model parameters , 2002 .

[22]  M. Bradbury,et al.  A mechanistic description of Ni and Zn sorption on Na-montmorillonite Part II: modelling , 1997 .

[23]  Robert M. Smith,et al.  NIST Critically Selected Stability Constants of Metal Complexes Database , 2004 .

[24]  A. Maes,et al.  Cesium ion equilibria in illite clay , 1982 .

[25]  J. Zachara,et al.  Kinetics of uranium(VI) desorption from contaminated sediments: effect of geochemical conditions and model evaluation. , 2009, Environmental science & technology.

[26]  S. Gorelick,et al.  Multiple‐Rate Mass Transfer for Modeling Diffusion and Surface Reactions in Media with Pore‐Scale Heterogeneity , 1995 .

[27]  H. Lilliefors.  On the Kolmogorov-Smirnov Test for Normality with Mean and Variance Unknown , 1967 .

[28]  C. Tiedeman,et al.  Effective Groundwater Model Calibration , 2007 .

[29]  Mary C. Hill,et al.  UCODE_2005 and six other computer codes for universal sensitivity analysis, calibration, and uncertainty evaluation constructed using the JUPITER API , 2006 .

[30]  J. Kruschke.  Doing Bayesian Data Analysis: A Tutorial with R and BUGS , 2010 .

[31]  J. Zachara,et al.  Effect of grain size on uranium(VI) surface complexation kinetics and adsorption additivity. , 2011, Environmental science & technology.

[32]  James A. Davis,et al.  Comparison of in situ uranium KD values with a laboratory determined surface complexation model , 2004 .

[33]  J. Zachara,et al.  Characterization of the intragranular water regime within subsurface sediments: Pore volume, surface area, and mass transfer limitations , 2011 .

[34]  A. Brix.  Bayesian Data Analysis, 2nd edn , 2005 .

[35]  D. Kent,et al.  Quantifying Differences in the Impact of Variable Chemistry on Equilibrium Uranium(VI) Adsorption Properties of Aquifer Sediments , 2011, Environmental science & technology.

[36]  S. P. Neuman,et al.  Estimation of spatial covariance structures by adjoint state maximum likelihood cross validation: 3. Application to hydrochemical and isotopic data , 1989 .

[37]  S. P. Neuman,et al.  estimation of spatial covariance structures by adjoint state maximum likelihood cross validation: 1. Theory , 1989 .

[38]  R. Guillaumont,et al.  Update on the chemical thermodynamics of uranium, neptunium, plutonium, americium and technetium , 2003 .

[39]  A. Raftery,et al.  How Many Iterations in the Gibbs Sampler , 1991 .

[40]  Chongxuan Liu,et al.  Kinetic desorption and sorption of U(VI) during reactive transport in a contaminated Hanford sediment. , 2005, Environmental science & technology.

[41]  Andrew Gelman,et al.  Inference from Simulations and Monitoring Convergence , 2011 .

[42]  M. Bradbury,et al.  A mechanistic description of Ni and Zn sorption on Na-montmorillonite Part I: Titration and sorption measurements , 1997 .

[43]  J. Zachara,et al.  Multispecies diffusion models: A study of uranyl species diffusion , 2011 .

[44]  Javier Samper,et al.  Reactive transport model and apparent Kd of Ni in the near field of a HLW repository in granite , 2012, Comput. Geosci..

[45]  S. Heald,et al.  Dissolution of uranyl microprecipitates in subsurface sediments at Hanford Site, USA , 2004 .

[46]  S. Brooks,et al.  Determination of the formation constants of ternary complexes of uranyl and carbonate with alkaline earth metals (Mg2+, Ca2+, Sr2+, and Ba2+) using anion exchange method. , 2006, Environmental science & technology.

[47]  C. Tiedeman,et al.  Effective Groundwater Model Calibration: With Analysis of Data, Sensitivities, Predictions, and Uncertainty , 2007 .

[48]  Cajo J. F. ter Braak,et al.  A Markov Chain Monte Carlo version of the genetic algorithm Differential Evolution: easy Bayesian computing for real parameter spaces , 2006, Stat. Comput..

[49]  Henning Prommer,et al.  A field‐scale reactive transport model for U(VI) migration influenced by coupled multirate mass transfer and surface complexation reactions , 2010 .