Formaldehyde and hydroperoxide distribution around the Arabian Peninsula – evaluation of EMAC model results with ship-based measurements

. Formaldehyde (HCHO), hydrogen peroxide (H 2 O 2 ) and organic hydroperoxides (ROOH) play a key role in atmospheric oxidation processes. They act as sources and sinks for HO x radicals (OH + HO 2 ), with OH as the primary oxidant that governs the atmospheric self-cleaning capacity. Measurements of these species allow for evaluation of chemistry-transport models which need to account for multifarious source distributions, transport, complex photochemical reaction pathways and deposition processes of these species. HCHO is an intermediate during the oxidation of volatile organic compounds (VOCs) and is an indicator of photochemical activity and combustion-related emissions. In this study, we use in situ observations of HCHO, H 2 O 2 and ROOH in the marine boundary layer (MBL) to evaluate results of the general circulation model EMAC (ECHAM5


Introduction
The effects of anthropogenic emission of greenhouse gases and aerosols and their increasing impact on climate and air quality represent a global threat.Industrialization enabled the economic development to the modern society, which is characterized by urbanization and immense population growth.Large shares of agriculture and industry are coupled to the utilization of fossil fuels, and thus, emission controls and the characterization of air quality and its health impacts are of increasing importance.Globally, fossil-fuel-related emissions account for about 65 % of the excess mortality and 70 % of the climate cooling by anthropogenic aerosols (Lelieveld et al., 2019).Worldwide some of the largest hydrocarbon reservoirs are being exploited by the oil and gas industry in the Middle East, increasing also the local ship traffic drastically.The emission of volatile organic compounds (VOCs) and nitrogen oxides (NO x = NO + NO 2 ) by combustion processes make this region a hotspot of air pollution and favor the production of tropospheric ozone (O 3 ).The Arabian Peninsula is overall characterized by unique atmospheric conditions (e.g., high temperatures and intense solar irradiation, accompanied by aridity, low cloudiness and occasional dust storms), which classifies the region as a unique environment to study the abundance of atmospheric pollutants and their processing through photochemical oxidation.
The oxidation capacity of the atmosphere determines its self-cleaning ability and is mainly controlled by hydroxyl (OH) radicals in the gas phase.OH oxidizes methane (CH 4 ) and other VOCs so that these gases are efficiently removed from the atmosphere, e.g., by transition into the particle phase and subsequent rainout.Lelieveld et al. (2016) showed that global OH concentrations are buffered with a mean recycling probability of 67 %, indicating that OH is not very sensitive to perturbations by natural or anthropogenic emission changes.This buffering mechanism is based on complementary primary and secondary production of OH, e.g., through photo dissociation of ozone (O 3 ), reservoir species and radical recycling mechanisms (Lelieveld et al., 2016).OH recycling is generally dominated by the reaction of peroxy radicals with NO, which is referred to as the NO x recycling mechanism of OH.
Besides hydroxyl radicals, peroxides are a main contributor to the oxidation capacity of the atmosphere, especially in the liquid phase.Further, H 2 O 2 plays a key role in atmospheric sulfate formation and acts as a temporary reservoir for OH.With its lifetime of several hours, H 2 O 2 enables horizontal and vertical transport of HO x by, e.g., advection/convection of air masses (Nussbaumer et al., 2021a).However, H 2 O 2 also transitions readily into the liquid phase and thus also acts as a net sink for HO x radicals via its dry and wet deposition.To understand the H 2 O 2 budget and its diurnal variability, it is necessary to consider all physical and chemical processes within the atmosphere.Besides the net photochemical production (production minus loss) and deposition, horizontal and vertical transport have to be considered.The variation in the H 2 O 2 mixing ratio in the absence of clouds during the day can be described by Eq. (1) (Fischer et al., 2019): with P chem as the sum of all photochemical production terms and L chem as the sum of photochemical losses.The third term describes vertical transport in the well-mixed boundary layer, which is determined by entrainment and deposition.ω e represents the entrainment velocity with [H 2 O 2 ] as the concentration difference between the boundary layer and the free troposphere.The deposition is determined by the deposition velocity (V dep ) and the boundary layer height (h BL ).The last term describes the effect of horizontal transport on the H 2 O 2 budget due to a gradient in H 2 O 2 mixing ratios The dominant photochemical source of H 2 O 2 is the recombination of HO 2 radicals which involves a collision partner (M) usually nitrogen (N 2 ), oxygen (O 2 ) or water vapor (H 2 O): The production of H 2 O 2 via Reaction (R1) competes with the reaction of nitrogen monoxide (NO) and HO 2 in Reaction (R2), which is one of the most important reactions in the troposphere to recycle OH radicals (Lelieveld et al., 2016).
The photochemical formation of peroxides therefore depends to a large extent on the abundance of NO x , as elevated mixing ratios of NO (∼ 100 pptv and higher) substantially suppress peroxide formation (Lee et al., 2000).Photochemical loss reactions of H 2 O 2 are the conversion by OH radicals to HO 2 radicals (Reaction R3) and photolysis as a source of OH (Reaction R4).
Note that Reactions (R3) and (R4) regenerate HO x , and thus only physical removal of H 2 O 2 from the atmosphere establishes a net loss of HO x .According to its relatively high Henry's law coefficient (∼ 10 5 mol L −1 atm −1 ), H 2 O 2 is highly soluble and thus efficiently removed by rain or fog (Klippel et al., 2011;Fischer et al., 2019).Dry deposition also contributes significantly to the removal of H 2 O 2 in the boundary layer with typical deposition velocities of 0.1-5 cm s −1 (Stickler et al., 2007;Nguyen et al., 2015), which leads to a local maximum of H 2 O 2 mixing ratios above the boundary layer (Stickler et al., 2007;Klippel et al., 2011).
In the marine boundary layer (MBL), H 2 O 2 concentration gradients are small so that horizontal transport becomes unimportant.Additionally, the MBL height is relatively constant with no significant diel variation, and thus vertical transport is weak, except close to convective clouds (Nussbaumer et al., 2021a;Fischer et al., 2015).Therefore, the H 2 O 2 distribution in the MBL depends largely on net photochemical tendencies and deposition processes (Fischer et al., 2015).The H 2 O 2 budget in the continental boundary layer is more complex, since all terms in Eq. ( 1) contribute significantly to the H 2 O 2 budget and the boundary layer height follows a relatively strong diel variation.In situ observations in various locations enable (together with meteorological and boundary layer height information) assessment of the role of H 2 O 2 for the oxidizing capacity of the atmosphere (Fischer et al., 2019).Various measurement techniques have been developed to determine its vertical and geographical distribution and understand its budget and its response to natural and anthropogenic perturbations (Hottmann et al., 2020;Fischer et al., 2015Fischer et al., , 2019;;Bozem et al., 2017;Klippel et al., 2011;Snow et al., 2007;Lee et al., 2000;Sauer et al., 1997).
Similar to H 2 O 2 , organic peroxides (ROOH) impact the oxidative potential of the atmosphere significantly, and they also act as HO x reservoirs (Lee et al., 2000).Methyl hydroperoxide (MHP, CH 3 OOH) is generally the most abundant gaseous, organic hydroperoxide, which is produced by the reaction of HO 2 with methylperoxy radicals (CH 3 O 2 ) formed, e.g., during the photochemical oxidation of methane (CH 4 ) (Reactions R5, R6) or by reactions of acetyl peroxy radicals (CH 3 C(O)O 2 ) with HO 2 and NO, which can dominate the production of CH 3 O 2 (Crowley et al., 2018).Note that the production of MHP competes with the production of formaldehyde (HCHO, Reaction R7) from the methylperoxy radical (Nussbaumer et al., 2021b).Besides the photochemical pathways, H 2 O 2 and MHP have also been observed to be directly released from biomass burning (Lee et al., 1997).
The main loss reactions of MHP are its photolysis (∼ 5 × 10 −6 s −1 at sun peak) (Reaction R8), the reaction with OH (Reaction R9) (with a lifetime of ∼ 15 h for most regions during AQABA; Air Quality and Climate Change in the Arabian Basin) and physical deposition processes, although it is ∼ 2 orders of magnitude less soluble than H 2 O 2 (O'Sullivan et al., 1996;Klippel et al., 2011).
H 2 O 2 and MHP can be found in comparable concentrations in many parts of the atmosphere, with the highest variations in the boundary layer (Reeves and Penkett, 2003;Klippel et al., 2011).Besides MHP, peracetic acid (PAA) is another abundant organic hydroperoxide in the troposphere.PAA production rates depend on HO 2 and the acetyl peroxy radical (CH 3 C(O)O 2 ), which is considered one of the four most abundant organic peroxy radicals (Tyndall et al., 2001;Crowley et al., 2018).Acetyl peroxy radicals also react rapidly with NO x ; thus the highest concentrations of PAA are expected in regions which are impacted by biogenic emissions in which HO 2 levels are high enough to compete with NO x (Berasategui et al., 2020;Phillips et al., 2013).Further organic peroxides are formed in the oxidation of isoprene and other volatile organic compounds (VOCs) (Wennberg et al., 2018;St. Clair et al., 2016;Reeves and Penkett, 2003;Sauer et al., 1999;O'Sullivan et al., 1996).Recent studies also indicate the oxidative potential of isoprene hydroxyl hydroperoxides (ISOPOOH) for sulfate formation in cloud droplets, which could even surpass that of H 2 O 2 in forested regions (Dovrou et al., 2019(Dovrou et al., , 2021)).Another major HO x reservoir is formaldehyde, which is a ubiquitous trace gas and the most abundant aldehyde in the troposphere.HCHO is highly reactive and acts as a major source of HO 2 via its photolysis to H and HCO radicals.It can be emitted directly from a variety of both biogenic and anthropogenic sources and is an intermediate during the oxidation of a large number of VOCs, making budget assessments highly complex.Previous studies designed to distinguish between secondary production and direct emissions of HCHO vary widely in their estimates and highlight the importance of local phenomena (Dienhart et al., 2021;Nussbaumer et al., 2021b;Luecken et al., 2018;Anderson et al., 2017;Wolfe et al., 2016;Stickler et al., 2006;Lee et al., 1997).
Photolysis of the NO 2 product (Reactions R2, R12) leads to tropospheric O 3 formation.Further sources of HCHO are the photochemical degradation of several VOCs, e.g., the ozonolysis of isoprene and other alkenes as well as the degradation of MHP, acetaldehyde, acetone and methanol (Nussbaumer et al., 2021b;Wennberg et al., 2018;Wolfe et al., 2016;Snow et al., 2007;Stickler et al., 2006).
Since the sources of HCHO are diverse, an alternative approach to the calculation of the HCHO budget is to derive the production rate of HCHO from measurements of OH reactivity towards VOCs, as demonstrated for the Air Quality and Climate Change in the Arabian Basin (AQABA) campaign data by Dienhart et al. (2021).Photochemical losses of HCHO are the reaction with OH and its photolysis (Reactions R14-R16) (e.g., Heikes et al., 2001).
Anthropogenic release of HCHO by the oil and gas industry, biomass burning, and secondary production can significantly enhance local HO 2 production (Parrish et al., 2012;Klippel et al., 2011;Lee et al., 1997).Since the atmospheric lifetime of HCHO is at least several hours and it is released during the photochemical oxidation of numerous VOCs, it is a suitable tracer for localized emissions from anthropogenic activity and combustion processes including biomass burning.The budget of HCHO can be described similarly to H 2 O 2 via Eq.( 1): its photochemical production pathways depend strongly on the abundance and the composition of VOCs.
In the free troposphere, the main sources of HCHO are the photochemical degradation of methane, methanol and MHP (Stickler et al., 2006), whereas in the boundary layer the oxidation of alkenes (e.g., isoprene, ethene) and alkanes and the photochemical degradation of, e.g., acetaldehyde, acetone, peroxyacetyl nitrate (PAN) and dimethyl sulfide (DMS) become more significant (Crowley et al., 2018;Nussbaumer et al., 2021b).Wolfe et al. (2016) showed that the link between HCHO and isoprene oxidation is a strong, nonlinear function of NO x .Primary emissions of HCHO are dominated by combustion processes, with the combustion of fossil fuels in industrialized areas (Williams et al., 2009;Wert et al., 2003) and biomass burning as a strong local source (Kluge et al., 2020;Coggon et al., 2019).Wet and dry deposition are major loss processes of HCHO, even though it is less soluble than H 2 O 2 .The highest mixing ratios of HCHO are typically found in the boundary layer and decrease with altitude in the free troposphere (Zhu et al., 2020;Anderson et al., 2017;Stickler et al., 2007).In clean MBL conditions, HCHO mixing ratios mainly depend on the abundance of HO x and are therefore rather homogenously distributed, whereby horizontal transport is not significant.In more polluted conditions, horizontal transport can significantly influence HCHO mixing ratios on a regional scale.Vertical transport of HCHO is often limited to within the MBL, as the boundary layer height is almost constant, except close to convective clouds where elevated mixing ratios of HCHO can be used as an indicator for recent convection (Anderson et al., 2017).
In this study we present the first ship-based measurements in the marine boundary layer of the Arabian Gulf and around the Arabian Peninsula.H 2 O 2 , organic peroxides and HCHO mixing ratios were evaluated during AQABA in summer 2017 and compared to results of the 3-D general circulation model EMAC (ECHAM5/MESSy2 Atmospheric Chemistry; European Center HAMburg, Modular Earth Submodel System).Dry deposition rates of H 2 O 2 and HCHO were determined at night using the method of Shepson et al. (1992).Photochemical equilibrium concentrations of H 2 O 2 were evaluated with measured OH, HO 2 and actinic flux measurements.

AQABA campaign
The Air Quality and Climate Change in the Arabian Basin (AQABA) measurement expedition took place from 25 June until 3 September 2017.Instrumentation of the ship (Kommandor Iona) was performed in La Seyne-sur-Mer (near Toulon, France), from where the first leg of the cruise started through the Mediterranean, the Suez Canal and the Red Sea to the first stop in Jeddah.The expedition continued 2 d later via the Gulf of Aden, the Indian Ocean, the Gulf of Oman and the Arabian Gulf (also Persian Gulf) to Kuwait.On the second leg, the ship returned with the same route (Fig. 1), without stopping in Jeddah, to end the expedition at Stromboli volcano.The Kommandor Iona was equipped with a weather station and five laboratory containers on the front deck with instrumentation for in situ and offline monitoring of a large variety of gaseous species, particles and radicals.Details about the measurements performed during AQABA can be found in a number of previous publications (Dienhart et al., 2021;Friedrich et al., 2021;Paris et al., 2021;Celik et al., 2020;Tadic et al., 2020;Wang et al., 2020;Bourtsoukidis et al., 2019;Pfannerstil et al., 2019;Eger et al., 2019) HCHO and hydroperoxides were measured using modified commercial Aero-Laser instruments (AL2021, AL4021, Aero-Laser, Germany), which were placed in a temperaturecontrolled container.With the exception of the aerosol and radical measurements (OH and HO 2 ), air was sampled from a high-flow (10 m 3 min −1 ) cylindrical stainless steel inlet (sampling height: 5.5 m above deck, diameter: 0.2 m), placed between the containers on the front deck of the ship.Air was drawn from the center of the high flow inlet into the air-conditioned laboratory containers using PFA (perfluoroalkoxy alkane) tubing.The 4.2 m long 1/2 in.PFA bypass was insulated to prevent condensation and was used with a flow rate of 12 L min −1 , which resulted in a residence time of ∼ 9 s for both instruments.This setup ensured no vessel contamination while sampling against the wind direction and minimized sampling artifacts, e.g., by preventing condensation.The sampling bypass was exchanged in Kuwait before the second leg.

HCHO measurements
HCHO measurements were performed based on the fluorometric Hantzsch reagent method (AL4021 is therefore called the "Hantzsch monitor") following the principle of Dasgupta et al. (1988) and the design of Kelly and Fortune (1994).In a first step, HCHO is stripped from an airflow of 1 L min −1 into 0.025 M H 2 SO 4 (sulfuric acid for analysis, 96 %, Acros Organics) with a flow of 0.55 mL min −1 at 10  Lazrus et al. (1985).Ambient air is collected through a bypass with 2.3 L min −1 and consequently passed through a glass coil together with a buffered (potassium hydrogen phthalate for analysis, PanReac; NaOH, 1 mol L −1 , Fluka) stripping solution (0.55 mL min −1 , pH 5.8).Hydroperoxides dissolve in the stripping solution with a stripping efficiency depending on their Henry's law constant (O'Sullivan et al., 1996).Typically, H 2 O 2 is dissolved quantitatively; CH 3 OOH (methyl hydroperoxide, MHP) is the smallest organic hydroperoxide, with a stripping efficiency of ∼ 60 % (Hottmann et al., 2020;Klippel et al., 2011).As the instrument does not differentiate between different organic hydroperoxides and as solvation is a critical step for quantification, the AL2021 delivers a lower estimate of the total organic hydroperoxide mixing ratios.The dissolved hydroperoxides are separated into two channels and subsequently detected via reaction to a fluorescent dye with horseradish peroxidase (HRP, Sigma-Aldrich) and 4-hydroxyphenylacetic acid (POPHA, Sigma-Aldrich).The dimer of POPHA, 6,6'dihydroxy-3,3'-biphenyldiacetic acid, is formed stoichiometrically and detected by fluorescence spectroscopy via excitation with a cadmium Pen-Ray ® lamp at 326 nm.Detection of the fluorescence radiation 90 • off axis is performed between 400-420 nm with a photomultiplier tube (model H957-01, Hamamatsu Photonics) for both channels.The enzyme catalase (Sigma-Aldrich) is injected into the reaction coil of channel B, prior to the reaction with HRP and POPHA, to selectively destroy H 2 O 2 .This technique allows for quantification of H 2 O 2 by calculation of the difference between channel A, which delivers the total mixing ratio of ROOH and H 2 O 2 , and channel B, which delivers the total mixing ratio of ROOH.Since this principle is dependent on the catalase efficiency, it is determined for every liquid calibration and was in the range of 95 %-100 % during AQABA.
In addition to the AL2021, we also operated an instrument for the detection of different organic peroxides separated by HPLC (high-performance liquid chromatography).Similar to the AL2021, it utilizes the selective dual-enzyme technique by post-column derivatization, and thus the HPLC enables quantification of separated organic hydroperoxides and H 2 O 2 in low levels of parts per trillion by volume.When the sea was rough, the movement of the ship interfered with the instrument, causing drifts of the baseline, which may have been caused by pressure variations within the constant-flow eluent pumps.Therefore, quantification of the organic hydroperoxides was not possible, and we only used the chromatograms for qualitative identification of the more abundant species.
In addition, gaseous standards were measured to calculate the inlet efficiency of the PFA bypass.Calibration gas flows were generated using permeation devices in temperaturecontrolled glass flasks, which were flushed at a constant flow rate of 80 sccm (standard cubic centimeters per minute) with zero air (zero-air generator CAP 60, Infiltec, Germany).HCHO calibration gas was created from a paraformaldehyde container (VICI AG, Switzerland) which was heated to 60 • C. The gaseous H 2 O 2 standard was generated from a permeation source built with a 15 cm long 1/8 in.polyethylene (PE) tube, which was filled with the 30 % H 2 O 2 solution, closed with PFA fittings (Swagelok, USA) and heated to 35 • C. The highly concentrated flow was then diluted with additional zero air.The permeation rates of both sources were measured based on the chromotropic acid reaction (Altshuller et al., 1961) and the reaction of H 2 O 2 with TiCl 4 described in Pilz and Johann (1974).Note that the AL2021 has known O 3 and NO interferences, which were accounted for in the final dataset.We found an interference of 36 pptv H 2 O 2 equivalents per 100 ppbv O 3 and 12 pptv H 2 O 2 per 100 ppbv NO.We did not find a significant O 3 interference in lab experiments for the AL4021.
Zero-gas measurements were performed every 3.5 h for 30 min to account for baseline drifts and to determine the instrument's stability.For this purpose, we used a bypass via a three-way valve with a silica gel cartridge (SiO 2 with orange indicator, Roth) to dry the sampled air followed by a scrubber cartridge containing hopcalite (MnO 2 /CuO, IAC-330, Infiltec, Germany) and platinum (Pt/Al 2 O 3 , IAC-114, Infiltec, Germany) as catalysts to destroy the remaining hydroperoxides, HCHO, other OVOCs (oxygenated VOCs) and O 3 .
Both instruments log data on a custom-built computing unit (V25) with a 3 s time resolution, but the data shown in this paper were at least averaged to the so-called effective time resolution, which was determined as the response time of the instrument (10 % to 90 % of the signal intensity during the injection of liquid standards).The limit of detection (LOD) was calculated as the 2σ deviation of all zero-air measurements during AQABA at the effective time resolution of 180 s.The precision (P ) was calculated by the 1σ deviation of the liquid standard calibrations throughout the whole measurement campaign; therefore it contains also the pipetting error during the preparation of the standards.The total measurement uncertainty (TMU) was calculated according to Gaussian error propagation.In this equation, S is the uncertainty of the standard, IE is the inlet efficiency and OI is the O 3 interference.
2.6 Further measurements OH and HO 2 were performed with the Hydroxyl Radical Measurement Unit based on Fluorescence Spectroscopy (HORUS) instrument and sampled from a separated inlet closer to the LIF (laser-induced fluorescence) instrument to achieve residence times as low as possible in the sampling.
The instrument utilizes LIF of the OH radical at 308 nm, which is created by a pulsed tunable laser system (Nd:YAG, neodymium-doped yttrium aluminum garnet) operated at a pulse frequency of 3 kHz.OH radicals are excited in a lowpressure detection cell (White cell setup) with a flow rate of 10 L min −1 .The detection of HO 2 is achieved by chemical conversion via the addition of NO downstream of the OH detection.The resulting sum of ambient plus chemically converted OH is measured in a second detection axis (Hens et al., 2014;Marno et al., 2020).HORUS achieved an instrumental LOD of 0.03-0.15pptv for OH and 0.22-2.01pptv for HO 2 with a TMU of 17 % (OH) and 20 % (HO 2 ), respectively.Wavelength-resolved actinic flux was measured with a spectral photometer (CCD spectroradiometer 85237; chargecoupled device) to calculate photolysis frequencies (j values).The radiometer was installed about 5 m above the front deck level, and it was cleaned every morning to remove sea salt and dust particles.Decreases in sensitivity due to sensor contamination were corrected with a linear interpolation between the cleaning events.The j values are not corrected for upwelling actinic flux from the sea surface, and therefore the TMU was estimated with > 10 %, depending on the reaction.More details about the setup and calibration of CCD spectroradiometers can be found in Bohn and Lohse (2017).Temperature, pressure, and wind direction and speed were measured with the European Common Automatic Weather Station (EUCAWS), a weather station specifically designed for ships.

Global atmospheric chemistry model EMAC
The ECHAM5/MESSy2 Atmospheric Chemistry (EMAC) model is a numerical chemistry general circulation model (CGCM), which describes tropospheric and middleatmosphere processes.EMAC is based on the fifth generation of the European Center HAMburg (ECHAM5) general circulation model (Roeckner et al., 2006) and uses the second version of the Modular Earth Submodel System (MESSy2) to link multi-institutional submodels (Joeckel et al., 2010).
Here we use EMAC with the chemistry mechanism MOM (Mainz Organics Mechanism) implemented with the submodel MECCA (Module Efficiently Calculating the Chemistry of the Atmosphere), which contains not only the basic HO x , NO x and CH 4 chemistry but also halogens, sulfur and mercury (Sander et al., 2019;Lelieveld et al., 2016).Development of MOM also included a variety of nonmethane hydrocarbons (NMHCs), aromatics and OVOCs including isoprene and terpene oxidation (Sander et al., 2019), and, re- Note that for the LOD of the AL4021, we found a significant change in the background noise, while operating the instrument when the sea was rough with strong waves.Excluding times of rough sea yields a LOD (2σ ) of 80 pptv.
cently, the model has been thoroughly evaluated with this chemical mechanism (Pozzer et al., 2022).Therefore, it is ideal to test the model with complex photo-oxidation during AQABA, especially in the Arabian Gulf, where a lot of the oil and gas industry is operating and which model results already identified as a hotspot of tropospheric O 3 (Lelieveld et al., 2009).The model simulations were carried out in the T106L31 resolution, which correspond to a grid of 1.1 • × 1.1 • (∼ 110 km) with 31 vertical pressure layers and a time resolution of 10 min.The EMAC data were interpolated bi-linearly along the GPS track of the ship with the S4D submodel (sampling in 4 dimensions; Jöckel et al., 2010).The model was initialized from a previous evaluated simulation (Pozzer et al., 2022) and started on 1 June 2017, covering the entire campaign.The dynamics have been weakly nudged (Jeuken et al., 1996;Jöckel et al., 2006) towards the ERA-Interim data (Berrisford et al., 2011) of the European Centre for Medium-Range Weather Forecasts (ECMWF) to reproduce the actual day-to-day meteorology in the troposphere.The anthropogenic emissions are based on the Emissions Database for Global Atmospheric Research (EDGAR v4.3.2).Further details are presented in Pozzer et al. (2022).Previous results of airborne and shipborne expeditions have been compared to EMAC (Fischer et al., 2015;Klippel et al., 2011); also the AQABA datasets of NO x , O 3 and VOCs during AQABA have been published (Tadic et al., 2020;Wang et al., 2020).

Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model
Back trajectories of air parcels encountered during the AQABA campaign have been calculated with the Hybrid Single-Particle Lagrangian Integrated Trajectory model, HYSPLIT (version 4, 2014), developed by the Air Resources Laboratory of the National Oceanic and Atmospheric Administration (NOAA).HYSPLIT is a hybrid between a Lagrangian and Eulerian model for the simulation of the movement and dispersion of small imaginary air parcels within the atmosphere (Draxler et al., 1998).The model can be accessed at https://www.arl.noaa.gov/hysplit/(last access: 15 November 2022).Back trajectories for the AQABA campaign were calculated with a start height of 200 m above sea level at the location of the Kommandor Iona with a resolution of 1 h (Figs.S17, S18).For the purpose of this study, the movement of air parcels was tracked for 48 h backwards.

Methods
In the MBL, the production of peroxides is dominated by the recombination of peroxyl radicals; thus the in situ measurements of HO 2 enable calculation of the production rates P (H 2 O 2 ) with Eq. ( 4).
Photochemical loss reactions are the photolysis and the reaction with OH, besides the deposition, which is the dominant loss at night.
Since the measurements were performed in the MBL, the water dependency of the HO 2 recombination becomes significant and was calculated with the relative humidity (RH) measurements via Eqs.( 10)-( 12) (https://iupac-aeris.ipsl.fr/test/front-office/datasheets/pdf/HOx14.pdf,21 July 2022).et al. (1992).In this calculation, we assume that the initial mixing ratio [X] 0 (j NO 2 < 10 −3 ) represents the mixing ratio on top of the nocturnal boundary layer (h BL ).If species X is homogeneously distributed within the BL and it follows exponential decay at night, the first-order decay is given by Eq. ( 11).Thus, the first-order decay plot (ln[X] t /[X] 0 versus the time) yields the deposition velocity of species X with a known boundary layer height.
Additionally, the deposition rate can be calculated with Eq. ( 12), assuming that the boundary layer is well mixed.
3 Results 3.1 Regional distribution of HCHO, H 2 O 2 and organic peroxides around the Arabian Peninsula The cruise track of the Kommandor Iona is shown in Fig. 1, subdivided into eight regions identified by different colors: the Mediterranean Sea (MS), Suez Canal and the Gulf of Suez (SU), Red Sea north (RN), Red Sea south (RS), Gulf of Aden (GA), Arabian Sea (AS), Gulf of Oman (GO), and Arabian Gulf (AG).The AG (also known as the Persian Gulf) and the SU are well known for their oil and gas industry and intensive ship traffic, respectively; hence primary emissions of NO x , CO and to a lesser extent HCHO were expected to affect mixing ratios of these trace gases.Nearby ship plumes (including the Kommandor Iona exhaust) and other point sources were identified with the use of NO x , CO, SO 2 and wind direction data and excluded from the dataset used in this study (Celik et al., 2020).The measurements were affected by aft winds in particular during the first leg, resulting in a contamination from the ships exhaust, thus limiting the amount of data available (Fig. S1 in the Supplement).In general, we did not find elevated mixing ratios of peroxides in ship plumes (presumably a result of the high NO x levels) and thus decided to keep the data, but corrections for known NO interferences were applied (Sect.2.4).Contrary to H 2 O 2 , HCHO is affected by ship exhaust plumes (Celik et al., 2020), and thus contaminated data have been filtered out.
In Fig. 2 we present mixing ratios of HCHO (upper panels), H 2 O 2 (middle panels) and ROOH (lower panels) colorcoded along the ship cruise track for the first (left panels) and second leg (right panels).Time series of these species can also be found in the Supplement (Figs.S1, S3 and S5).Box-and-whisker plots of the mixing ratios for HCHO, H 2 O 2 and ROOH for the eight regions are shown in Figs. 3 and 4. Numerical values are listed in Table S1 in the Supplement.
Mixing ratios of HCHO (upper panels of Figs. 2, S1) exhibit a large variability.The highest mixing ratios (12.6 ppbv) were measured in the center of the Arabian Gulf with northwesterly winds originating from Iraq/Kuwait during the first leg (Figs. 2, S18).Lower mixing ratios were detected in this area during the second leg, when the wind originated from the northeast, coming from Iran (Fig. 2, upper-right panel).The lowest HCHO median mixing ratios were measured in the RS (0.37 ppbv) during the second leg, in unpolluted air mass originating from Eritrea (Fig. S17).Low HCHO was also found over the GA (0.50 ppbv), the MS (0.77 ppbv) and the AS (0.86 ppbv).In general, low mixing ratios of HCHO are associated with low NO x , low O 3 (Tadic et al., 2021), low VOCs (Bourtsoukidis et al., 2019), low OH reactivity (Pfannerstill et al., 2019) and in particular low OH reactivity towards VOCs (Dienhart et al., 2021), while high mixing ratios of HCHO are associated with elevated values for these species.
To the best of our knowledge there are no shipborne measurements of HCHO available in the Red Sea and the Arabian Gulf to be compared to our data.In general, the measured mean mixing ratios during AQABA are generally larger compared to previous studies in the MBL.Wagner et al. (2001) performed shipborne measurements during the Indian Ocean Experiment (INDOEX) campaign in the central Indian Ocean with HCHO mixing ratios between 0.2-0.5 ppbv, with the lowest mixing ratios in the clean maritime background of the Southern Hemisphere and about 0.5 ppbv HCHO in continentally influenced air masses.Weller et al. (2000) reported ship-based HCHO measurements in the Atlantic, which reached a broad maximum with values of 1.0-1.2ppbv in the tropical Atlantic, decreasing towards the poles with values below 0.8 ppbv.These air masses represented pristine MBL conditions with average daytime NO of 3.1 pptv.During AQABA we did not encounter such very clean conditions with the lowest median NO x of 0.19 ppbv for the AS and 0.25 ppbv for the MS.This was likely due to sailing on major ship traffic routes which may have also led to enhanced background HCHO compared to remote MBL conditions.
The lowest median mixing ratios of H 2 O 2 were found in the GO and the GA (0.12 ppbv), followed by the AS (0.15 ppbv), while higher mixing ratios of H 2 O 2 were found in the RS, SU (0.25 ppbv) and MS (0.26 ppbv).Altogether, the H 2 O 2 measurements exhibit lower variation around the Arabian Peninsula compared to HCHO.Higher variability in H 2 O 2 was found in the SU and AG, although less than 25 % of the data exceed 0.50 ppbv, with the highest H 2 O 2 mixing ratios observed in the AG (0.92 ppbv) in the harbor of Kuwait.The AG, SU, MS and RN also show the strongest diurnal variations of up to ∼ 300 pptv (Fig. S3).
The organic peroxides showed higher variability compared to H 2 O 2 (Fig. 4), with the lowest median value in the AS (0.06 ppbv), followed by the GO (0.07 ppbv) and GA (0.10 ppbv).We found the lowest variability (whisker intervals) in the AS, which represents the cleanest conditions and the lowest variability in O 3 and NO x (Tadic et al., 2020).Higher levels of organic peroxides were detected in the SU (0.26 ppbv), AG (0.23 ppbv), MS (0.22 ppbv) and RN (0.20 ppbv) with the highest mixing ratios in the center of the Arabian Gulf during the first leg (2.26 ppbv).
The chromatograms of the HPLC-based instrument indicate significant abundances of four distinct inorganic and organic hydroperoxides in the AG (Fig. S10), which were identified as H 2 O 2 , MHP, PAA (peracetic acid) and EHP (ethyl hydroperoxide) based on their retention times and gaseous injections of PAA with a diffusion source.In addition to the continuous HPLC measurements, we also injected en-riched samples with sampling times varying between 12-36 h during various times along the ship track.Although these samples have a very limited time resolution, they were used for a qualitative assessment of the abundance of further organic hydroperoxides.Significantly enhanced amounts of EHP were only detected over the Arabian Gulf, although small amounts of EHP were also detected in the enriched samples of the MS (Fig. S10), where we detected aged air masses originating in Europe (Fig. S17).
The highest amounts of photochemical air pollution were detected over the AG, which is confirmed by the highest mixing ratios of HCHO and ROOH in this region when observing winds from the western coastline and Kuwait (Fig. S19).Less air pollution was observed during the second leg, when we were sampling air masses originating from Iran.In this region we also observed the strongest radiation and the highest temperatures during AQABA.The SU and RN also show enhanced contributions of VOCs (Wang et al., 2020) and elevated OH reactivity (Pfannerstill et al., 2019), mainly while passing oil rigs and on the way through the Suez Canal.Tadic et al. (2020) found the cleanest conditions, from both a NO x and O 3 perspective, for the AS and RS.This can be confirmed by the rather decreased mixing ratios of HCHO, ROOH and H 2 O 2 mixing ratios, which reflect low levels of HO x .Air masses transported from Eritrea also contained suppressed mixing ratios of ROOH during the first leg.In general, EMAC reproduces the regional trends of HCHO quite well (Fig. 3, left panel), although the whisker intervals (10 %-90 %) demonstrate that the model computes less variability, most likely due to missing localized pollution events, which is to be expected for a global model at relatively low spatial resolution.The median values of the measurement-model ratio vary between 0.3 and 1.6 (Fig. 3, right panel).EMAC overestimates HCHO in the cleaner regions of the MS, RS and GA, while it underestimates HCHO in the AS and the heavily polluted AG (Fig. S1).
Elevated NO x and O 3 measurements classified the AG and SU as the most polluted regions, followed by the RN and GO, which both were influenced to a higher extent by anthropogenic pollution (Tadic et al., 2020).The elevated HCHO during the high-ozone and VOC conditions on the first leg in the AG was not reproduced by EMAC.Here the model clearly underestimates HCHO (Fig. S1) and ROOH (Fig. S5).Since the elevated OVOCs correlated well with CO and O 3 (Wang et al., 2020), we can assume that we probed a highly polluted and photochemically active air mass, with both effective photochemical production and primary emissions of HCHO.Paris et al. (2021) identified natural-gas flaring as a major source of the elevated VOCs over the Arabian Gulf.During this event, maximum HCHO and ROOH mixing ratios were measured during AQABA, with values up to 12.6 ppbv HCHO and 2.26 ppbv ROOH in the center of the gulf (Figs. 2, S9, S19).The event was less pronounced in EMAC with up to 3.31 ppbv HCHO leading to an underestimation over the AG by about a factor of 4. EMAC does not simulate significantly elevated values of ROOH peaking at 0.49 ppbv with elevated contributions of PAA and EHP.Even though the model underestimates ROOH also by about a factor of 4 (Fig. S9).Wang et al. (2020) showed that EMAC simulates enhanced acetone and methyl ethyl ketone (MEK) during this event, although the model shows no significant increase in acetaldehyde.
In the Suez Canal and the Gulf of Suez (SU), the second most polluted region, the model also underestimates HCHO mixing ratios, even though EMAC simulates a significant increase in HCHO compared to the MS and RN.Wang et al. (2020) showed that these air masses have been influenced by biomass burning and increased anthropogenic emissions, e.g., by gas flaring similar to the Arabian Gulf.
EMAC overestimates HCHO in the less polluted regions of the MS, RS and GA.Tadic et al. (2020) also found that the simulations overestimate NO x and O 3 in these areas.EMAC significantly underestimates HCHO for the AS, especially at night (Fig. S1).Here, a well-pronounced diurnal cycle is simulated, while the observations indicate only a distinct diurnal variation in the eastern part of the AS, which got stronger in the GA (Fig. S1).Slightly elevated mixing ratios were observed in the AS, compared to clean MBL conditions, e.g., during INDOEX (Wagner et al., 2001), and are most likely caused by other primary sources and oxidation of further VOCs, not by methane oxidation only.Previous observations in the remote MBL showed significantly lower HCHO in the range of 0.2-0.4ppbv (Wagner et al., 2001).Thus, we assume that the air masses encountered in the AS during AQABA were still influenced by anthropogenic pollution, which is supported by the elevated NO x (Tadic et al., 2020) and acrolein as the main contributor to OH reactivity in the AS and GA (Pfannerstill et al., 2019).An additional source for HCHO was the ozonolysis of ethene, which reached maximum values of 0.24 ppbv with a median of 0.08 ppbv in the AS (Bourtsoukidis et al., 2019).Wang et al. (2020) also show strongly enhanced acetaldehyde in the AS relative to model simulations, indicating a missing oceanic source in the model.Tripathi et al. ( 2020) also performed VOC and seawater measurements of phytoplankton species, which demonstrate the high biological activity in the region.They determined elevated ethene (8.92 ± 3.50 ppbv) and propene (3.38±1.30ppbv) in marine air originating from the Arabian Sea.Just recently, Tegtmeier et al. (2022) highlighted the complexity of the air phase composition over the Indian Ocean, with the major differences between the Indian summer and winter monsoon.
Altogether, EMAC reproduces observed HCHO mixing ratios on average within a factor of 2. However, the model clearly underestimates air pollution over the Arabian Gulf, which leads to the assumption of missing sources in EMAC and may also be related to the limited resolution of EMAC.The model simulates the AS cleaner than it was observed, but the relatively low levels of HCHO over the RS with winds from Eritrea (Figs.S2, S17) were overestimated.Analysis of the air mass origin showed that air masses over the Arabian Sea represented clean and aged conditions transported from the center of the Indian Ocean (Fig. S18).Given the multitude of potential HCHO sources both from direct emissions and a large variety of photochemical precursors and the limited resolution of EMAC, the agreement within a factor of 2 is satisfactory.The comparison of simulated HCHO based on a more highly resolved model (WRF-Chem; Weather Research and Forecasting) did not improve the accuracy of HCHO and was thus not included in the paper.
The measurement-model comparison for H 2 O 2 reveals that EMAC systematically overestimates H 2 O 2 mixing ratios by up to an order of magnitude.With the exception of the Arabian Sea (0.7 ppbv), the model predicts H 2 O 2 with median mixing ratios in excess of 1.3 ppbv, with the highest median values of 2.0 and 1.8 ppbv for the SU and GO, respectively.While the measurements cover a whisker range (10 % to 90 % of the data) of only 0.1 to 0.4 ppbv, this consistent overestimation by the model indicates either a significant overestimation of H 2 O 2 sources or missing sinks in the model or a combination of both.
In order to compare the observations of ROOH with model results, we summed up individual simulated organic hydroperoxide species, which were identified in the qualitative HPLC measurements (Fig. S4): methyl hydroperoxide (MHP), peracetic acid (PAA) and ethyl hydroperoxide https://doi.org/10.5194/acp-23-119-2023Atmos.Chem.Phys., 23, 119-142, 2023  (EHP).Please note that measured ROOH is a lower limit of the sum of organic hydroperoxides, since different sampling efficiencies for the individual species, which depend on the Henry's law constants, are not accounted for.It can be assumed that MHP is the dominant contributor to the total organic hydroperoxides in the clean MBL, which has a sampling efficiency of only 60 % (Fischer et al., 2015).In remote areas, this would lead to an underestimation of measured ROOH by a factor of 0.6.With significant contributions of higher organic hydroperoxides, which are generally more soluble, this underestimation tends to be smaller.EMAC also tends to overestimate the organic hydroperoxides, with the lowest median value of 0.26 ppbv in the MS (observations of 0.22 ppbv) and the highest of 0.42 ppbv in the AS (observations of 0.06 ppbv) (Fig. 4 and Table S1).For the whole dataset, the simulated ROOH cover a whisker range (10 % to 90 % of the data) between 0.19 and 0.50 ppbv, while the observations span 0.04 to 0.68 ppbv.Although the measurements can be reproduced within the 25 %-75 % box range in some regions, median values differ between about a factor of 1 to 7 between the simulations and observations.Please note that the AG was the only region where we measured four separated hydroperoxides in the in situ results of the HPLC, with the largest contribution of MHP and EHP.This enhancement was also found in the EMAC results, as EHP and PAA mixing ratios increased in the AG, especially during the high-pollution events of the first leg (Fig. S10).Although the simulations of ROOH match the observations better than H 2 O 2 , EMAC overestimates the organic peroxides, especially in the clean regions while cruising close to the coast.The AG shows the highest variability in ROOH, which is to be expected due to the complex photochemistry of VOCs (Bourtsoukidis et al., 2019;Pfannerstill et al., 2019).The model simulates strong diel cycles (∼ 300 pptv) throughout the whole dataset, while the observations only indicate comparable variations for the RN, SU, MS and AG.A distinct decline in ROOH mixing ratios was observed for the AS and GO, which is not reproduced by EMAC.
Possible explanations for the systematic overestimation of both H 2 O 2 and organic hydroperoxides by EMAC can be an overestimation of photochemical sources or an underestimation of loss processes in the model or due to the unknown sampling efficiencies of the organic peroxides.To investigate photochemical misrepresentations, we compared the observed and simulated OH and HO 2 daytime values (j NO 2 ≥ 10 −3 s −1 ; Figs. 5, S7, S8).Since the source term of H 2 O 2 depends quadratically on HO 2 concentrations (Eq.3), simulations of H 2 O 2 are highly sensitive to HO 2 , while its photochemical loss scales linearly with OH (Eq.4).The model overestimates both OH and HO 2 throughout the whole campaign.The highest observed HO 2 median values were found in the SU (19.3 pptv) and RN (16.3 pptv), followed by the RS (14.8 pptv) and MS (14.4 pptv).The remaining regions show significantly less HO 2 with the smallest median value for the GO (4.5 pptv).The smallest whisker ranges of the dataset in the AS and GO demonstrate suppressed HO 2 in these regions.Surprisingly small mixing ratios were also detected over the polluted Arabian Gulf (6.8 pptv).
According to the observations, EMAC simulates the highest HO 2 median values for the SU (27.9 pptv) and RN (26.1 pptv), while mixing ratios in the MS (19.4 pptv) and RS (11.2 pptv) are smaller.Significant overestimation of HO 2 was found for the GO (23.1 pptv) and AG (19.0 pptv), showing similar enhancements of HO 2 .In the rather clean regions during AQABA, e.g., GA (16.7 pptv) and AS (19.6 pptv), the model generates significantly enhanced daytime HO 2 compared to the observations, while diurnal variation in HO 2 matches the observations in the RS and MS (Fig. S8).Altogether, we examined an overestimation of HO 2 by about a factor of 2 (EMAC dataset was aligned to the observations), with average daytime mixing ratios of 11.3 pptv for the observations and 19.7 pptv for the EMAC simulation, respectively.
The highest median OH values were observed for the SU (0.13 pptv), RN (0.13 pptv), MS (0.12 pptv) and GO (0.11 pptv).Slightly less OH was detected in the RS (0.07 pptv), AS (0.05 pptv) and AG (0.05 pptv).The boxand-whisker plot ranges indicate the highest variations in OH in the GA, where we detected the highest OH mixing ratios of 0.6 pptv close to Bab-el-Mandeb on 16 August 2017 (Fig. S7).EMAC simulates the highest OH mixing ratios in the SU (0.49 pptv), GO (0.47 pptv) and RN (0.44 pptv), while the lowest median values of OH were simulated in the RS (0.13 pptv), followed by the MS (0.26 pptv) and AS (0.26 pptv).
The EMAC results for OH follow similar regional trends compared to HO 2 , although they reflect a stronger pronounced overestimation, since the measurements do not reflect a substantial increase in OH in the MS, SU, RN and RS.This results in overestimated daily median values within a factor of 2 to 5.

Photochemical production and loss of H 2 O 2
The comparison to EMAC in the previous section showed that overestimations of HO x by EMAC affect simulations of HCHO, H 2 O 2 and ROOH.How misrepresentation of HO x will affect HCHO is complex due to the many HCHO sources and the fact that both sources and sinks are strongly related to OH concentrations.Therefore, we will concentrate on H 2 O 2 in the following calculation of its photochemical production and loss terms in order to evaluate the discrepancy between modeled and measured H 2 O 2 .H 2 O 2 is highly sensitive towards deviations of HO x between the observations and the model, as its production depends quadratically on HO 2 (Eq.3) but its loss only linearly on OH (Eq.4).Thus, H 2 O 2 can be used to evaluate the discrepancy between measured and modeled HO x .
The daytime production rates of H 2 O 2 (j NO 2 ≥ 10 −3 s −1 ) are displayed in Fig. 6 in order to compare the results of the observations and EMAC, with the corresponding timelines presented in Fig. S11 and a scatterplot of k HO 2 +HO 2 in Fig. S12.Based on the observations, the highest production rates were found in the SU with a median production of 202 pptv h −1 and the largest whisker range of up to ∼ 480 pptv h −1 .A similar range was simulated by EMAC with an increased median value of 298 pptv h −1 .The MS and RN demonstrate comparable results, where the box ranges of the observations and the model agree to some extent, although the overestimated HO 2 by EMAC outweighs the slightly smaller reaction constant of k HO 2 +HO 2 (Fig. S12).Overall, EMAC tends to overestimate P (H 2 O 2 ), except in the RS, where the model correctly simulates the lowest H 2 O 2 production rates of the dataset with a median value of 63 pptv h −1 -about a factor of 5 lower than in the SU.However, the observations display reduced P (H 2 O 2 ) in the remaining regions, which resulted in a stronger pronounced https://doi.org/10.5194/acp-23-119-2023Atmos.Chem.Phys., 23, 119-142, 2023 discrepancy with at least a factor of 5 for the GA, AS, GO and AG.
Besides dry deposition, photochemical losses of H 2 O 2 are the reaction with OH and photolysis, which were calculated according to Eq. ( 4) and displayed without multiplication of the H 2 O 2 mixing ratio for the sake of comparability (Fig. 6).The simulated photolysis rate constants (j values) demonstrate good agreement with the observations with a factor of 1.2, during the rarely cloudy conditions of AQABA.EMAC only overestimates j H 2 O 2 for air masses very close to the coastline, e.g., in the SU (up to ∼ 1.5×10 −6 s −1 ) and to a lesser extent in the GO, while the photolysis rate was underestimated in the RS (Fig. S11).Loss of H 2 O 2 due to photolysis was less important than the reaction with OH for most regions, while photolysis prevailed for the AS, RS and AG.In contrast to the observations, photochemical losses of H 2 O 2 were dominated by the reaction with OH and were overestimated within a factor of ∼ 2-5 by the model with the best agreement in the RS.Loss due to dry deposition could not be determined during the day, but nighttime deposition velocities are calculated in Sect.3.4.
To further put these results into perspective, Fig. 7 presents the net photochemical production of H 2 O 2 (P (H 2 O 2 )-L(H 2 O 2 )), whereby the loss rate constants were multiplied with the H 2 O 2 mixing ratio.Please note that in this term neither physical loss processes (e.g., deposition) nor transport are represented, and thus it only reflects the effect of local photochemistry on the H 2 O 2 mixing ratio.The slightly overestimated photochemical production in the MS, SU and RN by EMAC is compensated by elevated losses via reaction with OH in these regions so that both datasets agree well in the MS and RN and demonstrate a strongly pronounced diurnal variation peaking at ∼ 580 and ∼ 500 pptv h −1 , respectively (Fig. 8).Net photochemical production of H 2 O 2 outweighs the model results in the SU and RS, especially at noon (Fig. 8).The GA, AS, GO and AG remain overestimated by EMAC so that the elevated losses in the model do not compensate for the enhanced production rates due to the quadratic dependence on the HO 2 concentration.The observations demonstrate less net photochemical production in the GA, AS, GO and AG with noontime values below 200 pptv h −1 due to the decreased HO 2 .Pfannerstill et al. (2019) reported the highest OH reactivity in the AG (11.6 s −1 ) and SU (10.4-10.8s −1 ), comparable results for the GO (8.4 s −1 ) and GA (8.0 s −1 ), and the lowest OH reactivity for the AS (4.9 s −1 ).Air masses in the AG demonstrated by far the highest contribution of reactivity towards OVOCs (∼ 40 %), alkanes and alkenes (together ∼ 14 %).Air masses in the GO showed slightly higher contributions of reactions with NO x compared to the AG, while the AS represents the cleanest conditions during AQABA with respect to NO x (Tadic et al., 2020).A potential explanation for the surprisingly low HO 2 mixing ratios in the AG could be suppressed OH recycling by means of enhanced organic peroxy radicals (RO 2 ) and in general a large contribution of reactions with OVOCs, alkanes, alkenes and aromatics, as OH recycling through these reactions is slower compared to NO x recycling (via HO 2 + NO).Enhanced ROOH indicates a larger contribution of RO x reactions (= OH + HO 2 + RO 2 ), which would also slow down OH recycling.There are no measurements of organic peroxy radicals available, but Tadic et al. ( 2020) calculated noontime estimates of HO 2 + RO 2 , with the highest noontime median values of ∼ 75 pptv in the AG (see Fig. 7 of Tadic et al., 2020).The remaining regions of AQABA show noontime median values in the range of ∼ 10-35 pptv.Elevated RO 2 in the AG is also supported by the enhanced mixing ratios of ROOH, the only region where we detected MHP, PAA and EHP in the in situ measurements of the HPLC (Fig. S10).
The decreased observations of net photochemical production of H 2 O 2 in the AS, GO and GA are generally caused by lower HO x mixing ratios, which are not reproduced by EMAC.The box range indicates similar values of net photochemical production in the MS, which agree mostly within a factor of 2. However, the H 2 O 2 observations display an average diurnal variation of ∼ 0.2 ppbv with the highest mean Figure 6.Box plots of H 2 O 2 production rates (pptv h −1 ) for the observations (grey) and the EMAC model (yellow) during daytime (j NO 2 ≥ 10 −3 s −1 ) and the corresponding loss rate constants (s −1 ) for photolysis (observations in grey and EMAC in yellow) and reaction with OH (green, orange).The box represents 25 % to 75 % of the data, and the whiskers are from 10 % to 90 % with the median as the black line.The corresponding timelines are presented in Fig. S11.Note that the deposition term from Eq. ( 4) is not included in the calculation of the losses.4) is not included in this calculation.
mixing ratios of ∼ 0.4 ppbv at 12:00 UTC, while the EMAC results indicate a variation of ∼ 1.2 ppbv with the highest mixing ratios of ∼ 2.1 ppbv also at 12:00 UTC (Fig. S13).In general, the observations only display weak diurnal variations compared to the large diurnal variation in net photochemical production, which implies that other loss processes (e.g., deposition, wash out due to sea spray) contributed significantly to the diurnal variability in H 2 O 2 in the marine boundary layer.

Dry deposition
An alternative reason for deviations between H 2 O 2 observations and EMAC predictions could lie in the physical processes of deposition to the ocean surface and entrainment through the top of the MBL.Both processes are related to the absolute value and the diurnal variability in the boundary layer height.Due to the coarse resolution of EMAC grid cells (approx.120 km), the boundary layer height (h BL ) in the model is often affected by diurnal variation due to neighboring continental cells, especially close to the coast.Figure 9 shows EMAC simulations of h BL compared to ERA5 data (fifth-generation ECMWF reanalysis).ERA5 resolves the global atmosphere in hourly intervals for 30 km grids at 137 vertical levels up to 0.01 hPa, and thus its horizontal resolution is a factor of 4 higher than that of EMAC.The ERA5 dataset is available within the Copernicus Climate Change Service (https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/era5; last access: 27 February 2021).The BL simulated by EMAC is very shallow at night and increases rapidly in height after sunrise, which may reflect continental influence in the EMAC grid boxes.The only regions where EMAC shows no clear continental influence are the AS and the MS, even though the model clearly underestimated the h BL over the Arabian Sea and even more over the Mediterranean Sea.Local maxima on 27 and 30 August are prominent while passing Crete and the Strait of Messina.
This misrepresentation of the MBL height and its diurnal variation by EMAC has two consequences.First, according to Eq. ( 12), the deposition loss k dep for a given deposition velocity V dep is inversely proportional to the boundary layer height h BL .Overestimations of h BL by EMAC in particular during the day would thus lead to an underestimation of the deposition sink, while it would lead to an overestimation at night.Additionally, diurnal variations in h BL lead to entrainment of free-tropospheric air into the MBL, in particular during the early morning (Fischer et al., 2015(Fischer et al., , 2019)).While vertical profiles of HCHO mixing ratios decrease with height (Anderson et al., 2017;Klippel et al., 2011;Heikes et al., 2001), H 2 O 2 and MHP mixing ratios increase up to a local maximum above the boundary layer (Allen et al., 2022a;Nussbaumer et al., 2021a;Klippel et al., 2011).Therefore, intrusion of air masses from the lower troposphere will most https://doi.org/10.5194/acp-23-119-2023Atmos.Chem.Phys., 23, 119-142, 2023  4) is not included in this calculation.likely result in a decrease in HCHO in the MBL, while peroxide mixing ratios would likely increase as shown in Fischer et al. (2015).
To further evaluate the influence of deposition on H 2 O 2 and HCHO levels, deposition velocities were derived from nighttime observations (j NO 2 < 10 −3 s −1 ) of their loss rates following the method of Shepson et al. (1992) (Eq. 13).Here we use the exponential decays of the HCHO and H 2 O 2 mixing ratios versus time to deduce nighttime loss rates over the Arabian Sea and the Mediterranean Sea, where the EMAC simulation of h BL was most accurate (Figs. 9, 10).The slope of the linear regression yields the respective deposition rate constant (k dep ) at night, assuming negligible nighttime chemistry, i.e., nonsignificant production of H 2 O 2 and HCHO due to ozonolysis and nighttime oxidants (NO 3 , Cl) and their respective losses.We only used the data of the second leg in the AS (less data coverage on the first leg due to contamination), which represent high-humidity conditions (87.6 ± 3.1 %) and strong headwinds (10.3±2.3 m s −1 ), while winds over the Mediterranean Sea were slower (6.1 ± 2.2 m s −1 , 73.6 ± 7.0 %).In general, deduced loss rates for H 2 O 2 and HCHO show higher variability in the AS than in the MS, with the highest values for both species on 12 August and the lowest values on 9 August in the eastern part of the AS after a partly cloudy day.
A mean deposition loss k dep (H 2 O 2 ) of 1.83 ± 0.93 × 10 −5 s −1 was determined for the Arabian Sea, similar to the results over the Mediterranean Sea (1.51 ± 0.85 × 10 −5 s −1 ).We thereby determined a minimum of 0.37±0.15×10−5 s −1 and a maximum deposition loss rate of 2.87±1.15×10−5 s −1 for H 2 O 2 with a variability within an order of magnitude.Dry deposition rates of HCHO are comparable to H 2 O 2 with k dep (HCHO) = 1.34±0.46×10 −5 s −1 over the AS and 1.91 ± 0.93 × 10 −5 s −1 over the MS.The deposition losses of HCHO cover a similar range of 0.33 ± 0.13 × 10 −5 s −1 to 2.74 ± 1.10 × 10 −5 s −1 .The results for the MS cover a smaller number of nights, as we experienced local increases in HCHO from 24 to 26 August 2017.These enhancements were associated with slightly elevated NO 2 , indicating local pollution events, and thus we excluded these nights.
Based on the deposition losses we calculated deposition velocities according to Eq. ( 12) using mean values of h BL from ERA5 data for the corresponding time frame.For the resulting V dep we assume an uncertainty of at least 40 % (Fig. 11, Table S2).V dep (HCHO) values cover a range of 0.23-2.22cm s −1 , with the highest values during the night of 27-28 August 2017 in the Mediterranean Sea.Mean values (±1σ ) of 0.77 ± 0.29 and 1.49 ± 0.76 were determined for the AS and MS, respectively.The deposition velocities of H 2 O 2 cover a similar range of 0.26-2.34cm s −1 , also with the highest values during the night of 27-28 August 2017.This resulted in mean values of 1.03 ± 0.52 cm s −1 for the AS and 1.21 ± 0.69 cm s −1 for the MS.
Averaged values were compared to the V dep used by EMAC in Fig. 11.In general, observation-based V dep and model values for both species are of similar magnitude for the AS (with the exception of the very low values derived during the night of 9-10 August 17), while values of V dep are underestimated by at least a factor of 2 for the Mediterranean Sea.Additionally, EMAC simulates reduced variability compared to the observations with mean values (±1σ ) of 0.78 ± 0.16 cm s −1 for the AS and 0.32 ± 0.16 cm s −1 for the MS, respectively.Deposition velocities of H 2 O 2 show enhanced values compared to HCHO due to the larger Henry's law coefficient of H 2 O 2 , which resulted in mean values of 1.05 ± 0.27 cm s −1 for the AS and 0.37 ± 0.21 cm s −1 for the MS.
Striking similarities were found for both species, as the calculated k dep and also the derived V dep seem to follow the same trend in the Arabian Sea.A linear correlation coefficient of R 2 = 0.77 was found for a linear fit of V dep (HCHO) against V dep (H 2 O 2 ) for the Arabian Sea.The simulated deposition velocity in EMAC depends linearly on the wind speed (Fischer et al., 2015), which explains the higher values derived during the period with strong headwinds over the Arabian Sea.The observations confirm larger deposition velocities of H 2 O 2 for the AS, while we determined enhanced values of V dep (HCHO) for the MS.Both species show lower deposition velocities close to the coast on 30 July 2017 and in general a higher variability than the simulations by EMAC.
The observed nighttime values of V dep (H 2 O 2 ) match previously derived values within the literature: Allen et al. (2022b) found similar values of V dep (H 2 O 2 ) in the MBL with a range of 1.00-1.32cm s −1 , which corresponds to a loss of 5 %-10 % HO x in the marine boundary layer.Stickler et al. (2007) determined a mean V dep (H 2 O 2 ) of 1.3 cm s −1 in the MBL with a range of 0.1 to 1.8 cm s −1 depending on the entrainment rate.Fischer et al. (2019) calculated nighttime deposition velocities in the continental boundary layer for five different campaigns in Europe and determined values in the range of 0.16 to 0.60 cm s −1 at night and 0.56 to 6.04 cm s −1 during the day.Nguyen et al. (2015) performed flux measurements and derived a diel cycle of V dep (H 2 O 2 ) with values ≤ 1 cm s −1 at night and a maximum of ∼ 6 cm s −1 at noon.In sum, our observations agree with previous measurements in the continental and marine boundary layer at night.
Deposition velocities of HCHO are generally expected to be lower than that of H 2 O 2 due to its higher uptake resistance, which is related to, e.g., solubility and hydrolysis of the trace gas (Stickler et al., 2007;Ganzeveld and Lelieveld, 1995).This can be confirmed with our results in the AS, where V dep (HCHO) is on average a factor of 0.8 smaller than V dep (H 2 O 2 ).The results for the MS do not confirm this expectation, as we determined a factor of 1.2 higher deposition velocities of HCHO.To the best of our knowledge, reports of the deposition velocity of HCHO in the MBL are sparse.Nussbaumer et al. (2021b) derived a nighttime V dep (HCHO) of 0.47 cm s −1 in the continental boundary layer.Sumner et al. (2001) calculated a V dep (HCHO) of 0.65 ± 0.36 cm s −1 at night at a mixed deciduous-coniferous forest site, while Stickler et al. (2007) suggest a constant value of 0.36 cm s −1 over the ocean based on single-column model results.In comparison to the evaluated V dep in the continental boundary layer, our measurements indicate more efficient deposition over the open ocean.This may be due to a more efficient near-surface transport, e.g., due to high wind speeds and turbulence, as it is expected that high wind speeds lead to a more efficient deposition.Altogether, EMAC simulates relatively accurate dry deposition velocities of HCHO and H 2 O 2 for the AS with a deviation of less than 5 % for the derived mean values, while V dep was underestimated by at least a factor of 2 for the MS.Additionally, the observations demonstrate a higher variability than the model.Please note that deposition velocities for the observations were only calculated at night, and thus deposition losses of H 2 O 2 and HCHO during daytime remain uncertain, although previous observations indicate a stronger deposition loss during the day (e.g., Fischer   et al., 2019).Please note that our observations might be affected by additional loss processes, e.g., due to interactions with sea spray or in general heterogeneous chemistry.

Conclusions
In situ measurements of HCHO, H 2 O 2 , ROOH, OH, HO 2 and actinic flux were carried out in the MBL around the Arabian Peninsula during the AQABA campaign in summer 2017 (1 July-1 September 2017).Mixing ratios of HCHO varied greatly within a range of 0.05 to 12.63 ppbv, with the highest mixing ratios over the polluted Arabian Gulf.Elevated ROOH mixing ratios (up to 2.26 ppbv) also reflected the high level of photochemical air pollution in the AG, while mixing ratios in other regions were below 0.30 ppbv.H 2 O 2 was the least abundant HO x reservoir with mixing ratios mostly within the range of 0.10-0.40ppbv, and no significant enhancements were detected during the pollution event in the center of the AG.The highest daytime median values of HO 2 were detected for the area around the Suez Canal (15.90 pptv), while significantly lower values were found for the AG (4.65 pptv).In the AG, we also encountered the lowest daytime median values of OH (0.04 pptv), while we determined a daytime median of 0.09 pptv for the whole dataset.
In general, the elevated values of HCHO and ROOH in the Arabian Gulf are consistent with the increased mixing ratios of O 3 (Tadic et al., 2020), VOCs (Bourtsoukidis et al., 2019) and OVOCs (Wang et al., 2020) (Tadic et al., 2020), although we rarely encountered really clean MBL conditions during AQABA (Celik et al., 2020).
The comparison of measurements with results from the general circulation model EMAC indicates that the model simulates HCHO mixing ratios on average within a factor of 2, while it clearly overestimates H 2 O 2 and to a lesser extent ROOH.This effect can to a large degree be explained by the overestimated mixing ratios of HO x .Production and loss of HCHO scale linearly with OH, while the photochemical production of H 2 O 2 depends quadratically on HO 2 and only linearly on OH.The calculation of net photochemical production of H 2 O 2 (P (H 2 O 2 )-L(H 2 O 2 )) revealed an overestimation by the model, with the largest deviations over the Arabian Sea, the Gulf of Oman and the Arabian Gulf.These regions were characterized by reduced mixing ratios of HO 2 , which are not reproduced by EMAC.Even though net photochemical production matched over the MS, simulated H 2 O 2 was overestimated, which may be related to transport or underestimated deposition during the day.
Loss of HCHO and H 2 O 2 due to dry deposition, their predominant loss process at night, was studied for the regions far away from the coastline so that the diurnal variation in the boundary layer height in EMAC interfered less with the derived results.The analysis revealed mean values (±1σ ) of V dep (HCHO) with 0.77 ± 0.29 cm s −1 for the Arabian Sea and 1.49 ± 0.76 cm s −1 for the Mediterranean Sea.Enhanced values of V dep (H 2 O 2 ) compared to V dep (HCHO) were determined for the AS (1.03 ± 0.52 cm s −1 ), while the deposition of H 2 O 2 to the MS was on average slower than that of HCHO with 1.21±0.69cm s −1 .EMAC generally simulates accurate deposition velocities of HCHO and H 2 O 2 over the Arabian Sea, while the deposition of both species was underestimated by at least a factor of 2 over the Mediterranean Sea.Further, the simulated V dep shows less variability compared to the observations.Losses due to deposition could not be compared during the day and remain uncertain.
The analysis of the deposition velocity revealed underestimated losses due to dry deposition in the MBL, which leads to overestimated HO x precursors in EMAC.An updated scheme for the simulation of dry deposition of landatmosphere exchange in EMAC was developed by Emmerichs et al. (2021).They include, e.g., improved stomatal uptake and dry deposition dependent on soil moisture, which decreased HCHO mixing ratios during boreal summer by up to 25 % at the ground level.We emphasize that the ocean surface might be an additional surface where dry deposition of trace gases might need to be adjusted.This effect may be most important during the day and on a local scale but certainly could have an effect on the HO x budget in EMAC.Additionally, EMAC is limited by its coarse spatial resolution (∼ 110 km), which leads to spurious diurnal variation in the boundary layer height when cruising close to the coast.
The overestimated HO x in EMAC most likely results from the overall enhanced VOC oxidation and regionally overestimated O 3 .Moreover, the model was shown to be unable to reproduce the local phenomena encountered, e.g., the air pollution event in the AG.This leads to the assumption that the model may reproduce HCHO within a factor of 2 but possibly due to the wrong reasons, i.e., overestimated HO x which compensates for missing sources in the model.This assumption is supported by the observations of several other OVOCs which were not matched by EMAC (Wang et al., 2020).
The systematic overestimation of H 2 O 2 is at least partly explained by the overrated HO 2 , although net photochemical production of H 2 O 2 revealed that the model matches the observations in some regions well, as the overall overestimated OH compensates partially for too high values of HO 2 .The decreased concentrations of HO 2 observed over the GA, AS, GO and AG are not matched by EMAC.The reduced HO x encountered in air masses over the Arabian Gulf can most likely be attributed to elevated mixing ratios of RO 2 , which are reflected in the enhanced ROOH and the estimates of RO 2 + HO 2 by Tadic et al. (2020).Despite matching results for the net photochemical production in some regions, the observations of H 2 O 2 reflect less diurnal variation and overall lower mixing ratios than simulated by EMAC.This implies that further loss processes, e.g., the deposition during daytime, and in general heterogeneous chemistry remain a major uncertainty in the photochemical budget of H 2 O 2 in the MBL.We therefore emphasize the importance of H 2 O 2 and organic peroxide in situ measurements, which were valuable to evaluate simulated deposition velocities and the accuracy of HO x simulations.the Max Planck Digital Library (https://keeper.mpdl.mpg.de/library/708bfb36-1cf8-4d6a-b5f8-acec03ed5e5e/2017_AQABA/; Lelieveld et al., 2022) and are available to all scientists, who agree to the AQABA data protocol.
Author contributions.DD, BB, HF and JL designed and supervised the study and the AQABA campaign.DD performed the HCHO and hydroperoxide measurements during the second leg of the campaign and evaluated the HCHO dataset.BB performed the HCHO and hydroperoxide measurements during the first leg and provided the hydroperoxide dataset.HH, MM, ST and RR performed the LIF OH and HO 2 measurements during AQABA.JS and JNC performed the actinic flux measurements and calculated photolysis rates.PGE and JNC carried out the O 3 measurements during the cruise.AP performed the EMAC model runs.DW calculated the HYSPLIT trajectories.
Competing interests.The contact author has declared that none of the authors has any competing interests.
Disclaimer.Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Acknowledgements.We graciously acknowledge cooperation with the Cyprus Institute (CyI), the King Abdullah University of Science and Technology (KAUST) and the Kuwait Institute for Scientific Research (KISR).We thank Hays Ships Ltd, Captain Pavel Kirzner and the crew of Kommandor Iona for the great support during all weather or wavy conditions and for an unforgettable time on board.We would like to especially thank Marcel Dorf and Claus Koeppel for the organization of the campaign and Hartwig Harder for the management on board.Last but not least we are grateful to the whole AQABA community and for a successful campaign.
Financial support.The article processing charges for this openaccess publication were covered by the Max Planck Society.
Review statement.This paper was edited by Tanja Schuck and reviewed by two anonymous referees.

Figure 1 .
Figure 1.The ship track (black) of the Kommandor Iona during the AQABA cruise subdivided according to different chemical regimes into eight regions: Mediterranean Sea (MS: yellow), Suez Canal (SU: pink), Red Sea north (RN: orange), Red Sea south (RS: green), Gulf of Aden (GA: blue), Arabian Sea (AS: turquoise), Gulf of Oman (GO: brown) and Arabian Gulf (AG: red).HYSPLIT trajectories for the respective regions are presented in the Supplement (Figs.S17, S18).The map was created with data provided by the Natural Earth website (http://www.naturalearthdata.com, last access: 3 June 2021).

Figure 2 .
Figure 2. Overview and data coverage of HCHO, H 2 O 2 and organic hydroperoxide measurements during both legs of the AQABA ship campaign (graphs on the left represent the first leg).Contaminated HCHO data (e.g., by ship exhausts) were removed from the dataset with a stack filter (based on the NO, CO and SO 2 observations); therefore there is less HCHO data coverage during the first leg in the Arabian Sea.

3. 2
EMAC model comparisonThe observations of HCHO, H 2 O 2 and ROOH were compared to numerical results of the model EMAC.The highly complex photochemistry around the Arabian Peninsula is well suited to evaluate the MOM photochemistry mechanism.The high pollution levels, e.g., in the AG contrast with the rather clean regions like the AS and RS, which represent mostly aged air masses with less anthropogenic influence, although clean MBL conditions (NO x < 50 pptv) were rarely encountered during AQABA.Here we use simulations from the lowest vertical level of EMAC (∼ 30 m) at a temporal resolution of 10 min.Time series (Figs.S1, S3 and S5) and scatterplots (Figs.S2, S4 and S6) are shown in the Supplement.Regional variations are highlighted in Figs.3 and 4by box-and-whisker plots and by the measurement-model ratio (Figs.3 and 4); the EMAC data were adapted to the measurements with a time resolution of 10 min for these plots.Numerical values are listed in TableS1.Additionally, frequency distribution of mixing ratios for observations and model simulations for the individual regions are presented in the Supplement (Figs.S14, S15, S16).

Figure 3 .
Figure 3. Formaldehyde observations (green) and EMAC simulations (yellow) divided into the eight regions during AQABA: Mediterranean Sea (MS), Suez Canal (SU), Red Sea north (RN), Red Sea south (RS), Gulf of Aden (GA), Arabian Sea (AS), Gulf of Oman (GO) and Arabian Gulf (AG).We used the lowermost model results (∼ 30 m) for the comparison in 10 min averages (a).The box represents 25 % to 75 % of the data, and the whiskers are from 10 % to 90 % with the median as the black line.Panel (b) shows the ratio between the observations and the model simulations.

Figure 4 .
Figure 4. H 2 O 2 (blue) and organic hydroperoxide (ROOH, orange) observations compared to the EMAC simulations (yellow).We used the lowermost model results (∼ 30 m) for the comparison.The box represents 25 % to 75 % of the data, and the whiskers are from 10 % to 90 % with the median as the black line.The right panel shows the ratio between the observations and the simulations.

Figure 5 .
Figure 5. OH and HO 2 daytime values (j NO 2 ≥ 10 −3 s −1 ) of the observations (grey) and the EMAC simulations (yellow).The box represents 25 % to 75 % of the data, and the whiskers are from 10 % to 90 % with the median as the black line.The EMAC data were adapted to the measurements with a time resolution of 10 min so that the diurnal variations are reflected accurately.

Figure 7 .
Figure 7. Box plot of the net photochemical production of H 2 O 2 (P (H 2 O 2 )-L(H 2 O 2 )) of the observations and the EMAC model results (yellow) during the day (j NO 2 ≥ 10 −3 s −1 ).The box represents 25 % to 75 % of the data, and the whiskers are from 10 % to 90 % with the median as the black line.The corresponding timelines are presented in Fig. S11.Note that the deposition term from Eq. (4) is not included in this calculation.

Figure 8 .
Figure 8.Diurnal variation in the net photochemical production of H 2 O 2 (P (H 2 O 2 )-L(H 2 O 2 )) of the observations (grey) and the EMAC model results (orange) during the day (j NO 2 ≥ 10 −3 s −1 ).The corresponding timelines are presented in Fig. S11.Note that the deposition term from Eq. (4) is not included in this calculation.

Figure 9 .
Figure 9.Comparison of the simulated h BL in EMAC with ERA5 (fifth-generation ECMWF reanalysis) data at a resolution 4 times higher (31 km grid) than EMAC (∼ 110 km grid) and hourly averaged data.EMAC's h BL indicates continental influence by the strong diurnal variation while cruising close to the coast (SU, RN, GA, GO and AG), with the highest values determined for the ports of Jeddah (10-13 July), Djibouti (16 July) and Kuwait (1-3 August).The measured photolysis frequency j NO 2 serves as a reference of sunlight intensity during AQABA.

Figure 10 .
Figure 10.Determination of the deposition rates k dep for H 2 O 2 and HCHO in 10 min averages at night (j NO 2 < 10 −3 s −1 ).

Figure 11 .
Figure 11.Comparison of the calculated deposition velocity V dep for HCHO and H 2 O 2 in the Arabian Sea (AS) and the Mediterranean Sea (MS).

Table 1 .
Instrument characteristics of the HCHO and hydroperoxide measurements during the AQABA campaign.
Dienhart et al.:Formaldehyde and hydroperoxide distribution around the Arabian Peninsula and HCHO in clean air masses at night is dominated by deposition.With the assumption of a constant, horizontally homogenous boundary layer and a linear concentration profile within the BL, the exponential decay can be used to estimate the deposition velocity (V dep ) with the method of Shepson At night, the photochemical production and loss reactions due to OH can be neglected; therefore, the decay of H 2 O 2 https://doi.org/10.5194/acp-23-119-2023Atmos.Chem.Phys., 23, 119-142, 2023 126 D.
Pfannerstill et al. (2019)reactivity during AQABA, wherePfannerstill et al. (2019)derived the highest median loss rate of OH with 18.8 s −1 .The lowest mixing ratios of HCHO, H 2 O 2 and ROOH were found for air masses in the Arabian Sea and the southern Red Sea during the second leg of the campaign, where the lowest concentrations of NO x and O 3 were present