Caroline Tebes-Stevens

Multicomponent transport with coupled geochemical and microbiological reactions: Model description and example simulations

A reactive transport code (FEREACT) has been developed to examine the coupled effects of two-dimensional steady-state groundwater flow, equilibrium aqueous speciation reactions, and kinetically-controlled interphase reactions. The model uses an iterative two-step (SIA-1) solution algorithm to incorporate the effects of the geochemical and microbial reaction processes in the governing equation for solute transport in the subsurface. This SIA-1 method improves upon the convergence behavior of the traditional sequential iterative approach (SIA) through the inclusion of an additional first-order term from the Taylor Series expansion of the kinetic reaction rate expressions. The ability of FEREACT to simulate coupled reactive processes was demonstrated by modeling the transport of a radionuclide (cobalt, 60Co2+) and an organic ligand (ethylenediaminetetraacetate, EDTA4−) through a column packed with an iron oxide-coated sand. The reaction processes considered in this analysis included equilibrium aqueous speciation reactions and three types of kinetic reactions: adsorption, surface dissolution, and biodegradation.

Calculation of reaction parameter sensitivity coefficients in multicomponent subsurface transport models

A significant source of uncertainty in multicomponent transport models is due to reaction parameter estimation error. In this paper, two codes are developed to compute reaction parameter sensitivity coefficients, which provide a quantitative measure of the impact of this uncertainty. One of the codes solves the full coupled system of sensitivity equations, while the other applies the operator splitting approach to decouple the sensitivity equations for each chemical component. CPU time comparisons demonstrate that operator splitting significantly reduces the simulation time for solving both the state and sensitivity equations. These comparisons also confirm that solution of the sensitivity equations is more efficient than a traditional direct perturbation sensitivity analysis. Sensitivity calculations for an example transport problem with cobalt and nitrilotriacetate (NTA) are used to gain insight into the relative significance of reaction processes and to rank individual reaction parameters in terms of importance. In this example, speciation reactions cause an indirect linking of the NTA biodegradation and the cobalt sorption reactions, making the NTA concentration nearly as sensitive to the sorption parameters as to the degradation parameters.

In Vitro Metabolism of the Fungicide and Environmental Contaminant trans-bromuconazole and Implications for Risk Assessment

trans-Bromuconazole is a chiral chemical representative of a class of triazole derivatives known to inhibit specific fungal cytochrome P-450 (CYP) reactions. Kinetic measurements and delineation of metabolic pathways for triazole chemicals within in vitro hepatic microsomes are needed for accurate risk assessment and predictive in vivo physiological modeling. The studies described here were conducted with rat liver microsomes to determine Michaelis–Menten saturation kinetic parameters (V max and K M) for trans-bromuconazole using both substrate depletion and product formation reaction velocities. Kinetic parameters determined for trans-bromuconazole depletion at varying protein levels incubated at physiological temperature 37°C resulted in a K M value of 1.69 μM and a V max value of 1398 pmol/min/mg protein. The concomitant linear formation of two metabolites identified using liquid chromatography/time-of-flight mass spectrometry (LC/MS-TOF) and LC-MS/MS indicated hydroxylation of the trans-bromuconazole dichlorophenyl ring moiety. K M values determined for the hydroxylated metabolites were 0.87 and 1.03 μM, with V max values of 449 and 694 pmol/min/mg protein, respectively. Chemical inhibition assays and studies conducted with individual purified human recombinant enzymes indicated the CYP3A subfamily was primarily responsible for biotransformation of the parent substrate. Additionally, trans-bromuconazole was found to undergo stereoselective metabolism as evidenced by a change in the enantiomeric ratio (trans−/trans +) with respect to time.

Evaluating the sensitivity of a subsurface multicomponent reactive transport model with respect to transport and reaction parameters.

