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.

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 C10-C30 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.

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.

Petroleum dynamics in the sea and influence of subsea dispersant injection during Deepwater Horizon.

During the Deepwater Horizon disaster, a substantial fraction of the 600,000-900,000 tons of released petroleum liquid and natural gas became entrapped below the sea surface, but the quantity entrapped and the sequestration mechanisms have remained unclear. We modeled the buoyant jet of petroleum liquid droplets, gas bubbles, and entrained seawater, using 279 simulated chemical components, for a representative day (June 8, 2010) of the period after the sunken platform's riser pipe was pared at the wellhead (June 4-July 15). The model predicts that 27% of the released mass of petroleum fluids dissolved into the sea during ascent from the pared wellhead (1,505 m depth) to the sea surface, thereby matching observed volatile organic compound (VOC) emissions to the atmosphere. Based on combined results from model simulation and water column measurements, 24% of released petroleum fluid mass became channeled into a stable deep-water intrusion at 900- to 1,300-m depth, as aqueously dissolved compounds (∼23%) and suspended petroleum liquid microdroplets (∼0.8%). Dispersant injection at the wellhead decreased the median initial diameters of simulated petroleum liquid droplets and gas bubbles by 3.2-fold and 3.4-fold, respectively, which increased dissolution of ascending petroleum fluids by 25%. Faster dissolution increased the simulated flows of water-soluble compounds into biologically sparse deep water by 55%, while decreasing the flows of several harmful compounds into biologically rich surface water. Dispersant injection also decreased the simulated emissions of VOCs to the atmosphere by 28%, including a 2,000-fold decrease in emissions of benzene, which lowered health risks for response workers.

Gas-Phase Ozone Reactions with a Structurally Diverse Set of Molecules: Barrier Heights and Reaction Energies Evaluated by Coupled Cluster and Density Functional Theory Calculations.

Reactions with ozone transform organic and inorganic molecules in water treatment systems as well as in atmospheric chemistry, either in the aqueous phase, at gas/particle interfaces, or in the gas phase. Computed thermokinetic data can be used to estimate the reactivities of molecules toward ozone in cases where no experimental data are available. Although the gas-phase reactivity of olefins with ozone has been characterized extensively in the literature, this is not the case for the richer chemistry of ozone with polar molecules, which occurs in the aqueous phase or in microhydrated environments. Here, we selected a number of model reactions with small molecules (ethene, ethyne, hydrogen cyanide, hydrogen chloride, ammonia, bromide, and trimethylamine) to study the accuracy of different quantum chemical methods for describing the reactivities of these molecules with ozone. We calculated benchmark electronic energies of gas-phase reactions of these systems with single-reference coupled cluster (CC) theory. These benchmark results for the binding energy in the van der Waals complex, the energy of the transition structure, and the reaction energy were estimated to be accurate within 1-2 kcal mol-1. Singlet oxygen (1O2) is a common product of ozone reactions. Coupled cluster calculations with up to perturbative quadruples (CCSDT( Q)) were needed to obtain reaction energies accurate within 1 kcal mol-1 when this species was involved. In (micro)hydrated environments or at interfaces, coupled cluster methods are prohibitively expensive in most cases. We tested the suitability of some contemporary density functional theory (DFT) methods to reproduce the benchmark electronic energy differences. Range-separated functionals were found to be promising candidates to estimate forward barrier heights, with LC-ωPBE rivaling the accuracy of CCSD( T). For energies of reaction, however, DFT methods exhibited large systematic errors, depending on their fraction of orbital exchange. This was found to worsen when 1O2 is a product, and no safe recommendation can be given for DFT reaction energies in such cases.

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.

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 (pKa) values of 1.5 ± 1 for NH3Cl+, 0.8 ± 1 for NH3Br+, 11.8 ± 1 for NHCl2, 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.

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.

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.

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.

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.

Chloramines, bromamines, and bromochloramines are halogen-containing oxidants that arise from the reaction of hypohalous acids with ammonia in water. Although relevant to both water disinfection chemistry and biochemistry, these molecules are difficult to study in the laboratory, and their thermochemical properties remain poorly established. We developed a benchmark level ab initio calculation protocol, termed TA14, adapted from the Weizmann theory and Feller-Peterson-Dixon approaches to determine the molecular structures and thermochemical properties of these compounds. We find that the halamine molecules are bound largely, and in some cases entirely, by electron correlation forces. This presumably explains their high reactivity as electrophilic oxidants. We provide computed heats of formation at 0 K (ΔfH(0)(0 K)) and at 298 K (ΔfH(0)(298 K)) and Gibbs free energies of formation at 298 K (ΔfG(0)(298 K)) for the 9 inorganic chloramines, bromamines, bromochloramines in gas phase. Based on comparisons to previous theoretical and experimental data for a set of 11 small molecules containing N, O, H, Cl, and Br, we propose uncertainties ranging from 1 to 3 kJ mol(-1) for computed thermodynamic properties of the halamines. Reported thermochemical data enable the determination of equilibrium constants for reactions involving halamines, opening possibilities for more quantitative studies of the chemistry of these poorly understood compounds.

Controlling factors in the rates of oxidation of anilines and phenols by triplet methylene blue in aqueous solution.

Anilines and phenols are structurally similar compound classes that both are susceptible to oxidation by excited state triplet sensitizers but undergo oxidation by different mechanisms. To gain an understanding of the factors that control the rate of oxidation of anilines and phenols by triplet excited states, a kinetic study was performed on the oxidation of substituted anilines and phenols by methylene blue. The rate constants of one-electron transfer from anilines to triplet state methylene blue and their dependence on the reaction free energy are well fit to a Sandros-Boltzmann model. The observed rate constants are also well modeled when aniline oxidation potentials derived computationally are used. For phenols, the proton-coupled electron transfer rate constants were found to correlate primarily with O-H bond dissociation free energy and secondarily with phenol pKa. Rate constants for phenols could be modeled using computed bond dissociation free energies. These results provide a basis for predicting aniline and phenol oxidation rates, which could be valuable, for example, in assessing the likely persistence and fate of aniline- and phenol-based aqueous environmental pollutants.

Composition and fate of gas and oil released to the water column during the Deepwater Horizon oil spill.

Quantitative information regarding the endmember composition of the gas and oil that flowed from the Macondo well during the Deepwater Horizon oil spill is essential for determining the oil flow rate, total oil volume released, and trajectories and fates of hydrocarbon components in the marine environment. Using isobaric gas-tight samplers, we collected discrete samples directly above the Macondo well on June 21, 2010, and analyzed the gas and oil. We found that the fluids flowing from the Macondo well had a gas-to-oil ratio of 1,600 standard cubic feet per petroleum barrel. Based on the measured endmember gas-to-oil ratio and the Federally estimated net liquid oil release of 4.1 million barrels, the total amount of C(1)-C(5) hydrocarbons released to the water column was 1.7 10(11) g. The endmember gas and oil compositions then enabled us to study the fractionation of petroleum hydrocarbons in discrete water samples collected in June 2010 within a southwest trending hydrocarbon-enriched plume of neutrally buoyant water at a water depth of 1,100 m. The most abundant petroleum hydrocarbons larger than C(1)-C(5) were benzene, toluene, ethylbenzene, and total xylenes at concentrations up to 78 µg L(-1). Comparison of the endmember gas and oil composition with the composition of water column samples showed that the plume was preferentially enriched with water-soluble components, indicating that aqueous dissolution played a major role in plume formation, whereas the fates of relatively insoluble petroleum components were initially controlled by other processes.

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.