Articles | Volume 13, issue 12
Atmos. Meas. Tech., 13, 6915–6931, 2020
Atmos. Meas. Tech., 13, 6915–6931, 2020

Research article 18 Dec 2020

Research article | 18 Dec 2020

Ozone deposition to a coastal sea: comparison of eddy covariance observations with reactive air–sea exchange models

Ozone deposition to a coastal sea: comparison of eddy covariance observations with reactive air–sea exchange models
David C. Loades1, Mingxi Yang2, Thomas G. Bell2, Adam R. Vaughan1, Ryan J. Pound1, Stefan Metzger3,4, James D. Lee1,5, and Lucy J. Carpenter1 David C. Loades et al.
  • 1Wolfson Atmospheric Chemistry Laboratories, Department of Chemistry, University of York, University Road, York, YO10 5DD, UK
  • 2Plymouth Marine Laboratory, Prospect Place, Plymouth, PL1 3DH, UK
  • 3National Ecological Observatory Network Program, Battelle, 1685 38th Street, Boulder, CO 80301, USA
  • 4Department of Atmospheric and Oceanic Sciences, University of Wisconsin-Madison, 1225 West Dayton Street, Madison, WI 53706, USA
  • 5National Centre for Atmospheric Science, University of York, University Road, York, YO10 5DD, UK

