An instrument for in situ measurement of total ozone reactivity

Abstract

We present an instrument for the measurement of total ozone reactivity – the reciprocal of the chemical lifetime of ozone (${\mathrm{O}}_{\mathrm{3}}$) – in the troposphere. The Total Ozone Reactivity System (TORS) was developed with the objective to study the role of biogenic volatile organic compounds (BVOCs) as chemical sinks of tropospheric ozone. The instrument was extensively characterized and tested in the laboratory using individual BVOCs and small plants (lemon thyme, Thymus citriodorus) in a Teflon bag and proved able to measure reactivities corresponding to $>\mathrm{4.5}\times {\mathrm{10}}^{-\mathrm{5}}$${\mathrm{s}}^{-\mathrm{1}}$ (at 5 min averaging time), with an estimated total uncertainty of $\sim$32%. Such reactivities correspond to $>$20 ppb of $\mathit{\alpha }$-pinene or $>$150 ppb of isoprene in isolation – larger than typical ambient levels but observable in environmental chamber and enclosure experiments as well as in BVOC-rich environments. The functionality of TORS was demonstrated in quasi-ambient conditions with a deployment in a horticultural glasshouse containing a range of aromatic plants. The measurements of total ozone reactivity made in the glasshouse showed a clear diurnal pattern, following the emissions of BVOCs, and are consistent with mixing ratios of tens of parts per billion of monoterpenes and several parts per billion of sesquiterpenes.

1 Introduction

Ozone (${\mathrm{O}}_{\mathrm{3}}$) is a key component of the troposphere: it is known to be damaging for human health and vegetation and to reduce crop yields, and it is an important greenhouse gas . Ozone is also a primary source of the $\mathrm{OH}$ radical, the main atmospheric oxidant, and acts as an oxidant itself . Because of its importance to tropospheric chemistry, the ozone budget has long been a subject of considerable interest. Ozone is not directly emitted but is formed in the troposphere via photolysis of nitrogen dioxide (${\mathrm{NO}}_{\mathrm{2}}$), followed by reaction of atomic oxygen with molecular oxygen (Reactions  and ). $$\begin{array}{}\text{R1}& {\displaystyle }& {\displaystyle }{\mathrm{NO}}_{\mathrm{2}}+\mathrm{h}\mathit{\nu }\to \mathrm{O}{(}^{\mathrm{3}}\mathrm{P})+\mathrm{NO}(\mathit{\lambda }<\mathrm{420}\mathrm{nm})\text{R2}& {\displaystyle }& {\displaystyle }\mathrm{O}{(}^{\mathrm{3}}\mathrm{P})+{\mathrm{O}}_{\mathrm{2}}+\mathrm{M}\to {\mathrm{O}}_{\mathrm{3}}+\mathrm{M}\end{array}$$

Ozone is destroyed in the troposphere via a series of processes, both physical (e.g. dry and wet deposition) and chemical . The latter involve photolysis (Reactions  and ) and reactions with a range of inorganic molecules and unsaturated volatile organic compounds (Reactions –). $$\begin{array}{}\text{R3}& {\displaystyle }& {\displaystyle }{\mathrm{O}}_{\mathrm{3}}+\mathrm{h}\mathit{\nu }\to \mathrm{O}{(}^{\mathrm{3}}\mathrm{P})+{\mathrm{O}}_{\mathrm{2}}(\mathit{\lambda }<\mathrm{340}\mathrm{nm})\text{R4}& {\displaystyle }& {\displaystyle }{\mathrm{O}}_{\mathrm{3}}+\mathrm{h}\mathit{\nu }\to \mathrm{O}{(}^{\mathrm{1}}\mathrm{D})+{\mathrm{O}}_{\mathrm{2}}(\mathit{\lambda }<\mathrm{340}\mathrm{nm})\text{R5}& {\displaystyle }& {\displaystyle }{\mathrm{O}}_{\mathrm{3}}+\mathrm{NO}\to {\mathrm{NO}}_{\mathrm{2}}+{\mathrm{O}}_{\mathrm{2}}\text{R6}& {\displaystyle }& {\displaystyle }{\mathrm{O}}_{\mathrm{3}}+{\mathrm{NO}}_{\mathrm{2}}\to {\mathrm{NO}}_{\mathrm{3}}+{\mathrm{O}}_{\mathrm{2}}\text{R7}& {\displaystyle }& {\displaystyle }{\mathrm{O}}_{\mathrm{3}}+\mathrm{OH}\to {\mathrm{HO}}_{\mathrm{2}}+{\mathrm{O}}_{\mathrm{2}}\text{R8}& {\displaystyle }& {\displaystyle }{\mathrm{O}}_{\mathrm{3}}+{\mathrm{HO}}_{\mathrm{2}}\to \mathrm{OH}+\mathrm{2}{\mathrm{O}}_{\mathrm{2}}\text{R9}& {\displaystyle }& {\displaystyle }{\mathrm{O}}_{\mathrm{3}}+\text{unsaturated VOCs}\to \text{products}\end{array}$$

Ozone photolysis forms atomic oxygen in ground ($\mathrm{O}{(}^{\mathrm{3}}\mathrm{P})$) and excited ($\mathrm{O}{(}^{\mathrm{1}}\mathrm{D})$) states, which primarily react with molecular oxygen in the atmosphere to reform ozone in a null cycle (Reactions , ). However, a small fraction ($\sim \mathrm{10}$ % in the lower troposphere) of $\mathrm{O}{(}^{\mathrm{1}}\mathrm{D})$ reacts with water vapour to form $\mathrm{OH}$ radicals (Reaction ), a key process for the oxidative capacity of the atmosphere. $$\begin{array}{}\text{R10}& {\displaystyle }& {\displaystyle }\mathrm{O}{(}^{\mathrm{1}}\mathrm{D})+{\mathrm{O}}_{\mathrm{2}}+\mathrm{M}\to {\mathrm{O}}_{\mathrm{3}}+\mathrm{M}\text{R11}& {\displaystyle }& {\displaystyle }\mathrm{O}{(}^{\mathrm{1}}\mathrm{D})+{\mathrm{H}}_{\mathrm{2}}\mathrm{O}\to \mathrm{2}\mathrm{OH}\end{array}$$

The chemical losses of ozone are one of the least understood parts of the tropospheric ozone budget, in particular because of the many unknowns related to the abundance of volatile organic compounds (VOCs) and their reactivity with ozone . The organic compounds that react with ozone contain one or more double carbon bonds (e.g. alkenes and dialkenes), and many of these species are emitted by plants during their metabolic processes. These biogenic VOCs (BVOCs) constitute a large fraction of the carbon loading of the atmosphere: estimates suggest that the total biogenic sources of VOCs can be 8–10 times larger than the total anthropogenic sources . Isoprene is by far the most important reactive BVOC, accounting for $\sim \mathrm{50}$ % of the total natural emissions of non-methane hydrocarbons by mass, followed by monoterpenes (15 %), methanol (9 %), $\mathrm{CO}$ (7 %), and a range of other organic compounds which include acetone (4 %) and sesquiterpenes (3 %) .

However, due to the limitations of the analytical techniques used to measure VOCs, it is very likely that both the number and the mass of BVOCs in the troposphere are severely underestimated . Estimates vary depending on the approach used, but it is thought that between 20 % and 60 % (and possibly more) of the total organic carbon pool in the troposphere is currently unaccounted for . A significant fraction of these unmeasured organic compounds is constituted by biogenic VOCs: besides isoprene, only a few monoterpenes (e.g. $\mathit{\alpha }$-pinene, $\mathit{\beta }$-pinene, limonene) and even fewer sesquiterpenes (e.g. $\mathit{\beta }$-caryophyllene) are routinely measured in the atmosphere . Sesquiterpenes are particularly challenging to measure, due to their low vapour pressure, and therefore their ambient levels are most likely underestimated .

Most BVOCs are reactive with $\mathrm{OH}$ and many of them, such as isoprene, monoterpenes, and sesquiterpenes, include double carbon bonds and therefore also react with ${\mathrm{O}}_{\mathrm{3}}$: the existence of a significant pool of unknown and unmeasured BVOCs has important consequences for the quantification of the ozone budget, which are crucial to understanding the environmental and societal impacts of ozone pollution. The oxidation of BVOCs by ozone is especially important because it forms additional pollutants, such as secondary organic aerosol, and impacts key chemical processes, such as the conversion of $\mathrm{NO}$ to ${\mathrm{NO}}_{\mathrm{2}}$ and therefore the formation of ozone via Reactions () and (), as well as the radical budget . Missing a large fraction of BVOCs means that all these processes remain potentially underestimated.

One way to address this problem is to expand the number of VOCs that can be measured. This approach has achieved some success, thanks to comprehensive two-dimensional gas chromatographic techniques . However, the chemical complexity of the composition of ambient air makes it difficult, if not impossible, to completely quantify the VOC loading of the atmosphere. An alternative approach is to measure an integrated chemical property of all VOCs, such as the chemical lifetime, which includes all the reactions that remove a given species, in this case ${\mathrm{O}}_{\mathrm{3}}$ (Reactions –). Instruments to directly measure the total ozone reactivity have been demonstrated by and and used for laboratory studies of gas-phase ozonolysis reactions .

This paper presents an instrument designed to measure total ozone reactivity under ambient, environmental chamber, and branch enclosure conditions. The Total Ozone Reactivity System (TORS) was developed from the instrument described by with several modifications, as described below. Section  describes the theoretical and operating principles of TORS, while Sects.  and  describe the design and the laboratory characterization of the TORS instrument, respectively. In Sect.  the TORS instrument is evaluated with three types of experiments, increasingly approaching ambient conditions: laboratory measurements with known concentrations of selected BVOCs, laboratory measurements with small plants, and quasi-ambient measurements in a horticultural glasshouse.

2 Ozone reactivity

