Melanie A. Mayes

Impact of Uranyl−Calcium−Carbonato Complexes on Uranium(VI) Adsorption to Synthetic and Natural Sediments

Adsorption on soil and sediment solids may decrease aqueous uranium concentrations and limit its propensity for migration in natural and contaminated settings. Uranium adsorption will be controlled in large part by its aqueous speciation, with a particular dependence on the presence of dissolved calcium and carbonate. Here we quantify the impact of uranyl speciation on adsorption to both goethite and sediments from the Hanford Clastic Dike and Oak Ridge Melton Branch Ridgetop formations. Hanford sediments were preconditioned with sodium acetate and acetic acid to remove carbonate grains, and Ca and carbonate were reintroduced at defined levels to provide a range of aqueous uranyl species. U(VI) adsorption is directly linked to UO(2)(2+) speciation, with the extent of retention decreasing with formation of ternary uranyl-calcium-carbonato species. Adsorption isotherms under the conditions studied are linear, and K(d) values decrease from 48 to 17 L kg(-1) for goethite, from 64 to 29 L kg (-1) for Hanford sediments, and from 95 to 51 L kg(-1) for Melton Branch sediments as the Ca concentration increases from 0 to 1 mM at pH 7. Our observations reveal that, in carbonate-bearing waters, neutral to slightly acidic pH values ( approximately 5) and limited dissolved calcium are optimal for uranium adsorption.

Development of microbial-enzyme-mediated decomposition model parameters through steady-state and dynamic analyses

We developed a microbial-enzyme-mediated decomposition (MEND) model, based on the Michaelis-Menten kinetics, that describes the dynamics of physically defined pools of soil organic matter (SOC). These include particulate, mineral-associated, dissolved organic matter (POC, MOC, and DOC, respectively), microbial biomass, and associated exoenzymes. The ranges and/or distributions of parameters were determined by both analytical steady-state and dynamic analyses with SOC data from the literature. We used an improved multi-objective parameter sensitivity analysis (MOPSA) to identify the most important parameters for the full model: maintenance of microbial biomass, turnover and synthesis of enzymes, and carbon use efficiency (CUE). The model predicted that an increase of 2 degrees C (baseline temperature 12 degrees C) caused the pools of POC-cellulose, MOC, and total SOC to increase with dynamic CUE and decrease with constant CUE, as indicated by the 50% confidence intervals. Regardless of dynamic or constant CUE, the changes in pool size of POC, MOC, and total SOC varied from -8% to 8% under +2 degrees C. The scenario analysis using a single parameter set indicates that higher temperature with dynamic CUE might result in greater net increases in both POC-cellulose and MOC pools. Different dynamics of various SOC pools reflected the catalytic functions of specific enzymes targeting specific substrates and the interactions between microbes, enzymes, and SOC. With the feasible parameter values estimated in this study, models incorporating fundamental principles of microbial-enzyme dynamics can lead to simulation results qualitatively different from traditional models with fast/slow/passive pools.

Multispecies transport of metal-EDTA complexes and chromate through undisturbed columns of weathered fractured saprolite