Correspondence: David C. Loades (


A fast-response (10 Hz) chemiluminescence detector for ozone (O3) was used to determine O3 fluxes using the eddy covariance technique at the Penlee Point Atmospheric Observatory (PPAO) on the south coast of the UK during April and May 2018. The median O3 flux was 0.132 mg m−2 h−1 (0.018 ppbv m s−1), corresponding to a deposition velocity of 0.037 cm s−1 (interquartile range 0.017–0.065 cm s−1) – similar to the higher values previously reported for open-ocean flux measurements but not as high as some other coastal results. We demonstrate that a typical single flux observation was above the 2σ limit of detection but had considerable uncertainty. The median 2σ uncertainty of deposition velocity was 0.031 cm s−1 for each 20 min period, which reduces with the square root of the sample size. Eddy covariance footprint analysis of the site indicates that the flux footprint was predominantly over water (> 96 %), varying with atmospheric stability and, to a lesser extent, with the tide. At very low wind speeds when the atmosphere was typically unstable, the observed ozone deposition velocity was elevated, most likely because the footprint contracted to include a greater land contribution in these conditions. At moderate to high wind speeds when atmospheric stability was near-neutral, the ozone deposition velocity increased with wind speed and showed a linear dependence with friction velocity. This observed dependence on friction velocity (and therefore also wind speed) is consistent with the predictions from the one-layer model of Fairall et al. (2007), which parameterises the oceanic deposition of ozone from the fundamental conservation equation, accounting for both ocean turbulence and near-surface chemical destruction, while assuming that chemical O3 destruction by iodide is distributed over depth. In contrast to our observations, the deposition velocity predicted by the recently developed two-layer model of Luhar et al. (2018) (which considers iodide reactivity in both layers but with molecular diffusivity dominating over turbulent diffusivity in the first layer) shows no major dependence of deposition velocity on wind speed and underestimates the measured deposition velocities. These results call for further investigation into the mechanisms and control of oceanic O3 deposition.

1 Introduction

Tropospheric ozone is important due to its considerable effects on human health (Medina-Ramón et al., 2006), agricultural yields (Heck et al., 1982), and global warming (Stevenson et al., 2013). Dry deposition is a major sink of tropospheric ozone, comprising as much as 25 % of total loss from the troposphere (Ganzeveld et al., 2009; Lelieveld and Dentener, 2000; Pound et al., 2020). Deposition to the sea surface is the greatest source of uncertainty in global estimates of total ozone dry deposition (Hardacre et al., 2015) due to deposition occurring at a slow and highly uncertain rate but over a vast area.

Ozone deposition flux is commonly parameterised according to Eq. (1) (Pacyna, 2008):

(1) F = - v d O 3 ,

where F is the flux in mol cm−2 s−1, vd is the deposition velocity in cm s−1, and [O3] is the ozone concentration in mol cm−3. In models, the deposition velocity is commonly calculated using a series of resistance terms, each defining barriers to deposition in separate layers (Wesely and Hicks, 2000):

(2) v d = R a + R b + R c - 1 .

Ra is the aerodynamic resistance, independent of the species being considered. Rb represents the resistance through the quasi-laminar thin layer of air in contact with a surface – this varies with the species' diffusivity. Lastly Rc is the surface resistance, which is typically the largest barrier to deposition for insoluble gases – roughly 95 % of total resistance in the case of ozone (Chang et al., 2004; Lenschow et al., 1982).

There are few reported observations of ozone deposition to the sea surface. Early work to determine oceanic O3 deposition velocity was either laboratory-based (Garland et al., 1980; McKay et al., 1992) or used box enclosure loss rate experiments in the field (Aldaz, 1969; Galbally and Roy, 1980). Such experiments are valuable in determining surface resistance (describing the affinity of a surface for absorbing a given gas) for ozone deposition. However, these experiments are limited in their ability to represent real-world physical processes such as turbulence at the air–sea interface. More recent flux measurements have been made with the eddy covariance method, which is the best way of observing fluxes in the atmospheric surface layer without perturbing it. Eddy covariance measurements have been made from coastal towers (Gallagher et al., 2001; Whitehead et al., 2009; McVeigh et al., 2010), aircraft (Lenschow et al., 1982; Kawa and Pearson, 1989), and ships (Bariteau et al., 2010; Helmig et al., 2012). The deposition velocities reported in the few eddy covariance studies over saltwater vary greatly: 0.01–0.15 cm s−1, with wind speed dependencies evident in some measurements and not in others.

The reported eddy covariance measurements use two different techniques to measure ozone at high frequency, both utilising chemiluminescent reactions of ozone. In the instruments used for tower-based measurements (Gallagher et al., 2001; McVeigh et al., 2010; Whitehead et al., 2009), ozone is reacted with a coumarin-based dye on the surface of a silica gel disk. Aircraft (Kawa and Pearson, 1989; Lenschow et al., 1982) and ship-borne (Bariteau et al., 2010; Helmig et al., 2012) instruments have instead utilised the reaction between ozone and gas-phase nitric oxide.

Ozone deposition to the ocean depends both upon physical exchange, facilitated by diffusion and turbulence, and chemical reaction at the water surface (Chang et al., 2004; Fairall et al., 2007; Luhar et al., 2018). Iodide in seawater has been identified as a key reactant (Garland et al., 1980). There has been considerable recent progress in understanding the global distribution of oceanic surface iodide (Chance et al., 2014, 2019; Macdonald et al., 2014; Sherwen et al., 2019). However, there has only been one report of the dependence of the iodide–ozone rate constant with temperature (Magi et al., 1997), and this remains a considerable uncertainty in global models. Dissolved organic material (DOM) has been suggested to be of similar importance for ozone deposition as iodide (Martino et al., 2012; Shaw and Carpenter, 2013), especially given its enrichment in the sea-surface microlayer (SML) (Zhou and Mopper, 1997). The complex and variable composition of DOM makes assessing its global reactivity with ozone a challenge.

Early work by Garland et al. (1980) formulated a description of ozone loss to seawater based on surface properties:

(3) v dw = a D ,

where a is the reactivity of iodide with ozone, D is the diffusivity of ozone in water, and vdw is the waterside deposition velocity, related to surface resistance (Rc) by

(4) R c = 1 α v dw ,

where α is the dimensionless solubility (liquid/gas) of ozone in water. This interpretation incorporates the chemical properties of the reaction but neglects turbulent diffusion and underestimates the deposition velocity in cold water. Fairall et al. (2007) allowed deposition velocity to vary with oceanic turbulence by considering the O3–iodide reaction beyond the molecular sublayer, obtaining the dependence

(5) v dw = a D K 1 ξ 0 K 0 ξ 0 .

K0 and K1 are modified Bessel functions of the second kind, of order 0 and 1 respectively, and

(6) ξ 0 = 2 κ u w a D ,

where κ is the von Kármán constant ( 0.4), and u∗w is the waterside friction velocity. This is sometimes referred to as a one-layer model, due to the assumption that reactivity is uniform with depth. This one-layer approach has been reported to match observations better than a using a fixed-surface resistance term but overestimates deposition velocity by a factor of 2–3 in colder waters, where the rate of reaction between ozone and iodide is slower.

An alternative, two-layer scheme is explored by Fairall et al. (2007) and expanded upon by Luhar et al. (2017). The authors consider an enhancement in reactivity in a very thin layer (reaction-diffusion sublayer) at the surface, while the water beneath has only very minor background reactivity. In a revision of the two-layer scheme, Luhar et al. (2018) assumed turbulent transfer to be negligible compared with chemical removal of ozone within the reaction-diffusion sublayer but with both turbulence and chemistry accounted for in the layer beneath, defining the waterside deposition velocity:

(7) v dw = a D ψ K 1 ξ δ cosh ( λ ) + ψ K 0 ξ δ sinh ( λ ) ψ K 1 ξ δ sinh ( λ ) + ψ K 0 ξ δ cosh ( λ ) .

The terms ψ, ξδ, and λ in Eq. (7) all vary according to the reaction-diffusion sublayer depth, δm:


Equations (7)–(10) describe the two-layer scheme that will be discussed in this work. The method of assigning a value to δm is discussed by Luhar et al. (2018), who found that a fixed depth of 3 µm was a good fit to the data of Helmig et al. (2012). When a variable reaction-diffusion sublayer depth was considered as being proportional to the reaction-diffusion length scale (lm=D/a), Luhar et al. (2018) found it necessary to multiply lm by a factor of 0.7 to obtain a δm value that fitted reasonably with observations. Pound et al. (2020) were however able to obtain a good fit to observational data without this factor by using the oceanic iodide parameterisation of Sherwen et al. (2019) in place of that of Macdonald et al. (2014). Pound et al. (2020) define the reaction-diffusion layer depth according to Eq. (11).

(11) δ m = D a

The dependence of deposition velocity with wind speed (or friction velocity, u, which scales linearly with wind speed over the ocean) within the Fairall et al. (2007) and Luhar et al. (2018) models is markedly different, and it is not clear which is a better fit to existing observations. The deposition velocity estimated by the one-layer model of Fairall et al. (2007) increases linearly with friction velocity and compares favourably with the TexAQS06 and GOMECC07 cruises (Helmig et al., 2012). However, observations made during other cruises discussed by Helmig et al. (2012) show no dependence on friction velocity. The two-layer model of Luhar et al. (2018) predicts almost no influence of friction velocity on deposition velocity, except at very low (< 2 m s−1) wind speeds.

Better characterisation of the effects of wind speed and the chemical composition of the surface water on ozone deposition velocity to the sea surface would significantly improve our understanding of the global tropospheric O3 budget (Ganzeveld et al., 2009; Pound et al., 2020). Here we present coastal ozone flux measurements made at Penlee Point Atmospheric Observatory (PPAO;, last access: 10 December 2020) on the south-west coast of the UK using a fast-response gas-phase chemiluminescence detector (CLD). Factors affecting the variation and uncertainty in the observed deposition velocity are discussed, including the effects of changing relative contributions from sea and land within the flux footprint.

2 Materials and methods

2.1 Measurement location

The PPAO is situated on a headland just south-west of Plymouth, UK (5019.08 N, 411.35 W). The observatory is located 11 m a.m.s.l., with an extendable mast on the roof. It lies 30–60 m away from the sea, depending on the tide, with the intervening land predominantly bare rock with some grass immediately surrounding the tower. For the work presented here, the top of the tower was extended to 19 m a.m.s.l. The dominant wind directions are from the south-west, followed by the north-east (Fig. 1). The focus of this work is the south-west (180–240) wind sector, which brings in air from the Atlantic Ocean and English Channel to the site (Yang et al., 2016).

Figure 1Wind directions and speeds at the PPAO during the study period (left). Radial percentage values indicate the portion of all observed wind that fell within a given sector. Local geography of the PPAO (right). © Google Earth.

2.2 Experimental set-up

The ozone chemiluminescence detector was adapted from an Eco Physics® CLD 886 NOx detector, working on the same principle as the instrument used by Helmig et al. (2012). A supply of excess NO is introduced to the sample, which reacts with O3 to generate NO2 in an excited state. The relaxation process leads to emission of a photon that is amplified and detected using a photomultiplier tube (PMT). In order to maintain a low number of dark counts, the PMT is cooled to 5 C by a Peltier cooler. Clean dry air is continuously pumped over the PMT to avoid the build-up of water (Fig. 2).

Figure 2Schematic of the ozone chemiluminescence detector.


Sample air was drawn from the top of the tower through  10 m of 3/8′′ PFA tubing by a vacuum pump at 13.5 SLPM. This maintained a turbulent flow in the main sampling line (Reynolds number  3000). A flow of 300 mL min−1 was drawn from this sample line through 1/8′′ PFA tubing and into the analyser using an internal vacuum pump (Fig. 2), limited by a critical orifice. Before entering the analyser, the sample air was first passed through a dryer consisting of 60 cm of Nafion™ tubing coiled in a container of desiccant (indicating Drierite) to reduce humidity. A three-way solenoid valve allowed for a sample of indoor air passed through a charcoal filter to remove O3 to record an instrument zero. A 50 mL min−1 flow of 2 % NO in N2 was supplied separately to the analyser at a pressure of 4 bar through approximately 1.5 m of 1/8′′ PFA tubing. The NO and O3 were then mixed immediately before the reaction chamber (at  26 mbar pressure) and the resulting chemiluminescence detected by the PMT.

The CLD counts were logged at 10 Hz and converted into ozone mixing ratios using the signal from a co-located, recently calibrated 2B Model 205 Dual Beam Ozone Monitor. The CLD sensitivity was determined to be 240 counts s−1 ppbv−1 and showed no obvious dependence on ambient humidity (Fig. S1 in the Supplement), providing evidence for the efficacy of the dryer. Instrument dark counts were 480 ± 40 count s−1, leading to a 10 Hz signal-to-noise ratio of 33 for the average 46 ppbv O3 measured during this work.

Three-dimensional wind data were obtained from a Gill WindMaster Pro 3D sonic anemometer at 10 Hz. Humidity, air pressure, and temperature data were logged at 0.25 Hz from a Gill MetPak Pro. Vertical wind data were adjusted by +16.6 % and +28.9 % in magnitude for positive and negative values, respectively, in line with the corrections recommended for a reported firmware bug in the Gill WindMaster instruments (, last access: 10 December 2020).

2.3 Pre-flux processing

The eddy covariance method (EC) relies on the simultaneous measurement of vertical wind speed (w) and the relevant scalar (in this case, ozone dry mixing ratio). These values were determined at 10 Hz in order to resolve the full range of eddies responsible for vertical ozone transport. It is necessary to calculate eddy covariance fluxes over a suitable averaging interval to reduce random noise and capture transport from large eddies whilst avoiding too long a period such that non-turbulent transport and non-stationarity become more important. An averaging time of around 30 min is often recommended (Foken, 2008). Previous measurements of O3 flux have used averaging intervals from 10 min (Helmig et al., 2012) to 1 h (Gallagher et al., 2001), and a 20 min period was chosen for this work. Prior to the flux calculation, data were despiked using a median filter despiking method (Brock, 1986; Starkenburg et al., 2016) using an order of N= 4 (9 points in a window). This involves binning the differences from the normalised data into exponentially more bins until bins exist within the range of the histogram that have zero values. Difference values beyond these empty bins are then identified as spikes and removed. For the flux calculation, data were linearly detrended to determine deviation from the mean within the averaging interval. A double rotation was applied to the wind data in each averaging interval to align the u axis with the mean wind and remove any tilt in the wind vector, resulting in a mean vertical wind of zero. A planar fit method (Wilczak et al., 2001) was considered as an alternative to double rotation, but a single set of planar fit coordinates was found to be inappropriate for the Penlee site. Instead, an approach defining separate planar fit coefficients for each 10 sector (e.g. Mammarella et al., 2007; Yuan et al., 2011) was used, resulting in a median 7 % increase in flux compared with the double rotation method. This sector-wise approach does, however, introduce discontinuous adjustments at the boundaries of the somewhat arbitrarily chosen sectors. A possible solution is to define the tilt angle as a continuous function of the wind direction (Ross and Grant, 2015), but given the minor difference between the fluxes resulting from the sector planar fit and double rotation methods, the latter was chosen for this work.

Due to the Nafion™ dryer and the fixed temperature and pressure of the reaction chamber, density corrections known as WPL (Webb–Pearman–Leuning) corrections (Webb et al., 1980) were unnecessary for determining an accurate ozone mixing ratio. However, the presence of water vapour was taken into account for the determination of ancillary parameters such as the Obukhov length used in footprint modelling. It should be noted that in addition to its effect on mixing ratio, water vapour also quenches the chemiluminescence of the reaction of NO with O3. This can be dealt with either by determining the instrument sensitivity over a range of water vapour conditions (at the cost of some sensitivity) and applying a correction or by sufficiently drying the sample air. The latter approach was taken here. Despite a range of humidity (2.8 × 10−5–1.8 × 10−2 mol mol−1; Fig. S1) over the 42 d observation period, the two instruments compare well when using a fixed sensitivity for the CLD. The sensitivity value of 240 ppbv s−1 also compares favourably to 213 ppbv s−1, which was estimated using a supply of known ozone in the absence of water vapour (supplied from a calibrated Thermo model 49i-PS ozone primary standard) during lab tests prior to deployment. These results suggest that the dryer removed any major water vapour effect on the detection of ozone concentration and flux.

The sample air must travel to the detector through the inlet tubing, which introduces a time lag relative to the instantaneously measured wind data. The two datasets must therefore be realigned in order to calculate the covariance. A cross-correlation function (CCF) was calculated at different time lags, with a high-pass Butterworth filter applied to the input values. The presence of a negative peak in the resulting CCF spectrum indicated a strong anticorrelation between ozone concentration and vertical wind, characteristic of deposition. Individual CCF plots were noisy and gave scattered lag values, with a high density around 4 s. Daily average CCF plots indicated clear peaks in all but one case and drifted from 3.9 to 4.1 s over the course of the experiment (e.g. Fig. 3). This is likely a consequence of slight particulate build-up in the sample line filters over the course of the measurements. Individual 20 min flux interval lags were accepted if they fell between 3.5 and 4.5 s to allow for some variability in conditions (e.g. atmospheric pressure), vacuum pump strength, etc. Lags that fell outside of these boundaries were then set to a value determined by a linear fit of the accepted data (Fig. S2). Simply setting the lag to 4 s in all instances was found to decrease the flux by 5 % relative to the method used here (CCF lag determination maximises the flux magnitude). The expected lag was also estimated from the inlet set-up: a 13.5 L min−1 flow rate through 10 m of 3/8′′ tubing plus a 300 mL min−1 sample flow through 2 m of 1/8′′ tubing yields a calculated lag of 4.2 s, similar to the CCF-determined values.

Figure 3Example cross correlation function (CCF) for ozone and vertical wind on 10 April. The negative peak minimum indicates that ozone data lag 3.9 s behind the wind data. Dashed blue lines denote the 95 % significance threshold.


Following these steps, the ozone flux was calculated on a 20 min basis using eddy4R (Metzger et al., 2017), with a workflow customised for our measurements. Flux values were then used to determine the deposition velocity according to Eq. (1). The molar flux was calculated using the instantaneous vertical wind, ozone mixing ratio, and density of dry air. Similarly, the ozone concentration used in Eq. (1) was calculated for dry air using the mean ozone mixing ratio for the averaging interval to avoid introducing a dependence on water vapour to the deposition velocity.

2.4 Data selection

A series of selection criteria were applied to the calculated 20 min flux data. Firstly, periods with more than 10 % missing data were excluded. Missing data were most commonly caused by periods of maintenance or when heavy rain disrupted the sonic anemometer readings. Data were also selected by wind direction – only data between the true wind direction of 180 and 240 were accepted to avoid observing deposition on the headland to the north-west.

A selection criterion based on ozone variation, as used by Bariteau et al. (2010), was introduced to avoid periods of non-stationarity i.e. significantly different conditions within an averaging interval (such as a sudden change in the air mass passing by the sensor or a change in wind direction). Data were excluded if the ozone concentration drifted significantly (> 6 ppbv in 20 min) or if the standard deviation in ozone was above 2 ppbv. Data with a standard deviation in wind direction of > 10 were also removed to avoid non-stationarity of wind, as performed by Yang et al. (2016) for the same site.

Flux footprint analysis was used to investigate the potential for land influence within the footprint area. Land influence may increase as the footprint contracts during the unstable conditions coinciding most frequently with low wind speeds. Using the flux footprint parameterisation of Kljun et al. (2015), footprints were calculated for each averaging interval. These were defined using tide-adjusted measurement height, roughness length, friction velocity, wind speed (and direction), crosswind variability, and stability conditions, then aggregated into 1 m s−1 wind speed bins. Using these aggregated footprints, the percentage of land area contribution in the footprint area was estimated to increase from 1 %–2 % at high wind speeds, when atmospheric stability was predominantly neutral, to 15 % at winds below 2 m s−1, when the atmosphere was generally unstable (Fig. 4). It should be noted that the footprint model is designed for flat homogeneous terrain – not a heterogeneous coastal site. This will therefore introduce some additional uncertainty to footprint extent and land coverage, beyond that inherent to the parameterisation.

Figure 4Land cover percentage within the average flux footprint for 1 m s−1 wind speed bins as calculated with the Kljun et al. (2015) flux footprint parameterisation. The presence of land within the footprint area was greater during periods of low wind speed and atmospheric instability.


Roughness lengths (z0), derived from eddy covariance measurements using the logarithmic wind profile and Eqs. (12)–(15), were also elevated at low wind speeds (Fig. 5).

(12) z 0 = z / e k U u - Ψ m z L ,

where z0 is the roughness length in m, z is the measurement height in m, κ is the von Kármán constant, U is the wind speed in m s−1, u is the friction velocity in m s−1 (determined directly from the covariance of the fluctuations of horizontal and vertical wind components), and ΨmzL is the integral of the universal function (with dimensionless Obukhov stability zL calculated from observed heat flux and u), defined as (Businger et al., 1971; Högström, 1988)



(15) x = 1 - 19.3 z L 1 / 4 .

Roughness lengths at high wind speeds are scattered approximately around 0.0002 m, which is expected for an open sea fetch (World Meteorological Organisation, 2008), but a large increase can be seen at wind speeds < 3 m s−1 (Fig. 5). Roughness length can be slightly higher during very low wind speed, low u conditions (Vickers and Mahrt, 2006). However the scale of the increase at the PPAO is indicative of a surface with more roughness elements, such as the rocks and grass found on the headland. Greater inaccuracies in the double rotation method at low wind speeds can mean that the removal of horizontal wind from the rotated vertical component is incomplete, further contributing to the elevated surface roughness values. Additionally, higher deposition velocities were observed during periods of very low winds, contrasting with the trend of increasing deposition velocity with wind speed proposed by Chang et al. (2004) and observed during open-ocean cruises by Helmig et al. (2012). Yang et al. (2016, 2019a) observed a similar enhancement in CO2 transfer at low wind speeds and chose to filter out low wind speed data. The above discussion indicates the need for a filter to exclude land-influenced flux data. A wind speed filter of > 3 m s−1 was used in this work, whereby median fluxes and deposition velocities are reported for the whole dataset (or model work), though filters on the basis of z0 could also be used to a similar effect.

Figure 5Roughness length for each averaging interval (black dots) with a smoothed local regression (LOESS) line (solid red, 95 % confidence interval shaded). Points left of the 3 m s−1 filter threshold (dashed red) are not used in subsequent discussions of oceanic deposition velocity. The y axis is limited for clarity, with 17 points < 10−9 m not shown.


Previous eddy covariance work on CO2 flux over land has applied filters on the basis of friction velocity (e.g. Barr et al., 2013) to avoid underestimation of flux during periods of poorly developed turbulence, especially at night (Aubinet, 2008). However past measurements of oceanic ozone deposition velocity have not reported using such a filter (Gallagher et al., 2001; Helmig et al., 2012; McVeigh et al., 2010) because very low wind speeds and u are uncommon over the ocean. For our data, removing data with u< 0.1 cm s−1 in addition to the criteria in Table 1 made no difference to the observed median deposition velocity. Therefore, given that a wind speed filter was already applied, no additional friction velocity filter was included.

Table 1Selection criteria applied to calculated fluxes, with the number (and percent) of points remaining.

Download Print Version | Download XLSX

Longer averaging intervals than 20 min were also considered, but 60 min averaging caused a large loss of data to the selection criteria. Missing data, as well as the non-stationarity of wind and ozone, contributed to an overall 23 % reduction in total data accepted when using 60 min averaging compared with 20 min averaging. This shorter averaging time was therefore retained.

2.5 Flux uncertainty

Flux uncertainty can be estimated in a number of ways, and in this work we make use of an empirical method (Langford et al., 2015; based on Wienhold, 1995) and a theoretical method (Fairall et al., 2000). In the method of Langford et al. (2015), cross-correlation functions (discussed in Sect. 2.3) are calculated at a series of improbable lag times (150–180 s) for each averaging interval, and the root mean squared deviation of these values is taken to be representative of the random error of the flux measurement. Alternatively, the theoretical estimation of flux uncertainty of Fairall et al. (2000) can be made according to the expression

(16) Δ F χ = Δ w X σ w σ X T / τ wca ,

where ΔFX is the flux uncertainty, w' is the instantaneous vertical wind velocity fluctuation, X' is the instantaneous ozone fluctuation, σw is the standard deviation in vertical wind velocity, σX is the standard deviation in ozone concentration, T is the length of the averaging interval in seconds, and τwca is the integral timescale for the instantaneous covariance time series wX. A factor with a value of 1–2 is sometimes also included in the numerator of Eq. (16) to reflect uncertainty in this relationship (Blomquist et al., 2010). A factor of 1 is used in this work. The integral timescale τwca can either be determined from a flux co-spectrum peak frequency,

(17) τ wca = 1 2 π f max ,

or empirically according to

(18) τ wca = b z U ,

where z is the measurement height in metres, U is the mean wind speed, and b is a value that varies with atmospheric stability. The value of b has been reported variably as 0.3–3 for near-neutral conditions (Blomquist et al., 2010; Lenschow and Kristensen, 1985) and on the order of 10–12 for convective or unstable conditions (Blomquist et al., 2010; Fairall, 1984). The application of these methods to our data is discussed further in Sect. 3.5.

3 Results

3.1 Flux and deposition velocity values

From 10 April to 21 May 2018, the median O3 deposition velocity was 0.037 cm s−1 (interquartile range 0.017–0.063 cm s−1), with a median mass flux of 0.132 mg m−2 h−1 and a median ozone concentration of 48 ppbv (Figs. 6 and 7). The resulting distribution of vd values was compared to that obtained with the lag time set to 180 s and was significantly different from the results of the disjoined data (Kolmogorov–Smirnov test, p< 0.001; Fig. S3), rejecting the null hypothesis that the two sets of values could be taken by chance from the same distribution. This confirms that the experimental set-up used here has a sufficiently low limit of detection to discern the flux from noise over the whole duration of the measurements. The 2σ flux uncertainty was determined for each 20 min period (see Sect. 3.5), with a median uncertainty of 0.113 mg m−2 h−1, corresponding to a deposition velocity uncertainty of 0.031 cm s−1. A typical single flux observation is therefore above the 2σ limit of detection, albeit with considerable uncertainty, although this uncertainty reduces with the square root of the sample size where averaged results are presented.

Figure 6Ozone deposition velocity (a), mass flux (b), ozone concentration (c) and wind speed (d) histograms for all periods that passed the filtering criteria. Mean values are represented by blue lines and median values by red lines. Deposition velocity and mass flux are plotted in the range 0.25–0.50 cm s−1 and 1.0–1.0 mg m−2 h−1 respectively for clarity, with arrows indicating the number of points beyond these limits.


Figure 7Time series of ozone deposition velocity  (a), ozone mass flux (b), mean ozone concentration (c), and mean wind speed (d) from 10 April to 21 May 2018. Grey crosses represent 20 min values, with red dots for 6 h means with standard errors. All concentration and wind speed data are shown from 10 April to 21 May, with only deposition/flux values that passed filtering criteria shown in panels (a) and (b). Periods with an accepted wind direction (180–240) are shaded. Flux and deposition velocity data are thus only presented from these periods and when the wind speed was > 3 m s−1 (d). The y axis in panels (a) and (b) is limited to 0.1–0.2 cm s−1 and 0.8–0.4 mg m−2 h−1 respectively for clarity.


Previous eddy covariance ozone deposition velocity measurements have yielded values of 0.009–0.034 cm s−1 over five open-ocean cruises (Helmig et al., 2012), with higher values typically corresponding to warmer oceans. Additionally, tower-based measurements have reported deposition velocities at coastal locations to be 0.025 cm s−1 (McVeigh et al., 2010), 0.030 cm s−1 (Whitehead et al., 2009), and 0.13 cm s−1 (Gallagher et al., 2001). These measurements were carried out at Mace Head (west Ireland), Weybourne (east UK), and Roscoff (north-west France) respectively. Our median vd of 0.037 cm s−1 is towards the upper end of previous work, though much lower than Gallagher et al. (2001).

3.2 Wind speed dependence

Reports on the dependence of vd on wind speed and friction velocity (u) have varied considerably; the cruise observations discussed by Helmig et al. (2012) vary from strong to zero dependence, while both McVeigh et al. (2010) and Gallagher et al. (2001) observed tentative relationships. We examine this relationship for our data in Fig. 8. Individual values that passed the filtering criteria exhibited a large degree of scatter and are therefore presented alongside median values within wind speed bins of 1 m s−1. Note that vd values removed by the wind speed filter (Sect. 2.4) are shown in the shaded region of Fig. 8 and demonstrate the elevated vd at low wind speeds. Outside of the excluded low wind speed region, vd values are relatively constant up to 10 m s−1. Above 10 m s−1, vd begins to increase, though data are sparse above 14 m s−1.

Figure 8Deposition velocity dependence on wind speed. 20 min values are shown in grey, with bin-averaged medians (1 m s−1) and interquartile ranges shown as red dots with bars. A second-order polynomial fit is plotted as a dotted red line, with a 95 % confidence interval (red shaded area). The grey region below 3 m s−1 indicates values removed by the wind speed filter (Sect. 2.4) that are not included in the fit.


The wind speed dependency of vd has been discussed in a number of other studies. Chang et al. (2004) reported a 5-fold increase in vd (0.0158–0.0775 cm s−1) from 0 to 20 m s −1, with vd near constant below 4 m s−1 and approximately doubling from 4–10 m s−1. Tower-based eddy covariance measurements by Gallagher et al. (2001) exhibited increasing ozone deposition velocity as wind speed increases, with vd tripling over the range u=0.05–0.5 m s−1. Using the same type of instrument, McVeigh et al. (2010) reported a similar trend, fitting an exponential curve to their data. Lastly, deposition velocity during two of the five cruises reported by Helmig et al. (2012) increased with increasing wind speeds. The dependence observed in our data is discussed further in Sect. 4.2.

3.3 Land influence

Aggregate flux footprint analysis of the PPAO site (as discussed in Sect. 2.4) shown in Fig. 9 suggests that the spatial contribution of land surfaces to our observed deposition velocity is approximately 3.9 %. However, deposition velocity to land is typically greater than to the ocean, amplifying the potential influence of land deposition on our data. If our observations were adjusted for 3.9 % spatial contribution of grassland (vd≈0.25 cm s−1, median land deposition value from datasets analysed by Hardacre et al., 2015), then our calculated median coastal water vd would be 0.028 cm s−1 (23 % lower than we measured). In reality the terrain is a mixture of grassland and rocky shoreline, varying in extent with the tide, so the land vd discussed above may be an overestimate. It should also be noted that the grassland deposition velocity value used here is itself prone to considerable uncertainty due to the variability of the datasets used in the model. Although there are insufficient data over the land to the north-west to reliably determine a vd value to the land around the PPAO, an estimate can be made by obtaining a least-squares solution using the land cover determined in Fig. 4 and the observed vd values in Fig. 8. Data from wind speeds > 14 m s−1 were not used (only four data points). Using all data from 2–13 m s−1 yielded values of 0.167 ± 0.080 and 0.034 ± 0.016 cm s−1 for land and sea respectively, suggesting a lesser effect from land than using the fixed value from Hardacre et al. (2015). Given that the land contribution in Fig. 4 does not stabilise until 9 m s−1, it is possible that constant vd between 4 and 10 m s−1 wind speeds (Fig. 8) may be a consequence of land influence and wind speed enhancement counteracting one another. Estimated water-only vd values, calculated by subtracting the product of the land fraction and the land vd value from the measured vd, are shown in Fig. 10.

Figure 9Flux footprint climatology for all 20 min data that passed the selection criteria output from the Kljun et al. (2015) footprint model. Each contour represents the area contributing 10 % of the observed flux, up to 90 % for the outermost contour. A binary land/sea classification estimated a mean land contribution of 3.9 %.

Figure 10Median deposition velocities in 1 m s−1 wind speed bins for combined land and sea surfaces as measured (red) and for sea only (blue). Sea only values were calculated by subtracting the land contribution, estimated from the land cover and land deposition determined by least-squares regression. Periods with wind speeds below 3 m s−1 were not included in the final results.


It is worth reiterating that the Kljun footprint model is designed for use in homogenous environments, which is not the case for our site. Furthermore, the double rotation applied to the wind data will result in varying pitch angles relative to the water surface, introducing a dependence of the footprint extent on this pitch angle. These limitations may be important for work relying on direct interpretations of the flux footprint, such as comparisons to emissions inventories (Squires et al., 2020; Vaughan et al., 2017). In contrast to an inventory comparison, we only use the flux footprint model to develop a strategy for robust data selection and generate an aggregate footprint from several individual footprints. This approach follows the works of Amiro (1998), Göckede et al. (2006, 2008), Kirby et al. (2008), Metzger (2018), and Xu et al. (2018), who have demonstrated the utility of aggregation for deriving robust footprint-based metrics in heterogeneous environments.

3.4 Tidal influence

The PPAO site flux footprint also experiences periodic variations associated with the tide, which alters the effective measurement height and changes the land type in the footprint when the shoreline is exposed. Whitehead et al. (2009) provide an extreme example of this, reporting vd increasing from 0.030 cm s−1 at high tide to 0.21 cm s−1 at low tide during the day. This large variation in their work was a consequence of a 9 m tidal range exposing the sea floor up to 3 km from the shore. At Penlee, the tide also causes periodic movement of the river plume around the headland, altering the salinity and composition of the surface water (Yang et al., 2016). This altered composition could affect the reactivity of ozone at the sea surface. Such effects will be examined in future work. Tower height above the water was determined for all flux calculations using tidal data from the British Oceanographic Data Centre (BODC), measured approximately 6 km upstream. Periodograms were also used to look for periodic variation in deposition velocity from exposed shoreline or riverine water, but none could be identified above the variability in the data. We note that previous measurements of air–sea exchange of momentum (Yang et al., 2016), CO2 (Yang et al., 2019a), and sea spray (Yang et al., 2019b) at the PPAO were also unable to identify tidal cycles in the data. Gallagher et al. (2001) report a tentative (though statistically insignificant) diurnal cycle for coastal water during observations made at Weybourne in East Anglia, UK. However, no diurnal variability was observed in the PPAO O3 flux data (as might be expected due to deposition to land), again implying minimal land influence in our filtered observations.

3.5 Measurement uncertainty

To understand the variability in our vd observations, a flux limit of detection was obtained empirically according to the method of Langford et al. (2015) (Sect. 2.5). Limits of detection were calculated for each averaging interval due to its dependence on wind speed and atmospheric stability, giving a median 2σ flux limit of detection of 0.113 mg m−2 h−1. At the average ozone concentration of 48 ppbv, this equates to a deposition velocity of 0.031 cm s−1, with 305 of the 491 averaging intervals exceeding their individually determined 2σ limit of detection.

To determine a theoretical uncertainty using Eq. (16), the peak frequency of the co-spectrum shown in Fig. 11 (0.07 Hz) was used to determine τwca as approximately 2.2 s during near-neutral conditions and wind speeds of 12.1 m s−1. Using Eqs. (17) and (18), this corresponds to a value for b of 1.5, similar to the literature values for near-neutral conditions (Blomquist et al., 2010; Lenschow and Kristensen, 1985). Since individual 20 min co-spectra were too noisy, this b value was used with Eq. (18) to determine τwca for each 20 min period. It should be noted that the value of b is stability-dependent. However, since stability was near-neutral for most periods (z/L=-0.39 to 0.15, 20th–80th percentile), the effects of varying stability on b are expected to be small.

Figure 11Average ozone flux co-spectrum for the 17 April, normalised to area = 1, shown in blue with a smoothed local regression (LOESS, dashed line) and 95 % confidence interval (blue shading). Wind speeds were 10.3–12.3 m s−1, and dimensionless Obukhov lengths were 0.14–0.17, representing near-neutral, slightly stable conditions. Expected co-spectral shape predicted by Kaimal et al. (1972) shown in black.


Using these integral timescales, a theoretical flux uncertainty can be calculated for each averaging interval using Eq. (16). The theoretical values obtained were much higher than those found empirically – the median theoretical 2σ limit of detection was 0.241 mg m−2 h−1 compared with the empirical value of 0.113 mg m−2 h−1. We note however that this is an approximation, derived from the work of Lenschow and Kristensen (1985), who multiplied the right-hand side of Eq. (16) by 2 to derive an upper limit on flux uncertainty.

Equation (16) demonstrates how the variability of ozone and of vertical wind within averaging intervals are directly related to uncertainty in the measured flux. White noise in the wind measurement is expected to be very small, whereas random noise in the ozone instrument likely represents a significant contribution to the total variance of ozone observed at 10 Hz. Given the relatively low sensitivity of the instrument used in this work (240 counts ppbv −1 s−1 compared to 2800 counts ppbv−1 s−1 reported by Helmig et al., 2012), autocovariances were calculated for each averaging interval using the 10 Hz ozone data to examine the extent to which variance in ozone concentration is caused by instrument white noise. White noise only correlates with itself at zero lag time, so it can be estimated from the difference between the first and second points in an autocovariance plot (Blomquist et al., 2010). Instrument white noise derived using this approach was found to contribute 45 %–98 % to the total ozone variance (10th–90th percentile), with a median σnoise of 1.4 ppbv. A more sensitive ozone instrument could therefore significantly improve the flux uncertainty at a 20 min averaging interval.

Besides the random uncertainty discussed above, systematic errors are also worthy of some consideration, specifically, whether the highest and lowest frequencies of turbulence have been adequately observed. High-frequency information can be lost if measurements are made too infrequently or if the sample is attenuated significantly in the sample tubing. Measurements at 10 Hz, as performed here, are widely considered sufficient to observe this high-frequency structure. Sensor separation was minimised by locating the sample inlet directly beneath the sonic anemometer ( 20 cm below). Laminar flow was also avoided through the length of the sample line (Reynolds number = 3000). As a result, the co-spectrum in Fig. 11 shows no major loss of high-frequency information compared to theory. Since fluxes were calculated over 20 min averaging intervals using linear detrending, there is also a chance that low-frequency information may not be fully observed. Firstly, using a simple block average in place of linear detrending had little effect on the median flux observed (+1.7 %), implying that linear detrending is not causing much low-frequency information loss. Using an averaging interval of 1 h instead of 20 min gave a slightly larger magnitude flux (+4.1 %) as well. However, the longer period led to much greater data loss (22 %) to the selection criteria in Sect. 2.4; hence the 20 min average was used for this work. This suggests that any low-frequency loss is approximately 5 % of the total flux – a small amount relative to the calculated 2σ random uncertainty (85 %).

4 Discussion

4.1 Model comparison

For the average meteorological conditions observed during this work, the one-layer model of Fairall et al. (2007) predicts a deposition velocity of 0.037 cm s−1, assuming reaction of ozone with iodide only. Here, one layer refers to considering the water column to have uniform reactivity to ozone with depth. This is not the same as considering the chemical reaction only in the reaction-diffusion sublayer and both chemical reaction and turbulent transfer in the layer beneath (the two-layer model). By contrast, the revised two-layer model of Luhar et al. (2018) predicts a deposition velocity of 0.018 cm −1 for the same conditions using a reaction-diffusion sublayer (δm) of 4.2 µm, parameterised using Eq. (11). An iodide concentration of  600 nmol dm−3 would be necessary to yield the observed deposition velocity – much higher than a typical oceanic value of  80 nmol dm−3 (Chance et al., 2014). However, DOM (Shaw and Carpenter, 2013), chlorophyll (Clifford et al., 2008), and surfactants (McKay et al., 1992) have also been shown to increase ozone deposition velocity. Therefore the effective pseudo-first-order rate constant for the reaction of ozone with water, a, is likely to be higher than accounted for by iodide alone. Chang et al. (2004) defined this total reactivity as

(19) a = i k i C i ,

where a is the effective pseudo-first-order rate constant for the reaction of ozone with water, and ki and Ci are the second-order rate constant and concentration of species i, respectively. We include an estimate of the effects of DOM reactivity using a typical oceanic DOM concentration of 52 µmol dm−3 (Massicotte et al., 2017) and a rate constant of 3.7 × 10−6 dm3 mol−1 s−1 (average of the values reported by Sarwar et al., 2016, and Coleman et al., 2010). Doing so increases a from 544 to 737 s−1 and leads to average deposition velocities for our field campaign of 0.048 and 0.028 cm s−1 for the models of Fairall and Luhar, respectively.

The magnitude of the effect of DOM on O3 deposition velocity remains highly uncertain due to the uncertainties in how O3 interacts with DOM and surfactants, variability in the sea-surface microlayer (SML) composition, and the effect of temperature. The coastal waters near the PPAO experience large phytoplankton growth during the spring bloom (Cushing, 1959; Smayda, 1998), and the organic content and composition of the SML could be very different compared to the open ocean. The seasonal and spatial variations in these O3-reactive substances could in turn drive differences in ozone deposition velocity. For example, Bariteau et al. (2010) reported vd increasing from 0.034 to 0.065 cm s−1 as the waters changed from open ocean to coastal during the TexAQS-2006 cruise. It is unclear how much of the observed gradient is a result of SML composition or of terrestrial influence. Similarly, the model of Ganzeveld et al. (2009) underestimated coastal ozone deposition velocities when DOM reactivity was omitted, suggesting that this may be a particularly important factor in coastal environments. While the model of Fairall et al. (2007) appears to match our observed vd well, it is possible that this is a consequence of some missing reactivity. Inclusion of DOM causes the one-layer model to overestimate vd, as reported by Luhar et al. (2018).

4.2 Wind speed dependence

In their discussion on wind speed dependence, Helmig et al. (2012) found that their data fit reasonably well with the parameterisation of Fairall et al. (2007):

(20) v d α a D c + α 6 κ u w ,

where α is the dimensionless solubility of ozone in water, a is the effective rate constant for the reaction of ozone with molecules in the surface water in s−1, Dc is the molecular diffusion coefficient of ozone in water in m2 s−1, κ is the von Kármán constant (0.4), and u∗w is the waterside friction velocity in m s−1. The fit shown in blue in Fig. 12 was determined using parameter values relevant to the experiment at the PPAO, with u∗w derived from u, assuming atmospheric surface stress to be equal to the waterside surface stress (Luhar et al., 2017):

(21) u w = ρ air ρ water u ,

where ρair and ρwater are the densities of air and water respectively. α, a, and Dc were determined empirically according to Eq. (22) (Morris, 1988), Eq. (23) (Magi et al., 1997), and Eq. (24) (Johnson and Davis, 1996):


where Ts is the sea-surface temperature (in K), and [I] is the aqueous iodide concentration in mol dm−3. We note that Eq. (23) only accounts for the reactivity of ozone with iodide in the sea surface. Other species present in the SML have also been shown to react with ozone (Martino et al., 2009; Shaw and Carpenter, 2013), but given the uncertainty surrounding their reactivity and any temperature dependence, they have been omitted here. Fixed Ts (284  K) and [I] (85 nmol dm−3) values from April–May 2018 and representative of the footprint of PPAO (Sherwen et al., 2019) were used to determine α, a, and Dc and thus vd (cm s−1) using Eq. (20) (dashed blue line in Fig. 12). This can be simplified to


In comparison, the linear fit (dashed red line in Fig. 12) of our experimental 20 min vd values against u (with standard errors) is


Our results therefore show comparable but slightly lower dependence on friction velocity (and therefore also wind speed) than predicted by the parameterisation of Fairall et al. (2007). Comparison of our data to this parameterisation yielded a root mean square error (RMSE) of 0.0522 cm −1 and a mean bias of 0.0020 cm s−1 (a positive bias here denoting observations greater than the model). Given the assumptions of the simplified model (Eq. 20) and the uncertainties in various parameters, not least the rate constant for the reaction of O3 with I (e.g. Moreno and Baeza-Romero, 2019), this agreement is perhaps surprising. The two-layer model of Luhar et al. (2018) for the same data is shown in black in Fig. 12. Considering only iodide reactivity (i.e. omitting any enhancement in the reaction rate due to the presence of organic material in both models), this model appears to under-predict deposition velocity compared with the one-layer model of Fairall et al. (2007) and lacks any major dependence on wind speed except during very calm conditions. Comparison of our data to the two-layer model gave a higher RMSE and mean bias (0.0584 and 0.0247 cm s−1 respectively).

Figure 12Deposition velocity dependence on friction velocity. 20 min values are shown in grey. Bin-averaged median fluxes (0.05 m s−1 bins) are presented with interquartile ranges in red. Dependence of O3 deposition velocity on friction velocity is presented with a linear fit in red (95 % confidence interval shaded), with the dependence predicted by Fairall et al. (2007) in blue and that predicted by Luhar et al. (2018) in black.


The two-layer model is set up to account for ozone reactions with chemical species other than iodide. Inclusion of these additional reactions would increase the predicted deposition velocity to be more similar to our observations. However, the two-layer model also predicts that vd does not strongly depend upon variations in wind speed, which is in contrast with our observations.

5 Summary and conclusions

An ozone chemiluminescence detector adapted from an Eco Physics® CLD 886 NOx detector was used to measure the ozone deposition velocity to the sea surface at a coastal site near Plymouth, on the south-west coast of the UK. The median observed deposition velocity was 0.037 cm s−1, comparable with previous values from tower-based measurements of 0.025 cm s−1 (McVeigh et al., 2010) and 0.030 cm s−1 (Whitehead et al., 2009). Furthermore, our data are at the upper end of the values obtained by Helmig et al. (2012) during ship-based, open-ocean measurements (0.009–0.034 cm−1). Cross-covariance was used to empirically determine a 2σ limit of detection for the O3 flux for each averaging interval. This limit of detection had a median value of 0.113 mg m−2 h−1 and was exceeded in 305 out of 491 flux intervals. Autocovariance of high-frequency ozone data indicated that instrument noise was a significant component in the observed ozone variability (45 %–98 %), and lowering the noise level would reduce the flux uncertainty.

In moderate to high winds, the observed deposition velocity showed a linear dependence on friction velocity in the mean. This is comparable to that predicted by the one-layer model of Fairall et al. (2007) considering only ozone–iodide reaction. However, including estimated (but unverified) contributions from ozone–DOM reactions causes the one-layer model to overpredict the observations.

For the conditions encountered during the campaign, the two-layer model of Luhar et al. (2018) yields a vd of 0.018 cm s−1 with iodide reaction only and 0.026 cm s−1 with reactions of both iodide and contributions from DOM. While the latter value is close to our median observation, the two-layer model does not reproduce the observed wind speed dependence in vd.

Elevated deposition velocities were observed at low wind speeds, contrary to predictions (Chang et al., 2004) and to previous observations (Helmig et al., 2012). We attribute this observation to a contribution to vd from land within the footprint during periods of low wind. Periods with wind speeds > 3 m s−1 (corresponding to approximately < 10 % land cover in the footprint) were used to evaluate vd. However, the possibility of land influence could not be completely removed, with our oceanic vd estimates potentially overestimated by 8 %, even after wind speed filtering. The potential for tidal effects on vd (exposing shoreline and input of river water with different chemical composition) was also examined, though no clear periodicity could be observed, either at the tidal frequency or on a diurnal timescale.

Future work will link the properties of the sea-surface microlayer in the footprint area to observed O3 fluxes. A longer time series with more observations of microlayer chemical composition may help to elucidate the influence of biogeochemical parameters, seasonal variation, and wind speed dependence, which have not been definitively characterised to date.

Code and data availability

The public code portion of eddy4R is maintained at (Durden et al., 2019;, last access: 10 December 2020). 20 min data have been submitted to the Centre for Environmental Data Analysis (CEDA) at (Loades et al., 2020). The corresponding author can be contacted directly for the full high-frequency data.


The supplement related to this article is available online at:

Author contributions

Experimental work was carried out by DCL, TGB, and MY. DCL also conducted the formal analysis and visualisation of the results, with relevant supervision from TGB and MY. SM developed the eddy4R codebase, with ARV providing modification for its use here. RJP provided software for instrumentation and validation of model applications to the data. JDL and LJC supervised the interpretation of the results. The work was proposed by LJC, who also acquired the necessary funding. DCL prepared the paper, with all authors contributing to the editing process.

Competing interests

The authors declare that they have no conflict of interest.


Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author and do not necessarily reflect the views of the National Science Foundation.


Lucy J. Carpenter and David C. Loades acknowledge funding from the Natural Environment Research Council (NERC), UK, through the grant “Iodide in the ocean: distribution and impact on iodine flux and ozone loss” (NE/N009983/1). David C. Loades also thanks NERC for the funding of his PhD project (NERC SPHERES DTP NE/L002574/1). Lucy J. Carpenter acknowledges funding from the European Research Council (ERC) under the European Union's Horizon 2020 programme (project O3-SML; grant agreement no. 833290). The National Ecological Observatory Network is a project sponsored by the National Science Foundation and managed under co-operative agreement by Battelle. This material is based upon work supported by the National Science Foundation (grant DBI-0752017). Trinity House (, last access: 10 December 2020) owns the Penlee Point Atmospheric Observatory (PPAO) site, who allows Plymouth Marine Laboratory (PML) to use the building to house instrumentation. Access to the site is arranged thanks to Mount Edgcumbe Estate (, last access: 10 December 2020). PPAO research (including the contributions of Thomas G. Bell and Mingxi Yang to this paper) is supported by NERC via the national capability ACSIS project (grant no. NE/N018044/1). We thank Frances Hopkins (PML), Daniel Philips (University of East Anglia), and Oban Jones (PML) for assistance at the field site. This work is contribution number 8 from the PPAO.

Financial support

This research has been supported by the Natural Environment Research Council (grant no. NE/L002574/1), Natural Environment Research Council (grant “Iodide in the ocean: distribution and impact on iodine flux and ozone loss”, grant no. NE/N009983/1), Natural Environment Research Council (ACSIS project, grant no. NE/N018044/1), and the European Research Council (ERC) under the European Union's Horizon 2020 programme (project O3-SML; grant agreement no. 833290).

Review statement

This paper was edited by Christof Ammann and reviewed by Christopher Fairall, Ivan Mammarella, and one anonymous referee.


Aldaz, L.: Flux Measurements of Atmospheric Ozone Over Land Water, J. Geophys. Res., 74, 6943–6946, 1969. 

Amiro, B. D.: Footprint climatologies for evapotranspiration in a boreal catchment, Agr. Forest Meteorol., 90, 195–201,, 1998. 

Aubinet, M.: Eddy Covariance CO2 Flux Measurements in Nocturnal Conditions: an Analysis of the Problem, Ecol. Appl., 18, 1368–1378, 2008. 

Bariteau, L., Helmig, D., Fairall, C. W., Hare, J. E., Hueber, J., and Lang, E. K.: Determination of oceanic ozone deposition by ship-borne eddy covariance flux measurements, Atmos. Meas. Tech., 3, 441–455,, 2010. 

Barr, A. G., Richardson, A. D., Hollinger, D. Y., Papale, D., Arain, M. A., Black, T. A., Bohrer, G., Dragoni, D., Fischer, M. L., Gu, L., Law, B. E., Margolis, H. A., Mccaughey, J. H., Munger, J. W., Oechel, W., and Schaeffer, K.: Use of change-point detection for friction-velocity threshold evaluation in eddy-covariance studies, Agr. Forest Meteorol., 171–172, 31–45,, 2013. 

Blomquist, B. W., Huebert, B. J., Fairall, C. W., and Faloona, I. C.: Determining the sea-air flux of dimethylsulfide by eddy correlation using mass spectrometry, Atmos. Meas. Tech., 3, 1–20,, 2010. 

Brock, F. V.: A Nonlinear Filter to Remove Impulse Noise from Meteorological Data, J. Atmos. Ocean. Tech., 3, 51–58,<0051:anftri>;2, 1986. 

Businger, J. A., Wyngaard, J. C., Izumi, Y., and Bradley, E. F.: Flux- profile relationships in the atmospheric surface layer, J. Atmos. Sci., 28, 181–189, 1971. 

Chance, R., Baker, A. R., Carpenter, L., and Jickells, T. D.: The distribution of iodide at the sea surface, Environ. Sci. Process. Impacts, 16, 1841–1859,, 2014. 

Chance, R. J., Tinel, L., Sherwen, T., Baker, A. R., Bell, T., Brindle, J., Campos, M. L. A. M., Croot, P., Ducklow, H., Peng, H., Hopkins, F., Hoogakker, B., Hughes, C., Jickells, T. D., Loades, D., Macaya, D. A. R., Mahajan, A. S., Malin, G., Phillips, D., Roberts, I., Roy, R., Sarkar, A., Sinha, A. K., Song, X., Winkelbauer, H., Wuttig, K., Yang, M., Peng, Z., and Carpenter, L. J.: Global sea-surface iodide observations, 1967–2018, Sci. Data, 6, 286,, 2019. 

Chang, W., Heikes, B. G., and Lee, M.: Ozone deposition to the sea surface: Chemical enhancement and wind speed dependence, Atmos. Environ., 38, 1053–1059,, 2004. 

Clifford, D., Donaldson, D. J., Brigante, M., D'Anna, B., and George, C.: Reactive uptake of ozone by chlorophyll at aqueous surfaces, Environ. Sci. Technol., 42, 1138–1143,, 2008. 

Coleman, L., Varghese, S., Tripathi, O. P., Jennings, S. G., and O'Dowd, C. D.: Regional-Scale Ozone Deposition to North-East Atlantic Waters, Adv. Meteorol., 2010, 1–16,, 2010. 

Cushing, D. H.: The seasonal variation in oceanic production as a problem in population dynamics, ICES J. Mar. Sci., 24, 455–464,, 1959. 

Durden, D., Pingintha-Durden, N., Metzger, S., Sturtevant, C., Xu, K., Luo, H., Drysdale, W. S., Vaughan, A. R., and Roberti, J. A.: NEONScience/eddy4R: eddy4R-Docker 1.0.1 (Version 1.0.1), Zenodo,, 2019. 

Fairall, C. W.: Interpretation of eddy-correlation measurements of particulate deposition and aerosol flux, Atmos. Environ., 18, 1329–1337, 1984. 

Fairall, C. W., Hare, J. E., Edson, J. B., and McGillis, W.: Parameterization and Micrometeorological Measurement of Air–Sea Gas Transfer, Bound.-Lay. Meteorol., 96, 63–106,, 2000. 

Fairall, C. W., Helmig, D., Ganzeveld, L., and Hare, J.: Water-side turbulence enhancement of ozone deposition to the ocean, Atmos. Chem. Phys., 7, 443–451,, 2007. 

Foken, T.: Micrometeorology, 1st edn., Springer-Verlag, Berlin, Heidelberg, Germany, 2008. 

Galbally, I. E. and Roy, C. R.: Destruction of ozone at the earth's surface, Q. J. Roy. Meteor. Soc., 106, 599–620,, 1980. 

Gallagher, M. W., Beswick, K. M., and Coe, H.: Ozone deposition to coastal waters, Q. J. Roy. Meteor. Soc., 127, 539–558,, 2001. 

Ganzeveld, L., Helmig, D., Fairall, C. W., Hare, J., and Pozzer, A.: Atmosphere-ocean ozone exchange: A global modeling study of biogeochemical, atmospheric, and waterside turbulence dependencies, Global Biogeochem. Cy., 23, 1–16,, 2009. 

Garland, J. A., Elzerman, A. W., and Penkett, A. A.: The mechanism for dry deposition of ozone to seawater surface, J. Geophys. Res., 85, 7488–7492, 1980. 

Göckede, M., Markkanen, T., Hasager, C. B., and Foken, T.: Update of a footprint-based approach for the characterisation of complex measurement sites, Bound.-Lay. Meteorol., 118, 635–655,, 2006. 

Göckede, M., Foken, T., Aubinet, M., Aurela, M., Banza, J., Bernhofer, C., Bonnefond, J. M., Brunet, Y., Carrara, A., Clement, R., Dellwik, E., Elbers, J., Eugster, W., Fuhrer, J., Granier, A., Grünwald, T., Heinesch, B., Janssens, I. A., Knohl, A., Koeble, R., Laurila, T., Longdoz, B., Manca, G., Marek, M., Markkanen, T., Mateus, J., Matteucci, G., Mauder, M., Migliavacca, M., Minerbi, S., Moncrieff, J., Montagnani, L., Moors, E., Ourcival, J.-M., Papale, D., Pereira, J., Pilegaard, K., Pita, G., Rambal, S., Rebmann, C., Rodrigues, A., Rotenberg, E., Sanz, M. J., Sedlak, P., Seufert, G., Siebicke, L., Soussana, J. F., Valentini, R., Vesala, T., Verbeeck, H., and Yakir, D.: Quality control of CarboEurope flux data – Part 1: Coupling footprint analyses with flux data quality assessment to evaluate sites in forest ecosystems, Biogeosciences, 5, 433–450,, 2008. 

Hardacre, C., Wild, O., and Emberson, L.: An evaluation of ozone dry deposition in global scale chemistry climate models, Atmos. Chem. Phys., 15, 6419–6436,, 2015. 

Heck, W. W., Taylor, O. C., Adams, R., Bingham, G., Miller, J., Preston, E., and Weinstein, L.: Assessment of Crop Loss from Ozone, J. Air Pollut. Control Assoc., 32, 353–361,, 1982. 

Helmig, D., Lang, E. K., Bariteau, L., Boylan, P., Fairall, C. W., Ganzeveld, L., Hare, J. E., Hueber, J., and Pallandt, M.: Atmosphere-ocean ozone fluxes during the TexAQS 2006, STRATUS 2006, GOMECC 2007, GasEx 2008, and AMMA 2008 cruises, J. Geophys. Res.-Atmos., 117, 1–15,, 2012. 

Högström, U.: Non-dimensional wind and temperature profiles in the atmospheric surface layer: A re-evaluation, Bound.-Lay. Meteorol., 42, 55–78,, 1988. 

Johnson, P. N. and Davis, R. A.: Diffusivity of ozone in water, J. Chem. Eng. Data, 41, 1485–1487,, 1996. 

Kaimal, J. C., Wyngaard, J. C., Izumi, Y., and Coté, O. R.: Spectral characteristics of surface-layer turbulence, Q. J. Roy. Meteor. Soc., 98, 563–589,, 1972. 

Kawa, S. R. and Pearson, R.: Ozone budgets from the dynamics and chemistry of marine stratocumulus experiment, J. Geophys. Res., 94, 9809,, 1989. 

Kirby, S., Dobosy, R., Williamson, D., and Dumas, E.: An aircraft-based data analysis method for discerning individual fluxes in a heterogeneous agricultural landscape, Agr. Forest Meteorol., 148, 481–489,, 2008. 

Kljun, N., Calanca, P., Rotach, M. W., and Schmid, H. P.: A simple two-dimensional parameterisation for Flux Footprint Prediction (FFP), Geosci. Model Dev., 8, 3695–3713,, 2015. 

Langford, B., Acton, W., Ammann, C., Valach, A., and Nemitz, E.: Eddy-covariance data with low signal-to-noise ratio: time-lag determination, uncertainties and limit of detection, Atmos. Meas. Tech., 8, 4197–4213,, 2015. 

Lelieveld, J. and Dentener, F. J.: What controls tropospheric ozone?, J. Geophys. Res., 105, 3531,, 2000. 

Lenschow, D. H. and Kristensen, L.: Uncorrelated noise in turbulence measurements, J. Atmos. Ocean. Tech., 2, 68–81,<0068:UNITM>2.0.CO;2, 1985. 

Lenschow, D. H., Pearson, R., and Stankov, B. B.: Measurements of ozone vertical flux to ocean and forest, J. Geophys. Res.-Oceans, 87, 8833–8837,, 1982. 

Loades, D., Bell, T., and Yang, M.: Ozone eddy covariance measurements at Penlee Point Atmospheric Observatory, April–May 2018, Centre for Environmental Data Analysis,, 2020. 

Luhar, A. K., Galbally, I. E., Woodhouse, M. T., and Thatcher, M.: An improved parameterisation of ozone dry deposition to the ocean and its impact in a global climate–chemistry model, Atmos. Chem. Phys., 17, 3749–3767,, 2017. 

Luhar, A. K., Woodhouse, M. T., and Galbally, I. E.: A revised global ozone dry deposition estimate based on a new two-layer parameterisation for air–sea exchange and the multi-year MACC composition reanalysis, Atmos. Chem. Phys., 18, 4329–4348,, 2018. 

MacDonald, S. M., Gómez Martín, J. C., Chance, R., Warriner, S., Saiz-Lopez, A., Carpenter, L. J., and Plane, J. M. C.: A laboratory characterisation of inorganic iodine emissions from the sea surface: dependence on oceanic variables and parameterisation for global modelling, Atmos. Chem. Phys., 14, 5841–5852,, 2014. 

Magi, L., Schweitzer, F., Pallares, C., Cherif, S., Mirabel, P., and George, C.: Investigation of the Uptake Rate of Ozone and Methyl Hydroperoxide by Water Surfaces, J. Phys. Chem. A, 101, 4943–4949,, 1997. 

Mammarella, I., Kolari, P., Rinne, J., Keronen, P., Pumpanen, J., and Vesala, T.: Determining the contribution of vertical advection to the net ecosystem exchange at Hyytiäläforest, Finland, Tellus B, 59, 900–909,, 2007. 

Martino, M., Mills, G. P., Woeltjen, J., and Liss, P. S.: A new source of volatile organoiodine compounds in surface seawater, Geophys. Res. Lett., 36, 2–6,, 2009. 

Martino, M., Lézé, B., Baker, A. R., and Liss, P. S.: Chemical controls on ozone deposition to water, Geophys. Res. Lett., 39, 39–43,, 2012. 

Massicotte, P., Asmala, E., Stedmon, C., and Markager, S.: Global distribution of dissolved organic matter along the aquatic continuum: Across rivers, lakes and oceans, Sci. Total Environ., 609, 180–191,, 2017. 

McKay, W. A., Stephens, B. A., and Dollard, G. J.: Laboratory Measurements of Ozone Deposition To Sea-Water and Other Saline Solutions, Atmos. Environ. A-Gen., 26, 3105–3110, 1992. 

McVeigh, P., O'Dowd, C., and Berresheim, H.: Eddy Correlation Measurements of Ozone Fluxes over Coastal Waters West of Ireland, Adv. Meteorol., 2010, 1–7,, 2010. 

Medina-Ramón, M., Zanobetti, A., and Schwartz, J.: The effect of ozone and PM10 on hospital admissions for pneumonia and chronic obstructive pulmonary disease: A national multicity study, Am. J. Epidemiol., 163, 579–588,, 2006. 

Metzger, S.: Surface-atmosphere exchange in a box: Making the control volume a suitable representation for in-situ observations, Agr. Forest Meteorol., 255, 68–80,, 2018. 

Metzger, S., Durden, D., Sturtevant, C., Luo, H., Pingintha-Durden, N., Sachs, T., Serafimovich, A., Hartmann, J., Li, J., Xu, K., and Desai, A. R.: eddy4R 0.2.0: a DevOps model for community-extensible processing and analysis of eddy-covariance data based on R, Git, Docker, and HDF5, Geosci. Model Dev., 10, 3189–3206,, 2017. 

Moreno, C. and Baeza-Romero, M. T.: A kinetic model for ozone uptake by solutions and aqueous particles containing I- and Br-, including seawater and sea-salt aerosol, Phys. Chem. Chem. Phys., 19, 19835–19856,, 2019. 

Morris, J. C.: The aqueous solubility of ozone – a review, Ozone News, 1, 14–16, available at: (last access: 10 December 2020), 1988. 

Pacyna, J. M.: Atmospheric Deposition, in: Encyclopedia of Ecology, Five-Volume Set, Elsevier Science, Amsterdam, the Netherlands, 275–285, 2008. 

Pound, R. J., Sherwen, T., Helmig, D., Carpenter, L. J., and Evans, M. J.: Influences of oceanic ozone deposition on tropospheric photochemistry, Atmos. Chem. Phys., 20, 4227–4239,, 2020. 

Ross, A. N. and Grant, E. R.: A new continuous planar fit method for calculating fluxes in complex, forested terrain, Atmos. Sci. Lett., 16, 445–452,, 2015. 

Sarwar, G., Kang, D., Foley, K., Schwede, D., Gantt, B., and Mathur, R.: Technical note: Examining ozone deposition over seawater, Atmos. Environ., 141, 255–262,, 2016. 

Shaw, M. D. and Carpenter, L. J.: Modification of ozone deposition and I2 emissions at the air-aqueous interface by dissolved organic carbon of marine origin, Environ. Sci. Technol., 47, 10947–10954,, 2013. 

Sherwen, T., Chance, R. J., Tinel, L., Ellis, D., Evans, M. J., and Carpenter, L. J.: A machine-learning-based global sea-surface iodide distribution, Earth Syst. Sci. Data, 11, 1239–1262,, 2019. 

Smayda, T. J.: Patterns of variability characterizing marine phytoplankton, with examples from Narragansett Bay, ICES J. Mar. Sci., 55, 562–573,, 1998. 

Squires, F. A., Nemitz, E., Langford, B., Wild, O., Drysdale, W. S., Acton, W. J. F., Fu, P., Grimmond, C. S. B., Hamilton, J. F., Hewitt, C. N., Hollaway, M., Kotthaus, S., Lee, J., Metzger, S., Pingintha-Durden, N., Shaw, M., Vaughan, A. R., Wang, X., Wu, R., Zhang, Q., and Zhang, Y.: Measurements of traffic-dominated pollutant emissions in a Chinese megacity, Atmos. Chem. Phys., 20, 8737–8761,, 2020. 

Starkenburg, D., Metzger, S., Fochesatto, G. J., Alfieri, J. G., Gens, R., Prakash, A., and Cristóbal, J.: Assessment of despiking methods for turbulence data in micrometeorology, J. Atmos. Ocean. Tech., 33, 2001–2013,, 2016. 

Stevenson, D. S., Young, P. J., Naik, V., Lamarque, J.-F., Shindell, D. T., Voulgarakis, A., Skeie, R. B., Dalsoren, S. B., Myhre, G., Berntsen, T. K., Folberth, G. A., Rumbold, S. T., Collins, W. J., MacKenzie, I. A., Doherty, R. M., Zeng, G., van Noije, T. P. C., Strunk, A., Bergmann, D., Cameron-Smith, P., Plummer, D. A., Strode, S. A., Horowitz, L., Lee, Y. H., Szopa, S., Sudo, K., Nagashima, T., Josse, B., Cionni, I., Righi, M., Eyring, V., Conley, A., Bowman, K. W., Wild, O., and Archibald, A.: Tropospheric ozone changes, radiative forcing and attribution to emissions in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), Atmos. Chem. Phys., 13, 3063–3085,, 2013. 