2.1 Reactivity measurements

The atmospheric lifetime ($\mathit{\tau }$) of a generic species $A$ is defined as the inverse of its total chemical loss rate, i.e. of its chemical reactivity ($R$):

1$${\mathit{\tau }}_{\left(A\right)}={\displaystyle \frac{\mathrm{1}}{R_A}}={\displaystyle \frac{\mathrm{1}}{{\mathrm{\Sigma }}_i{k}_i[\mathrm{X}{]}_i}},$$ where $[\mathrm{X}{]}_i$ is the concentration of a molecule reacting with $A$ and $k_i$ is the corresponding bimolecular rate coefficient. The chemical lifetime of a species has long been considered a useful diagnostic parameter, especially to investigate the loss terms of key atmospheric oxidants, such as the $\mathrm{OH}$ radical, which are generally less well known than the respective production terms .

The comparison between the directly measured lifetime and the lifetime calculated from an independent determination of $k_i$ and $[\mathrm{X}{]}_i$ (Eq. ) allows us to understand whether all the loss terms for a given species have been accounted for. This approach has been used with success to investigate the budget of the $\mathrm{OH}$ radical , of the ${\mathrm{HO}}_{\mathrm{2}}$ radical , and of the ${\mathrm{NO}}_{\mathrm{3}}$ radical . Measurements of $\mathrm{OH}$ reactivity have also helped to discover previously unknown chemistry in terms of recycling and losses of the $\mathrm{OH}$ radical .

The same principle can be used to investigate the chemical loss of ozone. From Reactions ()–(), the rate of production and loss of tropospheric ozone can be calculated as 2$$\begin{array}{rl}{\displaystyle \frac{\mathrm{d}\left[{\mathrm{O}}_{\mathrm{3}}\right]}{\mathrm{d}t}}=& -k_{\mathrm{NO}}\left[\mathrm{NO}\right]\left[{\mathrm{O}}_{\mathrm{3}}\right]-k_{{\mathrm{NO}}_{\mathrm{2}}}\left[{\mathrm{NO}}_{\mathrm{2}}\right]\left[{\mathrm{O}}_{\mathrm{3}}\right]-k_{\mathrm{OH}}\left[\mathrm{OH}\right]\left[{\mathrm{O}}_{\mathrm{3}}\right]\\ & -k_{{\mathrm{HO}}_{\mathrm{2}}}\left[{\mathrm{HO}}_{\mathrm{2}}\right]\left[{\mathrm{O}}_{\mathrm{3}}\right]-{\mathrm{\Sigma }}_i{k}_i\left[\text{VOC}{]}_i\right[{\mathrm{O}}_{\mathrm{3}}]-\mathit{\alpha }{j}_{{\mathrm{O}}_{\mathrm{3}}}[{\mathrm{O}}_{\mathrm{3}}]\\ & +j_{{\mathrm{NO}}_{\mathrm{2}}}\left[{\mathrm{NO}}_{\mathrm{2}}\right],\end{array}$$ where $\mathit{\alpha }$ is the fraction of atomic oxygen that does not reform ${\mathrm{O}}_{\mathrm{3}}$ via reaction with ${\mathrm{O}}_{\mathrm{2}}$ (Reaction ). Halogens ($\mathrm{X}$ = $\mathrm{Cl}$, $\mathrm{Br}$, $\mathrm{I}$) can also react with ozone to form halogen oxides ($\mathrm{XO}$), although this process is likely minor in continental environments far from the main halogen sources . Under the typical unpolluted conditions that can be encountered in a forested environment ($\mathrm{NO}=\mathrm{50}$ ppt, ${\mathrm{NO}}_{\mathrm{2}}=\mathrm{500}$ ppt, $\mathrm{OH}=\mathrm{5}\times {\mathrm{10}}^{\mathrm{6}}$ $\mathrm{molecule}{\mathrm{cm}}^{-\mathrm{3}}$, ${\mathrm{HO}}_{\mathrm{2}}=\mathrm{1}\times {\mathrm{10}}^{\mathrm{8}}$ $\mathrm{molecule}{\mathrm{cm}}^{-\mathrm{3}}$; ), the loss of ozone to ${\mathrm{NO}}_{\mathrm{2}}$, $\mathrm{OH}$, and ${\mathrm{HO}}_{\mathrm{2}}$ is a factor of 50–100 times slower than the loss of ozone to $\mathrm{NO}$ (Table ). Under more polluted conditions, ${\mathrm{NO}}_{\mathrm{2}}$, $\mathrm{OH}$ and ${\mathrm{HO}}_{\mathrm{2}}$ reactions are even less important as ozone sinks compared to $\mathrm{NO}$.

Table 1

Rate coefficients of the reactions with ozone ($k_{{\mathrm{O}}_{\mathrm{3}}}$, ${\mathrm{cm}}^{\mathrm{3}}\mathrm{molec}{.}^{-\mathrm{1}}{\mathrm{s}}^{-\mathrm{1}}$), chemical lifetimes for $\left[{\mathrm{O}}_{\mathrm{3}}\right]=\mathrm{120}$ ppb ($\mathit{\tau }$, min), absolute ozone reactivities ($R_{{\mathrm{O}}_{\mathrm{3}}}$, ${\mathrm{s}}^{-\mathrm{1}}$), and relative ozone reactivities (with respect to $\mathit{\alpha }$-pinene) of selected BVOCs and inorganic compounds. The rate coefficients are from and calculated for standard conditions (298 K, 1 atm); the ozone reactivities are calculated for mixing ratios of 1 ppb of each species.

Class	Species	$k_{{\mathrm{O}}_{\mathrm{3}}}$ (${\mathrm{cm}}^{\mathrm{3}}\mathrm{molec}{.}^{-\mathrm{1}}{\mathrm{s}}^{-\mathrm{1}}$)	$\mathit{\tau }$ (min)	$R_{{\mathrm{O}}_{\mathrm{3}}}$ (${\mathrm{s}}^{-\mathrm{1}}$)	Relative $R_{{\mathrm{O}}_{\mathrm{3}}}$
Hemiterpenes	isoprene	$\mathrm{1.28}\times {\mathrm{10}}^{-\mathrm{17}}$	441	$\mathrm{3.15}\times {\mathrm{10}}^{-\mathrm{7}}$	0.13
Monoterpenes	$\mathit{\alpha }$-terpinene	$\mathrm{1.90}\times {\mathrm{10}}^{-\mathrm{14}}$	0.30	$\mathrm{4.68}\times {\mathrm{10}}^{-\mathrm{4}}$	198
	$\mathit{\beta }$-ocimene	$\mathrm{5.10}\times {\mathrm{10}}^{-\mathrm{16}}$	11.1	$\mathrm{1.25}\times {\mathrm{10}}^{-\mathrm{5}}$	5.31
	myrcene	$\mathrm{4.70}\times {\mathrm{10}}^{-\mathrm{16}}$	12.0	$\mathrm{1.16}\times {\mathrm{10}}^{-\mathrm{5}}$	4.90
	limonene	$\mathrm{2.20}\times {\mathrm{10}}^{-\mathrm{16}}$	25.7	$\mathrm{5.41}\times {\mathrm{10}}^{-\mathrm{6}}$	2.29
	$\mathit{\alpha }$-pinene	$\mathrm{9.60}\times {\mathrm{10}}^{-\mathrm{17}}$	58.8	$\mathrm{2.36}\times {\mathrm{10}}^{-\mathrm{6}}$	1.00
	sabinene	$\mathrm{8.30}\times {\mathrm{10}}^{-\mathrm{17}}$	67.9	$\mathrm{2.04}\times {\mathrm{10}}^{-\mathrm{6}}$	0.86
	3-carene	$\mathrm{4.90}\times {\mathrm{10}}^{-\mathrm{17}}$	115	$\mathrm{1.21}\times {\mathrm{10}}^{-\mathrm{6}}$	0.51
	$\mathit{\beta }$-pinene	$\mathrm{1.89}\times {\mathrm{10}}^{-\mathrm{17}}$	298	$\mathrm{4.67}\times {\mathrm{10}}^{-\mathrm{7}}$	0.20
	camphene	$\mathrm{5.02}\times {\mathrm{10}}^{-\mathrm{19}}$	11 229	$\mathrm{1.24}\times {\mathrm{10}}^{-\mathrm{8}}$	0.01
Sesquiterpenes	$\mathit{\beta }$-caryophyllene	$\mathrm{1.20}\times {\mathrm{10}}^{-\mathrm{14}}$	0.50	$\mathrm{2.96}\times {\mathrm{10}}^{-\mathrm{4}}$	125
	$\mathit{\alpha }$-farnesene	$\mathrm{5.88}\times {\mathrm{10}}^{-\mathrm{16}}$	9.60	$\mathrm{1.45}\times {\mathrm{10}}^{-\mathrm{5}}$	6.13
	$\mathit{\alpha }$-copaene	$\mathrm{1.50}\times {\mathrm{10}}^{-\mathrm{16}}$	37.6	$\mathrm{3.69}\times {\mathrm{10}}^{-\mathrm{6}}$	1.56
	longifolene	$\mathrm{5.00}\times {\mathrm{10}}^{-\mathrm{19}}$	11 280	$\mathrm{1.23}\times {\mathrm{10}}^{-\mathrm{8}}$	0.01
Inorganic	$\mathrm{NO}$	$\mathrm{1.89}\times {\mathrm{10}}^{-\mathrm{14}}$	0.30	$\mathrm{4.65}\times {\mathrm{10}}^{-\mathrm{4}}$	197
	${\mathrm{HO}}_{\mathrm{2}}$	$\mathrm{2.01}\times {\mathrm{10}}^{-\mathrm{15}}$	2.80	$\mathrm{4.96}\times {\mathrm{10}}^{-\mathrm{5}}$	21.0
	${\mathrm{NO}}_{\mathrm{2}}$	$\mathrm{3.52}\times {\mathrm{10}}^{-\mathrm{17}}$	160	$\mathrm{8.67}\times {\mathrm{10}}^{-\mathrm{7}}$	0.37

