Atmospheric Chemistry of Toxic Contaminants. 3. Unsaturated Aliphatics: Acrolein,Acrylonitrile, Maleic Anhydride

Detailed mechanisms are outlined for the chemical reactions that contribute to In-situ formation and atmospheric removal of the unsaturated aliphatic contaminants acrolein, acrylonitrile, and maleic anhydride. In-situ formation of small amounts of acrolein and maleic anhydride may Involve the reaction of OH (and O₃) with 1,3-dienes and the reaction of OH with aromatic hydrocarbons, respectively. There is no known pathway for In-situ formation of acrylonitrile. Rapid removal of acrolein (half-life = less than one day) and of maleic anhydride (half-life = several hours) is expected from their rapid reactions with OH (major), O₃, and NO₃. These reactions lead to formaldehyde and glyoxal from acrolein and to dicarbonyls from maleic anhydride. Acrylonitrile is removed at a slower rate (half-life = 2–7 days) by reaction with OH, leading to formaldehyde and formyl cyanide.     

Atmospheric Chemistry of Toxic Contaminants. 2. Saturated Aliphatics: Acetaldehyde,Dioxane, Ethylene Glycol Ethers,Propylene Oxide

Detailed mechanisms are outlined for the chemical reactions that contribute to In-situ formation and atmospheric removal of the saturated aliphatic contaminants acetaldehyde, dioxane, ethylene glycol ethers (methyl, ethyl, n-butyl) and propylene oxide. In-situ formation Is of major Importance for acetaldehyde. In-situ removal Involves reaction with OH (all compounds) and, for acetaldehyde, photolysis and reaction with NO3. Acetaldehyde, dioxane, and the ethers are rapidly removed (half-lives of less than one day), leading to PAN (acetaldehyde) and to 2-oxodioxane and formaldehyde (dioxane). Reaction products of the glycol ethers include a large number of hydroxyesters, hydroxyacids, and hydroxycarbonyls. Propylene oxide reacts only slowly with OH, with an atmospheric half-life of 3-10 days, to yield formaldehyde, acetaldehyde, and PAN. Uncertainties in the reaction mechanisms for dioxane, the glycol ethers, and propylene oxide are discussed and include C-C vs C-0 bond scission in alkoxy radicals as well as alkoxy radical unimolecular decomposition vs reaction with oxygen.     

Atmospheric Chemistry of Toxic Contaminants 1. Reaction Rates and Atmospheric Persistence

Using structure-reactivity relationships between reaction rate constants and ionlzatlon potentials for structural homologues, estimates are presented for the rate constants of the reactions of ozone, the hydroxyl radical, and the nitrate radical with forty toxic air contaminants for which no or little data are available. These rate constants are in turn used to estimate the atmospheric persistence of saturated allphatics, unsaturated allphatics, and aromatic toxic organics. The corresponding atmospheric half-lives for removal by chemical reactions range from a few hours for the most reactive toxics (chloroprene, hexachlorocyclo-pentadiene, cresols, nitrosamines, maleic anhydride) to several months for the least reactive compounds (nitrobenzene, methyl bromide, phosgene).     

The Cl-initiated oxidation of CH3C(O)OCH=CH2, CH3C(O)OCH2CH=CH2, and CH2=CHC(O)O(CH2)3CH3 in the troposphere

Background, aim, and scope  

Unsaturated esters are emitted to the atmosphere from biogenic and anthropogenic sources, including those from the polymer industry. Little information exists concerning the atmospheric degradation of unsaturated esters, which are mainly initiated by OH radicals. Limited information is available on the degradation of alkenes by Cl atoms and almost no data exists for the reactions of unsaturated esters with Cl atoms. This data is necessary to assess the impact of such reactions in maritime environments where, under circumstances, OH radical- and Cl atom-initiated oxidation of the compounds can be important. Rate coefficients for the reactions of chlorine atoms with vinyl acetate, allyl acetate, and n-butyl acrylate have been determined at 298 ± 3 K and atmospheric pressure. The kinetic data have been used in combination with that for structurally similar compounds to infer the kinetic contributions from the possible reaction channels to the overall reaction rate.     