Abstract Laboratory-scale tracer experiments were conducted to investigate the geochemical and hydrological processes that govern the fate and transport of organically chelated radionuclides and toxic metals in undisturbed saturated columns of weathered, fractured shale saprolite. Three long-term, reactive contaminant injections were pulsed onto three separate soil columns, with the following influent mixtures: (1) 109CdEDTA2−, (2) 109CdEDTA2− and 57,58Co(II)EDTA2−, and (3) 109CdEDTA2−, 57Co(III)EDTA−, and H51CrO4−. Both single and multiple species experiments were conducted to determine the importance of interaction between the contaminants and competition for surface sites. Flow interruption was used to identify physical and chemical non-equilibrium (PNE and CNE) which were caused by multiple pore-region flow and rate-limited chemical reactions, respectively. Reactive contaminant transport through the fractured, weathered shale was affected by sorption, redox, and dissociation reactions, which were mediated by soil organic matter and surficial oxides of Fe, Mn, and Al. The transport of CdEDTA2− was significantly influenced by ligand-promoted dissolution of subsurface Fe and Al sources, resulting in the liberation of Cd2+, Al(III)EDTA− and Fe(III)EDTA−. Flow interruption confirmed that the surface-mediated dissociation reaction was time-dependent, with the stability of the CdEDTA2− complex dependent on its residence time within the soil. The migration of Co(II)EDTA2− was dominated by oxidization to the highly stable Co(III)EDTA− species, and elevated effluent Mn2+ suggested that surficial Mn(IV) oxides likely catalyzed the redox reaction, though Fe-oxides may have also contributed to the reaction. Dissociation (12%) of the Co(II)EDTA2− complex was first observed during flow interruption, indicating that rate-limited dissociation of the complex by Fe-oxides may be significant under equilibrium conditions. The transport of HCrO4− was significantly altered by the reduction of mobile Cr(VI) to irreversibly bound Cr(III). The reduction reaction was catalyzed by surface-bound natural organic matter and flow interruption confirmed that the reaction was time-dependent. There was little evidence of competitive effects between the various contaminants in the multispecies experiments, since each was influenced by a different geochemical process during transport through the soil. The results of this study further support research findings that suggest anionic toxic metals and radionuclide–organic complexes can be significantly influenced by soil geochemical processes that can both enhance and impede the subsurface migration of these contaminants.

Transport of Multiple Tracers in Variably Saturated Humid Region Structured Soils and Semi-arid Region Laminated Sediments

The processes governing physical nonequilibrium (PNE)—coupled preferential flow and matrix diffusion—are diverse between humid and semi-arid regions, and are directly related to climate and rock/sediment type, and indirectly related to subsequent soil profile development. The fate and transport of contaminants in these variably saturated undisturbed media is largely a function of the influence of PNE processes. Large cores of laminated silts and sands were collected from the US Department of Energy Pacific Northwest National Laboratory (PNNL) in semi-arid south central Washington. Additional cores of weathered, fractured interbedded limestone and shale saprolite were collected from the Oak Ridge National Laboratory (ORNL) in humid eastern Tennessee. PNNL cores were collected parallel (FBP) and perpendicular (FXB) to bedding, and the ORNL core was 30° to bedding. Saturated and unsaturated transport experiments were performed using multiple nonreactive tracers that had different diffusion coefficients (Br−, PFBA, and PIPES), in order to identify the influence of PNE on the fate and transport of solutes. In the ORNL structured saprolite, solute transport was governed by coupled preferential flow and matrix diffusion, as evidenced by tracer separation and highly asymmetric breakthrough curves (BTC). BTCs became more symmetric as preferential flowpaths became inactive during drainage. Tracer separation persisted during unsaturated flow suggesting the continued importance of nonequilibrium mass transfer between flowpaths and the immobile water that was held in the soil matrix. No evidence of PNE was observed under near-saturated conditions in the semi-arid region (PNNL) laminated silts and sands. Unsaturated flow in cores with discontinuous layering resulted in preferential flow and the development of perched, immobile water as evidenced by early breakthrough and separation of tracers. Conversely, transport parallel to laterally continuous beds did not result in preferential flow, the development of perched water, or tracer separation regardless of water content. These observations suggested that desaturation had two effects: (1) grain size variations between individual beds resulted in different antecedent water contents, and (2) the exchange of water and solutes between individual sedimentary beds was subsequently inhibited. Under unsaturated conditions, these effects may promote either stable lateral flow, or unstable vertical finger flow coupled with the development of perched, immobile water.

CXTFIT/Excel-A modular adaptable code for parameter estimation, sensitivity analysis and uncertainty analysis for laboratory or field tracer experiments

