Direct and indirect photodegradation in aquatic systems mitigates photosensitized toxicity in screening-level substance risk assessments of selected petrochemical structures

Photochemical processes are typically not incorporated in screening-level substance risk assessments due to the complexity of modeling sunlight co-exposures and resulting interactions on environmental fate and effects. However, for many substances, sunlight exerts a profound influence on environmental degradation rates and ecotoxicities. Recent modeling advances provide an improved technical basis for estimating the effect of sunlight in modulating both substance exposure and toxicity in the aquatic environment. Screening model simulations were performed for 25 petrochemical structures with varied uses and environmental fate properties. Model predictions were evaluated by comparing the ratios of predicted exposure concentrations with and without light to the corresponding ratios of toxicity thresholds under the same conditions. The relative ratios of exposure and hazard in light vs. dark were then used to evaluate how inclusion of light modulates substance risk analysis. Results indicated that inclusion of light reduced PECs by factors ranging from 1.1- to 63-fold as a result of photodegradation, while reducing PNECs by factors ranging from 1- to 49-fold due to photoenhanced toxicity caused by photosensitization. Consequently, the presence of light altered risk quotients by factors that ranged from 0.1- to 17-fold, since the predicted increase in substance hazard was mitigated by the reduction in exposure. For many structures, indirect photodegradation decreases environmental exposures independently of the direct photolysis pathway which is associated with enhanced phototoxicity. For most of the scenarios and chemicals in the present work, photosensitization appears to be mitigated by direct and indirect degradation from sunlight exposure.

Disentangling the History of Deep Ocean Disposal for DDT and Other Industrial Waste Off Southern California