Products and mechanisms of the gas-phase reactions of OH radicals and O3 with 2-methyl-3-buten-2-ol

Products of the gas-phase reactions of OH radicals (in the presence of NO) and O₃ with the biogenic organic compound 2-methyl-3-buten-2-ol have been investigated using gas chromatography with flame ionization detection (GC-FID), combined gas chromatography–mass spectrometry (GC-MS), gas chromatography with Fourier transform infrared detection (GC-FTIR), in situ FT-IR spectroscopy and in situ atmospheric pressure ionization tandem mass spectrometry (API-MS/MS). Formaldehyde, 2-hydroxy-2-methylpropanal and acetone were identified from both the OH radical and O₃ reactions, glycolaldehyde and organic nitrate (s) were also observed from the OH radical reaction, and the OH radical formation yield from the O₃ reaction was measured. The formaldehyde, 2-hydroxy-2-methylpropanal, glycolaldehyde, acetone and organic nitrate yields from the OH radical reaction were 0.29±0.03, 0.19±0.07, 0.61±0.09, 0.58±0.04 and 0.05±0.02, respectively, and the formaldehyde, 2-hydroxy-2-methylpropanal and OH radical formation yields from the O₃ reaction were 0.29±0.03, 0.30±0.06 (0.47 from FT-IR measurements) and 0.19 (uncertain to a factor of 1.5), respectively. Acetone was also observed from the O₃ reaction, but appeared to be formed from secondary reactions. Reaction mechanisms are presented and discussed.     

Atmospheric Oxidation of Methyl Chloride Methylene Chloride,and Chloroform

Chlorinated atmospheric pollutants are presently receiving much attention because of expected chlorine-ozone interactions in the stratosphere.^1,2 The fully halogenated pollutants,such as CCI₃F, CCl₂F₂, and CCl₄, have no known removal processes that operate in the troposphere. These compounds are accumulating on a global scale, their atmospheric mixing ratios having reached about 2 X 10^-10 for CCl₂F₂,1.2 X 10^-10 for CCl₃F, and 9 X 10^-11 for CCl4.^3-5     

Atmospheric Chemistry of Toxic Contaminants. 6. Nitrosamines: Dialkyl Nitrosamines and Nitrosomorpholine

Detailed mechanisms are outlined for the reactions that contribute to in-sltu formation and atmospheric removal of dlmethylnitrosamine, diethylnitrosamlne, methyl-ethylnltrosamine, and nitrosomorphollne. In-sltu formation involves the rapid reaction of amines with the hydroxyl radical, leading to nltrosamlnes, nltramlnes, amides, and aldehydes. Nitrosamlne photolysis accounts for their rapid daytime removal, leading to amlno radicals whose atmospheric reactions are also discussed.     

Atmospheric chemistry of 1-methyl-2-pyrrolidinone

Rate constants for the atmospheric reactions of 1-methyl-2-pyrrolidinone with OH radicals, NO₃ radicals and O₃ have been measured at 296±2 K and atmospheric pressure of air, and the products of the OH radical and NO₃ radical reactions investigated. Using relative rate techniques, rate constants for the gas-phase reactions of OH and NO₃ radicals with 1-methyl-2-pyrrolidinone of (2.15±0.36)×10^-11 cm³ molecule^-1 s^-1 and (1.26±0.40)×10^-13 cm³ molecule^-1 s^-1, respectively, were measured, where the indicated errors include the estimated overall uncertainties in the rate constants for the reference compounds. An upper limit to the rate constant for the O₃ reaction of <1×10^-19 cm³ molecule^-1 s^-1 was also determined. These kinetic data lead to a calculated tropospheric lifetime of 1-methyl-2-pyrrolidinone of a few hours, with both the daytime OH radical reaction and the nighttime NO₃ radical reaction being important loss processes. Products of the OH radical and NO₃ radical reactions were analyzed by gas chromatography with flame ionization detection and combined gas chromatography–mass spectrometry. N-methylsuccinimide and (tentatively) 1-formyl-2-pyrrolidinone were identified as products of both of these reactions. The measured formation yields of N-methylsuccinimide and 1-formyl-2-pyrrolidinone were 44±12% and 41±12%, respectively, from the OH radical reaction and 59±16% and ∼4%, respectively, from the NO₃ radical reaction. Reaction mechanisms consistent with formation of these products are presented.     