The input variables for a numerical model of reactive solute transport in groundwater include both transport parameters, such as hydraulic conductivity and infiltration, and reaction parameters that describe the important chemical and biological processes in the system. These parameters are subject to uncertainty due to measurement error and due to the spatial variability of properties in the subsurface environment. This paper compares the relative effects of uncertainty in the transport and reaction parameters on the results of a solute transport model. This question is addressed by comparing the magnitudes of the local sensitivity coefficients for transport and reaction parameters. General sensitivity equations are presented for transport parameters, reaction parameters, and the initial (background) concentrations in the problem domain. Parameter sensitivity coefficients are then calculated for an example problem in which uranium(VI) hydrolysis species are transported through a two-dimensional domain with a spatially variable pattern of surface complexation sites. In this example, the reaction model includes equilibrium speciation reactions and mass transfer-limited non-electrostatic surface complexation reactions. The set of parameters to which the model is most sensitive includes the initial concentration of one of the surface sites, the formation constant (Kf) of one of the surface complexes and the hydraulic conductivity within the reactive zone. For this example problem, the sensitivity analysis demonstrates that transport and reaction parameters are equally important in terms of how their variability affects the model results.

Runoff water quality during drought in a zero-order georgia piedmont pasture: nitrogen and total organic carbon.

Approximately 11% of the Southern Piedmont (1.8 million ha) is used for pasture and hay production, mostly under low-input management. Few studies have investigated in the region long-term nitrogen and carbon losses in surface runoff, which can be significant. We present 1999 to 2009 hydrologic and water quality data from a rotationally grazed, 7.8-ha, zero-order pasture (W1) near Watkinsville in the Georgia Piedmont. Annual rainfall was 176 to 463 mm below the long-term average (1240 mm) in 7 of the 11 yr. There were 20 runoff events during 86 mo of below-average rainfall (deficit period), compared with 54 events during 46 mo of nondeficit period. Mean event flow-weighted concentration (in mg L) was 0.96 for nitrate-nitrogen (NO-N), 0.97 for ammonium-nitrogen (NH-N), 3.70 for total nitrogen (TN), and 9.12 for total organic carbon (TOC) ( = 43-47; limited due to instrument problem). Nutrient loads (in kg ha per event) averaged 0.04 for NO-N, 0.03 for NH-N, 0.19 for TN, and 0.54 for TOC. Total loads for N and TOC were 6 to 11 times greater from nondeficit than from deficit periods. The observed N concentrations, while well below maximum drinking water standard limits, could pose risk for eutrophication, which can be stimulated at lower concentrations. However, the ability of headwater streams, such as the one downstream of W1, to reduce nutrient concentrations might partially alleviate this concern. The results of this study point to the need to use a long-term dataset that includes measurements made in drought and wet years when evaluating the efficacy of water quality standards.

H2 consumption during the microbial reductive dehalogenation of chlorinated phenols and tetrachloroethene.

Competition for molecular hydrogen exists amonghydrogen-utilizing microorganismsin anoxic environments, and evidence suggeststhat lower hydrogen concentrations areobserved with more energetically favorableelectron-accepting processes. The transferof electrons to organochlorines via reductivedehalogenation reactions plays an importantrole in hydrogen dynamics in impacted systems. Westudied the flux of aqueous hydrogenconcentrations in methanogenic sediment microcosmsprior to and during reductivedehalogenation of a variety of substituted chlorophenols(CP) and tetrachloroethene(perchloroethylene, PCE). Mean hydrogen concentrationsduring reductive dehalogenationof 2,4-CP, 2,3,4-CP, and PCP were 3.6 nM, 4.1 nM,and 0.34 nM, respectively. Sedimentmicrocosms that were not dosed with chlorophenolsyet were actively methanogenicmaintained a significantly higher mean hydrogenconcentration of 9.8 nM. Duringactive PCE dehalogenation, sediment microcosmsmaintained a mean hydrogenconcentration of 0.82 nM. These data indicate thatduring limiting hydrogen production,the threshold ecosystem hydrogen concentration iscontrolled by microbial populationsthat couple hydrogen oxidation to thermodynamicallyfavorable electron acceptingreactions, including reductive dehalogenationof chloroaromatic compounds. Wealso present revised estimates for the Gibbsfree energy available from the reductivedehalogenation of a variety of substitutedchlorophenols based on recently publishedvalues of vapor pressure, solubility, and pKafor these compounds.