During the night – when photolysis rates are zero, and the concentrations of $\mathrm{OH}$ and ${\mathrm{HO}}_{\mathrm{2}}$ are typically 2 orders of magnitude lower than during the day – $\mathrm{NO}$ is titrated to ${\mathrm{NO}}_{\mathrm{2}}$ (Reaction ) soon after sunset, leading to the formation of ${\mathrm{NO}}_{\mathrm{3}}$ radicals (Reaction ). The rate coefficient of ${\mathrm{O}}_{\mathrm{3}}+{\mathrm{NO}}_{\mathrm{2}}$ is small ($k=\mathrm{3.52}\times {\mathrm{10}}^{-\mathrm{17}}$ ${\mathrm{cm}}^{\mathrm{3}}\mathrm{molec}{.}^{-\mathrm{1}}{\mathrm{s}}^{-\mathrm{1}}$; Table ) compared to the rate coefficients of most monoterpenes and sesquiterpenes with ${\mathrm{O}}_{\mathrm{3}}$ (since its emissions are light dependent, isoprene is normally not present at night; ): for example, 1 ppb of ${\mathrm{NO}}_{\mathrm{2}}$ has the same ${\mathrm{O}}_{\mathrm{3}}$ reactivity as 0.37 ppb of $\mathit{\alpha }$-pinene and 0.16 ppb of limonene (at 298 K). This means that ${\mathrm{NO}}_{\mathrm{2}}$ can be a significant ozone loss during the night if its concentration is high – which is not usually the case in unpolluted forested environments. Equation () can thus be simplified to 3$$\begin{array}{rl}{\displaystyle \frac{\mathrm{d}\left[{\mathrm{O}}_{\mathrm{3}}\right]}{\mathrm{d}t}}=& -k_{\mathrm{NO}}\left[\mathrm{NO}\right]\left[{\mathrm{O}}_{\mathrm{3}}\right]-k_{{\mathrm{NO}}_{\mathrm{2}}}\left[{\mathrm{NO}}_{\mathrm{2}}\right]\left[{\mathrm{O}}_{\mathrm{3}}\right]\\ & -{\mathrm{\Sigma }}_i{k}_i\left[\mathrm{VOC}{]}_i\right[{\mathrm{O}}_{\mathrm{3}}]-\mathit{\alpha }{j}_{{\mathrm{O}}_{\mathrm{3}}}[{\mathrm{O}}_{\mathrm{3}}]+j_{{\mathrm{NO}}_{\mathrm{2}}}[{\mathrm{NO}}_{\mathrm{2}}],\end{array}$$ where the $\mathrm{NO}$ and photolysis terms are significant only during the day and the ${\mathrm{NO}}_{\mathrm{2}}$ term is significant only during the night (and under relatively high-${\mathrm{NO}}_x$ conditions). The concentrations of ${\mathrm{O}}_{\mathrm{3}}$, $\mathrm{NO}$, and ${\mathrm{NO}}_{\mathrm{2}}$, their rate coefficients with ${\mathrm{O}}_{\mathrm{3}}$ and the photolysis rates of ${\mathrm{O}}_{\mathrm{3}}$ and ${\mathrm{NO}}_{\mathrm{2}}$ can be measured and/or are well known. Therefore, Eq. () can be used to evaluate the contribution of the volatile organic compounds (${\mathrm{\Sigma }}_i{k}_i[\text{VOC}{]}_i$). As mentioned in Sect. , only a handful of the VOCs that react with ozone are routinely measured, which likely leads to underestimation of the VOC contribution to ozone loss. If the loss rate of ozone can be directly measured, it is possible to determine the total VOC loading and compare it with the measured VOCs, thus allowing an estimate of missing (i.e. unmeasured) VOCs.

2.2 TORS concept

The Total Ozone Reactivity System (TORS) is based on the technique developed by . At its core, the system is a dark flow tube reactor in which a known amount of ozone is reacted with a sample (e.g. $\mathrm{NO}$, unsaturated VOCs). If the change in the concentration of the co-reactants in the sample is small, the reaction follows pseudo-first-order kinetics and the change in ozone concentration is described by

4$$[{\mathrm{O}}_{\mathrm{3}}{]}_t=[{\mathrm{O}}_{\mathrm{3}}{]}_{\mathrm{0}}\times e^{-k^{\prime }t},$$ where $[{\mathrm{O}}_{\mathrm{3}}{]}_{\mathrm{0}}$ and $[{\mathrm{O}}_{\mathrm{3}}{]}_t$ are, respectively, the initial and final ozone concentrations, $k^{\prime }$ is the pseudo-first-order rate coefficient and $t$ is the reaction time, which corresponds to the reactor residence time and can be determined experimentally (Sect. ). Equation () can be solved analytically provided the reaction time and the initial and final ${\mathrm{O}}_{\mathrm{3}}$ concentrations are known: 5$$k^{\prime }={\displaystyle \frac{-\ln \left([{\mathrm{O}}_{\mathrm{3}}{]}_t/[{\mathrm{O}}_{\mathrm{3}}{]}_{\mathrm{0}}\right)}{t}}.$$

TORS provides a direct measurement of $k^{\prime }$, which includes the chemical reactions inside the reactor (described by Eq. ) plus other ${\mathrm{O}}_{\mathrm{3}}$ removal processes, such as the loss of ozone on the reactor wall. Since the reactor is completely dark, the photolysis terms in Eq. () can be neglected and the only contributors to the ${\mathrm{O}}_{\mathrm{3}}$ chemical loss are $\mathrm{NO}$ and VOCs. The focus of this study is on the reactivity of VOCs (Sect. ), and therefore the loss of ${\mathrm{O}}_{\mathrm{3}}$ due to $\mathrm{NO}$ and to the reactor wall needs to be subtracted. For the purpose of this work, we define ozone reactivity ($R_{{\mathrm{O}}_{\mathrm{3}}}$) as 6$$\begin{array}{rl}{R}_{{\mathrm{O}}_{\mathrm{3}}}& =k^{\prime }-R_{\mathrm{NO}}-R_{\mathrm{wall}}=k^{\prime }-k_{\mathrm{NO}}\left[\mathrm{NO}\right]-R_{\mathrm{wall}}\\ & ={\mathrm{\Sigma }}_i{k}_i[\text{VOC}{]}_i,\end{array}$$ where $R_{\mathrm{NO}}$ is the removal of ${\mathrm{O}}_{\mathrm{3}}$ by reaction with $\mathrm{NO}$ and $R_{\mathrm{wall}}$ is the loss of ${\mathrm{O}}_{\mathrm{3}}$ on the reactor wall. This definition of $R_{{\mathrm{O}}_{\mathrm{3}}}$ is directly comparable to the definition by , since those experiments were conducted in zero air (i.e. at $R_{\mathrm{NO}}=\mathrm{0}$). The ozone loss on the reactor wall ($R_{\mathrm{wall}}$) is a potentially significant parameter in the TORS technique and requires accurate and precise determination (Sect. ). Another potentially important factor is the effect of secondary reactions, which can increase the loss of ozone (Reactions –) causing overestimation of $R_{{\mathrm{O}}_{\mathrm{3}}}$ or can decrease the concentration of VOCs in the sample via reactions other than ozonolysis (e.g. if there is significant formation of ${\mathrm{NO}}_{\mathrm{3}}$ in the reactor via Reaction ), causing underestimation of $R_{{\mathrm{O}}_{\mathrm{3}}}$.

Particularly important can be the formation of $\mathrm{OH}$ radicals from VOC ozonolysis : to remove this interference, an $\mathrm{OH}$ scrubber, such as cyclohexane, can be added to the system. Cyclohexane does not react with ${\mathrm{O}}_{\mathrm{3}}$, but reacts quickly with $\mathrm{OH}$ ($k=\mathrm{6.95}\times {\mathrm{10}}^{-\mathrm{12}}$ ${\mathrm{cm}}^{\mathrm{3}}\mathrm{molec}{.}^{-\mathrm{1}}{\mathrm{s}}^{-\mathrm{1}}$) and forms organic peroxy radicals which combine with ${\mathrm{HO}}_{\mathrm{2}}$ to form products that do not react with ${\mathrm{O}}_{\mathrm{3}}$ . Therefore, cyclohexane acts as an efficient $\mathrm{OH}$ scrubber when present at high mixing ratios (parts per million level). The effects of the ozone wall loss and of the secondary chemistry of ozonolysis products were investigated with a box-model simulation of the TORS chemistry (Sect. ) and with laboratory characterization experiments (Sect. ).

2.3 Simulation of TORS chemistry

A box model was used to simulate the chemical reactions occurring in the TORS reactor. The main objective of the model was to assess the effect of the $\mathrm{OH}$ radicals generated by the ozonolysis of VOCs and the role of the $\mathrm{OH}$ scrubber, as well as the impact of potential interfering chemistry such as removal of ozone by $\mathrm{OH}$, ${\mathrm{HO}}_{\mathrm{2}}$, and ${\mathrm{NO}}_{\mathrm{2}}$ (Reactions –) and removal of VOCs by $\mathrm{OH}$ and ${\mathrm{NO}}_{\mathrm{3}}$ radicals.