Atmospheric degradation of alkylfurans with chlorine atoms: Product and mechanistic study

As part of a study on the oxidation mechanism of heterocyclic aromatic compounds, some aspects of the atmospheric chemistry of several alkyl derivatives of furan have been investigated. The aim of this work was to identify the products of the reactions of chlorine atoms with 2-methylfuran, 2-ethylfuran and 2,5-dimethylfuran. Experiments were performed in two different smog chambers at 296 ± 2 K and 1000 ± 20 mbar of synthetic air. The experimental investigation was carried out using in situ long-path FTIR absorption spectroscopy and both SPME-GC/FID-ECD and SPME-GC/MS as sampling and detection techniques. The major primary products from the addition reaction channel were 4-oxo-2-pentenoyl chloride and formaldehyde for the reactions of 2-methylfuran and 2,5-dimethylfuran; 4-oxo-2-hexenoyl chloride and acetaldehyde for the reaction of 2-ethylfuran and 5-chloro-2(5H)-furanone for the reactions of both 2-methylfuran and 2-ethylfuran. Other minor products were 4-oxo-2-pentenal, 4-oxo-2-hexenal and 3-hexene-2,5-dione for the 2-methylfuran, 2-ethylfuran and 2,5-dimethylfuran reactions, respectively. From the abstraction pathway, HCl, furfural, 2-acetylfuran, 5-methylfurfural, maleic anhydride and 5-hydroxy-2(5H)-furanone were detected. The formation of furfural, 2-acetylfuran and 5-methylfurfural confirmed the H-atom abstraction from the alkyl group of 2-methylfuran, 2-ethylfuran and 2,5-dimethylfuran, respectively. This mechanism was not observed in previous studies with OH and NO₃ radicals. A mechanism is proposed to explain the main reaction products observed. The observed products confirm that addition of Cl atoms to the double bond of the alkylfuran is the dominant reaction pathway.     

Condensed atmospheric photooxidation mechanisms for isoprene

Two condensed mechanisms for the atmospheric reactions of isoprene, which differ in the number of species used to represent isoprene's reactive products, have been developed for use in ambient air quality modehng. They are based on a detailed isoprene mechanism that has recently been developed and extensively evaluated against environmental chamber data. The new condensed mechanisms give very close predictions to those of the detailed mechanism for ozone, OH radicals, nitric acid, H₂O₂, formaldehyde, total PANS, and for incremental effects of isoprene on ozone formation in one day simulations. The effects of the condensations become somewhat greater in multi-day simulations, particularly in cases where NO₃ reactions are important at nighttime, but the ozone predictions are still very close. On the other hand, the SAPRC-90, RADM-2, and Carbon Bond IV isoprene mechanisms give quite different predictions of these quantities. It is recommended that the new mechanisms replace those currently used in airshed simulations where isoprene emissions are important.     

Atmospheric oxidation capacity in the summer of Houston 2006: Comparison with summer measurements in other metropolitan studies

