Abstract

Ozone fluxes, derived from gradient measurements in Northeast Atlantic coastal waters, were observed to depend on both tide height and solar radiation. Peak ozone fluxes of μg m−2 s−1 occurred during low-tide conditions when exposed microalgae fields contributed to the flux footprint. Additionally, at mid-to-high tide, when water surfaces contribute predominantly to the flux footprint, fluxes of the order of μg m−2 s−1 were observed. Considering only fluxes over water covered surfaces, and using an advanced ozone deposition model that accounts for surface-water chemistry enhancing the deposition sink, it is demonstrated that a photochemical enhancement reaction with dissolved organic carbon (DOC) is required to explain the enhanced ozone deposition during daylight hours. This sink amounts to an ozone loss rate of up to 0.6 ppb per hour under peak solar irradiance and points to a missing sink in the marine boundary layer ozone budget.

1. Introduction

Tropospheric ozone has negative impacts on human health and vegetation [1] and also acts as an important greenhouse gas. Recent trend analysis in background North Atlantic air at the Mace Head research station in western Ireland shows background ozone levels to be on the increase by 0.16 ppb per annum (1988 to 2007) [2]. Mace Head is a clean marine site, representing air transported over the North Atlantic towards Europe in the marine boundary layer. Consequently, rising ozone levels at Mace Head represent rising background ozone levels. The rise in background ozone concentrations at Mace Head cannot be explained by changes in anthropogenic emission patterns, and model studies have not been able to recreate observed trends [3, 4]; hence, there is a need to investigate processes influencing background ozone concentrations. Jenkin [5] identified gradual hemispheric ozone increase due to global-scale effects, influencing the ozone advected into Western Europe over the Atlantic as one of the three major influences on locally measured ozone data in the UK. The rise in background ozone concentration constitutes a significant contribution to regional ozone pollution, and so is of significant concern for pollution control policy makers. Hence, in order to assess the future trends in background ozone levels, the intercontinental transport of ozone should be well represented in assessment models as well as any relevant removal processes, one of which is the dry deposition of ozone to the ocean.

The deposition of ozone to surface waters via physical uptake (solubility, turbulent mixing) and chemical reactions in the surface layer [6] is an ozone removal process likely to influence surface ozone concentrations during intercontinental transport. From field measurements [7, 8], deposition velocity of ozone to the ocean is known to span the range of 0.53 to 1.0 mm s−1. Coleman et al. [9] describe the inclusion of the advanced ozone dry deposition parameterisation of Fairall et al. [10] in the regional climate model REMOTE. The Fairall ozone dry deposition parameterisation scheme accounts for turbulent and chemical enhancement of ozone deposition to water surfaces. Within their study, Coleman et al. [9] expanded the Fairall scheme to include ozone depletion via reactions with iodide and dissolved organic matter (DOM). The parameterisation was also adapted to account for variability of ozone diffusivity, solubility, and reactivity with sea surface temperature and to estimate emissions of volatile organoiodine (VOI) vapours following reactions of ozone with iodide and DOC (e.g., [11]). In the present study, the Fairall parameterisation was further developed to parameterise ozone-DOC reactions with a second-order reaction rate empirically derived from the laboratory work of [12]. Box model results of the updated Fairall parameterisation are compared to in situ ozone gradient fluxes measured at Mace Head.

2. Measurement Methods

Ozone gradient flux measurements were conducted on a 22 m tower at Mace Head atmospheric research facility [13] on the west coast of Ireland (53°20′N, 9°54′W). The shoreline to the westward side of the 22 m tower is inhomogeneous, rocky, and slanted with a tidal region extending 50 m to 150 m from the base of the tower. Two Thermo Scientific (model 49i) photometric O3 analysers with response time of 20 seconds sampled at 22 m and 10 m heights to obtain O3 concentration profiles at 1-minute intervals. Colocated with the O3 analysers, a Gill R3 3-dimensional ultrasonic anemometer sampling at 10 Hz was used to quantify micrometeorological parameters. The Monin-Obukhov length () was used in conjunction with measurement height () to calculate a dimensionless scaling parameter, , which defined the thermal stability state of the boundary layer. This parameter was used to ensure only that data pertaining to neutral or slightly unstable conditions were considered. Simultaneous filtering for data corresponding to the Mace Head clean marine sector (190° < wind direction < 350°) ensured the considered data corresponded to an oceanic footprint. The peak of the footprint is located approximately 1 km offshore, as in the study of Geever et al. [14]. The relative contributions to the flux footprints at midday on each measurement day are depicted in Figure 1.