The model was assembled using the AtChem2 modelling toolkit with the inorganic chemistry and the oxidation mechanisms of $\mathit{\alpha }$-pinene and cyclohexane taken from the Master Chemical Mechanism (MCM v3.3.1; ). The cyclohexane mechanism was updated to include a more accurate representation of the ring-opening path of the cyclohexoxy radical, following , although the model results were not substantially different from those of a model using the MCM standard cyclohexoxy radical scheme ($\le \mathrm{0.1}$ %). It must be noted that the model results depend on the VOC used, as the $\mathrm{OH}$ yield from ozone $+$ alkene reactions can vary from 0.16 for ethene to 0.90 for 2,3-dimethyl-2-butene . The choice of $\mathit{\alpha }$-pinene for the model, as well as for the laboratory experiments (Sect. ), is due to the fact that $\mathit{\alpha }$-pinene is one of the most abundant BVOCs . It also has a high $\mathrm{OH}$ yield (0.8; ), allowing an estimate of the upper bounds of the potential interferences caused by $\mathrm{OH}$ chemistry.

The model was initialized with a range of $\mathit{\alpha }$-pinene mixing ratios (0.1–50 ppb) and, for each, with a range of cyclohexane mixing ratios (0–20 ppm). The initial mixing ratios of $\mathrm{NO}$ and ${\mathrm{NO}}_{\mathrm{2}}$ were 50 and 500 ppt, respectively, representative of an unpolluted forested environment . The initial ${\mathrm{O}}_{\mathrm{3}}$ mixing ratio was set to 120 ppb, as used during the experimental work (Sect. ). The model runtime was 300 s, covering the range of potential instrument residence times (Sect. ). A summary of the model results is shown in Fig. .

Figure 1

(a, b) Modelled removal of $\mathit{\alpha }$-pinene and ${\mathrm{O}}_{\mathrm{3}}$ (relative to their initial concentrations) as a function of cyclohexane level. (c) Average modelled concentrations of $\mathrm{OH}$ as a function of cyclohexane level. (d) Sensitivity of modelled $R_{{\mathrm{O}}_{\mathrm{3}}}$ to cyclohexane level as a function of $\mathit{\alpha }$-pinene mixing ratio.

[Figure omitted. See PDF]

The removal of $\mathit{\alpha }$-pinene during the residence time in the reactor was $\sim \mathrm{7}$ % and 4 %–5 % for initial $\mathit{\alpha }$-pinene mixing ratios of 0.1 ppb and $\ge \mathrm{0.5}$ ppb, respectively (Fig. a). In the absence of cyclohexane, the removal of $\mathit{\alpha }$-pinene was slightly higher (1–2 percentage points) due to reaction with $\mathrm{OH}$ radicals. These numbers indicate that the consumption of $\mathit{\alpha }$-pinene inside the reactor was over an order of magnitude smaller than its initial concentration, and, therefore, the chemical system approached the pseudo-first-order conditions required for the TORS method (Sect. ). The error in the determination of ozone reactivity caused by the assumption of pseudo-first-order conditions can be estimated at $<\mathrm{4}$ % for $\mathit{\alpha }$-pinene mixing ratios $>\mathrm{10}$ ppb. As expected , the reactions with oxidation products of cyclohexane were not significant sinks for ozone: apart from the wall loss, which is not included in the model, the loss of ${\mathrm{O}}_{\mathrm{3}}$ inside the reactor was less than 1.5 % for $\mathit{\alpha }$-pinene mixing ratios up to 50 ppb and was independent of the concentration of cyclohexane (Fig. b).

The average modelled concentrations of the $\mathrm{OH}$ radical at different levels of cyclohexane are shown in Fig. c. In the absence of cyclohexane, the model calculated $\mathrm{OH}$ concentrations between $\mathrm{1.3}\times {\mathrm{10}}^{\mathrm{6}}$ and $\mathrm{4}\times {\mathrm{10}}^{\mathrm{6}}$ $\mathrm{molecule}{\mathrm{cm}}^{-\mathrm{3}}$ for $\mathit{\alpha }$-pinene mixing ratios of 0.1 and 50 ppb, respectively. With addition of the $\mathrm{OH}$ scrubber, the simulated concentration of $\mathrm{OH}$ in the reactor decreased by 2 orders of magnitude at mixing ratios of cyclohexane between 1 ppm (for $\mathit{\alpha }$-pinene $<\mathrm{5}$ ppb) and 10 ppm (for $\mathit{\alpha }$-pinene $=\mathrm{50}$ ppb). Increasing the level of cyclohexane above 10 ppm did not cause further significant decrease in the calculated concentration of $\mathrm{OH}$ or a reduction in the loss of ozone and $\mathit{\alpha }$-pinene (Fig. a–c), at least within the range of $\mathit{\alpha }$-pinene concentrations explored by the model. Figure d shows that the ozone reactivities determined with 1 and 20 ppm of cyclohexane were essentially the same for $\mathit{\alpha }$-pinene initial mixing ratios up to 50 ppb. Moreover, the model results show that the differences between the ozone reactivity calculated with and without cyclohexane were between $+\mathrm{1}$ % and $-\mathrm{6}$ %, depending on the $\mathit{\alpha }$-pinene level (Fig. d). This demonstrates that $\mathrm{OH}$ chemistry has a small overall impact on the determination of total ozone reactivity, a conclusion that is supported by the laboratory experiments (Sect. ).

Model-calculated ${\mathrm{HO}}_{\mathrm{2}}$ concentrations were less than $\mathrm{1}\times {\mathrm{10}}^{\mathrm{8}}$ $\mathrm{molecule}{\mathrm{cm}}^{-\mathrm{3}}$, meaning that ozone reactivity with ${\mathrm{HO}}_{\mathrm{2}}$ was 2 orders of magnitude lower than ozone reactivity with $\mathit{\alpha }$-pinene (at 10 ppb of $\mathit{\alpha }$-pinene). Only at very low concentrations of $\mathit{\alpha }$-pinene ($<\mathrm{0.1}$ ppb) was ${\mathrm{HO}}_{\mathrm{2}}$ a significant sink for ozone. The presence of ${\mathrm{NO}}_{\mathrm{3}}$ radicals in the reactor is a potentially important interference for the TORS technique, both because its formation consumes ${\mathrm{O}}_{\mathrm{3}}$ (Reaction ) and because ${\mathrm{NO}}_{\mathrm{3}}$ reacts quickly with $\mathit{\alpha }$-pinene ($k=\mathrm{6.2}\times {\mathrm{10}}^{-\mathrm{12}}$ ${\mathrm{cm}}^{\mathrm{3}}\mathrm{molec}{.}^{-\mathrm{1}}{\mathrm{s}}^{-\mathrm{1}}$). Ambient ${\mathrm{NO}}_{\mathrm{3}}$ is likely lost in the inlet before the reactor, since the transmission of ${\mathrm{NO}}_{\mathrm{3}}$ through the inlet – a 6 mm diameter, 5 m long Teflon tube with a residence time of $\sim \mathrm{4}$ s – is poor . However, ${\mathrm{NO}}_{\mathrm{3}}$ can be formed inside the reactor via Reaction () and the model calculated ${\mathrm{NO}}_{\mathrm{3}}$ formation of the order of ${\mathrm{10}}^{\mathrm{6}}$ $\mathrm{molecule}{\mathrm{cm}}^{-\mathrm{3}}$ for $\mathit{\alpha }$-pinene mixing ratios $>\mathrm{5}$ ppb. Although the rate coefficient of $\mathit{\alpha }$-pinene $+{\mathrm{NO}}_{\mathrm{3}}$ is 5 orders of magnitude larger than the rate coefficient of $\mathit{\alpha }$-pinene $+{\mathrm{O}}_{\mathrm{3}}$, the ozone concentration in the reactor is 6–7 orders of magnitude higher than the concentration of ${\mathrm{NO}}_{\mathrm{3}}$. Therefore, the reactivity of $\mathit{\alpha }$-pinene with ${\mathrm{O}}_{\mathrm{3}}$ was 1–2 orders of magnitude larger than its reactivity with ${\mathrm{NO}}_{\mathrm{3}}$ (Table ). It must also be noted that ${\mathrm{NO}}_{\mathrm{3}}$ formation in the reactor is only an issue for ambient measurements under moderate- or high-${\mathrm{NO}}_x$ conditions, not for laboratory, enclosure, and environmental chamber experiments under low- or zero-${\mathrm{NO}}_x$ conditions.

To summarize, the model of the TORS reactor suggests that, under the typical operating conditions described in Sect. , the concentrations of ${\mathrm{HO}}_{\mathrm{2}}$ and ${\mathrm{NO}}_{\mathrm{3}}$ are too small to compete with BVOCs for reaction with ${\mathrm{O}}_{\mathrm{3}}$. Additionally, the model provides no indication that the products of the oxidation of cyclohexane, when used as $\mathrm{OH}$ scrubber, can significantly affect the determination of total ozone reactivity. While parts-per-million levels of cyclohexane effectively eliminate the $\mathrm{OH}$ radicals formed by BVOC ozonolysis reactions, the model suggests that ozone reactivities determined with and without an $\mathrm{OH}$ scrubber differ by $<\mathrm{6}$ %. The model results are in agreement with the discussion in Sect. , where it was concluded that the decay of ${\mathrm{O}}_{\mathrm{3}}$ in the TORS reactor is predominantly due to the reactions with $\mathrm{NO}$ and VOCs, alongside loss on the reactor wall (Eq. ). It is important to note, however, that the conclusions drawn from the model simulations may vary depending on the chemical conditions in the reactor, as several factors affect the chemistry inside the TORS reactors: the type and mixture of VOCs in the sample, their $\mathrm{OH}$ yields, the ambient concentrations of $\mathrm{NO}$ and ${\mathrm{NO}}_{\mathrm{2}}$, and, to a lesser extent, ambient temperature and pressure (which influence the rate coefficients of chemical reactions). To date, the effect of temperature and pressure is negligible, as the system has been operated under near-ambient conditions.

3 Instrumentation

3.1 Description of TORS