Both similarities and differences in summertime atmospheric photochemical oxidation appear in the comparison of four field studies: TEXAQS2000 (Houston, 2000), NYC2001 (New York City, 2001), MCMA2003 (Mexico City, 2003), and TRAMP2006 (Houston, 2006). The compared photochemical indicators are OH and HO₂ abundances, OH reactivity (the inverse of the OH lifetime), HO_x budget, OH chain length (ratio of OH cycling to OH loss), calculated ozone production, and ozone sensitivity. In terms of photochemical activity, Houston is much more like Mexico City than New York City. These relationships result from the ratio of volatile organic compounds (VOCs) to nitrogen oxides (NO_x), which are comparable in Houston and Mexico City, but much lower in New York City. Compared to New York City, Houston and Mexico City also have higher levels of OH and HO₂, longer OH chain lengths, a smaller contribution of reactions with NO_x to the OH reactivity, and NO_x-sensitivity for ozone production during the day. In all four studies, the photolysis of nitrous acid (HONO) and formaldehyde (HCHO) are significant, if not dominant, HO_x sources. A problematic result in all four studies is the greater OH production than OH loss during morning rush hour, even though OH production and loss are expected to always be in balance because of the short OH lifetime. The cause of this discrepancy is not understood, but may be related to the under-predicted HO₂ in high NO_x conditions, which could have implications for ozone production. Three photochemical indicators show particularly high photochemical activity in Houston during the TRAMP2006 study: the long portion of the day for which ozone production was NO_x-sensitive, the calculated ozone production rate that was second only to Mexico City's, and the OH chain length that was twice that of any other location. These results on photochemical activity provide additional support for regulatory actions to reduce reactive VOCs in Houston in order to reduce ozone and other pollutants.     

Atmospheric fate of acrylic acid and acrylonitrile: Rate constants with Cl atoms and OH radicals in the gas phase

Rate coefficients for the reactions of hydroxyl radicals and chlorine atoms with acrylic acid and acrylonitrile have been determined at 298 K and atmospheric pressure. The decay of the organics was followed using a gas chromatograph with a flame ionization detector (GC-FID) and the rate constants were determined using a relative rate method with different reference compounds. Room temperature rate constants are found to be (in cm³ molecule⁻¹ s⁻¹): k₁(OH+CH₂CHC(O)OH)=(1.75±0.47)×10⁻¹¹, k₂(Cl+CH₂CHC(O)OH)=(3.99±0.84)×10⁻¹⁰, k₃(OH+CH₂CHCN)=(1.11±0.33)×10⁻¹¹ and k₄(Cl+CH₂CHCN)=(1.11±0.23)×10⁻¹⁰ with uncertainties representing ±2σ. This is the first kinetic study for these reactions under atmospheric pressure. The rate coefficients are compared with previous determinations taking into account the effect of pressure on the rate constants. The effect of substituent atoms or groups on the overall rate constants is analyzed in comparison with other unsaturated compounds in the literature. In addition, atmospheric lifetimes based on the homogeneous sinks of acrylic acid and acrylonitrile are estimated and compared with other tropospheric sinks for these compounds.     

Kinetics of OH-initiated oxidation of oxygenated organic compounds in the aqueous phase: new rate constants, structure–activity relationships and atmospheric implications