Vaughan, A. R., Lee, J. D., Shaw, M. D., Misztal, P. K., Metzger, S., Vieno, M., Davison, B., Karl, T. G., Carpenter, L. J., Lewis, A. C., Purvis, R. M., Goldstein, A. H., and Hewitt, C. N.: VOC emission rates over London and South East England obtained by airborne eddy covariance, Faraday Discuss., 200, 599–620,, 2017. 

Vickers, D. and Mahrt, L.: Evaluation of the air-sea bulk formula and sea-surface temperature variability from observations, J. Geophys. Res.-Oceans, 111, C05002,, 2006. 

Webb, E. K., Pearman, G. I., and Leuning, R.: Correction of flux measurements for density effects due to heat and water vapour transfer, Q. J. Roy. Meteor. Soc., 106, 85–100,, 1980. 

Wesely, M. L. and Hicks, B. B.: A review of the current status of knowledge on dry deposition, Atmos. Environ., 34, 2261–2282,, 2000. 

Whitehead, J. D., Mcfiggans, G. B., Gallagher, M. W., and Flynn, M. J.: Direct linkage between tidally driven coastal ozone deposition fluxes, particle emission fluxes, and subsequent CCN formation, Geophys. Res. Lett., 36, 1–5,, 2009.  

Wienhold, F.: Micrometeorological measurement and source region analysis of nitrous oxide fluxes from an agricultural soil, Atmos. Environ., 29, 2219–2227,, 1995. 