The operating principles of TORS are described in Sect.  and , and a diagram of the TORS instrument is shown in Fig. . The reactor is a 1 m long polytetrafluoroethylene (PTFE) tube with an external diameter of 90 mm and an internal diameter of 87.33 mm. Several different materials and geometries for the reactor were tested during the instrument development phase. This design was found to allow a residence time inside the reactor sufficient for the ozonolysis reactions to take place to a suitably measurable extent, while minimizing the consumption of VOCs – in order to maintain pseudo-first-order conditions – and the loss of ozone on the reactor wall (Sect. ).

Figure 2

Diagram of the TORS instrument with typical flow settings (Sect. ).

[Figure omitted. See PDF]

An ozone flow is generated by irradiating a flow of zero air with a UV mercury lamp (UVP Ltd., UK); a zero-air flow is added downstream of the mercury lamp to control the concentration of ozone (Fig. ). The ozone flow is mixed with the sample flow just before the reactor, and the initial ozone concentration ($[{\mathrm{O}}_{\mathrm{3}}{]}_{\mathrm{0}}$) is measured at this point, while the final ozone concentration ($[{\mathrm{O}}_{\mathrm{3}}{]}_t$) is measured at the exit point of the reactor (Fig. ). Depending on the instrument settings, this setup produces an ${\mathrm{O}}_{\mathrm{3}}$ mixing ratio in the reactor of 100–140 ppb. Ozone concentrations are measured using two identical UV photometric ${\mathrm{O}}_{\mathrm{3}}$ monitors (model 49i, Thermo Fisher Scientific, USA). The model 49i ${\mathrm{O}}_{\mathrm{3}}$ monitor has a stated detection limit of 1 ppb and a precision of 0.25 ppb at 1 min averaging time. The reactor can be bypassed using two three-way Teflon valves, so that the two ozone monitors can simultaneously measure the ${\mathrm{O}}_{\mathrm{3}}$ concentration before it enters the reactor, thus allowing the ozone measurements to be corrected for any difference between the two monitors (Sect. ). In addition, a temperature ($T$) and relative humidity (RH) probe (HMP110, Vaisala Oyj, Finland) is inserted in the reactor to monitor temperature and relative humidity. All the flows in the TORS instrument are controlled with mass flow controllers (Brooks Instrument LLC, USA) using a custom-built control box (IGI Systems Ltd, UK). The signals from the ozone monitors, the $T$–RH probe, and the mass flow controllers are logged on a laptop and processed with bespoke software in the R programming language.

A potentially important factor for TORS is the stability of the ozone source: highly variable levels of ${\mathrm{O}}_{\mathrm{3}}$ in the ozone flow (Fig. ) can affect the determination of the ozone reactivity and increase the signal-to-noise ratio of the instrument. The ozone mixing ratio generated by the mercury lamp was found to vary, on average, by 0.4–0.6 ppb (5 min, $\mathrm{2}\mathit{\sigma }$), i.e. less than 1 %. The ozone reactivity measurements reported in Sect.  were all averaged to 5 min.

3.2 Supporting VOC measurements

A proton-transfer-reaction quadrupole-interface time-of-flight mass spectrometer (PTR-QiTOF-MS, Ionicon Analytik GmbH, Austria) was used to measure VOC concentrations during the laboratory experiments (Sect.  and ). The instrument was operated according to the standard operating conditions recommended by the manufacturer (drift pressure $=\mathrm{3.8}$ mbar, drift tube temperature $=\mathrm{80}$ ${}^{\circ }$C and $E/N=\mathrm{129}$ Td), using ${\mathrm{H}}_{\mathrm{3}}{\mathrm{O}}^{+}$ as the reagent ion. Calibration was performed using a TO-14A aromatics standard mixture (Airgas Inc., USA). This mixture does not contain biogenic compounds, so a mass transmission curve calculated using the calibration gas was used for quantification. Recent work by showed that a PTR-MS operated under standard conditions is able to accurately measure concentration of uncalibrated compounds (to within 30 %) using a mass transmission curve, if these compounds have high proton affinity and do not undergo unknown fragmentation. This assumption is likely valid for BVOCs, such as isoprene and monoterpenes , which were the primary focus of this work. The limits of detection determined using zero air for calibrated compounds were on the order of 20–80 ppt.

Measurements of BVOCs can be problematic for monoterpenes and sesquiterpenes, as the PTR technique cannot distinguish between isomers . Typically, when measured by PTR-MS the main fragment ions for monoterpenes and sesquiterpenes are at $m/z$ 81.07 and 149.1, respectively, and this fragmentation is independent of the structure of the isomers . Therefore total monoterpene and sesquiterpene concentrations were estimated using the abundances of the protonated parent ions ($m/z$ 137 and 205, respectively) and of the main fragment ions ($m/z$ 81 and 149, respectively). Other compounds associated with biogenic emissions (e.g. substituted monoterpene alcohols, after loss of a neutral ${\mathrm{H}}_{\mathrm{2}}\mathrm{O}$) can also be detected at these fragment ions and thus the estimated concentrations shown here may be considered upper limits.

The PTR-MS instrument was not available for the experiments outside the laboratory (Sect. ). Instead, samples were taken with adsorption tubes, which were desorbed using a TD UNITY 2 thermal desorption unit (Markes Int., UK) and subsequently analysed with a gas chromatograph (GC 7890B, Agilent Technologies, USA) interfaced with a BenchTOF-Select time-of-flight mass spectrometer with tandem ionization (Markes Int., UK). Further information on the analytical technique and the GC protocol can be found in . The gas chromatography–mass spectrometry (GC–MS) analysis was only qualitative due to unavailability of appropriate calibration standards, and the data were used to identify the VOCs in the samples.

4 Characterization of TORS

4.1 Residence time

The residence time is one of the key parameters in Eq. () and therefore needs to be determined as accurately as possible. The geometry of the reactor and the total flow (sample $+$ ozone) affect the residence time. Three methods were used to determine the residence time: (1) direct measurement using a double injection of acetone, (2) indirect measurement via determination of $\mathrm{NO}$ reactivity with ${\mathrm{O}}_{\mathrm{3}}$, and (3) calculation using the flow rate and the internal volume of reactor (5990 ${\mathrm{cm}}^{\mathrm{3}}$).

Method 1 is illustrated in Fig. : the reactor was connected to the PTR-MS and aliquots of 0.05 $\mathrm{\mu }\mathrm{L}$ of acetone, in a flow of zero air, were injected simultaneously at the entrance and at the exit of the reactor. Acetone was used as a tracer for these experiments because it is easily ionized by ${\mathrm{H}}_{\mathrm{3}}{\mathrm{O}}^{+}$ and does not undergo fragmentation, and consequently it is straightforward to detect by PTR-MS at its protonated mass ($m/z$ 59). The simplest way to determine the residence time is to measure the time lag between the detection of the first and the second acetone peak signals by the PTR-MS (“maxima” calculations in Fig. ). However, the flow dynamics inside the reactor are complex and, as a consequence, there is no single value of residence time, but a distribution . Hence, the acetone signal from each injection can also be used to determine the mean residence time in the reactor (“means” calculations in Fig. ): the residence times estimated using the means calculations are about 40 % larger than those estimated using the maxima calculations.

Figure 3

Two acetone double-injection experiments used to estimate the residence time in the PTFE reactor for a flow rate of 1600 sccm. The “maxima” calculation uses the differences in time of the peak acetone signals; the “means” calculation uses the differences between the mean elapsed time of the two acetone signals.

[Figure omitted. See PDF]

Method 2 uses the TORS technique to measure the reactivity of ${\mathrm{O}}_{\mathrm{3}}$ ($\sim \mathrm{20}$ ppb) with $\mathrm{NO}$: the sample flow contained only $\sim \mathrm{100}$ ppb of $\mathrm{NO}$ (diluted from a certified gas cylinder, $\mathrm{4.90}\pm \mathrm{0.25}$ ppm in ${\mathrm{N}}_{\mathrm{2}}$, by BOC UK) and, since the rate coefficient of $\mathrm{NO}+{\mathrm{O}}_{\mathrm{3}}$ is known ($\mathrm{1.89}\times {\mathrm{10}}^{-\mathrm{14}}$ ${\mathrm{cm}}^{\mathrm{3}}\mathrm{molec}{.}^{-\mathrm{1}}{\mathrm{s}}^{-\mathrm{1}}$ at 298 K, with an uncertainty of 17 %; ), the only unknown variable in Eq. () was the reaction time $t$. Since under these conditions the chemical system deviates from pseudo-first-order conditions, the mean concentration of $\mathrm{NO}$ inside the reactor was used to analyse the experimental results. The averages of several experiments conducted at different flow rates are shown in Fig. a. Method 3 assumes perfect instant mixing and plug flow in the reactor: the calculation for a range of flow rates is shown in Fig. b.

Figure 4

(a) $\mathrm{NO}$ reactivity experiments analysed using the mean $\mathrm{NO}$ concentration inside the reactor. One experiment analysed using the initial $\mathrm{NO}$ concentration is also shown, for reference. (b) Residence times as a function of reactor flow – determined by three different methods – and polynomial fit to the three methods (black dashed line). The results of the acetone-injection method are taken from Fig. . The red star indicates the residence time used in this work (140 s for a reactor flow of 2470 sccm).

[Figure omitted. See PDF]

All three methods offer internally consistent results, within the respective uncertainties (Fig. b). Method 3 agrees well with the acetone-injection means calculation but overestimates the residence times determined with method 2 by $\sim \mathrm{18}$ %. Method 2 relies on well-known kinetic parameters and implicitly takes into account the real distribution of flow paths through the reactor; it also provides a simple test of the TORS functionality using $\mathrm{NO}$ instead of a BVOC. In the experiments described in Sect. , we used a value of 140 s for the residence time, determined by fitting a second-degree polynomial to all three methods, as shown in Fig. b for a reactor flow of 2470 sccm (standard cubic centimetres per minute).

4.2 Ozone wall loss