The kinetics of OH oxidation of several organic compounds of atmospheric relevance were measured in the aqueous phase. Relative kinetics were performed using various organic references and OH sources. After validation of the protocol, temperature-dependent rate constants for the reactions of OH radical with ethyl ter-butyl ether (, E_a/R=580 (±560) K), n-butyl acetate ( (±0.4)×10⁹ M⁻¹ s⁻¹, E_a/R=1000 (±200) K), acetone ( (±0.05)×10⁹ M⁻¹ s⁻¹, E_a/R=1400 (±500) K), methyl ethyl ketone (, E_a/R=1200 (±200) K), methyl iso-butyl ketone (, E_a/R=1200 (±300) K) and methylglyoxal (, E_a/R=1100 (±300) K) were determined. A non-Arrhenius behavior was found for phenol, in good agreement with the contribution of an OH addition to the mechanism, which also includes H-abstraction by OH radicals. Global rate constants of acetaldehyde, propionaldehyde, butyraldehyde and valeraldehyde were studied at 298 K only, as these compounds partly hydrate in the aqueous phase. All the obtained data (except those of phenol) complemented by literature data were used to investigate three methods to estimate rate constants for H-abstraction reactions of OH radicals in aqueous solutions when measured data were not available: Evans-Polanyi-type correlations, comparisons with gas-phase data, structure activity relationships (SAR). The results show that the SAR method is promising; however, the data set is currently too small to extend this method to temperatures other than 298 K. The atmospheric impact of aqueous phase OH oxidation of water-soluble organic compounds is discussed with the determination of their global atmospheric lifetimes, taking into account both gas- and aqueous-phase reactivities. The results show that atmospheric droplets can act as powerful photoreactors to eliminate soluble organic compounds from the atmosphere.     

Atmospheric alcohols and aldehydes concentrations measured in Osaka,Japan and in Sao Paulo,Brazil

The use of alcohol fuel has received much attention since 1980s. In Brazil, ethanol-fueled vehicles have been currently used on a large scale. This paper reports the atmospheric methanol, ethanol and isopropanol concentrations which were measured from May to December 1997, in Osaka, Japan, where alcohol fuel was not used, and from 3 to 9 February 1998, in Sao Paulo, Brazil, where ethanol fuel was used. The alcohols were determined by the alkyl nitrite formation reaction using gas chromatography (GC-ECD) analysis. The concentration of atmospheric alcohols, especially ethanol, measured in Sao Paulo were significantly higher than those in Osaka. In Osaka, the average concentrations of atmospheric methanol, ethanol, and isopropanol were 5.8±3.8, 8.2±4.6, and 7.2±5.9 ppbv, respectively. The average ambient levels of methanol, ethanol, and isopropanol measured in Sao Paulo were 34.1±9.2, 176.3.±38.1, and 44.2±13.7 ppbv, respectively. The ambient levels of aldehydes, which were expected to be high due to the use of alcohol fuel, were also measured at these sampling sites. The atmospheric formaldehyde average concentration measured in Osaka was 1.9±0.9 ppbv, and the average acetaldehyde concentration was 1.5±0.8 ppbv. The atmospheric formaldehyde and acetaldehyde average concentrations measured in Sao Paulo were 5.0±2.8 and 5.4±2.8 ppbv, respectively. The C₂H₅OH/CH₃OH and CH₃CHO/HCHO were compared between the two measurement sites and elsewhere in the world, which have already been reported in the literature. Due to the use of ethanol-fueled vehicles, these ratios, especially C₂H₅OH/CH₃OH, are much higher in Brazil than these measured elsewhere in the world.     

Relative rate measurements of reactions of unsaturated alcohols with atomic chlorine as a function of temperature

The kinetics of two structurally similar unsaturated alcohols, 3-butene-2-ol and 2-methyl-3-butene-2-ol (MBO232), with Cl atoms have been investigated for the first time, as a function of temperature using a relative method. As far as we know, the present work also provides the first value for 3-buten-2-ol. The coefficient at room temperature was also obtained for 2-propene-1-ol (allyl alcohol). The reactions were investigated using a 400 L Teflon reaction chamber coupled with gas chromatograph-coupled with flame-ionization detection (GC-FID) detection. The experiments were performed at atmospheric pressure and at temperatures between 256 and 298 K in air or nitrogen as the bath gas. The obtained kinetic data were used to derive the Arrhenius expressions, k_MBO232=(2.83±2.50)×10⁻¹⁴ exp (2670±249)/T, k_3-buten-2-ol=(0.65±1.60)×10⁻¹⁵ exp (3656±695)/T (in units of cm³ molecule⁻¹ s⁻¹). Finally, results and atmospheric implications are discussed and compared with the reactivity with OH and NO₃ radicals.     