Ocean disposal of industrial waste from technical DDT [mainly 1,1'-(2,2,2-trichloroethane-1,1-diyl)bis(4-chlorobenzene), or 4,4'-DDT] manufacture occurred historically in the Southern California Bight. However, the paucity of historical records highlights uncertainties as to the mode, location, and timing of disposal or ongoing ecological effects of these wastes. This study combines sampling, chemical analysis, and numerical modeling of deep San Pedro Basin sediments revealing substantial DDT contamination that extends at least 25 km from the mainland. These findings narrate bulk DDT waste disposal to the offshore that peaked in the 1950s, prior to the onset of formal regulations; was agnostic to later-designated disposal sites; and has experienced sluggish transformation. Our findings further indicate an attenuating secondary source for the DDT daughter product, 1-chloro-4-[2,2-dichloro-1-(4-chlorophenyl)ethenyl]benzene (4,4'-DDE), which still deposits into deep San Pedro Basin sediments. While demonstrating the severity of DDT contamination to the region, these findings further define the burial potential of DDT wastes and inform the past, present, and future contamination potential that is needed to understand and predict ecological consequences. This work also points firmly to bulk, not containerized, disposal of DDT waste and to potential alternative contents of collocated waste.

Comprehensive Two-Dimensional Gas Chromatography with Peak Tracking for Screening of Constituent Biodegradation in Petroleum UVCB Substances

Petroleum substances, as archetypical UVCBs (substances of unknown or variable composition, complex reaction products, or biological substances), pose a challenge for chemical risk assessment as they contain hundreds to thousands of individual constituents. It is particularly challenging to determine the biodegradability of petroleum substances since each constituent behaves differently. Testing the whole substance provides an average biodegradation, but it would be effectively impossible to obtain all constituents and test them individually. To overcome this challenge, comprehensive two-dimensional gas chromatography (GC × GC) in combination with advanced data-handling algorithms was applied to track and calculate degradation half-times (DT50s) of individual constituents in two dispersed middle distillate gas oils in seawater. By tracking >1000 peaks (representing ∼53-54% of the total mass across the entire chromatographic area), known biodegradation patterns of oil constituents were confirmed and extended to include many hundreds not currently investigated by traditional one-dimensional GC methods. Approximately 95% of the total tracked peak mass biodegraded after 64 days. By tracking the microbial community evolution, a correlation between the presence of functional microbial communities and the observed progression of DT50s between chemical classes was demonstrated. This approach could be used to screen the persistence of GC × GC-amenable constituents of petroleum substance UVCBs.

Relative sensitivity of hydrodynamic, thermodynamic, and chemical processes for simulating the buoyant multiphase plume and surfacing flows of an oil and gas blowout

Deepwater hydrocarbon releases experience complex chemical and physical processes. To assess simplifications of these processes on model predictions, we present a sensitivity analysis using simulations for the Deepwater Horizon oil spill. We compare the buoyant multiphase plume metrics (trap height, rise time etc), the hydrocarbon mass flowrates at the near-field plume termination and their mass fractions dissolved in the water column and reaching the water surface. The baseline simulation utilizes a 19-component hydrocarbon model, live-fluid state equations, hydrate dynamics, and heat and mass transfer. Other simulations turn-off each of these processes, with the simplest one using inert oil and methane gas. Plume metrics are the least sensitive to the modeled processes and can be matched by adjusting the release buoyancy flux. The mass flowrate metrics are more sensitive. Both liquid- and gas-phase mass transfer should be modeled for accurate tracking of soluble components (e.g. C₁ - C₇ hydrocarbons) in the environment.

Modeling the GCxGC Elution Patterns of a Hydrocarbon Structure Library To Innovate Environmental Risk Assessments of Petroleum Substances

Comprehensive two-dimensional gas chromatography (GCxGC) offers unrivaled separation of petroleum substances, which can contain thousands of constituents or more. However, interpreting substance compositions from GCxGC data is costly and requires expertise. To facilitate environmental risk assessments, industries provide aggregated compositional information known as "hydrocarbon blocks" (HCBs), but these proprietary methods do not transparently associate the HCBs with GCxGC chromatogram data. These obstacles frustrate efforts to study the environmental risks of petroleum substances and associated environmental samples. To address this problem, we developed a GCxGC elution model for user-defined petroleum substance compositions. We calibrated the elution model to experimental GCxGC retention times of 56 known hydrocarbons by fitting three tunable model parameters to two candidate instrument methods. With the calibrated model, we simulated retention times for a library of 15,447-15,455 hydrocarbon structures (plus 40-48 predicted as chromatographically unretained) spanning 11 classes of petroleum substance constituents in the C₁₀-C₃₀ range. The resulting simulation data reveal that GCxGC retention times are quantitatively associated with hydrocarbon class and carbon number information throughout the GCxGC chromatogram. These innovations enable the development of transparent and efficient technical methods to investigate the chemical compositions and environmental properties of petroleum substances, including in environmental and lab-weathered samples.

Development of an offshore response guidance tool for determining the impact of SSDI on released gas and benzene from artificial subsea oil well blowout simulations

We present an analysis of 2225 simulations of artificial oil well blowouts in nearshore and offshore waters of Newfoundland, Canada. In the simulations, we coupled the VDROP-J and TAMOC models to simulate the fate and transport of oil and gas from the release to the sea surface. Simulations were conducted with and without subsea dispersant injection. We analyzed the simulation database to quantify the mass fraction of oil and gas that surfaces, the mass fraction of released benzene that surfaces, and the horizontal offset to the surfacing zone. These data are also synthesized to yield empirical correlations to predict these output metrics from key input parameters. These correlations are summarized in an excel spreadsheet that allows rapid evaluation of spill dynamics with minimal initial knowledge of spill details. We call this tool an offshore response guidance table, which allows exploration of spill dynamics under diverse spill and response options.

Deepwater Horizon 2010: Subsea dispersants protected responders from VOC exposure

Controversy remains on the use of Sub-Sea Dispersant Injection (SSDI) during the Deepwater Horizon (DWH) spill to minimize the exposure of responders on surface vessels to volatile organic compounds (VOC). Here, we use extensive evidence (>90,000 VOC measurements) collected near the oil well MC252 site during the DWH spill and demonstrate at a high level of statistical confidence that SSDI enhanced the safety and health conditions of the responders at the water surface through the reduction of airborne VOC concentrations in a dose-dependent manner. VOC levels on ships' decks were clearly diminished (p < 0.001) during subsea dispersant use, and incidents of peak concentrations (>50 ppm VOC) that could have been an immediate concern to worker health were reduced by a factor of ~6 to 19 when dispersants were delivered at the intended rate. SSDI thus played an important role in minimizing potential exposure to VOC, and should be embedded in guidelines and regulations for dispersant use.

Dynamics of Live Oil Droplets and Natural Gas Bubbles in Deep Water

Explaining the dynamics of gas-saturated live petroleum in deep water remains a challenge. Recently, Pesch et al. [ Environ. Eng. Sci. 2018, 35 (4), 289-299] reported laboratory experiments on methane-saturated oil droplets under emulated deep-water conditions, providing an opportunity to elucidate the underlying dynamical processes. We explain these observations with the Texas A&M Oil spill/Outfall Calculator (TAMOC), which models the pressure-, temperature-, and composition-dependent interactions between oil-gas phase transfer; aqueous dissolution; and densities and volumes of liquid oil droplets, gas bubbles, and two-phase droplet-bubble pairs. TAMOC reveals that aqueous dissolution removed >95% of the methane from ∼3.5 mm live oil droplets within 14.5 min, prior to gas bubble formation, during the experiments of Pesch et al. Additional simulations indicate that aqueous dissolution, fluid density changes, and gas-oil phase transitions (ebullition, condensation) may all contribute to the fates of live oil and gas in deep water, depending on the release conditions. Illustrative model scenarios suggest that 5 mm diameter gas bubbles released at a <470 m water depth can transport methane, ethane, and propane to the water surface. Ethane and propane can reach the water surface from much deeper releases of 5 mm diameter live oil droplets, during which ebullition occurs at water depths of <70 m.

Prediction of aqueous free energies of solvation using coupled QM and MM explicit solvent simulations

A method based on molecular dynamics simulations which employ two distinct levels of theory is proposed and tested for the prediction of Gibbs free energies of solvation for non-ionic solutes in water. The method consists of two additive contributions: (i) an evaluation of the free energy of solvation predicted by a computationally efficient molecular mechanics (MM) method; and (ii) an evaluation of the free energy difference between the potential energy surface of the MM method and that of a more computationally intensive first-principles quantum-mechanical (QM) method. The latter is computed by a thermodynamic integration method based on a series of shorter molecular dynamics simulations that employ weighted averages of the QM and MM force evaluations. The combined computational approach is tested against the experimental free energies of aqueous solvation for four solutes. For solute-solvent interactions that are found to be described qualitatively well by the MM method, the QM correction makes a modest improvement in the predicted free energy of aqueous solvation. However, for solutes that are found to not be adequately described by the MM method, the QM correction does not improve agreement with experiment. These preliminary results provide valuable insights into the novel concept of implementing thermodynamic integration between two model chemistries, suggesting that it is possible to use QM methods to improve upon the MM predictions of free energies of aqueous solvation.

Predicting Partitioning and Diffusion Properties of Nonpolar Chemicals in Biotic Media and Passive Sampler Phases by GC × GC

The chemical parameters needed to explain and predict bioavailability, biodynamics, and baseline toxicity are not readily available for most nonpolar chemicals detected in the environment. Here, we demonstrate that comprehensive two-dimensional gas chromatography (GC × GC) retention times can be used to predict 26 relevant properties for nonpolar chemicals, specifically: partition coefficients for diverse biotic media and passive sampler phases; aquatic baseline toxicity; and relevant diffusion coefficients. The considered biotic and passive sampler phases include membrane and storage lipids, serum and muscle proteins, carbohydrates, algae, mussels, polydimethylsiloxane, polyethylene, polyoxymethylene, polyacrylate, polyurethane, and semipermeable membrane devices. GC × GC-based chemical property predictions are validated with a compilation of 1038 experimental property data collected from the literature. As an example application, we overlay a map of baseline toxicity to fathead minnows onto the separated analyte signal of a polychlorinated alkanes (chlorinated paraffins) technical mixture that contains 7820 congeners. In a second application, GC × GC-estimated properties are used to parametrize multiphase partitioning models for mammalian tissues and organs. In a third example, we estimate chemical depuration kinetics for mussels. Finally, we illustrate an approach to screen the GC × GC chromatogram for nonpolar chemicals of potentially high concern, defined based on their GC × GC-estimated biopartitioning properties, diffusion properties, and baseline toxicity.

Formation and reactivity of inorganic and organic chloramines and bromamines during oxidative water treatment

The formation and further reactions of halamines during oxidative water treatment can be relevant for water quality. In this study, we investigated the formation and reactivity of several inorganic and organic halamines (monochloramine, N-chloromethylamine, N-chlorodimethylamine, monobromamine, dibromamine, N-bromomethylamine, N,N-dibromomethylamine, and N-bromodimethylamine) by kinetic experiments, transformation product analysis, and quantum chemical computations. Kinetic model simulations were conducted to evaluate the relevance of halamines for various water treatment scenarios. Halamines were quickly formed from the reaction of chlorine and bromine with ammonia or organic amines. Species-specific second-order rate constants for the reaction of chlorine and bromine with ammonia, methyl- and dimethylamine were in the order of 10⁶-10⁸ M^-1s^-1. The formed halamines were found to be reactive towards phenolic compounds, forming halogenated phenols via electrophilic aromatic substitution (phenol and resorcinol) or quinones via electron transfer (catechol and hydroquinone). At near neutral pH, apparent second-order rate constants for these reactions were in the order of 10^-4-10^-1 M^-1s^-1 for chloramines and 10¹-10² M^-1s^-1 for bromamines. Quantum chemical computations were used to determine previously unknown aqueous pK_a values, gas phase bond dissociation energies (BDE) and partial atomic charges of the halamines, allowing a better understanding of their reactivities. Kinetic model simulations, based on the results of this study, showed that during chlorination inorganic and organic chloramines are the main halamines formed. However, their further reactions with organic matter are outcompeted kinetically by chlorine. During ozonation, mainly inorganic bromamines are formed, since ozone quickly oxidizes organic amines. The further reactions of bromamine are typically outcompeted by ozone and thus generally of minor importance. The use of peracetic acid for saline ballast water treatment can result in the formation of substantial amounts of bromamines, which can react with dissolved organic matter and contribute to the formation of brominated products.

Elevated Concentrations of 4-Bromobiphenyl and 1,3,5-Tribromobenzene Found in Deep Water of Lake Geneva Based on GC×GC-ENCI-TOFMS and GC×GC-μECD

We quantified the concentrations of two little-studied brominated pollutants, 1,3,5-tribromobenzene (TBB) and 4-bromobiphenyl (4BBP), in the deep water column and sediments of Lake Geneva. We found aqueous concentrations of 625 ± 68 pg L^-1 for TBB and 668 ± 86 pg L^-1 for 4BBP over a depth range of 70-191.5 m (near-bottom depth), based on duplicate measurements taken at five depths during three separate 1 month sampling periods at our sampling site near Vidy Bay. These levels of TBB and 4BBP were 1 or 2 orders of magnitude higher than the quantified aqueous concentrations of the components of the pentabrominated biphenyl ether technical mixture, which is a flame retardant product that had a high production volume in Europe before 2001. We observed statistically significant vertical concentration trends for both TBB and 2,2',4,4',6-pentabromobiphenyl ether in the deep water column, which indicates that transport and/or degradation processes affect these compounds. These measurements were enabled by application of a comprehensive two-dimensional gas chromatograph coupled to an electron capture negative chemical ionization time-of-flight mass spectrometer (GC×GC-ENCI-TOFMS) and to a micro-electron capture detector (GC×GC-μECD). GC×GC-ENCI-TOFMS and GC×GC-μECD were found to be >10× more sensitive toward brominated pollutants than conventional GC×GC-EI-TOFMS (with an electron impact (EI) ionization source), the latter of which had insufficient sensitivity to detect these emerging brominated pollutants in the analyzed samples. GC×GC also enabled the estimation of several environmentally relevant partitioning properties of TBB and 4BBP, further confirming previous evidence that these pollutants are bioaccumulative and have long-range transport potential.

Equilibria and Speciation of Chloramines, Bromamines, and Bromochloramines in Water

The stabilities and speciation of the halamines in water are difficult to characterize experimentally. We provide theoretical estimates of aqueous standard free energies of formation for inorganic chloramines, bromamines, and bromochloramines, based on high-accuracy theoretical standard free energies of formation in gas phase combined with quantum chemical estimates of Henry's law constant. Based on comparisons between several theoretical and experimental datasets, we assign an error of 1.1-1.2 log unit for equilibrium constants of several reactions leading to halamines in water. The reactions of ammonia with HOCl or HOBr that lead to dichloramine, trichloramine, and tribromamine are found to be thermodynamically more favorable than was previously believed. The newly reported equilibrium data also allow us to propose rate constant values for some hydrolysis and disproportionation reactions of dichloramine, monobromamine, and bromochloramine. Finally, theoretical results indicate aqueous acid dissociation constant (pK_a) values of 1.5 ± 1 for NH₃Cl⁺, 0.8 ± 1 for NH₃Br⁺, 11.8 ± 1 for NHCl₂, and 12.5 ± 1 for NHBrCl. The present report provides a comprehensive data set describing the free energies of the neutral inorganic halamines, the anionic conjugate base species, and the cationic conjugate acid species, with approximately uniform uncertainty bounds assigned throughout.

Correction: Benchmark thermochemistry of chloramines, bromamines, and bromochloramines: halogen oxidants stabilized by electron correlation

Correction for 'Benchmark thermochemistry of chloramines, bromamines, and bromochloramines: halogen oxidants stabilized by electron correlation' by Daniela Trogolo et al., Phys. Chem. Chem. Phys., 2015, 17, 3584-3598.

Simulating Gas-Liquid-Water Partitioning and Fluid Properties of Petroleum under Pressure: Implications for Deep-Sea Blowouts

With the expansion of offshore petroleum extraction, validated models are needed to simulate the behaviors of petroleum compounds released in deep (>100 m) waters. We present a thermodynamic model of the densities, viscosities, and gas-liquid-water partitioning of petroleum mixtures with varying pressure, temperature, and composition based on the Peng-Robinson equation-of-state and the modified Henry's law (Krychevsky-Kasarnovsky equation). The model is applied to Macondo reservoir fluid released during the Deepwater Horizon disaster, represented with 279-280 pseudocomponents, including 131-132 individual compounds. We define >n-C8 pseudocomponents based on comprehensive two-dimensional gas chromatography (GC × GC) measurements, which enable the modeling of aqueous partitioning for n-C8 to n-C26 fractions not quantified individually. Thermodynamic model predictions are tested against available laboratory data on petroleum liquid densities, gas/liquid volume fractions, and liquid viscosities. We find that the emitted petroleum mixture was ∼29-44% gas and ∼56-71% liquid, after cooling to local conditions near the broken Macondo riser stub (∼153 atm and 4.3 °C). High pressure conditions dramatically favor the aqueous dissolution of C1-C4 hydrocarbons and also influence the buoyancies of bubbles and droplets. Additionally, the simulated densities of emitted petroleum fluids affect previous estimates of the volumetric flow rate of dead oil from the emission source.

Impacts of Polar Changes on the UV-induced Mineralization of Terrigenous Dissolved Organic Matter

Local climates in the Northern and Southern Hemisphere are influenced by Arctic Amplification and by interactions of the Antarctic ozone hole with climate change, respectively. Polar changes may affect hydroclimatic conditions in temperate regions, for example, by increasing the length and intensity of precipitation events at Northern Hemisphere midlatitudes. Additionally, global warming has led to the thawing of ancient permafrost soils, particularly in Arctic regions, due to Arctic Amplification. Both heavy precipitation events and thawing of permafrost are increasing the net transfer of terrestrially derived dissolved organic matter (DOM) from land to surface waters. In aquatic ecosystems, UV-induced oxidation of terrigenous DOM (tDOM) produces atmospheric CO2 and this process is one of several mechanisms by which natural organic matter in aquatic and soil environments may play an important role in climate feedbacks. The Arctic is particularly affected by these processes: for example, melting of Arctic sea ice allows solar UV radiation to penetrate into the ice-free Arctic Ocean and to cause photochemical reactions that result in bleaching and mineralization of tDOM. Open questions, in addition to those shown in the Graphical Abstract, remain regarding the resulting contributions of tDOM photomineralization to CO2 production and global warming.

Immediate ecotoxicological effects of short-lived oil spills on marine biota

Marine environments are frequently exposed to oil spills as a result of transportation, oil drilling or fuel usage. Whereas large oil spills and their effects have been widely documented, more common and recurrent small spills typically escape attention. To fill this important gap in the assessment of oil-spill effects, we performed two independent supervised full sea releases of 5 m(3) of crude oil, complemented by on-board mesocosm studies and sampling of accidentally encountered slicks. Using rapid on-board biological assays, we detect high bioavailability and toxicity of dissolved and dispersed oil within 24 h after the spills, occurring fairly deep (8 m) below the slicks. Selective decline of marine plankton is observed, equally relevant for early stages of larger spills. Our results demonstrate that, contrary to common thinking, even small spills have immediate adverse biological effects and their recurrent nature is likely to affect marine ecosystem functioning.

Explicit solvent simulations of the aqueous oxidation potential and reorganization energy for neutral molecules: gas phase, linear solvent response, and non-linear response contributions

First principles simulations were used to predict aqueous one-electron oxidation potentials (Eox) and associated half-cell reorganization energies (λaq) for aniline, phenol, methoxybenzene, imidazole, and dimethylsulfide. We employed quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations of the oxidized and reduced species in an explicit aqueous solvent, followed by EOM-IP-CCSD computations with effective fragment potentials for diabatic energy gaps of solvated clusters, and finally thermodynamic integration of the non-linear solvent response contribution using classical MD. A priori predicted Eox and λaq values exhibit mean absolute errors of 0.17 V and 0.06 eV, respectively, compared to experiment. We also disaggregate Eox into several well-defined free energy properties, including the gas phase adiabatic free energy of ionization (7.73 to 8.82 eV), the solvent-induced shift in the free energy of ionization due to linear solvent response (-2.01 to -2.73 eV), and the contribution from non-linear solvent response (-0.07 to -0.14 eV). The linear solvent response component is further apportioned into contributions from the solvent-induced shift in vertical ionization energy of the reduced species (ΔVIEaq) and the solvent-induced shift in negative vertical electron affinity of the ionized species (ΔNVEAaq). The simulated ΔVIEaq and ΔNVEAaq are found to contribute the principal sources of uncertainty in computational estimates of Eox and λaq. Trends in the magnitudes of disaggregated solvation properties are found to correlate with trends in structural and electronic features of the solute. Finally, conflicting approaches for evaluating the aqueous reorganization energy are contrasted and discussed, and concluding recommendations are given.

Accurate DFT Descriptions for Weak Interactions of Molecules Containing Sulfur

Dispersion corrected atom centered potentials (DCACPs) have been shown to significantly improve the density functional theory (DFT) description of weak interactions. In this work, we have calibrated a DCACP for sulfur in combination with the widely used Generalized Gradient Approximation (GGA) BLYP, thereby augmenting the existing library of DCACPs for the first- and second-row elements H, C, N, O, and rare gases. Three weakly bound complexes as well as elemental (orthorhombic) sulfur are used as test cases to evaluate the transferability of the DCACP to different chemical environments. It is found that the sulfur DCACP systematically improves the agreement of DFT-calculated weak interactions with respect to MP2 and CCSD(T) level results.

Development of Prediction Models for the Reactivity of Organic Compounds with Ozone in Aqueous Solution by Quantum Chemical Calculations: The Role of Delocalized and Localized Molecular Orbitals

Second-order rate constants (kO3) for the reaction of ozone with micropollutants are essential parameters for the assessment of micropollutant elimination efficiency during ozonation in water and wastewater treatment. Prediction models for kO3 were developed for aromatic compounds, olefins, and amines by quantum chemical molecular orbital calculations employing ab initio Hartree-Fock (HF) and density functional theory (B3LYP) methods. The kO3 values for aromatic compounds correlated well with the energy of a delocalized molecular orbital first appearing on an aromatic ring (i.e., the highest occupied molecular orbital (HOMO) or HOMO-n (n ≥ 0) when the HOMO is not located on the aromatic ring); the number of compounds tested (N) was 112, and the correlation coefficient (R(2)) values were 0.82-1.00. The kO3 values for olefins and amines correlated well with the energy of a localized molecular orbital (i.e., the natural bond orbital (NBO)) energy of the carbon-carbon π bond of olefins (N = 45, R(2) values of 0.82-0.85) and the NBO energy of the nitrogen lone-pair electrons of amines (N = 59, R(2) values of 0.81-0.83), respectively. Considering the performance of the kO3 prediction model and the computational costs, the HF/6-31G method is recommended for all aromatic groups and olefins investigated herein, whereas the HF/MIDI!, HF/6-31G*, or HF/6-311++G** methods are recommended for amines. Based on their mean absolute errors, the above models could predict kO3 within a factor of 4, on average, relative to the experimentally determined values. Overall, good correlations were also observed (R(2) values of 0.77-0.96) between kO3 predictions by quantum molecular orbital descriptors in this study and by the Hammett (σ) and Taft (σ*) constants from previously developed quantitative structure-activity relationship (QSAR) models. Hence, the quantum molecular orbital descriptors are an alternative to σ and σ*-values in QSAR applications and can also be utilized to estimate unknown σ or σ*-values.  .

GC×GC Quantification of Priority and Emerging Nonpolar Halogenated Micropollutants in All Types of Wastewater Matrices: Analysis Methodology, Chemical Occurrence, and Partitioning

We report the development and validation of a method to detect and quantify diverse nonpolar halogenated micropollutants in wastewater treatment plant (WWTP) influent, effluent, primary sludge, and secondary sludge matrices (including both the liquid and particle phases) by comprehensive two-dimensional gas chromatography (GC×GC) coupled to micro- electron capture detector (μECD). The 59 target analytes included toxaphenes, polychlorinated naphthalenes, organochlorine pesticides, polychlorinated biphenyls, polybrominated diphenyl ethers, and emerging persistent and bioaccumulative chemicals. The method is robust for a wide range of nonpolar halogenated micropollutants in all matrices. For most analytes, recoveries fell between 70% and 130% in all matrix types. GC×GC-μECD detections of several target analytes were confirmed qualitatively by further analysis with GC×GC coupled to electron capture negative chemical ionization-time-of-flight mass spectrometry (ENCI-TOFMS). We then quantified the concentrations and apparent organic solid-water partition coefficients (Kp) of target micropollutants in samples from a municipal WWTP in Switzerland. Several analyzed pollutants exhibited a high frequency of occurrence in WWTP stream samples, including octachloronaphthalene, PCB-44, PCB-52, PCB-153, PCB-180, several organochlorine pesticides, PBDE-10, PBDE-28, PBDE-116, musk tibetene, and pentachloronitrobenzene. Our results suggest that sorption to dissolved organic carbon (DOC) can contribute substantially to the apparent solids-liquid distribution of hydrophobic micropollutants in WWTP streams.

Molecular mechanism of NDMA formation from N,N-dimethylsulfamide during ozonation: quantum chemical insights into a bromide-catalyzed pathway

During ozonation of drinking water, the fungicide metabolite N,N-dimethylsulfamide (DMS) can be transformed into a highly toxic product, N-nitrosodimethylamine (NDMA). We used quantum chemical computations and stopped-flow experiments to evaluate a chemical mechanism proposed previously to describe this transformation. Stopped-flow experiments indicate a pK(a) = 10.4 for DMS. Experiments show that hypobromous acid (HOBr), generated by ozone oxidation of naturally occurring bromide, brominates the deprotonated DMS(-) anion with a near-diffusion controlled rate constant (7.1 ± 0.6 × 10(8) M(-1) s(-1)), forming Br-DMS(-) anion. According to quantum chemical calculations, Br-DMS has a pK(a) ∼ 9.0 and thus remains partially deprotonated at neutral pH. The anionic Br-DMS(-) bromamine can react with ozone with a high rate constant (10(5 ± 2.5) M(-1) s(-1)), forming the reaction intermediate (BrNO)(SO2)N(CH3)2(-). This intermediate resembles a loosely bound complex between an electrophilic nitrosyl bromide (BrNO) molecule and an electron-rich dimethylaminosulfinate ((SO2)N(CH3)2(-)) fragment, based on inspection of computed natural charges and geometric parameters. This fragile complex undergoes immediate (10(10 ± 2.5) s(-1)) reaction by two branches: an exothermic channel that produces NDMA, and an entropy-driven channel giving non-NDMA products. Computational results bring new insights into the electronic nature, chemical equilibria, and kinetics of the elementary reactions of this pathway, enabled by computed energies of structures that are not possible to access experimentally.

Benchmark thermochemistry of chloramines, bromamines, and bromochloramines: halogen oxidants stabilized by electron correlation

The free energy of the formation of NH₂Br at 298 K can be estimated by taking into account the total atomization energy of NH₂Br and the atomic and molecular contributions to the enthalpy and the entropy of formation of NH₂Br at 0 K and 298 K.

Exploring the aqueous vertical ionization of organic molecules by molecular simulation and liquid microjet photoelectron spectroscopy

To study the influence of aqueous solvent on the electronic energy levels of dissolved organic molecules, we conducted liquid microjet photoelectron spectroscopy (PES) measurements of the aqueous vertical ionization energies (VIEaq) of aniline (7.49 eV), veratrole alcohol (7.68 eV), and imidazole (8.51 eV). We also reanalyzed previously reported experimental PES data for phenol, phenolate, thymidine, and protonated imidazolium cation. We then simulated PE spectra by means of QM/MM molecular dynamics and EOM-IP-CCSD calculations with effective fragment potentials, used to describe the aqueous vertical ionization energies for six molecules, including aniline, phenol, veratrole alcohol, imidazole, methoxybenzene, and dimethylsulfide. Experimental and computational data enable us to decompose the VIEaq into elementary processes. For neutral compounds, the shift in VIE upon solvation, ΔVIEaq, was found to range from ≈-0.5 to -0.91 eV. The ΔVIEaq was further explained in terms of the influence of deforming the gas phase solute into its solution phase conformation, the influence of solute hydrogen-bond donor and acceptor interactions with proximate solvent molecules, and the polarization of about 3000 outerlying solvent molecules. Among the neutral compounds, variability in ΔVIEaq appeared largely controlled by differences in solute-solvent hydrogen-bonding interactions. Detailed computational analysis of the flexible molecule veratrole alcohol reveals that the VIE is strongly dependent on molecular conformation in both gas and aqueous phases. Finally, aqueous reorganization energies of the oxidation half-cell ionization reaction were determined from experimental data or estimated from simulation for the six compounds aniline, phenol, phenolate, veratrole alcohol, dimethylsulfide, and methoxybenzene, revealing a surprising constancy of 2.06 to 2.35 eV.

First day of an oil spill on the open sea: early mass transfers of hydrocarbons to air and water

During the first hours after release of petroleum at sea, crude oil hydrocarbons partition rapidly into air and water. However, limited information is available about very early evaporation and dissolution processes. We report on the composition of the oil slick during the first day after a permitted, unrestrained 4.3 m(3) oil release conducted on the North Sea. Rapid mass transfers of volatile and soluble hydrocarbons were observed, with >50% of ≤C17 hydrocarbons disappearing within 25 h from this oil slick of <10 km(2) area and <10 μm thickness. For oil sheen, >50% losses of ≤C16 hydrocarbons were observed after 1 h. We developed a mass transfer model to describe the evolution of oil slick chemical composition and water column hydrocarbon concentrations. The model was parametrized based on environmental conditions and hydrocarbon partitioning properties estimated from comprehensive two-dimensional gas chromatography (GC×GC) retention data. The model correctly predicted the observed fractionation of petroleum hydrocarbons in the oil slick resulting from evaporation and dissolution. This is the first report on the broad-spectrum compositional changes in oil during the first day of a spill at the sea surface. Expected outcomes under other environmental conditions are discussed, as well as comparisons to other models.

Mapping environmental partitioning properties of nonpolar complex mixtures by use of GC × GC

Comprehensive two-dimensional gas chromatography (GC × GC) is effective for separating and quantifying nonpolar organic chemicals in complex mixtures. Here we present a model to estimate 11 environmental partitioning properties for nonpolar analytes based on GC × GC chromatogram retention time information. The considered partitioning properties span several phases including pure liquid, air, water, octanol, hexadecane, particle natural organic matter, dissolved organic matter, and organism lipids. The model training set and test sets are based on a literature compilation of 648 individual experimental partitioning property data. For a test set of 50 nonpolar environmental contaminants, predicted partition coefficients exhibit root-mean-squared errors ranging from 0.19 to 0.48 log unit, outperforming Abraham-type solvation models for the same chemical set. The approach is applicable to nonpolar organic chemicals containing C, H, F, Cl, Br, and I, having boiling points ≤402 °C. The presented model is calibrated, easy to apply, and requires the user only to identify a small set of known analytes that adapt the model to the GC × GC instrument program. The analyst can thus map partitioning property estimates onto GC × GC chromatograms of complex mixtures. For example, analyzed nonpolar chemicals can be screened for long-range transport potential, aquatic bioaccumulation potential, arctic contamination potential, and other characteristic partitioning behaviors.

Resolving biodegradation patterns of persistent saturated hydrocarbons in weathered oil samples from the Deepwater Horizon disaster

Biodegradation plays a major role in the natural attenuation of oil spills. However, limited information is available about biodegradation of different saturated hydrocarbon classes in surface environments, despite that oils are composed mostly of saturates, due to the limited ability of conventional gas chromatography (GC) to resolve this compound group. We studied eight weathered oil samples collected from four Gulf of Mexico beaches 12-19 months after the Deepwater Horizon disaster. Using comprehensive two-dimensional gas chromatography (GC × GC), we successfully separated, identified, and quantified several distinct saturates classes in these samples. We find that saturated hydrocarbons eluting after n-C22 dominate the GC-amenable fraction of these weathered samples. This compound group represented 8-10%, or 38-68 thousand metric tons, of the oil originally released from Macondo well. Saturates in the n-C22 to n-C29 elution range were found to be partly biodegraded, but to different relative extents, with ease of biodegradation decreasing in the following order: n-alkanes > methylalkanes and alkylcyclopentanes+alkylcyclohexanes > cyclic and acyclic isoprenoids. We developed a new quantitative index designed to characterize biodegradation of >n-C22 saturates. These results shed new light onto the environmental fate of these persistent, hydrophobic, and mostly overlooked compounds in the unresolved complex mixtures (UCM) of weathered oils.

Reactivity of BrCl, Br₂, BrOCl, Br₂O, and HOBr toward dimethenamid in solutions of bromide + aqueous free chlorine

HOBr, formed via oxidation of bromide by free available chlorine (FAC), is frequently assumed to be the sole species responsible for generating brominated disinfection byproducts (DBPs). Our studies reveal that BrCl, Br(2), BrOCl, and Br(2)O can also serve as brominating agents of the herbicide dimethenamid in solutions of bromide to which FAC was added. Conditions affecting bromine speciation (pH, total free bromine concentration ([HOBr](T)), [Cl(-)], and [FAC](o)) were systematically varied, and rates of dimethenamid bromination were measured. Reaction orders in [HOBr](T) ranged from 1.09 (±0.17) to 1.67 (±0.16), reaching a maximum near the pK(a) of HOBr. This complex dependence on [HOBr](T) implicates Br(2)O as an active brominating agent. That bromination rates increased with increasing [Cl(-)], [FAC](o) (at constant [HOBr](T)), and excess bromide (where [Br(-)](o)>[FAC](o)) implicate BrCl, BrOCl, and Br(2), respectively, as brominating agents. As equilibrium constants for the formation of Br(2)O and BrOCl (aq) have not been previously reported, we have calculated these values (and their gas-phase analogues) using benchmark-quality quantum chemical methods [CCSD(T) up to CCSDTQ calculations plus solvation effects]. The results allow us to compute bromine speciation and hence second-order rate constants. Intrinsic brominating reactivity increased in the order: HOBr ≪ Br(2)O < BrOCl ≈ Br(2) < BrCl. Our results indicate that species other than HOBr can influence bromination rates under conditions typical of drinking water and wastewater chlorination.

Robust algorithm for aligning two-dimensional chromatograms

Comprehensive two-dimensional gas chromatography (GC × GC) chromatograms typically exhibit run-to-run retention time variability. Chromatogram alignment is often a desirable step prior to further analysis of the data, for example, in studies of environmental forensics or weathering of complex mixtures. We present a new algorithm for aligning whole GC × GC chromatograms. This technique is based on alignment points that have locations indicated by the user both in a target chromatogram and in a reference chromatogram. We applied the algorithm to two sets of samples. First, we aligned the chromatograms of twelve compositionally distinct oil spill samples, all analyzed using the same instrument parameters. Second, we applied the algorithm to two compositionally distinct wastewater extracts analyzed using two different instrument temperature programs, thus involving larger retention time shifts than the first sample set. For both sample sets, the new algorithm performed favorably compared to two other available alignment algorithms: that of Pierce, K. M.; Wood, Lianna F.; Wright, B. W.; Synovec, R. E. Anal. Chem.2005, 77, 7735-7743 and 2-D COW from Zhang, D.; Huang, X.; Regnier, F. E.; Zhang, M. Anal. Chem.2008, 80, 2664-2671. The new algorithm achieves the best matches of retention times for test analytes, avoids some artifacts which result from the other alignment algorithms, and incurs the least modification of quantitative signal information.

Role of environment for catalysis of the DNA repair enzyme MutY

Control of the N-glycosylase reaction by the DNA repair enzyme, MutY, entails the organization of solvent molecules. Classical molecular dynamics and QM/MM simulations were used to investigate the solvent and environment effects contributing to catalysis. Our findings suggest that the entire reaction is an energetically neutral process, in which the first step is rate determining, requiring protonation of adenine (N(7)) to initiate cleavage, and the second step is strongly exothermic, involving hydrolysis of an oxacarbenium ion intermediate. Although water molecules are catalytically active in both steps, the first step requires an entirely different level of solvent organization compared to the second. Needed to secure protonation at N(7), a long-term solvation pattern is established which facilitates the involvement of three out of the five structured water molecules in the active site. This persistent arrangement coordinates the catalytically active water molecules into prime positions to assist the proton transfer: (i) a water molecule frequently bridges the catalytic residues and (ii) the bridging water molecule is assisted by 1-2 other 'supporting' water molecules. To maintain this configuration, MutY, surprisingly, uses hydrophobic residues in combination with hydrophilic residues to tune the microenvironment into a 'water trap'. Hydrophilic residues prolong solvent residence times by maintaining hydrogen-bonding networks, whereas the hydrophobic residues constrain the positioning of the catalytic water molecules that assist the proton-transfer event. In this way, the enzyme uses both entropic and enthalpic considerations to guide the water-assisted reaction.

Development of bioreporter assays for the detection of bioavailability of long-chain alkanes based on the marine bacterium Alcanivorax borkumensis strain SK2

Long-chain alkanes are a major component of crude oil and therefore potentially good indicators of hydrocarbon spills. Here we present a set of new bacterial bioreporters and assays that allow to detect long-chain alkanes. These reporters are based on the regulatory protein AlkS and the alkB1 promoter from Alcanivorax borkumensis SK2, a widespread alkane degrader in marine habitats. Escherichia coli cells with the reporter construct reacted strongly to octane in short-term (6 h) aqueous suspension assays but very slightly only to tetradecane, in line with what is expected from its low water solubility. In contrast, long-term assays (up to 5 days) with A. borkumensis bioreporters showed strong induction with tetradecane and crude oil. Gel-immobilized A. borkumensis reporter cells were used to demonstrate tetradecane and crude oil bioavailability at a distance from a source. Alcanivorax borkumensis bioreporters induced fivefold more rapid and more strongly when allowed physical contact with the oil phase in standing flask assays, suggesting a major contribution of adhered cells to the overall reporter signal. Using the flask assays we further demonstrated the effect of oleophilic nutrients and biosurfactants on oil availability and degradation by A. borkumensis. The fluorescence signal from flask assays could easily be captured with a normal digital camera, making such tests feasible to be carried out on, e.g. marine oil responder vessels in case of oil accidents.

Disentangling oil weathering at a marine seep using GC x GC: broad metabolic specificity accompanies subsurface petroleum biodegradation

Natural seeps contribute nearly half of the oil entering the coastal ocean. However, environmental fate studies generally monitor fewer than 5% of these petroleum compounds. Hence, the rates and relevance of physical, chemical, and biological weathering processes are unknown for the large majority of hydrocarbons, both released from natural seeps and also from human activities. To investigate the specific compositional changes occurring in petroleum during subsurface degradation and submarine seepage, we studied the natural oil seeps offshore Santa Barbara, California with comprehensive, two-dimensional gas chromatography (GC x GC). With this technique, we quantified changes in the molecular diversity and abundance of hydrocarbons between subsurface reservoirs, a proximal sea floor seep, and the sea surface overlying the seep. We also developed methods to apportion hydrocarbon mass losses due to biodegradation, dissolution, and evaporation, for hundreds of tracked compounds that ascended from the subsurface to the sea floor to the sea surface. The results provide the first quantitative evidence of broad metabolic specificity for anaerobic hydrocarbon degradation in the subsurface and reveal new trends of rapid hydrocarbon evaporation at the sea surface. This study establishes GC x GC as a powerful technique for differentiating biological and physical weathering processes of complex mixtures at a molecular level.

Determination of biodiesel blending percentages using natural abundance radiocarbon analysis: testing the accuracy of retail biodiesel blends

Blends of biodiesel and petrodiesel are being used increasingly worldwide. Due to several factors, inaccurate blending of these two mixtures can occur. To test the accuracy of biodiesel blending, we developed and validated a radiocarbon-based method and then analyzed a variety of retail biodiesel blends. Error propagation analysis demonstrated that this method calculates absolute blend content with +/- 1% accuracy, even when real-world variability in the component biodiesel and petrodiesel sources is taken into account. We independently confirmed this accuracy using known endmembers and prepared mixtures. This is the only published method that directly quantifies the carbon of recent biological origin in biodiesel blends. Consequently, it robustly handles realistic chemical variability in biological source materials and provides unequivocal apportionment of renewable versus nonrenewable carbon in a sample fuel blend. Analysis of retail biodiesel blends acquired in 2006 in the United States revealed that inaccurate blending happens frequently. Only one out of ten retail samples passed the specifications that the United States Department of Defense requires for blends that are 20% biodiesel (v/v; referred to as B20).

Methyl-tert-hexyl ether and methyl-tert-octyl ether as gasoline oxygenates: anticipating widespread risks to community water supply wells

The widespread contamination of groundwater resources associated with methyl-tert-butyl ether (MtBE) use has prompted a search for replacement oxygenates in gasoline. Among the alternatives currently under development are higher methyl-tert-alkyl ethers, notably methyl-tert-hexyl ether (MtHxE) and methyl-tert-octyl ether (MtOcE). As was the case with MtBE, the introduction of these ethers into fuel supplies guarantees their migration into groundwater resources. In the present study, a screening-level risk assessment compared predicted well water concentrations of these ethers to concentrations that might cause adverse effects. A physicochemical model which has been successfully applied to the prediction of MtBE concentrations in community water supply wells (CSWs) was used to predict well water concentrations of MtHxE and MtOcE. The results indicate that these ethers are likely to contaminate water supply wells at slightly lower levels than MtBE as a result of migrating from leaking underground fuel tanks to CSWs. Because very little data is available on the physicochemical and environmental properties of MtHxE and MtOcE, estimation methods were employed in conjunction with the model to predict well water concentrations. Model calculations indicated that MtHxE and MtOcE will be present in many CSWs at concentrations approaching the concentrations that have caused widespread public health concern for MtBE. Based on these results and the possibility that MtHxE and MtOcE are potential carcinogens, testing of the toxicological properties of these ethers is recommended before they are used to replace MtBE in gasoline.

Disentangling oil weathering using GC x GC. 2. Mass transfer calculations

Hydrocarbon mass transfers to the atmosphere and water column drive the early weathering of oil spills and also control the chemical exposures of many coastal wildlife species. However, in the field, mass transfer rates of individual hydrocarbons to air and water are often uncertain. In the Part 1 companion to this paper, we used comprehensive two-dimensional gas chromatography (GC x GC) to identify distinct signatures of evaporation and dissolution encoded in the compositional evolution of weathered oils. In Part 2, we further investigate patterns of mass removal in GC x GC chromatograms using a mass transfer model. The model was tailored to conditions at a contaminated beach on Buzzards Bay, MA, after the 2003 Bouchard 120 oil spill. The model was applied to all resolved hydrocarbon compounds in the C11-C24 boiling range, based on their GC x GC-estimated vapor pressures and aqueous solubilities. With no fitted parameters, the model successfully predicted GC x GC chromatogram patterns of mass removal associated with evaporation, water-washing, and diffusion-limited transport. This enabled a critical field evaluation of the mass transfer model and also allowed mass apportionment estimates of hundreds of individual hydrocarbon compounds to air and water. Ultimately, this method should improve assessments of wildlife exposures to oil spill hydrocarbons.

Disentangling oil weathering using GC x GC. 1. chromatogram analysis

Historically, the thousands of compounds found in oils constituted an "unresolved complex mixture" that frustrated efforts to analyze oil weathering. Moreover, different weathering processes inflict rich and diverse signatures of compositional change in oil, and conventional methods do not effectively decode this elaborate record. Using comprehensive two-dimensional gas chromatography (GC x GC), we can separate thousands of hydrocarbon components and simultaneously estimate their chemical properties. We investigated 13 weathered field samples collected from the Bouchard 120 heavy fuel oil spill in Buzzards Bay, Massachusetts in 2003. We first mapped hydrocarbon vapor pressures and aqueous solubilities onto the compositional space explored by GC x GC chromatograms of weathered samples. Then we developed methods to quantitatively decouple mass loss patterns associated with evaporation and dissolution. The compositional complexity of oil, traditionally considered an obstacle, was now an advantage. We exploited the large inventory of chemical information encoded in oil to robustly differentiate signatures of mass transfer to air and water. With this new approach, we can evaluate mass transfer models (the Part 2 companion to this paper) and more properly account for evaporation, dissolution, and degradation of oil in the environment.

Biodegradation and environmental behavior of biodiesel mixtures in the sea: An initial study

Biodiesel, a mixture of fatty acid methyl esters (FAMEs) derived from animal fats or vegetable oils, is rapidly moving towards the mainstream as an alternative source of energy. However, the behavior of biodiesel, or blends of biodiesel with fossil diesel, in the marine environment have yet to be fully understood. Hence, we performed a series of initial laboratory experiments and simple calculations to evaluate the microbial and environmental fate of FAMEs. Aerobic seawater microcosms spiked with biodiesel or mixtures of biodiesel and fossil diesel revealed that the FAMEs were degraded at roughly the same rate as n-alkanes, and more rapidly than other hydrocarbon components. The residues extracted from these different microcosms became indistinguishable within weeks. Preliminary results from physical-chemical calculations suggest that FAMEs in biodiesel mixtures will not affect the evaporation rates of spilled petroleum hydrocarbons but may stabilize oil droplets in the water column and thereby facilitate transport.

Methyl tertiary hexyl ether and methyl tertiary octyl ether as gasoline oxygenates: assessing risks from atmospheric dispersion and deposition

Methyl tertiary hexyl ether (MtHxE) and methyl tertiary octyl ether (MtOcE) are currently being developed as replacement oxygenates for methyl tertiary butyl ether (MtBE) in gasoline. As was the case with MtBE, the introduction of these ethers into fuel supplies guarantees their introduction into the environment as well. In this study, a screening-level risk assessment was performed by comparing predicted environmental concentrations (PEC) of these ethers to concentrations that might cause adverse effects to humans or ecosystems. A simple box model that has successfully estimated urban air concentrations of MtBE was adapted to predict atmospheric concentrations of MtHxE and MtOcE. Expected atmospheric concentrations of these ethers were also estimated using the European Union System for the Evaluation of Substances (EUSES) multimedia fate model, which simultaneously calculates PECs in the various environmental compartments of air, water, soil, and sediment. Because little or no data are available on the physicochemical, environmental, and toxicological properties of MtHxE and MtOcE, estimation methods were used in conjunction with EUSES to predict both the PECs and the concentrations at which these ethers might pose a threat. The results suggest that these ethers would contaminate the air of a moderately sized U.S. city (Boston, MA) at levels similar to those found previously for MtBE. The risk assessment module in EUSES predicted risk characterization ratios of 10(-3) and 10(-2) for MtHxE and MtOcE, respectively, in Boston, and 10(-2) and 10(-1) in very large urban centers, suggesting that these ethers pose only a minimal threat to ecosystems at the anticipated environmental concentrations. The assessment also indicates that these compounds are possible human carcinogens and that they may be present in urban air at concentrations that pose an unacceptable cancer risk. Therefore, testing of the toxicological properties of these compounds is recommended before they replace MtBE in gasoline.