One of the most important parameters and key uncertainties of the TORS technique is the ozone wall loss ($R_{\mathrm{wall}}$), as discussed in Sect. . Other species can be lost on the reactor surface, but these are likely to have little or no impact on the determination of the ozone reactivity either because they do not react with ${\mathrm{O}}_{\mathrm{3}}$ (e.g. multifunctional products of VOC oxidation) or because they are present at very low concentrations.

Measured $R_{{\mathrm{O}}_{\mathrm{3}}}$ is obtained by subtracting the loss of ${\mathrm{O}}_{\mathrm{3}}$ on the reactor wall ($R_{\mathrm{wall}}$) and, if present, to $\mathrm{NO}$ ($R_{\mathrm{NO}}$) from the total loss rate of ${\mathrm{O}}_{\mathrm{3}}$ (Eq. ). Small changes in the ${\mathrm{O}}_{\mathrm{3}}$ wall loss could lead to significant variation in the final determination of ozone reactivity. It is therefore important to minimize this parameter in order to reduce the uncertainty of the measurement. In the design stage, several materials and sizes were tested under dry and humid conditions: (1) a glass cylinder ($\mathrm{5}\times \mathrm{70}$ cm), (2) a quartz cylinder ($\mathrm{9}\times \mathrm{100}$ cm), (3) a PFA coil ($\mathrm{1.9}\times \mathrm{1524}$ cm), (4) a PTFE cylinder ($\mathrm{9}\times \mathrm{100}$ cm). Based on these experiments, the large-diameter (o.d. $=\mathrm{9}$ cm, i.d. $=\mathrm{8.73}$ cm) PTFE cylinder showed the lowest ozone wall loss and minimal dependence on humidity.

The ozone wall loss was regularly determined during measurements by switching to a flow of zero air instead of the regular sample with a three-way Teflon valve (Fig. ): in this operating mode there are no reactants ($\mathrm{NO}$ or VOCs) in the reactor and the measured ozone reactivity is equal to $R_{\mathrm{wall}}$ (Eq. ). During each experiment and measurement period, multiple determinations of the ozone wall loss were made. Figure  shows the average $R_{\mathrm{wall}}$ determined during a series of laboratory experiments (Sect.  and ), as a function of measured humidity, temperature, and time. $R_{\mathrm{wall}}$ showed a weak dependence on the relative humidity inside the reactor ($R_{\mathrm{wall}}=\mathrm{9.6}\times {\mathrm{10}}^{-\mathrm{7}}\times \text{RH}+\mathrm{4.4}\times {\mathrm{10}}^{-\mathrm{5}}$, with $R^{\mathrm{2}}=\mathrm{0.198}$). Ambient humidity is usually higher than the range shown in Fig. , but because of the dilution of the sample flow by the dry ozone flow before entering the reactor (Fig. ), the relative humidity in the reactor is always less than 50 %. There was no clear dependence of $R_{\mathrm{wall}}$ on temperature, at least in the limited range experienced in the laboratory (14–24 ${}^{\circ }$C). The ozone wall loss can be expected to vary with time, as the surface of the reactor is passivated and exposed to ambient air, but there was no obvious temporal trend over the 9-month period of these experiments (Fig. ).

Figure 5

Mean and $\mathrm{2}\mathit{\sigma }$ standard deviation of the ozone wall loss ($R_{\mathrm{wall}}$) in the PTFE reactor as a function of humidity, temperature, and time. The data are taken from 15 experiments, with multiple determinations of $R_{\mathrm{wall}}$ per experiment, over a 9-month period.

[Figure omitted. See PDF]

Although there was no clear pattern with respect to the measured parameters, it is apparent that there was significant variability in the ozone wall loss and, therefore, it is necessary to measure $R_{\mathrm{wall}}$ often during an experiment and/or ambient measurements. The average standard deviation of $R_{\mathrm{wall}}$ was $\mathrm{2.4}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$ and the interquartile range was 0.5–$\mathrm{1.2}\times {\mathrm{10}}^{-\mathrm{4}}$ ${\mathrm{s}}^{-\mathrm{1}}$ (mean $=\mathrm{9.8}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$, median $=\mathrm{7.1}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$), which corresponds to the reactivity of 21–51 ppb of $\mathit{\alpha }$-pinene (mean $=\mathrm{41.5}$ ppb, median $=\mathrm{30}$ ppb). The range of the measured ozone wall losses suggests that the limit of detection of TORS is of the order of a few tens of parts per billion of $\mathit{\alpha }$-pinene or the equivalent concentration, in terms of reactivity, of other BVOCs (Table ); the limit of detection was quantified with laboratory experiments using known concentrations of $\mathit{\alpha }$-pinene, as described in Sect. .

4.3 TORS operation

The flows in the TORS instrument are constrained by several competing factors (Fig. ): first, the total flow (sample $+$ ozone) must be larger than the inlet flows of both ${\mathrm{O}}_{\mathrm{3}}$ monitors, which are fixed at $\sim \mathrm{1.4}$ slpm (standard litres per minute) each. Second, the residence time in the reactor must be long enough to allow the ozonolysis reactions to take place to a measurable extent, but short enough that the consumption of the reactants does not become significant, so that pseudo-first-order conditions are maintained (Sect. ).

During the design phase of the TORS instrument several combinations of flow settings were experimented with, eventually settling on a sample flow of $\sim \mathrm{2.3}$ slpm and an ozone flow – composed of a flow of zero air through the ozone lamp, plus a dilution flow to control the concentration of ozone, and a flow of cyclohexane as $\mathrm{OH}$ scrubber – of $\sim \mathrm{1.5}$ slpm (Fig. ). These settings result in a reactor flow rate of 2470 sccm, corresponding to a residence time of 140 s (Sect. ). Since the sample flow is mixed with the ozone flow (Fig. ), a correction factor needs to be applied to account for the dilution of the sample by the ozone flow. In the experiments and measurements described in Sect. , the mixing ratio of ozone at the entrance of the reactor was $\sim \mathrm{120}$ ppb. The instrument settings, such as the residence time and the ozone concentration, can in principle be varied to adjust the sensitivity of TORS for environments with different values of $R_{{\mathrm{O}}_{\mathrm{3}}}$.

Once the flows are set, the basic operation cycle of the TORS instrument consists of three main steps.

• 1.

  “Bypass mode” to check the agreement between the two O${}_{\mathit{3}}$ monitors. The sample flow is substituted with an equal flow of zero air and the reactor is bypassed for a period of approximately 15 min. In this mode, both monitors measure the ${\mathrm{O}}_{\mathrm{3}}$ concentration in zero air (i.e. with no reactants) before the entrance of the reactor, so that a correction factor can be derived if they differ. In this work, the difference between the ${\mathrm{O}}_{\mathrm{3}}$ monitors was checked every day at the start and at the end of each experiment/measurement period and was typically $\sim \mathrm{1}$ ppb.

• 2.

  “Wall loss mode” to determine R${}_{\mathit{wall}}$ (Eq. ). The sample flow is substituted with an equal flow of zero air inside the reactor for a period of approximately half an hour. The effect of the humidity change in the reactor within such a short period of time was found to be negligible. In this work, the wall loss was determined every 2–3 h.

• 3.

  “Sampling mode”. The main operation mode of the instrument has a sample flow containing BVOCs and/or $\mathrm{NO}$.

The rate coefficients required by Eq. () to retrieve the ozone reactivity are calculated using the actual temperature measured inside the reactor (Fig. ) and ambient pressure. This procedure was followed during all the experiments and measurements described in Sect. .

5 Evaluation of TORS

The TORS instrument was tested in a series of experiments to evaluate its functionality and potential. The experiments were designed and conducted in order of increasing complexity from individual species in laboratory conditions to the complex BVOC mixture of a horticultural glasshouse. The TORS instrument was at first tested in the laboratory using pure $\mathit{\alpha }$-pinene (Sect. ), followed by emissions from small plants (Sect. ). After these tests showed that the instrument was behaving as expected under controlled conditions, it was deployed in a glasshouse containing various aromatic plants to demonstrate that TORS can measure total ozone reactivity under quasi-ambient conditions (Sect. ).

It must be noted that the rate coefficients of ${\mathrm{O}}_{\mathrm{3}}$ with BVOCs span several orders of magnitude (Table ). Therefore, for individual species in isolation, a measured ozone reactivity of (for example) $\mathrm{2.36}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$ corresponds equally to 1.9 ppm of camphene, 10 ppb of $\mathit{\alpha }$-pinene, 75 ppb of isoprene, and 80 ppt of $\mathit{\beta }$-caryophyllene. This affects the interpretation of $R_{{\mathrm{O}}_{\mathrm{3}}}$ measurements by TORS, particularly if the composition of the sample is not known, as would be the case when taking ambient measurements.

5.1 Laboratory experiments

Several laboratory experiments were carried out using known concentrations of a selected biogenic compound. A thermostated diffusion tube and pure $\mathit{\alpha }$-pinene (98 %, from Sigma-Aldrich) in zero air were used to provide a constant source of BVOC. The diffusion rate was controlled by varying the temperature of the diffusion tube and determined by regularly weighing it with a precision balance over a period of several weeks. The concentration of $\mathit{\alpha }$-pinene was then calculated from the diffusion rate and confirmed via direct measurements by PTR-MS (Sect. ), with an agreement of $\sim \mathrm{14}$ %.

Figure 6

(a) Measured and calculated mean ozone reactivities with $\mathit{\alpha }$-pinene. The blue dashed line indicates the linear regression of the measured $R_{{\mathrm{O}}_{\mathrm{3}}}$ values for $\mathit{\alpha }$-pinene mixing ratios $>\mathrm{14}$ ppb. (b) Ozone reactivity with $\mathit{\alpha }$-pinene, measured with and without cyclohexane as $\mathrm{OH}$ scrubber.

[Figure omitted. See PDF]