We implemented the widely used CXTFIT code in Excel to provide flexibility and added sensitivity and uncertainty analysis functions to improve transport parameter estimation and to facilitate model discrimination for multi-tracer experiments on structured soils. Analytical solutions for one-dimensional equilibrium and nonequilibrium convection dispersion equations were coded as VBA functions so that they could be used as ordinary math functions in Excel for forward predictions. Macros with user-friendly interfaces were developed for optimization, sensitivity analysis, uncertainty analysis, error propagation, response surface calculation, and Monte Carlo analysis. As a result, any parameter with transformations (e.g., dimensionless, log-transformed, species-dependent reactions, etc.) could be estimated with uncertainty and sensitivity quantification for multiple tracer data at multiple locations and times. Prior information and observation errors could be incorporated into the weighted nonlinear least squares method with a penalty function. Users are able to change selected parameter values and view the results via embedded graphics, resulting in a flexible tool applicable to modeling transport processes and to teaching students about parameter estimation. The code was verified by comparing to a number of benchmarks with CXTFIT 2.0. It was applied to improve parameter estimation for four typical tracer experiment data sets in the literature using multi-model evaluation and comparison. Additional examples were included to illustrate the flexibilities and advantages of CXTFIT/Excel. The VBA macros were designed for general purpose and could be used for any parameter estimation/model calibration when the forward solution is implemented in Excel. A step-by-step tutorial, example Excel files and the code are provided as supplemental material.

Microbial dormancy improves development and experimental validation of ecosystem model

Climate feedbacks from soils can result from environmental change followed by response of plant and microbial communities, and/or associated changes in nutrient cycling. Explicit consideration of microbial life-history traits and functions may be necessary to predict climate feedbacks owing to changes in the physiology and community composition of microbes and their associated effect on carbon cycling. Here we developed the microbial enzyme-mediated decomposition (MEND) model by incorporating microbial dormancy and the ability to track multiple isotopes of carbon. We tested two versions of MEND, that is, MEND with dormancy (MEND) and MEND without dormancy (MEND_wod), against long-term (270 days) carbon decomposition data from laboratory incubations of four soils with isotopically labeled substrates. MEND_wod adequately fitted multiple observations (total C–CO2 and 14C–CO2 respiration, and dissolved organic carbon), but at the cost of significantly underestimating the total microbial biomass. MEND improved estimates of microbial biomass by 20–71% over MEND_wod. We also quantified uncertainties in parameters and model simulations using the Critical Objective Function Index method, which is based on a global stochastic optimization algorithm, as well as model complexity and observational data availability. Together our model extrapolations of the incubation study show that long-term soil incubations with experimental data for multiple carbon pools are conducive to estimate both decomposition and microbial parameters. These efforts should provide essential support to future field- and global-scale simulations, and enable more confident predictions of feedbacks between environmental change and carbon cycling.

Representation of Dormant and Active Microbial Dynamics for Ecosystem Modeling

Dormancy is an essential strategy for microorganisms to cope with environmental stress. However, global ecosystem models typically ignore microbial dormancy, resulting in notable model uncertainties. To facilitate the consideration of dormancy in these large-scale models, we propose a new microbial physiology component that works for a wide range of substrate availabilities. This new model is based on microbial physiological states and the major parameters are the maximum specific growth and maintenance rates of active microbes and the ratio of dormant to active maintenance rates. A major improvement of our model over extant models is that it can explain the low active microbial fractions commonly observed in undisturbed soils. Our new model shows that the exponentially-increasing respiration from substrate-induced respiration experiments can only be used to determine the maximum specific growth rate and initial active microbial biomass, while the respiration data representing both exponentially-increasing and non-exponentially-increasing phases can robustly determine a range of key parameters including the initial total live biomass, initial active fraction, the maximum specific growth and maintenance rates, and the half-saturation constant. Our new model can be incorporated into existing ecosystem models to account for dormancy in microbially-driven processes and to provide improved estimates of microbial activities.

Activation Energy of Extracellular Enzymes in Soils from Different Biomes