Demonstration of a consensus approach for the calculation of physicochemical properties required for environmental fate assessments

Eight software applications are compared for their performance in estimating the octanol-water partition coefficient (Kow), melting point, vapor pressure and water solubility for a dataset of polychlorinated biphenyls, polybrominated diphenyl ethers, polychlorinated dibenzodioxins, and polycyclic aromatic hydrocarbons. The predicted property values are compared against a curated dataset of measured property values compiled from the scientific literature with careful consideration given to the analytical methods used for property measurements of these hydrophobic chemicals. The variability in the predicted values from different calculators generally increases for higher values of Kow and melting point and for lower values of water solubility and vapor pressure. For each property, no individual calculator outperforms the others for all four of the chemical classes included in the analysis. Because calculator performance varies based on chemical class and property value, the geometric mean and the median of the calculated values from multiple calculators that use different estimation algorithms are recommended as more reliable estimates of the property value than the value from any single calculator.

Estimation of microbial reductive transformation rates for chlorinated benzenes and phenols using a quantitative structure-activity relationship approach.

A set of literature data was used to derive several quantitative structure-activity relationships (QSARs) to predict the rate constants for the microbial reductive dehalogenation of chlorinated aromatics. Dechlorination rate constants for 25 chloroaromatics were corrected for the effects of hydrophobic partitioning and adjusted for the observed distribution of product species. A number of physicochemical properties and molecular parameters were considered for inclusion in the QSARs. Multivariate statistical analyses were used to select the optimal set of descriptors to minimize multicollinearity between the descriptors, as well as to minimize the p-value of the regression coefficients. The final QSAR included four descriptors: The logarithm of the octanol-water partition coefficient (Kow), the summation of the Hammett sigma constants, and the sigma induction constants in the ortho and meta positions relative to the transformation reaction center. The predictive ability of this QSAR was evaluated using 24 site-specific rate constants that were measured in five separate studies and were not used to derive the expression. The peer-reviewed literature was screened carefully to ensure that all rate constant data were representative of environmentally relevant conditions.

Prediction of Hydrolysis Products of Organic Chemicals under Environmental pH Conditions

Cheminformatics-based software tools can predict the molecular structure of transformation products using a library of transformation reaction schemes. This paper presents the development of such a library for abiotic hydrolysis of organic chemicals under environmentally relevant conditions. The hydrolysis reaction schemes in the library encode the process science gathered from peer-reviewed literature and regulatory reports. Each scheme has been ranked on a scale of one to six based on the median half-life in a data set compiled from literature-reported hydrolysis rates. These ranks are used to predict the most likely transformation route when more than one structural fragment susceptible to hydrolysis is present in a molecule of interest. Separate rank assignments are established for pH 5, 7, and 9 to represent standard conditions in hydrolysis studies required for registration of pesticides in Organisation for Economic Co-operation and Development (OECD) member countries. The library is applied to predict the likely hydrolytic transformation products for two lists of chemicals, one representative of chemicals used in commerce and the other specific to pesticides, to evaluate which hydrolysis reaction pathways are most likely to be relevant for organic chemicals found in the natural environment.

The environmental fate of synthetic organic chemicals

Abstract This article focuses on the routes of transport and abiotic processes involved in the environmental transformation of synthetic organic chemicals and how molecular structure controls the products and lifetimes of several important classes of organic chemicals. The chapter also discusses the current methods to reliably determine the rates and products of degradation of new chemicals based on combinations of chemical structure and environmental processes as well as use of laboratory and field measurements. Methods are also discussed for use of structure activity relations for this purpose.