The values of $R_{{\mathrm{O}}_{\mathrm{3}}}$ measured during several experiments were compared with the values calculated using the known concentrations of $\mathit{\alpha }$-pinene in the sample (Eq. ). The agreement between the calculated and measured ozone reactivity for $\mathit{\alpha }$-pinene mixing ratios larger than 40 ppb was about 25 % – within the combined uncertainties of the instrument and of the $\mathit{\alpha }$-pinene $+{\mathrm{O}}_{\mathrm{3}}$ rate coefficient (41 %; ) (Fig. a). At mixing ratios below 10 ppb of $\mathit{\alpha }$-pinene, the measured reactivities cannot be statistically distinguished from each other and from zero; in fact, the corresponding reactivity ($\mathrm{2.36}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$) is of the same magnitude as the average standard deviation of $R_{\mathrm{wall}}$ (Sect. ). Measured $R_{{\mathrm{O}}_{\mathrm{3}}}$ corresponding to concentrations of $\mathit{\alpha }$-pinene $>\mathrm{14}$ ppb are linearly correlated ($r^{\mathrm{2}}=\mathrm{0.993}$) with a slope of $\mathrm{1.2}\times {\mathrm{10}}^{-\mathrm{6}}\pm \mathrm{1.0}\times {\mathrm{10}}^{-\mathrm{7}}$ ${\mathrm{s}}^{-\mathrm{1}}{\mathrm{ppb}}^{-\mathrm{1}}$, corresponding to the sensitivity of the instrument, and an intercept of $\mathrm{6.4}\times {\mathrm{10}}^{-\mathrm{5}}\pm \mathrm{8.5}\times {\mathrm{10}}^{-\mathrm{6}}$ ${\mathrm{s}}^{-\mathrm{1}}$ (Fig. a).

Based on these experiments, the TORS detection limit, for a residence time of 140 s, can be estimated between $\mathrm{4.5}\times {\mathrm{10}}^{-\mathrm{5}}$ and $\mathrm{9.0}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$, corresponding to ozone reactivities equivalent to 20–40 ppb of $\mathit{\alpha }$-pinene (Table ). These values are consistent with the estimates based on the range of measured $R_{\mathrm{wall}}$, as discussed in Sect. ; the actual detection limit for a given set of measurements depends on the magnitude of the ozone wall loss, which can vary significantly (Fig. ). These values are also comparable to the detection limit of the instrument described by of $\mathrm{1.4}\times {\mathrm{10}}^{-\mathrm{4}}$ ${\mathrm{s}}^{-\mathrm{1}}$ for a residence time of 57 s. From Eqs. () and (), the total uncertainty of TORS can be estimated by propagating the uncertainties in the determination of the residence time (related to the uncertainty in the concentration of the $\mathrm{NO}$ cylinder and in the rate coefficient of the $\mathrm{NO}+{\mathrm{O}}_{\mathrm{3}}$ reaction, Sect. ), the uncertainty of the ozone monitors (Sect. ), and the median variability of $R_{\mathrm{wall}}$ during individual experiments (Sect. ) at $\sim \mathrm{32}$ %.

Some experiments were conducted without adding cyclohexane to the ozone flow (Fig. ) to verify the effect of the $\mathrm{OH}$ scrubber, as discussed in Sect. . Ozonolysis of BVOCs is known to generate $\mathrm{OH}$ radicals with different yields , which may lead to consumption of BVOCs by $\mathrm{OH}$ in the reactor, thus causing underestimation of $R_{{\mathrm{O}}_{\mathrm{3}}}$. However, the experiments where the $\mathrm{OH}$ scrubber was used did not show substantially different results from those where it was not used (Fig. b). Measured $R_{{\mathrm{O}}_{\mathrm{3}}}$ with $\sim \mathrm{30}$ ppb of $\mathit{\alpha }$-pinene was $\mathrm{6.2}\times {\mathrm{10}}^{-\mathrm{5}}\pm \mathrm{2.2}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$ with cyclohexane and $\mathrm{5.8}\times {\mathrm{10}}^{-\mathrm{5}}\pm \mathrm{1.9}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$ without cyclohexane: the corresponding $p$ value was 0.394, indicating that the difference between the two measurements is not statistically significant. The difference between the ozone reactivities determined with and without $\mathrm{OH}$ scrubber was 6 %–7 %, in agreement with the modelling results (Sect. ). While this is less than the precision of the ozone monitors (Sect. ), we note that, in principle, the comparison of total ozone reactivity measurements with and without an $\mathrm{OH}$ scrubber can yield additional information on the speciation of the VOC mixture in the sample.

5.2 Plant experiments

Laboratory experiments were carried out using small aromatic plants in a controlled environment to test the TORS instrument under more realistic conditions. Three plants of lemon thyme (Thymus citriodorus) were enclosed in a Teflon bag (Adtech Polymer Engineering Ltd, UK) with an approximate volume of 0.1 ${\mathrm{m}}^{\mathrm{3}}$, filled with a continuous flow of zero air. A halogen lamp (ROK 120 W) was located over the bag and a temperature–humidity probe (Vaisala HMP110) was inserted into the bag, together with a small fan (RS Components) to ensure homogeneous conditions. The natural release of water vapour via transpiration and evaporation by the plants caused the humidity in the bag to rise over the course of the experiment, but relative humidity remained below 50 % inside the TORS reactor due to dilution with the dry ozone flow (Fig. ) and therefore it did not affect the loss of ${\mathrm{O}}_{\mathrm{3}}$ on the reactor wall, as discussed in Sect. . The PTR-MS was connected to the bag to identify and quantify the BVOCs that constitute the plant emissions (Sect. ). The TORS instrument and the PTR-MS sampled continuously from the Teflon bag during the experiment, which had a duration of about 48 h.

Figure  shows the ozone reactivity measurements of one lemon thyme experiment, together with the reactivity calculated using the BVOC measurements by PTR-MS. The interquartile range of measured $R_{{\mathrm{O}}_{\mathrm{3}}}$ was 3–$\mathrm{10}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$ for the first experiment and 6–$\mathrm{9}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$ for the second experiment, with mean values of $\mathrm{6.5}\times {\mathrm{10}}^{-\mathrm{5}}$ and $\mathrm{8.1}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$, respectively. Measured ozone reactivity increased by about a factor of 2 when the lamp was turned on, due to increased emissions of all BVOCs and, in particular, of the more reactive ones (i.e. monoterpenes and sesquiterpenes, Table ). This is because when the lamp was switched on, the temperature inside the bag increased (by $\sim \mathrm{10}$ ${}^{\circ }$C, Fig. ), as well as the light. Isoprene emissions are controlled by both light and temperature, but monoterpenes and sesquiterpenes emissions are mostly controlled by temperature and have an exponential response to temperature . Therefore, the emissions of these more reactive compounds tend to increase faster than those of isoprene when temperature rises quickly.

Figure 7

Measured ozone reactivity ($R_{{\mathrm{O}}_{\mathrm{3}}}$) from three plants of lemon thyme (Thymus citriodorus) in a Teflon bag, compared to the ozone reactivity calculated from BVOC measurements by PTR-MS. The temperature measured inside the bag is also shown. All data are averaged to 5 min.

[Figure omitted. See PDF]

To calculate $R_{{\mathrm{O}}_{\mathrm{3}}}$ from the PTR-MS measurements using Eq. (), a number of assumptions have to be made. The only BVOC that the proton-transfer-reaction technique can uniquely identify is isoprene. All monoterpenes and sesquiterpenes have the same molecular weight (136.24 and 204.36 $\mathrm{g}{\mathrm{mol}}^{-\mathrm{1}}$, respectively) and therefore are very difficult to distinguish from each other using a soft ionization technique . Thus, the PTR-MS instrument effectively reports the sum of monoterpenes and the sum of sesquiterpenes. To account for this problem, estimated low and high $R_{{\mathrm{O}}_{\mathrm{3}}}$ limits were calculated. Lemon thyme is an evergreen broadleaf plant, whose main emissions (besides isoprene) are $\mathit{\alpha }$-pinene, $\mathit{\beta }$-pinene, and $\mathit{\beta }$-ocimene (monoterpenes) and $\mathit{\beta }$-caryophyllene and $\mathit{\alpha }$-farnesene (sesquiterpenes) . The lower-limit $R_{{\mathrm{O}}_{\mathrm{3}}}$ estimate was calculated assuming that the measured monoterpene signal was solely due to $\mathit{\beta }$-pinene and that the measured sesquiterpene signal was solely due to $\mathit{\alpha }$-farnesene. The higher-limit $R_{{\mathrm{O}}_{\mathrm{3}}}$ estimate was calculated assuming that the measured monoterpene signal was solely due to $\mathit{\beta }$-ocimene and that the measured sesquiterpene signal was solely due to $\mathit{\beta }$-caryophyllene. This provides a range of $R_{{\mathrm{O}}_{\mathrm{3}}}$ which likely includes that of the particular BVOC mixture emitted by the lemon thyme plants (Table ). The calculated low and high $R_{{\mathrm{O}}_{\mathrm{3}}}$ limit estimates are compared to the TORS measurements in Fig. . The TORS measured reactivities were within the range of these estimates and followed the same pattern, with higher values when the light was on and the temperature higher.

5.3 Glasshouse experiments