Enzyme dynamics are being incorporated into soil carbon cycling models and accurate representation of enzyme kinetics is an important step in predicting belowground nutrient dynamics. A scarce number of studies have measured activation energy (Ea) in soils and fewer studies have measured Ea in arctic and tropical soils, or in subsurface soils. We determined the Ea for four typical lignocellulose degrading enzymes in the A and B horizons of seven soils covering six different soil orders. We also elucidated which soil properties predicted any measurable differences in Ea. β-glucosidase, cellobiohydrolase, phenol oxidase and peroxidase activities were measured at five temperatures, 4, 21, 30, 40, and 60°C. Ea was calculated using the Arrhenius equation. β-glucosidase and cellobiohydrolase Ea values for both A and B horizons in this study were similar to previously reported values, however we could not make a direct comparison for B horizon soils because of the lack of data. There was no consistent relationship between hydrolase enzyme Ea and the environmental variables we measured. Phenol oxidase was the only enzyme that had a consistent positive relationship between Ea and pH in both horizons. The Ea in the arctic and subarctic zones for peroxidase was lower than the hydrolases and phenol oxidase values, indicating peroxidase may be a rate limited enzyme in environments under warming conditions. By including these six soil types we have increased the number of soil oxidative enzyme Ea values reported in the literature by 50%. This study is a step towards better quantifying enzyme kinetics in different climate zones.

Selective Sorption of Dissolved Organic Carbon Compounds by Temperate Soils

Background Physico-chemical sorption onto soil minerals is one of the major processes of dissolved organic carbon (OC) stabilization in deeper soils. The interaction of DOC on soil solids is related to the reactivity of soil minerals, the chemistry of sorbate functional groups, and the stability of sorbate to microbial degradation. This study was conducted to examine the sorption of diverse OC compounds (D-glucose, L-alanine, oxalic acid, salicylic acid, and sinapyl alcohol) on temperate climate soil orders (Mollisols, Ultisols and Alfisols). Methodology Equilibrium batch experiments were conducted using 0–100 mg C L−1 at a solid-solution ratio of 1∶60 for 48 hrs on natural soils and on soils sterilized by γ-irradiation. The maximum sorption capacity, Qmax and binding coefficient, k were calculated by fitting to the Langmuir model. Results Ultisols appeared to sorb more glucose, alanine, and salicylic acid than did Alfisols or Mollisols and the isotherms followed a non-linear pattern (higher k). Sterile experiments revealed that glucose and alanine were both readily degraded and/or incorporated into microbial biomass because the observed Qmax under sterile conditions decreased by 22–46% for glucose and 17–77% for alanine as compared to non-sterile conditions. Mollisols, in contrast, more readily reacted with oxalic acid (Qmax of 886 mg kg−1) and sinapyl alcohol (Qmax of 2031 mg kg−1), and no degradation was observed. The reactivity of Alfisols to DOC was intermediate to that of Ultisols and Mollisols, and degradation followed similar patterns as for Ultisols. Conclusion This study demonstrated that three common temperate soil orders experienced differential sorption and degradation of simple OC compounds, indicating that sorbate chemistry plays a significant role in the sorptive stabilization of DOC.

Spatial arrangement of organic compounds on a model mineral surface: implications for soil organic matter stabilization.

The complexity of the mineral-organic carbon interface may influence the extent of stabilization of organic carbon compounds in soils, which is important for global climate futures. The nanoscale structure of a model interface was examined here by depositing films of organic carbon compounds of contrasting chemical character, hydrophilic glucose and amphiphilic stearic acid, onto a soil mineral analogue (Al2O3). Neutron reflectometry, a technique which provides depth-sensitive insight into the organization of the thin films, indicates that glucose molecules reside in a layer between Al2O3 and stearic acid, a result that was verified by water contact angle measurements. Molecular dynamics simulations reveal the thermodynamic driving force behind glucose partitioning on the mineral interface: The entropic penalty of confining the less mobile glucose on the mineral surface is lower than for stearic acid. The fundamental information obtained here helps rationalize how complex arrangements of organic carbon on soil mineral surfaces may arise.