Wilczak, J. M., Oncley, S. P., and Stage, S. A.: Sonic anemometer tilt correction algorithms, Bound.-Lay. Meteorol., 99, 127–150,, 2001. 

World Meteorological Organisation: Guide to meteorological instruments and methods of observation, 7th edn., Geneva, Switzerland, WMO-No. 8, 2008. 

Xu, K., Metzger, S., and Desai, A. R.: Surface-atmosphere exchange in a box: Space-time resolved storage and net vertical fluxes from tower-based eddy covariance, Agr. Forest Meteorol., 255, 81–91,, 2018. 

Yang, M., Bell, T. G., Hopkins, F. E., Kitidis, V., Cazenave, P. W., Nightingale, P. D., Yelland, M. J., Pascal, R. W., Prytherch, J., Brooks, I. M., and Smyth, T. J.: Air–sea fluxes of CO2 and CH4 from the Penlee Point Atmospheric Observatory on the south-west coast of the UK, Atmos. Chem. Phys., 16, 5745–5761,, 2016. 

Yang, M., Bell, T. G., Brown, I. J., Fishwick, J. R., Kitidis, V., Nightingale, P. D., Rees, A. P., and Smyth, T. J.: Insights from year-long measurements of air–water CH4 and CO2 exchange in a coastal environment, Biogeosciences, 16, 961–978,, 2019a. 

Yang, M., Norris, S. J., Bell, T. G., and Brooks, I. M.: Sea spray fluxes from the southwest coast of the United Kingdom – dependence on wind speed and wave height, Atmos. Chem. Phys., 19, 15271–15284,, 2019b. 

Yuan, R., Kang, M., Park, S. Bin, Hong, J., Lee, D., and Kim, J.: Expansion of the planar-fit method to estimate flux over complex terrain, Meteorol. Atmos. Phys., 110, 123–133,, 2011. 

Zhou, X. and Mopper, K.: Photochemical production of low-molecular-weight carbonyl compounds in seawater and surface microlayer and their air-sea exchange, Mar. Chem., 56, 201–213,, 1997. 

Short summary
The loss of ozone to the sea surface was measured from the south coast of the UK and was found to be more rapid than previous observations over the open ocean. This is likely a consequence of different chemistry and biology in coastal environments. Strong winds appeared to speed up the loss of ozone. A better understanding of what influences ozone loss over the sea will lead to better model estimates of total ozone in the troposphere.