Mechanism of OH-initiated atmospheric photooxidation of the organophosphorus insecticide (C2H5O)3PS

O,O,O-triethyl phosphorothioate ((C₂H₅O)₃PS, TEPT) is a widely used organophosphorus insecticide. TEPT may be released into the atmosphere where it can undergo transport and chemical transformations, which include reactions with OH radicals, NO₃ radicals and O₃. The mechanism of the atmospheric reactions of TEPT has not been fully understood due to the short-lifetime of its oxidized radical intermediates, and the extreme difficulty in detection of these species experimentally. In this work, we carried out molecular orbital theory calculations for the OH radical-initiated atmospheric photooxidation of TEPT. The profile of the potential energy surface was constructed, and the possible channels involved in the reaction are discussed. The theoretical study shows that OH addition to the PS bond and H abstractions from the CH₃CH₂O moiety are energetically favorable reaction pathways. The dominant products TEP and SO₂ arise from the secondary reactions, the reactions of OH-TEPT adducts with O₂. The experimentally uncertain dominant product with molecular weight 170 is mostly due to (C₂H₅O)₂P(S)OH and not (C₂H₅O)₂P(O)SH.     

Gas-phase ozonolysis of the monoterpenoids (S)-(+)-carvone, (R)-(−)-carvone, (−)-carveol, geraniol and citral

Biogenic emissions of volatile organic compounds (VOCs) play a fundamental role in atmospheric chemistry. Vegetation is the most abundant natural source of VOCs, while terpenoids, as limonene, α and β pinene and mircene, top the plants emission list. Much interest has been demonstrated in oxidation and photooxidation reactions of VOCs, particularly of monoterpenoids, owing to their diversity and to uncertainties regarding their mechanism of reaction. Quantification of primary carbonylic compounds, as well as of biradical reaction components, is highly relevant to the understanding of the major reactions. In this context, taking into account both structural factors and the fact that these compounds are found in the essential oils of plants typically found in Brazil and that they may be present in the atmosphere from emission by the plants, the monoterpenoids (S)-(+)-carvone, (R)-(−)-carvone, (−)-carveol, geraniol and citral (a mixture of the isomers geranial and neral) were selected for this study.The ozonolysis reactions of the monoterpenoids were carried out under dark conditions for all experiments, due to their photochemical reactivity. The analysis of the results lets us propose a mechanism by which these reactions occur. The observed results of the ozonolysis of S and R carvone suggest that the stereochemistry of asymmetric carbon does not affect either in the yields of both formaldehyde and of OH radicals produced in the reaction, or in the reactivity of these compounds, for which the rate constants were in the scale of 10⁻⁶ s⁻¹.We found that, in the (−)-carveol's cis and trans mixture, even though the hydroxyl in the axial position—in the case of trans-(C) and cis-(D′) isomers—favors the attack by the ozone molecule on the external double bond, thus increasing the mixture's reactivity , it affects the average production of formaldehyde. The presence of geraniol and citral led to the production of formaldehyde, propanone, glyoxal, methyl–glyoxal and cyclohexanone (OH radicals) as reaction products. The influence of an electron attractor group bonded to the carbon of the double bond, on the reactivity of the double bond, could not be observed in the case of citral, due to strong interference occurring in the instrument in all experiments with this monoterpenoid. For this reason, only the kinetics of geraniol was monitored .     

Principal component analysis to assess the efficiency and mechanism for enhanced coagulation of natural algae-laden water using a novel dual coagulant system