In order to ensure data quality, intercomparisons between the O3 analysers were performed at a single height on a regular basis. The mean difference from these intercomparison periods (0.491 ppb) was then subtracted from the difference between them to give a final gradient. Only final gradients that were larger than the mean added to standard deviation of the difference between both instruments during an intercomparison were deemed acceptable for further analysis. Such filtering of data ensured readings were not biased by instrumental drift and eliminated noise, ensuring only strong-signal gradient data were considered.

Gradient fluxes, based on -theory [15] and previously deployed at Mace Head, described by Ceburnis et al. [16], were computed using the vertical profiles of ozone concentration between 10 m and 22 m and the turbulent-transfer coefficient derived from the micrometeorological measurements. This calculation assumes that the eddy diffusivity coefficients for momentum and O3 are similar, which is not always the case. Analogous to Fick’s law of molecular diffusion, the O3 flux () can be calculated by multiplying the vertical concentration gradient with the eddy diffusivity coefficient, , as follows: From the calculated fluxes, ozone deposition velocity () to ocean surface was calculated (). Ideally, measured fluxes should be validated by assuring energy budget closure, but unfortunately the necessary data were not available for assessment of energy closure.

3. Model Development

The Fairall parameterisation, including chemical reactivity was further developed based on the laboratory studies [12, 17] in which the chemical enhancement of the ozone deposition velocity was quantified for iodide and DOC reactants. The reaction rate for ozone and DOC, , was calculated by comparing the results from ozone-iodide reactions with the Fairall dry deposition theory and known ozone-iodide reaction kinetics [18] and determining systematic constants for the experimental setup. Results from the ozone-iodide experiments could then be extrapolated back to the ozone-DOC experiments, assuming consistent experimental setup. Fairall calculated the O3 dry deposition in terms of the sum of the surface resistance () and aerodynamic resistance (). Jones [17] performed experiments measuring loss rates of ozone deposition to saturated iodide solutions in a similar laboratory setup. In the case of ozone deposition to a saturated iodide solution, can be assumed to be zero, and so resistance to dry deposition is due to only. Using data of Jones, was estimated as the inverse to deposition velocity and found to be 3130 s m−1. This was used as typical for laboratory conditions such as those in which these experiments were executed. Combining and the rate constant used in Magi et al. [18] for sea-surface O3-I reactions, we can define the experimental setup in terms of Fairall’s dry deposition velocity equation. A surface turbulence term (τ) is required to be combined with to replicate the Martino et al. [12] experiments. Consequently, the dry deposition velocity to seawater can be written as follows: where is the solubility of ozone, is the 2nd order rate constant of ozone and reactant is the molar constant of reactant in the water surface, and is the diffusivity of ozone in seawater. First, by fitting the right hand side version of the equation to the Martino data (Figure 2), we calculate , then we calculate the 2nd order rate constant for ozone and DOC reactions (). Oceanic DOC concentrations were estimated using the relationship between chlorophyll and DOC as described by Kepkay and Wells [19] based on their study in the North Atlantic: . Typical iodide concentrations for summer-period coastal waters off the west coast of Ireland were taken as 50 nM I [20] and a typical chlorophyll-a concentration of 1 mg/m3 derived from MODIS (http://oceancolor.gsfc.nasa.gov). Other variables in the model, including ozone solubility and diffusivity were set up as outlined in Coleman et al. [9].

4. Results and Discussion

The measured and modelled fluxes and corresponding deposition velocities are illustrated in Figure 3 for the duration of the short campaign study (days 222–225 or August 10–13, 2010). The measured fluxes were averaged over 30 minute intervals. Uncertainty in measurements is represented by error bars in Figure 3. The primary source of uncertainty is in the O3 analysers which are precise only to the nearest ppb. Periods of low tide influence the measurements, as was similarly observed by Whitehead et al. [21] and are depicted in the shaded areas. U22 values during this period were mostly below 10 m s−1, with predominant NW wind directions. Gradient fluxes were mostly negative, indicating a surface sink. It can be noted here that low tide occurred around midday on all four days, coinciding with the daily solar radiation maximum, as shown in Figure 3.

Both the modelled and measured fluxes are in accordance with those observed by Gallagher et al. [7] using the eddy covariance technique. Although the modelled fluxes are the same order of magnitude as the measured fluxes, variations in measured fluxes are not reproduced in initial box model results. The surface transfer represents the main bottleneck in ozone deposition to water surfaces [10]. Reactivity (), a first order loss rate representing chemical depletion of ozone in the sea surface layer, represents the biggest uncertainty in surface transfer, and comparison with the observations suggests that some important reactivity processes are either missing or unknown.

Referring to Figure 3, it is apparent that, excluding fluxes associated with low-tide conditions, episodes of large negative fluxes coincide with high solar irradiation. This suggests that there may be photoenhancement of ozone reactions. Gas-phase reactions above the water surface but below the lowest measurement level may partially contribute to the discrepancy between the model results and measurements; however, given that halogen concentrations are more than an order of magnitude lower than O3 [22], the contribution of such reactions to the gradient is thought to be negligible; consequently, the increased flux was attributed to water-phase ozone reactions. Although Küpper et al. [23] found reactions of ozone and iodide to be light-independent, a recent study by Reeser et al. [6] demonstrated photoenhancement of reactions between ozone and organic substances. A possible reason for the photoenhancement of ozone-organic reactions can be deduced from the study of Zepp et al. [24] which describes how UV radiation alters the availability of aquatic dissolved organic matter (DOM) to microorganisms and accelerates photodegradation of DOM which may increase ocean reactivity to ozone by exposing more reactive species. In an attempt to replicate the nontidal ozone flux trends associated with irradiance, we introduce a photochemical enhancement factor for the modelled ozone-DOC reactions, whereby ozone reactions are enhanced by a factor of 20 in the presence of maximum solar radiation, and a linear scaling was applied to the reactivity term. This scaling is a crude first order approximation, and the scaling limits were chosen to match the box-model results to the in situ data, results are depicted in Figure 3. While we are lacking quantitative information relating to the enhancing effect of light on ozone-DOC reactivity, it is clear that introducing this photochemical enhancement significantly improves the replication of the experimental data by the model. The inclusion of this new photochemical term results in an increase of the ozone deposition loss rate of up to 0.6 ppb per hour compared to 0.2 ppb per hour without the photochemistry.

Episodes of high downward ozone observed on days 223.5 and 224.5 are not reproduced by the model, but are most likely due to a combination of high wind speeds and low-tidal exposure of macroalgae. The macroalgae exposure during low tide enhances ozone deposition on both physical and chemical levels: physically, ozone deposition to the rougher macroalgae surface would be more significant than to a water surface, while chemically, exposure of macroalgae to ozone stimulates oxidative stress, inducing release of molecular iodine [25] which rapidly reacts with ambient ozone, causing further depletion and perhaps an apparent increase in the downward ozone flux. Referring to Figure 1, a significant portion of the flux footprint lies within the tidal zone, indicating an augmented tidal influence during these periods. Increased depositional fluxes observed at nocturnal high tide on days 223 and 224 are not reproduced by the model. These high fluxes occur at periods of low wind speeds, and so it is likely that the flux enhancement occurs due to variations in water phase chemical reactivity. The model parameterises oceanic chemical reactivity using typical oceanic iodide and DOC concentrations and so variations in ocean tide reactivity are not accounted for in these results.

5. Conclusion

Photoenhancement of ozone-DOC reactions must be considered in the ozone dry deposition parameterisation to yield results corresponding to in situ gradient flux measurements. Presuming turbulence is adequately accounted for in surface transfer in the Fairall scheme, the ozone reactivity must be solar dependent. Inclusion of the solar radiation enhancement of sea-surface ozone reactions in the ozone deposition parameterisation causes an increase in modelled ozone flux of 38.5%, a significant enhancement to an ozone marine loss mechanism. Laboratory studies are required to quantify the effect of solar radiation on ozone with DOC, but this work indicates that failure to consider photoenhancement of ozone reactions with DOC would result in underestimation of ozone lost to sea surface via dry deposition, a loss mechanism especially pertinent for modelling intercontinental transport of ozone. Coupled with the knowledge that global concentrations of ocean-dwelling phytoplankton are in decline [26, 27], the decrease in organic matter in the Atlantic, and associated deposition flux reduction, may be a contributing factor in the previously unexplained rise in background ozone concentration as found in Mace Head [2].

Acknowledgments

The EPA, under Contracts CCRP5.3 (EASI-AQCIS), 2006-AQ-MS-50 (Ozone Trends in Ireland), and CCRP-09-FS-4-2, and Air Quality-Climate Interactions are acknowledged for support of this work.