In order to evaluate the TORS technique under quasi-ambient conditions, the instrument was deployed in a horticultural glasshouse containing a range of aromatic plants. The glasshouse is a similar environment to ambient and was subject to a continuous inflow of ambient air, but, being a semi-enclosed system, the concentrations of BVOCs emitted from the plants are higher and the concentrations of $\mathrm{NO}$ are lower than the external environment, resulting in a stronger $R_{{\mathrm{O}}_{\mathrm{3}}}$ signal. The glasshouse is located at the Winterbourne House and Garden (https://www.winterbourne.org.uk/, last access: 1 March 2020), adjacent to the University of Birmingham campus, and has an approximate volume of 200 ${\mathrm{m}}^{\mathrm{3}}$. The following plants were inside the glasshouse during the sampling period: fringed “French” lavender (Lavandula dentata var. candicans), lemon verbena (Aloysia triphylla), scented-leaf geraniums (Pelargonium x citriodorum “Prince of Orange”, Pelargonium Radula), and several varieties of the Citrus genus (Citrus x limon, Citrus x latifolia “Tahiti”, Citrus reticulata “Clementine”).

The TORS instrument was set up in a similar way as in the plant experiments (Sect. ), with regular determinations of the ozone wall loss using a flow of zero air instead of the ambient flow. The ozone wall loss during the measurement period varied between $\mathrm{4.9}\times {\mathrm{10}}^{-\mathrm{5}}$ and $\mathrm{1.1}\times {\mathrm{10}}^{-\mathrm{4}}$ ${\mathrm{s}}^{-\mathrm{1}}$ (first and third quartiles), with mean values between $\mathrm{0.7}\times {\mathrm{10}}^{-\mathrm{4}}$ and $\mathrm{1.1}\times {\mathrm{10}}^{-\mathrm{4}}$ ${\mathrm{s}}^{-\mathrm{1}}$. The measurements were taken over a period of 2 weeks in early June 2018; during this period the weather was dry (mean RH $=\mathrm{57}$ %) with temperatures reaching a maximum of 39 ${}^{\circ }$C inside the glasshouse (mean $=\mathrm{15}$ ${}^{\circ }$C). Cyclohexane was used as $\mathrm{OH}$ scrubber only during the second week of measurements. The PTR-MS was not available at the glasshouse, but two air samples were taken on two different days using adsorption tubes and qualitatively analysed by GC–MS (Sect. ). The GC data were used to determine the most important monoterpenes and sesquiterpenes in the air inside the glasshouse, based on their relative abundance.

The total ozone reactivity measurements made in the glasshouse are shown in Fig. . For the period without cyclohexane the interquartile range was 1.9–$\mathrm{4.6}\times {\mathrm{10}}^{-\mathrm{4}}$ ${\mathrm{s}}^{-\mathrm{1}}$ with a mean value of $\mathrm{3.3}\times {\mathrm{10}}^{-\mathrm{4}}$ ${\mathrm{s}}^{-\mathrm{1}}$. For the period with cyclohexane the interquartile range was 2.0–$\mathrm{4.2}\times {\mathrm{10}}^{-\mathrm{4}}$ ${\mathrm{s}}^{-\mathrm{1}}$ with a mean value of $\mathrm{3.3}\times {\mathrm{10}}^{-\mathrm{4}}$ ${\mathrm{s}}^{-\mathrm{1}}$. Taking into account the natural variability of plant emissions, these numbers suggest that the use of an $\mathrm{OH}$ scrubber does not significantly change the TORS measurements, in keeping with the laboratory experiments (Sect. ) and the model results of the chemistry inside the reactor (Sect. ).

Figure 8

Total ozone reactivity ($R_{{\mathrm{O}}_{\mathrm{3}}}+R_{\mathrm{NO}}$) measured in the Winterbourne House and Garden glasshouse compared to the ozone reactivity calculated from estimated concentrations of BVOCs. The temperature measured inside the glasshouse is also shown. All data are averaged to 5 min.

[Figure omitted. See PDF]

In the absence of BVOC measurements, an estimate of the ozone reactivity was calculated using the qualitative information obtained from the GC–MS analysis of the adsorption tubes and from the emission factors by : broadleaf evergreen plants emit isoprene, monoterpenes, and sesquiterpenes in a proportion of approximately $\mathrm{1}:\mathrm{0.1}:\mathrm{0.02}$. However, found that Citrus plants, several types of which were present in the glasshouse, emit more monoterpenes than isoprene. Since most monoterpenes are more reactive with ozone than isoprene (Table ), the estimates of $R_{{\mathrm{O}}_{\mathrm{3}}}$ discussed below are relatively insensitive to the actual isoprene concentration. Analysis of the adsorption tubes showed that the most important monoterpenes were limonene, $\mathit{\beta }$-pinene, camphene, and myrcene, and the most important sesquiterpenes were longifolene and farnesene. Based on these results, a high $R_{{\mathrm{O}}_{\mathrm{3}}}$ value was estimated assuming 100 ppb of isoprene, 25 ppb of myrcene, and 5 ppb of $\mathit{\alpha }$-farnesene. A low $R_{{\mathrm{O}}_{\mathrm{3}}}$ was estimated assuming 50 ppb of isoprene, 5 ppb of camphene, and 1 ppb of longifolene. The estimated low and high $R_{{\mathrm{O}}_{\mathrm{3}}}$ values were of the same magnitude as the TORS measurements (Fig. ).

In contrast to the laboratory experiments (Sect.  and ), which were performed using zero air, the ozone reactivity measurements in the glasshouse were affected by ambient $\mathrm{NO}$ (Reaction ). Measurements of $\mathrm{NO}$ were not available inside the glasshouse, so it is not possible to definitively quantify the contribution of $R_{\mathrm{NO}}$ to the total ozone reactivity measurements shown in Fig. . The average total ozone reactivity showed a clear diel pattern, with maximum values of about $\mathrm{5.2}\times {\mathrm{10}}^{-\mathrm{4}}$ ${\mathrm{s}}^{-\mathrm{1}}$ observed around 06:00 GMT (approximately 1 h after dawn), and was anticorrelated with ambient temperature (Fig. ). BVOC emissions are driven by both light and temperature and are therefore higher during the day ; likewise, $\mathrm{NO}$ concentrations are higher during the day, due to traffic emissions. Therefore, it may be expected that measured total ozone reactivity (from both BVOCs and $\mathrm{NO}$) is higher during the daylight hours, which was in fact observed during the laboratory experiments with the lemon thyme plants (Fig. ). Under typical ambient conditions, BVOCs react during the day with $\mathrm{OH}$ radicals at a faster rate than they react with ${\mathrm{O}}_{\mathrm{3}}$ , and, as a result, ozone reactivity may be expected to peak in the early morning, when the $\mathrm{NO}$ and BVOC emissions start increasing, but the concentration of $\mathrm{OH}$ is still too low to compete with ${\mathrm{O}}_{\mathrm{3}}$ for BVOC removal . In the quasi-ambient conditions of the glasshouse, $\mathrm{OH}$ formation is also possible – but may be influenced by environment-specific factors, such as heterogeneous production of $\mathrm{HONO}$ – and may affect the diurnal variability of BVOCs in a complex manner including the behaviour shown in Fig. .

6 Summary and future work

An instrument to measure total ozone reactivity, the Total Ozone Reactivity System (TORS), was developed, characterized, and tested under controlled conditions in the laboratory; both individual compounds and small plants were used. The instrument was deployed inside a horticultural glasshouse containing a range of aromatic plants to evaluate its functionality under quasi-ambient conditions.

The TORS instrument was able to measure ${\mathrm{O}}_{\mathrm{3}}$ reactivities with BVOCs ($R_{{\mathrm{O}}_{\mathrm{3}}}$) of 4.5–$\mathrm{9.0}\times {\mathrm{10}}^{-\mathrm{5}}$ ${\mathrm{s}}^{-\mathrm{1}}$ or more – with a residence time of 140 s, an averaging time of 5 min, and an estimated total uncertainty of $\sim \mathrm{32}$ %. These values correspond to 20–40 ppb of $\mathit{\alpha }$-pinene, 150–300 ppb of isoprene, or 160–320 ppt of $\mathit{\beta }$-caryophyllene. These mixing ratios are larger than typical ambient levels, but they can be observed in BVOC-rich forested environments and in enclosure studies ; they can also be easily reproduced in laboratory and environmental chamber experiments. An $\mathrm{OH}$ scrubber (cyclohexane) was used to remove the $\mathrm{OH}$ radicals formed by the ozonolysis of BVOCs; however, simulations of the chemistry inside the TORS reactor using a Master Chemical Mechanism (MCM v3.3.1) box model found that the formation of $\mathrm{OH}$ from $\mathrm{BVOC}+{\mathrm{O}}_{\mathrm{3}}$ reactions affected the measurements of $R_{{\mathrm{O}}_{\mathrm{3}}}$ by $<\mathrm{6}$ %, under the conditions used during the experiments.

Further work will aim to improve the stability of the signal and to reduce the signal noise and the detection limit. This may require using ozone monitors with higher precision and/or a more stable ${\mathrm{O}}_{\mathrm{3}}$ generator, as well as a detailed exploration of the various parameters affecting the TORS technique: gas flows, residence time, relative humidity, $\mathrm{OH}$ scrubber levels, and ozone concentrations. Moreover, our experimental data indicate that accurate measurements of ${\mathrm{NO}}_x$ are always required to be able to interpret the TORS observations. With these improvements and proper supporting measurements, the detection limit and the uncertainty of TORS can be improved and the technique will be able to make measurements under a wider range of conditions, to advance understanding of the role of natural emissions in the ozone budget and the oxidative capacity of the atmosphere.

Data availability

The original data are available from the authors upon request.

Author contributions

RS, LJK, and WJB designed the instrument and the experiments. LRC and MSA provided the measurements of volatile organic compounds. All authors contributed to the manuscript.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

We would like to thank Daniel Blenkhorn, Dewi Komalasari, Vasileios Matthaios, Daniel Rooney, and Karn Vohra for their help with some of the experiments and the University of Birmingham Biosciences Workshop for their continuous assistance. Special thanks are due to the staff of the Winterbourne House and Garden for their support and to Jun Matsumoto for his advice and kind hospitality.

Financial support

This research has been supported by the Natural Environment Research Council (grant no. NE/P003524/1) and the UoB-Waseda Research Collaborationn Fund.

Review statement

This paper was edited by Hendrik Fuchs and reviewed by two anonymous referees.