A novel dual coagulant system of polyaluminum chloride sulfate (PACS) and polydiallyldimethylammonium chloride (PDADMAC) was used to treat natural algae-laden water from Meiliang Gulf, Lake Taihu. PACS (Al_n(OH)_mCl_3n-m-2k(SO₄)_k) has a mass ratio of 10 %, a SO₄ ^2?/Al₃ ⁺ mole ratio of 0.0664, and an OH/Al mole ratio of 2. The PDADMAC ([C₈H₁₆NCl]_m) has a MW which ranges from 5?×?10⁵ to 20?×?10⁵ Da. The variations of contaminants in water samples during treatments were estimated in the form of principal component analysis (PCA) factor scores and conventional variables (turbidity, DOC, etc.). Parallel factor analysis determined four chromophoric dissolved organic matters (CDOM) components, and PCA identified four integrated principle factors. PCA factor 1 had significant correlations with chlorophyll-a (r?=?0.718), protein-like CDOM C1 (0.689), and C2 (0.756). Factor 2 correlated with UV₂₅₄ (0.672), humic-like CDOM component C3 (0.716), and C4 (0.758). Factors 3 and 4 had correlations with NH₃-N (0.748) and T-P (0.769), respectively. The variations of PCA factors scores revealed that PACS contributed less aluminum dissolution than PAC to obtain equivalent removal efficiency of contaminants. This might be due to the high cationic charge and pre-hydrolyzation of PACS. Compared with PACS coagulation (20 mg L^?1), the removal of PCA factors 1, 2, and 4 increased 45, 33, and 12 %, respectively, in combined PACS–PDADMAC treatment (0.8 mg L^?1?+?20 mg L^?1). Since PAC contained more Al (0.053 g/1 g) than PACS (0.028 g/1 g), the results indicated that PACS contributed less Al dissolution into the water to obtain equivalent removal efficiency.     

Atmospheric oxidation mechanisms of polychlorinated dibenzo-p-dioxins are different from those of benzene and dibenzofuran: a theoretical prediction

The reaction mechanisms of dibenzo-p-dioxin (DD) and 2,3,7,8-TCDD with OH radical have been studied using density functional theory calculations. Under the atmospheric conditions, ca 42% of DD + OH reaction proceeds as formation of DD − OH-β adduct, which will react with O₂ slowly; while the rest will proceed as formation of DD − OH-γ adduct, which will decompose to the substituted phenoxy radical P1 by the fused-ring C-O bond cleavage. For 2,3,7,8-TCDD + OH, the reaction will predominantly form the substituted phenoxy radical P2. The reaction mechanisms are drastically different from the peroxy mechanism for the atmospheric oxidations of benzene and dibenzofuran.     

Modeling in situ ozonation for the remediation of nonvolatile PAH-contaminated unsaturated soils

Mathematical models were developed to investigate the characteristics of gaseous ozone transport under various soil conditions and the feasibility of in situ ozone venting for the remediation of unsaturated soils contaminated with phenanthrene. On the basis of assumptions for the mass transfer and reactions of ozone, three approaches were considered: equilibrium, kinetic, and lump models. Water-saturation-dependent reactions of gaseous ozone with soil organic matter (SOM) and phenanthrene were employed. The models were solved numerically by using the finite-difference method, and the model parameters were determined by using the experimental data of Hsu [The use of gaseous ozone to remediate the organic contaminants in the unsaturated soils, PhD Thesis, Michigan State Univ., East Lansing, MI, 1995]. The transport of gas-phase ozone is significantly retarded by ozone consumption due to reactions with SOM and phenanthrene, in addition to dissolution. An operation time of 156 h was required to completely remove phenanthrene in a 5-m natural soil column. In actual situations, however, the operation time is likely to be longer than the ideal time because of unknown factors including heterogeneity of the porous medium and the distribution of SOM and contaminant. The ozone transport front length was found to be very limited (< 1 m). The sensitivity analysis indicated that SOM is the single most important factor affecting in situ ozonation for the remediation of unsaturated soil contaminated with phenanthrene. Models were found to be insensitive to the reaction mechanisms of phenathrene with either gas-phase ozone or dissolved ozone. More study is required to quantify the effect of OH* formation on the removal of contaminant and on ozone transport in the subsurface.