GEOPHYSICAL RESEARCH LETTERS, VOL. 30, NO. 16, 1862, doi:10.1029/2003GL017933, 2003 Oxygenated volatile organic chemicals in the oceans: Inferences and implications based on atmospheric observations and air-sea exchange models H. B. Singh,1 A. Tabazadeh,1 M. J. Evans,2 B. D. Field,2 D. J. Jacob,2 G. Sachse,3 J. H. Crawford,3 R. Shetter,4 and W. H. Brune,5 Received 10 June 2003; revised 21 July 2003; accepted 28 July 2003; published 28 August 2003. [1] Airborne measurements of oxygenated volatile organic chemicals (OVOC), OH free radicals, and tracers of pollution were performed over the Pacific during Winter/ Spring of 2001. We interpret atmospheric observations of acetaldehyde, propanal, methanol, and acetone with the help of a global 3-D model and an air-sea exchange model to assess their oceanic budgets. We infer that surface waters of the Pacific are greatly supersaturated with acetaldehyde and propanal. Bulk surface seawater concentration of 7 nM (10À9 mol LÀ1) and 2 nM and net fluxes of 1.1  10À12 g cmÀ2 sÀ1 and 0.4  10À12 g cmÀ2 sÀ1 are calculated for acetaldehyde and propanal, respectively. Large surface seawater concentrations are also estimated for methanol (100 nM) and acetone (10 nM) corresponding to an undersaturation of 6% and 14%, and a deposition velocity of 0.08 cm sÀ1 and 0.10 cm sÀ1, respectively. These data imply a large oceanic source for acetaldehyde and propanal, and a modest sink for methanol and acetone. Assuming a 50– 100 meter mixed layer, an extremely large oceanic reservoir of OVOC, exceeding the atmospheric reservoir by an order of magnitude, can be inferred to be present. Available seawater data are both preliminary and extremely limited but indicate rather low bulk OVOC concentrations and provide no support for the existence of a large oceanic reservoir. We speculate on the causes and implications of these INDEX TERMS: 0312 Atmospheric Composition findings. and Structure: Air/sea constituent fluxes (3339, 4504); 0322 Atmospheric Composition and Structure: Constituent sources and sinks; 0365 Atmospheric Composition and Structure: Troposphere—composition and chemistry. Citation: Singh, H. B., A. Tabazadeh, M. J. Evans, B. D. Field, D. J. Jacob, G. Sachse, J. H. Crawford, R. Shetter, and W. H. Brune, Oxygenated volatile organic chemicals in the oceans: Inferences and implications based on atmospheric observations and air-sea exchange models, Geophys. Res. Lett., 30(16), 1862, doi:10.1029/ 2003GL017933, 2003. [Singh et al., 2001; Wisthaler et al., 2002 and references there in]. These chemicals are expected to play an important role in the chemistry of the atmosphere by sequestering nitrogen oxides (NOx) in the form of PAN and by providing HOx free radicals in critical regions of the atmosphere [Singh et al., 1995; Wennberg et al., 1998]. OVOC are believed to have large terrestrial sources (%500 TgC yÀ1) but our quantitative knowledge about these is rudimentary [Fall, 1999]. Attempts to reconcile atmospheric observations with models have led to the inference that oceans must provide a sizable source of acetone [de Laat et al., 2001; Jacob et al., 2002]. Galbally and Kirstine [2002] have suggested that an extremely large oceanic reservoir of methanol may be present. It is also known that sunlight initiated reactions can decompose organic matter in surface oceans to form a variety of oxygenated chemicals [Zhou and Mopper, 1997]. The role of oceans in the global budgets of OVOC remains largely unexplored and there is little reliable data available. Here we investigate and assess the oceanic composition and budget of acetaldehyde (ethanal, CH3CHO), propanal (propionaldehyde, C2H5CHO), methanol (CH3OH), and acetone (propanone, CH3COCH3) based on airborne measurements over the Pacific and models of air-sea exchange. These OVOC are selected because they are among the most abundant and there is reason to believe that they are globally ubiquitous. 2. Results and Discussion [3 ] During the TRACE-P Spring 2001 experiment, the NASA/DC-8 flying laboratory was used to measure a large number of OVOC in the Pacific troposphere to an altitude of 12 km (Latitude 10 – 45°N; Longitude 100 – 230°E). Complementing these were in-situ measurements of hydroxyl radicals (OH), a key constituent that controls the atmospheric removal rates of OVOC, and multiple tracers of pollution (e.g., CO, C2Cl4). Additionally, spectral radiometers were employed to determine photolytic loss rates of the OVOC of interest in the troposphere. More information about measurement methods, mission payload, and flight profiles can be found in the overview paper by Jacob et al. [2003] and references there in. [4] All of the data discussed here were collected over open oceans some 500 km or more away from the Asian continent. A ‘‘pollution filter’’ based on CO and C2Cl4 distribution was employed to mitigate the effects of pollution and is described in more detail in Singh et al. [2003]. This filter eliminated all major pollution influences and resulted in mean tropospheric mixing ratios of 102(±20) ppb/CO 1. Introduction [2] In recent years it has become evident that significant concentrations of a large number of oxygenated organic chemicals (OVOC) are present in the global troposphere 1 2 3 NASA Ames Research Center, Moffett Field, CA, USA. Harvard University, Cambridge, MA, USA. NASA Langley Research Center, Hampton, VA, USA. 4 National Center for Atmospheric Research, Boulder, CO, USA. 5 Pennsylvania State University, University Park, PA, USA. This paper is not subject to U.S. copyright. Published in 2003 by the American Geophysical Union. ASC 13 - 1 ASC 13 - 2 SINGH ET AL.: OXYGENATED VOLATILE ORGANIC CHEMICALS IN THE OCEANS Figure 1. Measured and modeled mixing ratios of selected OVOC. Shaded area in (a) shows the MBL mixing ratio range from many regions of the globe. In the model run, oceans provide (in Tg yÀ1) 13 acetone source, 0 aldehyde source, and 15 methanol sink. Data altitudes are shifted by À0.25 km for clarity. and 3(±1) ppt/C2Cl4, and is assumed to represent nearbackground conditions. All analysis in this paper is based on this filtered data set. The observed vertical distribution of the four selected OVOC under these clean conditions is shown in Figure 1. [5] It is evident from Figure 1 that measured mixing ratios of acetone and methanol are lower in the marine boundary layer (MBL) compared to the free troposphere, while the reverse is the case for aldehydes. In Table 1 we present mixing ratio data in the MBL (0 – 2 km, dMBL) and over it (2 – 4 km) for the OVOC of interest in this study. Figure 1 also shows these mixing ratios simulated along the TRACE-P flight tracks by the GEOS-CHEM global tropospheric model. More details about the GEOS-CHEM model can be found in Jacob et al. [2002]. An in-depth analysis of the atmospheric distribution and global budgets of the large suite of OVOC measured during TRACE-P is being presented elsewhere [Singh et al., 2003]. 2.1. Ocean Flux [6] Oceanic fluxes are estimated assuming a steady state behavior in which the entrainment at the top of the MBL, production and losses within the MBL, and net oceanic fluxes are in balance. The various flux terms are described in equations 1– 3. Ft ¼ Ve  ÁC ð1Þ FMBL ¼ FMBL ðproductionÞ À FMBL ðlossÞ ð2Þ Fo ¼ ÀFMBL À Ft ð3Þ [7] Where Ft is the flux at the top of the MBL, FMBL is the net flux associated with chemical production and loss of OVOC in the MBL, and Fo is the calculated net oceanic flux. A mean entrainment velocity (Ve) of 0.4 cm sÀ1 is used in all these calculations [Faloona et al., 2003]. To avoid biases due to scatter in the data, gradient across the top of the MBL (ÁC) are determined by using median concentration differences between the 0 – 2 km and 2 – 4 km layers. The calculated value of Ft for each OVOC is presented in Table 1. Table 1. OVOC Mixing Ratios and Their Oceanic Flux Under Pristine Conditions Parametersa Mixing ratio (0 – 2 km) Mixing ratio (2 – 4 km) ÁC (ppt) Loss rate in MBL (daysÀ1) Ft (g cmÀ2 sÀ1) Fmbl(g cmÀ2 sÀ1) Fo(g cmÀ2 sÀ1) Oceanic source (Tg yÀ1)b a b CH3CHO 205 (204 ± 40)c 171 (173 ± 45) À34 1.43 À0.0  10À12 À1.1  10À12 +1.1  10À12 +125 C2H5CHO 60 (68 ± 24) 54 (60 ± 21) À6 1.43 À0.0  10À12 À0.4  10À12 +0.4  10À12 +45 CH3OH 563 (575 ± 211) 744 (840 ± 258) +181 0.12 +1.0  10À13 À0.4  10À13 À0.6  10À13 À8 CH3COCH3 437 (466 ± 97) 636 (642 ± 207) +199 0.05 +1.9  10À13 À0.7  10À13 À1.2  10À13 À14 An entrainment velocity of 0.4 cm sÀ1 and a mean MBL OH of 1.5  106 molec. cmÀ3 is used. Assumes that calculated oceanic flux is globally applicable. c Median mixing ratio (Mean ± 1s). See text for details. SINGH ET AL.: OXYGENATED VOLATILE ORGANIC CHEMICALS IN THE OCEANS ASC 13 - 3 [8] Since OVOC are both produced and destroyed in the atmosphere, FMBL is a measure of net chemical flux in the MBL. Loss rates in the MBL are calculated assuming that chemical loss is due to reaction with OH and photolysis {FMBL (loss) = CMBL  dMBL  (ki [OH]MBL + ji)}. The OH reaction rate constants (ki) for all species in this study were taken from the compilations of Atkinson et al. [2002]. Diurnally averaged OHMBL concentrations of 1.5  106 molec. cmÀ3 were derived from the in-situ OH measurements on board the DC-8 and were also in good agreement with GEOS-CHEM model simulations. Photolytic loss rates ( ji) were directly calculated from spectral radiometric measurements. Combined first order loss rates are shown in Table 1. [9] The OVOC production in the MBL is primarily determined by the concentrations of specific precursors and OH radicals. In the case of acetone and methanol, the main photochemical formation pathways involve oxidation of propane/i-butane and methane, respectively. An acetone yield of 0.8 from the propane/OH reaction was used [Atkinson et al., 2002]. A methanol mean molar yield of 0.035 from methane oxidation was derived from the GEOSCHEM model for the 0 –2 km region. Median MBL mixing ratios of propane (+i-butane) and methane were 76 ppt and 1772 ppb, respectively. For aldehydes, which are intermediate products of many hydrocarbons, we have used the chemical production rates as calculated by the GEOSCHEM model. The net MBL flux values due to chemistry (FMBL) are summarized in Table 1. [10] The net oceanic flux estimates (Fo) are positive (source to the atmosphere) for aldehydes and negative (sink) for methanol and acetone (Table 1). It is also evident that aldehyde fluxes are extremely large in comparison to methanol and acetone. These flux estimates are valid for the Pacific region and are not likely to be representative of the global oceans. As a first approximation however, if Table 1 fluxes are extrapolated to global scales, an oceanic source/sink (Tg yÀ1) of 125, 45, À8, and À14 can be calculated for CH 3 CHO, C 2 H 5 CHO, CH 3 OH, and CH3COCH3, respectively. The finding of an oceanic sink for CH3OH and CH3COCH3 is contrary to previous suggestions of a sizable oceanic source [de Laat et al., 2001; Jacob et al., 2002; Heikes et al., 2002]. As can be seen from Figure 1d, a 13 Tg yÀ1 net oceanic source of CH3COCH3 in the GEOS-CHEM model greatly over predicts MBL mixing ratios. Similarly, an assumed methanol oceanic sink of 15 Tg yÀ1 provided a reasonable fit to the observations (Figure 1c). While these results must be considered preliminary at this stage, they provide the first indications of the important role that oceans could play in the global budgets of OVOC. [11] Here it is appropriate to add that CH3OH and CH3COCH3 are relatively long-lived (%2 weeks) and their global concentration fields are reasonably well established [Singh et al., 1995; 2001; de Laat et al., 2001; Wisthaler et al., 2002]. While similar data for short-lived (<1 day) aldehydes are not available, CH3CHO mixing ratios in the MBL of 70– 400 ppt, indicated by the shaded area in Figure 1a, have been reported from many regions of the globe in both hemispheres [Singh et al., 2003 and references there in]. No data for C2H5CHO are available, but its mixing ratios are found to be strongly correlated with CH3CHO (R2 = 0.9) and it is expected to behave similarly. The ensemble of these observations supports the view that substantial aldehyde concentrations are present throughout the globe. We note that the large oceanic emissions of aldehydes estimated in Table 1 are insufficient to account for their mixing ratios in the free troposphere [Singh et al., 2003]. 2.2. Sea Water Concentrations [12] To explore oceanic budgets of these OVOC, we use a two-film model of air-sea exchange described by Liss and Slater [1974] and further discussed by Donelan and Wanninkhof [2001]. Based on this model, equations 4 – 5 were used to derive surface seawater concentration and saturation (1 + S) necessary to maintain the oceanic flux, Fo. ÁÀ1 À À Á À Á Fo ¼ K1 C1 À Cg =H ¼ 1=ak1 þ 1=Hkg C1 À Cg =H ð4Þ S ¼ HC1 =Cg À 1 ¼ HFo =K1 Cg ð5Þ where Cl and Cg are concentrations in the bulk liquid and gas phases, H (Csg/Csl) is the intrinsic dimensionless Henry’s law constant, kl and kg are exchange constants for liquid and gas phases, and a is an enhancement factor due to species reaction in solution. For both aldehydes a hydration constant {[RCH(OH)2]/[RCHO]} of 1.4 was used [Zhou and Mopper, 1997 and references there in]. Dimensionless intrinsic Henry’s Law constants shown in Table 2 were obtained from the data compiled by Sander [1999] and corrected for the reduced solubility in sea water (Â1.2) and in case of aldehydes for hydration effects (Â2.4). A liquid phase exchange constant (kl) of 11 cm hÀ1 and a gas phase exchange constant (kg) of 6912 (18/MW)1/2 cm hÀ1 were used based on recommended values for a mean surface wind speed of 6 m sÀ1 [Donelan and Wanninkhof, 2001, Asher, 1997]. These results are relatively insensitive to a, which is calculated to be 2.0 for aldehydes and 1.0 for other OVOC according to Liss and Slater [1974]. [13] As can be seen from Table 2, both gas phase (rg = 1/Hkg) and liquid phase (rl = 1/akl) resistances are important to the air-sea exchange of OVOC. Air-sea exchange models have been predominantly validated for sparingly soluble and relatively inert species (such as O2, N2, N2O, and CO2) where the liquid phase resistance dominates. Reactive organic chemicals such as OVOC, as well as chemicals for which resistance to transfer resides principally in the gas phase have not been extensively studied. Given these caveats, equation 4 is used to calculate bulk surface seawater concentrations (10À9 mol LÀ1 or nM) of 7, 2, 100, 10 for CH3CHO, C2H5CHO, CH3OH, and CH3COCH3, respectively (Table 2). A comparably large concentration of the hydrated ion of CH3CHO and C2H5CHO would also exist in seawater. These results indicate a very large seawater super-saturation of aldehydes. In comparison, Pacific waters appear to be under-saturated in methanol (6%) and acetone (14%) and provide a modest net sink. This sink corresponds to a deposition velocity (SKl/H) of 0.08 cm sÀ1 for methanol and 0.10 cm sÀ1 for acetone. [14] Bulk surface concentration (Cl) typically represents the top 50– 100 m well-mixed region of the oceans [Kantha and Clayson, 2003]. Assuming a 50 m mixed layer, the oceanic loading of these species over the Pacific can be calculated to be 2– 25 times their atmospheric loading ASC 13 - 4 SINGH ET AL.: OXYGENATED VOLATILE ORGANIC CHEMICALS IN THE OCEANS Table 2. Seawater OVOC Concentrations and Reservoir Estimated by Using an Air-Sea Exchange Model Chemicals CH3CHO C2H5CHOg CH3OH CH3COCH3 a b g Ha 9.1 9.1 2.2 1.5  10  10À3  10À4  10À4 À3 Klb (cm hÀ1) 14.3 13.5 1.0 3.9 rg/(rg + rl)c (%) 36 39 91 65 Fo (g cmÀ2 sÀ1) +1.1  10 +0.4  10À12 À0.6  10À13 À1.2  10À13 À12 Cld (nM) 7 2 100 10 Saturation (%) 800 800 94 86 Re 7 4 25 2 O-reservoirf (Tg) 4 1 50 10 Dimensionless Henry’s Law constants corrected for salt and hydration effects. Total air-sea exchange transfer velocity. c Percent resistance in gas phase. d Surface seawater concentrations (nM = 10À9 mol LÀ1). e Ratio of the oceanic to the atmospheric burden over the Pacific assuming a 50 m mixed layer depth. f Global oceanic reservoir assuming R is maintained and mixing ratios in the SH are 80% of NH. g For these aldehydes, a hydrated ion concentration comparable to the free molecule will also be present in seawater. (R in Table 2). Adjusting for hemispheric differences and extrapolating to global scales, a very large oceanic reservoir of these species (%65 Tg) can be calculated to be present (Table 2). In their assessment of the global budget of methanol, Galbally and Kirstine [2002] assumed saturation in the top 80 m to calculate an oceanic methanol reservoir some 66 times larger than the atmospheric reservoir. Heikes et al. [2002], in an independent review of the methanol budget, estimate a methanol oceanic source of 0 – 80 Tg yÀ1 but did not imply a reservoir. [15] Is there evidence for the large oceanic concentrations and the implied reservoir of OVOC estimated in Table 2? Two studies provide some preliminary indications. In the first study, Zhou and Mopper [1997] used a wet chemical technique to measure CH3CHO and CH3COCH3 in coastal and open ocean areas east of Florida. They report a bulk seawater concentration range of 1 – 15 nM for CH3CHO and 3 – 15 nM for CH3COCH3. Bulk water concentrations in two open ocean samples collected from a location some 100 km east of Bahamas were 1 nM for CH3CHO and 3 nM for CH3COCH3. However, the concentrations in the surface microlayer (%102 mm) were 10– 20 times larger. In a more recent investigation, using a sensitive mass spectrometric technique, Williams et al. [2003] report bulk seawater concentrations of <1 nM for CH3OH and CH3COCH3 in the open tropical Atlantic ocean. They also find that these low concentrations are nearly uniform in the top 50 m of the ocean. The atmospheric mixing ratios measured by them were similar to those in Table 1. The available oceanic data at the moment are preliminary and too limited to draw any firm conclusions. However, they do indicate that bulk seawater concentrations in remote oceans for several of these OVOC may be much lower than predicted by current air-sea exchange models (Table 2). These limited measurements provide no direct support for the presence of a large OVOC oceanic reservoir. [16] The super-saturation of aldehydes and under-saturation of acetone are indicative of the fact that the OVOC are being produced and destroyed in the oceans. We know little about these mechanisms. Formation processes could involve photo or biodegradation of organic matter. It is also possible that high OVOC concentrations inferred in Table 2 are present only in surface microlayers over the ocean [Zhou and Mopper, 1997]. For species such as aldehydes, where oceans are a source, a steady flux can be maintained from deep reservoirs but thinner layers must be continuously replenished if they are not to be quickly depleted. This almost certainly requires a mechanism for near-continuous OVOC formation in seawater. If these processes are driven by the availability of sunlight, then strong diurnal variations should be expected. Similarly, loss processes in seawater could include rapid consumption by bacteria, typically present in the top few meters, as well as photodegradation. At the moment it is not possible to reconcile the expected results from air-sea exchange models with the extremely limited and somewhat preliminary available oceanic data. 3. Conclusions [17] In this study we have used atmospheric measurements of OVOC and a standard air-sea exchange model to explore their oceanic budgets. Current models suggest that oceans should contain a very large reservoir of OVOC and may provide a net source in some cases and a sink in others. However, the reliability of these models for reactive organic species remains untested. Available oceanic data are extremely limited and uncertain. Mechanisms for the formation and destruction of OVOC in seawater are largely unknown. There is clearly an urgent need for more research on the abundance and fate of OVOC in the oceans. The development of reliable air-sea exchange models for reactive and moderately soluble organic species is also a worthwhile long-term objective. [18] Acknowledgments. This research was funded by the NASA Tropospheric Chemistry Program and the NSF Atmospheric Chemistry Program. Constructive discussions with R. Wanninkhof of NOAA and L. Iraci of Ames are much appreciated. References Asher, N. V., The sea-surface microlayer and its effect on global air-sea gas transfer, in The sea surface and global change, edited by P. S. Liss and R. A. Duce, 251 – 285, Cambridge Univ. Press, Cambridge, 1997. Atkinson, R., et al., IUPAC evaluated kinetic data, http://www.iupac-kinetic. ch.cam.ac.uk/, 2002. de Laat, A. T. J., J. A. de Gouw, J. Lelieveld, and A. Hansel, Model analysis of trace gas measurements and pollution impact during INDOEX, J. Geophys. Res., 106, 28,469 – 28,480, 2001. Donelan, M. A., and R. Wanninkhof, Gas Transfer at Water Surfaces – concepts and issues, in Gas Transfer at Water Surfaces, edited by M. Donelan, W. Drennan, E. Saltzman, and R. Wanninkhof, 1 – 10, AGU, Geophysical Monograph 127, Washington, 2001. Faloona, I., et al., Observations of entrainment in eastern Pacific marine stratocumulus using three conserved scalars, J. Atmos. Sci., submitted, 2003. Galbally, I. E., and W. Kirstine, The production of methanol by flowering plants and the global cycle of methanol, J. Atmos. Chem., 43(3), 195 – 229, 2002. Heikes, B. G., et al., Atmospheric methanol budget and ocean implication, Global Biogeochem. Cycles, 16(4), 80-1 to 80-13, 2002. SINGH ET AL.: OXYGENATED VOLATILE ORGANIC CHEMICALS IN THE OCEANS Jacob, D. J., et al., Atmospheric budget of acetone, J. Geophys. Res., 107, doi:10.1029/2001JD000694, 2002. Jacob, D. J., et al., The Transport and Chemical Evolution over the Pacific (TRACE-P) mission: Design, execution and overview of results, J. Geophys. Res., in press, 2003. Kantha, L. H., and C. A. Clayson, Ocean mixed layer, in Encyclopedia of Atmospheric Sciences, edited by J. R. Holton, J. A. Curry, and J. A. Pyle, 291 – 298, Academic Press, London, 2003. Liss, P. S., and P. G. Slater, Flux of gases across the air-sea interface, Nature, 247, 181 – 184, 1974. Sander, R., Compilation of Henry’s law constants for inorganic and organic species of potential importance in environmental chemistry, http://www. mpch-mainz.mpg.de/~sander/res/henry.html, version 3, 1999. Singh, H. B., M. Kanakidou, P. J. Crutzen, and D. J. Jacob, High concentrations and photochemical fate of oxygenated hydrocarbons in the global troposphere, Nature, 378, 50 – 54, 1995. Singh, H. B., et al., Evidence from the Pacific troposphere for large global sources of oxygenated organic compounds, Nature, 410, 1078 – 1081, 2001. Singh, H. B., et al., Analysis of the atmospheric distribution, sources, and sinks of oxygenated volatile organic chemicals (OVOC) based on measurements over the Pacific during TRACE-P, J. Geophys. Res., submitted, 2003. ASC 13 - 5 Wennberg, P. O., et al., Hydrogen radicals, nitrogen radicals, and the production of ozone in the upper troposphere, Science, 279, 49 – 53, 1998. Williams, J., et al., Atmospheric and oceanic measurements of organic species from the tropical Atlantic ocean, unpublished data, Max-Planck Institute, Mainz, Germany, 2003. Wisthaler, A., A. Hansel, R. R. Dickerson, and P. J. Crutzen, Organic trace gas measurements by PTR-MS during INDOEX 1999, J. Geophys. Res., 107(D19), doi:10.1029/2001JD000576, 2002. Zhou, X., and K. Mopper, Photochemical production of low-molecularweight carbonyl compounds in seawater and surface microlayer and their air-sea exchange, Mar. Chem., 56, 201 – 213, 1997. ÀÀÀÀÀÀÀÀÀÀÀ ÀÀÀÀÀÀÀÀÀÀÀ H. B. Singh and A. Tabazadeh, NASA Ames Research Center, Moffett Field, CA 94035, USA. M. J. Evans, B. D. Field, and D. J. Jacob, Harvard University, Cambridge, MA 02138, USA. G. Sachse and J. H. Crawford, NASA Langley Research Center, Hampton, VA 23665, USA. R. Shetter, National Center for Atmospheric Research, Boulder, CO 80307, USA. W. H. Brune, Pennsylvania State University, University Park, PA 16802, USA.