Interactive comment on “ A broadband cavity-enhanced spectrometer for atmospheric trace gas measurements and Rayleigh scattering cross sections in the cyan region ( 470 – 540 nm )

This paper presented a LED based cavity enhanced absorption spectroscopy at 470540 nm. Different with previous CEAS instruments. The CEAS moved the window to cyan region (470-540 nm) and measure Rayleigh scattering cross-sections of several pure gases. This part is fundamental and very meaningful to the dataset. The authors applied the instrument to measure NO2 and I2 simultaneously. The measurement of ambient NO2 by CEAS was compared with CRDs and presented the feasibility. With respect to the I2 measurement by CEAS, which has been reported in the neighbor or the same wavelength region (e.g., Vaughan et al., 2008; Bahrini et al., 2018). Although

The accuracy of retrievals depends on knowledge of relevant absorption and scattering (i.e., Rayleigh) cross sections (which both contribute to the optical extinction) and of the mirror reflectivity.Because of its high sensitivity, CEAS in turn has been utilized to measure these parameters.For example, Axson et al. (2011) and Kahan et al. (2012) recently reported the absorption cross sections of O 3 and H 2 O 2 , and Thalman and Volkamer (2010) reported scattering cross sections of N 2 , O 2 , Ar, and air for several wavelength intervals (345-390, 435-490, 515-545, 560-600, and 600-700 nm).
In spite of the large number of CEAS instruments that have been constructed, the mid-visible region has received relatively little attention to date, other than the pioneering study by Vaughan et al. (2008), even though the absorption maxima of several key atmospheric traces gases, such as NO 2 and the iodine species I 2 , IO, and OIO, are located in this wavelength interval (Fig. 1).The study of iodine chemistry has been of considerable interest due to the potential effects on the formation of new particles and atmospheric oxidizing capacity through, for example, catalytic destruction of O 3 , altering the partitioning of NO x (=NO + NO 2 ) and HO x (=HO+HO 2 ), or the activation of chlorine and bromine from sea salt aerosol in the marine and polar boundary layer, near salt lakes and volcanoes, and in the stratosphere (Saiz-Lopez et al., 2012).Prior laboratory measurements by Vaughan et al. (2008 used a 150 W Xenon arc lamp; such light sources tend to flicker, i.e., exhibit intensity fluctuations, which add noise to the absorption spectrum.For this and other reasons (such as compactness and heat generation), light-emitting diodes (LEDs) are now commonly used to generate broadband radiation (Washenfelder et al., 2008;Min et al., 2016).
In this paper, we describe an LED-powered incoherent broadband cavity-enhanced absorption spectroscopy (IB-BCEAS) instrument operated in the cyan region of the electromagnetic spectrum.We report laboratory measurements of scattering cross sections for N 2 , O 2 , Ar, CO 2 , CH 4 , and air from 480 to 535 nm and demonstrate detection of NO 2 and I 2 in laboratory-generated air.The IBBCEAS instrument was operated during the Ozone-depleting Reactions in a Coastal Atmosphere (ORCA) field campaign, conducted on 8-31 July 2015 at the Amphitrite Point Observatory (APO) in Ucluelet on the west coast of Vancouver Island, British Columbia, Canada (Tokarek et al., 2017).Kelp forests are present along this coastline (Watson and Estes, 2011), which are expected to emit I 2 (Dixneuf et al., 2009;Nitschke et al., 2011Nitschke et al., , 2015)).Mixing ratios of NO 2 retrieved from the cyan IBBCEAS data are compared to those measured by a co-located blue diode laser cavity ring-down spectrometer.(Coheur et al., 2002), O 3 (Burkholder and Talukdar, 1994), NO 2 (Voigt et al., 2002), I 2 (Spietz et al., 2006), IO and OIO (Spietz et al., 2005), and OBrO (Knight et al., 2000).The literature cross sections for H 2 O (shown in blue) were convolved with the resolution of the IBBCEAS (shown in red).The absorption cross section of O 4 is shown in Fig. S4 in the Supplement.
The potential of the instrument for quantification of iodine species in laboratory and field experiments is assessed.

Cavity-enhanced spectroscopy
The principle of broadband trace gas measurements by IB-BCEAS has been described elsewhere (Fiedler et al., 2003;Washenfelder et al., 2008).Briefly, broadband radiation is continuously injected and trapped between a set of highly reflective mirrors, forming a stable resonant cavity.The integrated cavity output intensity represents the combined extinction by the mirrors and the intra-cavity medium.The absorption coefficient, α abs (λ), is given in terms of the transmission signal through the cavity (Washenfelder et al., 2008) by Eq. ( 1): (1) Here, R L (see Sect. 3.3) is the ratio of the cell length (d ≈ 102 cm) divided by the length occupied by the sample (d 0 ≈ 82 cm; Sect.3.3), R(λ) is the average wavelength-dependent mirror reflectivity, d is the distance between the two reflective surfaces (i.e. the cavity length), α Ray (λ) is the sum of all Rayleigh scattering sample constituents, I 0 (λ) is the intensity spectrum in the absence of absorbers in the cavity cell, and I (λ) is the intensity spectrum measured in the presence of absorbers.When measurements are carried out in variable pressure environments (e.g., aircraft altitude change), an additional term ( α Ray (λ)) is added to Eq. ( 1) to account for pressure fluctuations (Min et al., 2016).If multiple absorbers are present in the cavity, the absorption coefficient becomes the sum of all species expressed as in Eq. ( 2) (Washenfelder et al., 2008).
Here, σ i (λ, T , p) is the wavelength-, temperature-, and pressure-dependent absorption cross section, and N i is the number density of the ith absorbing gas.

Determination of Rayleigh scattering cross sections
Cavity-enhanced spectroscopic techniques allow accurate and precise measurements of Rayleigh scattering cross sections of pure gases over broad wavelength regions (Naus and Ubachs, 2000;Sneep and Ubachs, 2005;Axson et al., 2011;Kahan et al., 2012;Thalman et al., 2014).If the mirror reflectivity R(λ) is known, the scattering cross section of any gas can be determined by measuring its extinction and that of a reference gas (Thalman et al., 2014), e.g., for CO 2 and using He as a reference, Here, α Ray (λ) is the optical extinction (in cm −1 ) caused by the intra-cavity scattering medium.The scattering cross sections are then calculated through division by the number density of the gas (N gas ); i.e., Scattering cross sections can be predicted if the refractive index, n, of a gas is known (Naus and Ubachs, 2000;Sneep and Ubachs, 2005), hereafter referred to as an "n-based" cross section, σ n : Here, ν is the frequency in wavenumbers (cm −1 ), n(ν) is the frequency-dependent refractive index of the gas, and F k is the King depolarization ratio, which describes the effect of molecular anisotropy (King, 1923).

Description of the IBBCEAS instrument
The IBBCEAS instrument consists of an LED light source, collimating optics, a high finesse optical cavity, focusing optics, specialized fibre bundle, and a spectrograph equipped with a charge-coupled device (CCD) camera (Fig. 2).
The optical cavity consists of two dielectric-coated, planoconcave fused silica substrate mirrors (Advanced Thin Films, Boulder, CO, USA) with maximum reflectivity from 460 to 550 nm, 2.54 cm in diameter, 0.635 cm thickness, and 1 m radius of curvature.To span a broad wavelength range, the mirrors were coated with two highly reflective substrates, resulting in a double maximum in their reflectivity.The mirrors were mounted on both ends of a 102 cm long cell with a gas sampling region of ∼ 73 cm and housed in a custombuilt mount equipped with purge gas ports.A flow of 50 sccm (standard cubic centimetres per minute) of ultrapure air (zero grade, Praxair) set using 50 µm critical orifices (Lenox Laser, Glen Arm, MD, USA) and a gas regulator back pressure of 138 kPa was directed to each mirror to protect the optical surfaces from contamination.The mirror mounts were attached to adjustable kinematic mounts (Newport U200-A, Irvine, CA, USA), each equipped with three set screws for mirror alignment.The gas sampling region was enclosed using 1.9 cm outer diameter (o.d.) and 1.59 cm inner diameter (i.d.) fluorinated ethylene propylene (FEP) Teflon ™ tubing (Saint Gobain Plastics, Chemfluor series).The Teflon ™ tube was held in place by a custom enclosure constructed from aluminium.Gases entered and exited the sampling region through 1.9 cm o.d.perfluoroalkoxy alkane (PFA) Teflon ™ fittings (Entegris Fluid Handling, Billerica, MA, USA).
The LED (Thorlabs M505L3, Newton, NJ, USA; 400 mW nominal output power) output was collimated into the cavity by a single f/0.8 aspheric condenser lens (Thorlabs ACL2520U-A).The cavity output was collected and focused by a 2.5 cm diameter f/1 lens through a non-polarizing quartz beam splitter (Thorlabs BS025) onto a 0.5 cm (diameter) f/4 fibre adapter that couples the radiation into a 2 m long 0.22 numerical aperture (NA) fibre bundle (Thorlabs BFL200HS02).The beam splitter has 90 % transmission efficiency in the plane parallel to the cavity and 10 % in the orthogonal plane and was used to direct the emission lines of a Hg(Ne) spectral calibration lamp (Newport Oriel 6032) towards the detector for calibration of the spectrometer wavelength scale and line width.The fibre bundle consists of seven 200 µm optical fibres arranged with circular and linear configurations on the input and output, respectively.The fibre bundle output was oriented linearly along the spectrograph entrance slit to optimize coupling of the cavity output, similar to Min et al. (2016) and Washenfelder et al. (2016), and illuminated the full vertical dimension of the CCD.Spectra were acquired by a 150 cm focal length f/4 dual grating Czerny-Turner spectrograph (Princeton Instruments Acton SP2156, Trenton, NJ, USA), which back-illuminates a 16 bit 1340 pixel × 100 pixel CCD array (Princeton Instruments, PIXIS 100B), mounted at the focal plane of the spectrograph and Peltier-cooled to −80 • C to reduce thermal noise.The spectrograph was configured with a 1200 groove mm −1 grating blazed at 500 nm.This configuration along with the grating positioned at 500 nm central wavelength yielded spectral coverage from 446.9 to 563.2 nm, spanning a total of 1340 discrete wavelength points.The Acton SP2156 spectrograph shipped with a mechanical shutter, which was removed in lieu of software control of the integration time via Labview ™ (National Instruments, Austin, TX, USA).
A fixed entrance slit width of 50 µm yielded approximately Gaussian profiles of the atomic Hg emission lines with a full-width at half maximum (FWHM) of (0.330 ± 0.003) nm at 540.06 nm.Other spectral lines (e.g., 520.39 nm) showed slightly lower resolution of roughly (0.360 ± 0.016) nm (Fig. S1).An integration time of 5-6 s saturated the CCD pixels 80 %-90 % of their well depth near 500 nm.
All mechanical and optical components were mounted on an optical breadboard (Thorlabs MB2448) bolted onto aluminium structural rails (80/20, Columbia City, IN, USA).The net weight of the instrument including the breadboard and railing is < 90 kg, and the overall power consumption is below 300 W.
The IBBCEAS instrument samples air through a 47 mm diameter, circular, permeable polytetrafluoroethylene (PTFE) Teflon  NJ, USA).The sample flow was in the range 1.5-2.5 slpm, resulting in a sample residence time of 5.9-3.5 s.The gas temperature was measured using a K-type thermocouple (Omega, Laval, QC, Canada) attached to the aluminium sample cell enclosure.Pressure was measured using a pressure transducer (MKS Baratron capacitance manometer 722B) mounted on the exhaust portion of the IBBCEAS between the gas sample cell and the MFC.
The light source is a small footprint (1 mm × 1 mm) broadband LED (M505L3, Thorlabs) mounted on the end of a heat sink supplied by the manufacturer.We modified this assembly by mounting a Peltier module (CUI Inc.CP30238, Tualatin, Oregon, USA) between the LED and the heat sink for temperature control using a proportional-integral-derivative (PID) controller (Omega CNi3253) to 30.0 ± 0.1 • C. If not stabilized, the LED output red-shifted ∼ 0.1 nm • C −1 temperature change.The LED was operated just below the maximum current of 1000 mA (3.3 W of electrical power).
The LED output spectrum was characterized by an asymmetric shape, an emission maximum at 507.5 nm, and a peak width of 22.5 nm FWHM, and was a good match with the mirror reflectivity curve (Fig. 3).

Determination of the mirror reflectivity
The response of the IBBCEAS instrument depends on R(λ) and d 0 (Eq.1), which need to be accurately known.If an optical cavity is filled with a non-absorbing and inert gas, and aerosols are removed by filtering, the optical extinction is due to gas scattering only.In this case, R(λ) can be derived from the extinction caused by two individual gases with known and different scattering cross sections such as N 2 and He (Washenfelder et al., 2008): . (6) Here, I X (λ) and α X Ray (λ) are the intensities and extinction coefficients of N 2 and He, respectively.Other pairs (e.g., He and Ar) may be used as well.For this work, we chose scattering cross sections from Peck and Khanna (1966) for N 2 , Cuthbertson and Cuthbertson (1932) for He, and Peck and Fisher (1964) for Ar.Typically, combinations with He provide a higher signal-to-noise due to He having a much smaller scattering cross section than other gases (Thalman et al., 2014).The above approach assumes that scattering and absorption within the optical cavity are small (Washenfelder et al., 2008) and neglects absorption and scattering due to the mirror substrate.
To determine R(λ) in this work, the cavity output intensity was recorded when the cell was filled with high-purity N 2 (99.998 %) or with He (99.998 %) to atmospheric pressure (∼ 890 hPa).In each case, the gases were delivered from a compressed gas cylinder (Praxair) and continuously injected through the purge ports until all other sample cell constituents were displaced through the open inlet port.To corroborate the results, the cavity was also filled with Ar (99.998 %), which has slightly larger scattering cross sections than N 2 .
Figure 3a shows cavity-enhanced transmission spectra with the optical cavity filled with N 2 , He, and Ar, each averaged over a 15 min time period.The intensity difference due to scattering by Ar and N 2 was more pronounced in the case when the sample cell was filled with He (Fig. 3a, inset), as expected based on their relative cross sections (Thalman et al., 2014).
Figure 3b shows the reflectivity spectra calculated using Eq. ( 6).The reflectivity curves are essentially smooth continua, with a maximum reflectivity of ∼ 99.98 %.The relative difference between the reflectivity calculated from the ratios of Ar/He and N 2 /He was <1.6 × 10 −6 in the 480-530 nm range.Arbitrarily, the N 2 /He mirror curve was chosen to calculate the mirror reflectivity for the remaining data analysis in this work with the exception of the scattering cross sections of N 2 and CO 2 (see Sect. 4.1.1 and 4.1.5)for which the Ar/He pair was used.The observed reflectivity corresponds to a total path length (≈ d/(1 − R)) of 5.5 and 6.3 km near 485 and 520 nm, respectively.Intermittent mea-N.Jordan et al.: A broadband cavity-enhanced spectrometer for atmospheric trace gas measurements surements of mirror reflectivity showed that it remained unchanged throughout this work; in particular, no difference in reflectivity was found prior to and following the ORCA campaign.
The overall uncertainty in the measured mirror reflectivity was ±2.3 %, resulting from addition by quadrature of a ±1 % uncertainty in the N 2 scattering cross section (Peck and Khanna, 1966), a ±2 % uncertainty in the He cross section (Washenfelder et al., 2008;Min et al., 2016), ±0.4 % uncertainty in temperature, ±0.1 % uncertainty in pressure, and ±0.1 % random noise due to photon counting statistics.

Determination of the effective optical absorption path
When sampling ambient air, knowledge of the path length over which sampled air is present (d 0 ) and the associated correction factor R L is essential for accurate IBBCEAS measurements when purge gases are used to maintain clean mirrors as in this work.The effective cavity length, d 0 , was determined by sampling pure oxygen and monitoring the O 2 dimer (Duan et al., 2018).This experiment was performed in the 350-390 nm wavelength region using highly reflective mirrors and an LED (Thorlabs M365LP1) suitable for this wavelength region, though in principle it could have been performed in the cyan region as well since both spectral regions contain O 4 absorption lines.Following Duan et al. (2018), d 0 is given by where [O 2 ] on and [O 2 ] off are the [O 2 ] retrieved using cross sections by Thalman and Volkamer (2013) with or without the purge flows.Consistent results were obtained when a constant mixing ratio of NO 2 was sampled (and purge flows were turned on and off) and from the ratio of expected absorption by O 4 (calculated from the square of the number density of O 2 and the O 4 cross section) relative to the observed O 4 absorption, though the latter may be less accurate as the absorption changes non-linearly in the region where purge and sample gases mix and requires accurate knowledge of mirror reflectivity.At a sample flow rate of 2.0 slpm, a purge flow rate of 0.1 slpm, and cell pressure of 870 hPa, we determined R L = d/d 0 = 1.28 ± 0.05 (1σ precision stated), with an accuracy of ±5 %.This R L value is smaller than calculated from the ratio of the inter-mirror distance d to the distance between the inlet and outlet of the sample flow (∼ 1.39), which suggests that sampled air partially mixes into the purge volume.At a volumetric flow rate of 11 L min −1 as deployed during ORCA, we observed R L = 1.20 ± 0.07 (1σ precision).

Reference absorption cross sections
Figure 1 shows trace absorbers in the 480-530 nm window calculated as extinction (α i (λ) = σ i (λ)N ) for typical atmo-spheric mixing ratios.IBBCEAS analysis was performed using the high-resolution NO 2 absorption cross section of Voigt et al. (2002) and the cross sections of Spietz et al. (2005) for OIO and Spietz et al. (2006) for I 2 , respectively.The NO 2 cross section was convolved with the instrument function corresponding to a sharp peak line (540.06nm) in the emission spectrum of the Hg(Ne) calibration lamp (Fig. S2), which was also used to calibrate the wavelength scale of the spectrometer.
The IBBCEAS instrument's resolution was 0.33 nm in the 480-530 nm window, whereas the NO 2 literature spectrum was recorded at a higher resolution of better than 0.005 nm (Voigt et al., 2002).A convolved NO 2 reference spectrum was used for the least-squares analysis.The literature cross sections of I 2 and OIO (Spietz et al., 2005(Spietz et al., , 2006) at a resolution of 0.59 nm were not convolved because their resolution was lower than that of the IBBCEAS instrument.
The NO 2 absorption cross section is somewhat pressuredependent, which is important considering that the IB-BCEAS instrument was operated at ∼ 466 hPa in the field.Voigt et al. (2002) reported cross sections at 1333 and 133 hPa (in N 2 ).Following convolution to the IBBCEAS instrument, these differences were judged to be small (2 %-5 % in the 480-530 nm range).

Spectral fitting
The DOAS intelligent system (DOASIS; Kraus, 2003) was used to retrieve gas-phase concentrations.Optical extinction data collected in laboratory experiments were fitted using the literature absorption spectra (Sect.3.3) and a thirddegree polynomial to represent broadband processes (e.g., temperature fluctuation, Mie scattering, and extinction by aerosols).The "spectral shifting" setting in DOASIS was set to ±0.2 nm, and stretching was disallowed.
The same parameters were used to fit the data from the ORCA field campaign except that in an effort to minimize the fit residuals, the convolved absorption cross section of H 2 O based on the high-resolution data by Coheur et al. (2002) was also included in the fit.The contribution of water vapour to the optical extinction was small since the relative humidity of the air sampled during the background measurement was the same as during the measurement of NO 2 or I 2 .Further, the optical extinction (α) values, as per Eq.(1), were smoothed using a fourth-order polynomial filter (Savitzky and Golay, 1964).Smoothing improved the root mean square of the residuals and the fit uncertainty by a factor of ∼ 2 (Fig. S3).

Description of the CRDS instrument
The IBBCEAS instrument quantified NO 2 in parallel to a blue diode laser CRDS instrument, which has been described previously (Paul and Osthoff, 2010;Odame-Ankrah, 2015).Briefly, mixing ratios of NO 2 are quantified by optical absorption using a continuous wave (cw) laser whose emission is centred at 405 nm (Power Technology IQµ2A105, Little Rock, AR, USA).The radiation is square-wave (on/off) modulated at a repetition rate of 1 kHz (50 % duty cycle) and enters a 111 cm long stable resonator formed by two highly reflective mirrors (2.54 cm diameter; Advanced Thin Films, Boulder, CO, USA) in a concentric configuration.The cavity output is coupled into a 200 µm diameter, 0.22 numerical aperture, multimode optical fibre (Thorlabs M25L01) illuminating a photomultiplier tube (Hamamatsu Photonics H9433-03MOD, Japan).The ring-down decay traces are co-added and fitted to an exponential using the discrete sums algorithm described by Everest and Atkinson (2008).
The CRDS NO 2 channel was operated in parallel to a CRDS NO x channel, to which a small flow of O 3 in O 2 was added to convert NO to NO 2 , similar to the method described by Fuchs et al. (2009).Following the O 3 addition point, the sampled air passed through a coiled reaction chamber with a total residence time of ∼ 7 s to ensure complete titration and a linear response up to ∼ 100 ppbv of NO (Odame-Ankrah, 2015).An identical chamber was added to the NO 2 channel to match the response times of the two channels.Mixing ratios of NO 2 and NO x measured by CRDS were accurate within ±10 % (Osthoff et al., 2018).

Generation and delivery of gases
Figure 2b shows the set-up for delivery and sampling of NO 2 .NO 2 was generated by mixing the output of a standard NO cylinder (100.2 ± 1.0 ppmv in N 2 ; Scott-Marrin, Riverside, CA, USA) delivered using a 20 sccm all-metal MFC (MKS Instruments 1479A) with O 3 produced by illuminating a flow of O 2 (Praxair) with a 254 nm Hg pen-ray lamp (Jelight, Irvine, CA, USA) and dilution with zero air.
Iodine was purchased in solid crystalline form (≥ 99.8 %, Sigma-Aldrich, Oakville, ON, Canada).Small amounts were placed in a Teflon ™ permeation tube (VICI Metronics, Poulsbo, WA, USA); these tubes were sealed on both ends with PTFE plugs compressed with stainless steel rings and were (individually) placed in a temperature-controlled permeation chamber (Vici Metronics, model 120-1) operated with a constant air flow of ∼ 0.3 L min −1 .The flow from the permeation chamber was diluted with zero air, yielding iodine concentrations in the range of ∼ 70 pptv up to 21 ppbv depending on chamber temperature and permeation tube dimensions (i.e., wall thickness and length).
Ar, CO 2 , CH 4 , N 2 , O 2 , air, and He were delivered from compressed gas cylinders (Praxair) in the same manner as described in Sect.3.2.

ORCA field campaign
The IBBCEAS instrument was operated alongside the CRDS during a 4-week-long field intensive conducted at the Amphitrite Point Observatory (APO) on the west coast of Vancouver Island, British Columbia, Canada (Tokarek et al., 2017).The APO was established as a marine background monitoring site but is sporadically impacted by NO x emissions mainly from ship and boat traffic in the region (McKendry et al., 2014).
The IBBCEAS and CRDS instruments sampled, in parallel, from a 5.5 m long common 0.63 cm o.d. and 0.48 cm i.d.FEP Teflon ™ inlet line, flooded with calibration gases and zero air at regular intervals.Zero air was generated using a custom-built zero air generation system (Odame-Ankrah, 2015), which delivered air at ambient relative humidity that was free of trace gases such as NO 2 , NO, I 2 , or O 3 .This approach minimized absorption by H 2 O from the extinction spectra which can create artefacts (Bahrini et al., 2018).A flow restriction was placed upstream of the inlet filter.When the IBBCEAS instrument was operated at a flow rate of 5 slpm, an average (±1σ ) sample cell pressure of 473 (±23) hPa and an inlet residence time of 2.3 s were achieved.The cell pressure varied over short time periods, with a standard deviation of ±2.4 hPa.A PAN gas chromatograph (Tokarek et al., 2014) and commercial NO/NO y and O 3 instruments (Thermo 42i and 49i), as well as a 662 nm diode laser N 2 O 5 channel (Osthoff et al., 2018), also sampled off this common inlet line, resulting in all instruments periodically sampling a variety of calibration gases.
The IBBCEAS transmission spectrum was integrated for 9 s, resulting in ∼ 90 % saturation of the CCD at 500 nm.Dark spectra (to characterize the spectrometer offset) were acquired daily with the same integration time (9 s).Wavelength calibration spectra of the Hg(Ne) lamp were collected daily over a period of 2 min and showed that the spectrometer wavelength scale remained relatively unchanged (±0.02 nm) over the region of interest.The mirror reflectivity was measured prior to and after the campaign and agreed within ±(1 × 10 −7 ).Data were averaged over 60 s post-campaign prior to analysis.

Measurements of cross sections of pure gases
The extinction cross sections of N 2 , CO 2 (based on mirror reflectivity determined from the Ar/He ratio of scattering coefficients), O 2 , air, Ar, and CH 4 (based on mirror reflectivity determined from the N 2 /He ratio of scattering coefficients) in the 480-535 nm wavelength range are shown in Fig. 4; values at selected wavelengths are summarized in Table 1.For N 2 , Ar, and CO 2 , the extinction is due to scattering only.In contrast, the spectra of O 2 , air, and CH 4 include optical losses due to absorption in addition to scattering.The systematic uncertainty of these data was ±2.5 %; the systematic errors arise mainly from uncertainties in the mirror reflectivity (±2.3 %, see Sect.3.2) and errors in the scattering cross section terms in Eq. ( 3).The precision of the data (calculated after subtracting a fourth-order polynomial from the observed cross sections) varies from ±2.5 % for the CO 2 data shown in Fig. 4e to ±0.9 % for the N 2 data shown in Fig. 4b.

Nitrogen
Figure 4b shows the measured scattering cross sections of N 2 at a temperature of 299.0 ± 0.1 K and pressure of 879.7 ± 0.7 hPa between 480 and 535 nm.Superimposed are n-based cross sections predicted using Eq. ( 5), data by Peck and Khanna (1966), and the King correction factor by Bates (1984).In addition, the single-wavelength CRDS measurement at 532.2 nm by Sneep and Ubachs (2005) and the nephelometer data of Shardanand and Rao (1977) are shown.
The IBBCEAS data are marginally larger than the nbased prediction, with the percent difference ranging from +0.2 % at 525.07 nm to +1.3 % at 485.03 nm (Table 1).These differences are well within the margin of measurement error (±2.5 %), set in this case by the accuracy of the scattering cross sections of He (±2 %) and Ar (< 1 %).The IBBCEAS data differ from the CRDS data of Sneep and Ubachs (2005) at 532.2 nm by +7.1 %, though the error bars overlap; the measurement uncertainty of the CRDS data alone is ±9.4 %.The IBBCEAS measurements also differ by +4.6 % at 488.0 nm and by +8.8 % at 514.5 nm from the data of Shardanand and Rao (1977) but are within the combined measurement error (±2.5 % and ±11 %, respectively).

Argon
The scattering cross sections of Ar at 299.0 ± 0.1 K and 879.7 ± 0.7 hPa are shown in Fig. 4b.Superimposed are the n-based cross sections calculated using data from Peck and Fisher (1964) and the King correction factor from Bates (1984).The single-wavelength measurement of Sneep and Ubachs (2005) and the nephelometer data of Shardanand and Rao (1977) are shown for comparison.Shardanand and Rao (1977).
The observed scattering cross sections are marginally larger than the n-based prediction, by +1.3 % at shorter wavelengths and +0.21 % at 515.06 nm (Table 1).These differences are within the measurement error of the IBBCEAS instrument (±2.5 %) inferred from the accuracy in the determination of the scattering cross sections of He (±2 %) and N 2 (±1 %) used in the determination of the mirror reflectivity.The IBBCEAS cross sections differ from the CRDS data by Sneep and Ubachs (2005) at 532.2 nm and are well within either measurement error bar (Fig. 4b).In contrast, the CEAS values differ from the nephelometer data by −10.8 % at 488.0 nm and by −5.1 % at 514.5 nm, though they are within the ±11 % uncertainty of the nephelometer.

Oxygen
The extinction cross sections of O 2 (at a temperature of 299 K and a pressure of 880 hPa) are shown in Fig. 4c next to the n-based predictions based on Bates (1984) and the nephelometer measurements of Shardanand and Rao (1977).The IBBCEAS data show two absorption bands due to the oxygen dimer, O 4 , centred at 477 and 532 nm (Thalman and Volka-mer, 2013).In the 490-515 nm wavelength interval, the contribution of these absorption bands to the total extinction is negligible (< 0.2 % at 515 nm); i.e., the extinction is dominated by scattering.In this region, the observed cross sections are slightly larger than the n-based data of Bates (1984): the differences range from +8.3 % at 495.08 nm to +2.5 % at 515.06 nm (Table 1).Relative to the nephelometer measurements, the IBBCEAS data differ by +8.7 % at 488 nm and by +15.3 % at 514.5 nm.
We calculated the absorption cross section of O 4 by subtracting the scattering cross section of Bates (1984), upshifted by +3.86 × 10 −28 cm 2 molecule −1 (i.e., by 6 %) to match the n-based data to the observed extinction from 487 to 516 nm, from the extinction shown in Fig. 4c, and dividing it by the number density of O 2 (2.13×10 19 molecules cm −3 ).The result is shown in Fig. S4.The error for this absorption cross section calculation is estimated at ±8.5 %, with the largest source of error being the correction to the scattering cross section extrapolated from the region where absorption is negligible, For the smaller band between 520 and 540 nm, the observed cross section of (1.08 ± 0.09) × 10 −46 cm 5 molecule −2 is in agreement with recent roomtemperature measurements by other groups (Supplement Table S1).For the larger band between 467 and 485 nm, the peak cross section of (6.2 ± 0.5) × 10 −46 cm 5 molecule −2 is smaller than that reported in recent measurements (Thalman and Volkamer, 2013;Sneep et al., 2006), though the shape of the peak is identical and the error bars overlap.

Air
Figure 4d shows the extinction spectrum of ultrapure air (zero grade, Praxair) at a temperature of (298.5 ± 0.1) K and pressure of 879.7±0.7 hPa.Superimposed is the n-based prediction from the refractive index data of Penndorf (1957) and the King correction factor from Bodhaine et al. (1999).
The sample cylinder contains a sufficiently high quantity of O 2 (∼ 19.5 % by volume, Praxair) for O 4 absorption bands to appear at 477.5 and 531.5 nm.In the 495-515 nm wavelength region, the cross sections observed by IBBCEAS are slightly larger than the n-based prediction, by +4.0 %, +3.0 %, and +6.8 % at 495, 505, and 515 nm, respectively (Table 1).

Carbon dioxide
Figure 4e shows the scattering cross sections of CO 2 at (299.0 ± 0.1) K and (879.7 ± 0.7 hPa for the 475-540 nm wavelength interval.Superimposed are the n-based prediction based on Bideau-Mehu et al. (1973), nephelometer measurements (Shardanand and Rao, 1977), and a CRDS measurement at 532.2 nm (Sneep and Ubachs, 2005).The observed scattering cross sections differ from the n-based prediction by −2.9 %, +2.4 %, and −2.8 % at 485, 505, and 525 nm, respectively (Table 1).In contrast, the CRDS mea-surement at 532.2 nm, which has a relatively high uncertainty of ±6.5 % at the 1σ level, is 11 % lower than the IBBCEAS measurement.The nephelometer measurements (Shardanand and Rao, 1977), on the other hand, are significantly larger than the other data shown and appear to be in error.

Methane
Figure 4f shows the extinction spectrum of CH 4 (research grade, 99.97 %) for the 475-550 nm interval at (298.5 ± 0.1) K and 879.7 ± 0.7 hPa.The CH 4 spectrum in this region contains a number of combination and overtone absorption bands near 484.7,486.1, 509.6, 522.0, and 543.4 nm that have been assigned by Giver (1978).These bands are relevant to studies simulating the atmospheres of Jovian planets (e.g., Uranus, Neptune) and Saturn's moon Titan that contain CH 4 (Karkoschka, 1994).Superimposed are CRDS data by Sneep and Ubachs (2005), the n-based scattering cross sections, calculated using the parameterization by Sneep and Ubachs, and refractive index data from Hohm (1993), assuming a depolarization ratio of unity and nephelometer data by Shardanand and Rao (1977) extrapolated using Eq.(5) (dashed line).In spite of the many absorption bands, there are small windows in which the optical extinction appears to be dominated by scattering and a comparison with literature is feasible.There is very little difference (within combined measurement errors) between IBBCEAS and interpolated nephelometer data at those wavelengths.At 492.06 and 527.28 nm, for example, the data are within −0.07 % and 0.17 %, respectively (Table 1).There is also reasonable agreement between the IB-BCEAS and CRDS data at 532.2 nm, which differs by 4.7 % but compares well with the total IBBCEAS uncertainty of 4.6 % for CH 4 .

Nitrogen dioxide
An example of a spectral retrieval for NO 2 in ambient air is shown on Fig. 5a.The figure shows the measured optical extinction (α) for a routine calibration period during the ORCA campaign on 22 July 2015 at 11:52:07 to 11:57:02 (UTC).Superimposed is a fit determined with DOASIS.An NO 2 mixing ratio of (11.6 ± 0.4) ppbv was retrieved, corresponding to an ambient air mixing ratio of (13.9 ± 1.0) ppbv when R L is factored in.The residual spectrum is lacking structure, indicating that other absorbers are not significant in this wavelength region.
During the ORCA campaign, the inlet of the IBBCEAS instrument (and of the CRDS instrument, which sampled in parallel) was overflowed every 30 min with a standard gas mixture of ∼ 20 ppbv NO x containing up to 16 ppbv of NO 2 in zero air and with ∼ 130 ppbv of NO added to ambient air.The zero air was generated using a scrubber constructed inhouse that delivered air with a similar moisture content as ambient air.A subset of these data (and the ambient air data sampled in between) is shown in Fig. 6a.
High concentrations of NO in air are prone to oxidation (by O 3 and, to a lesser degree, by O 2 ) to NO 2 (Atkinson et al., 2004).Because the CRDS had a longer inlet residence time (7 s) than the IBBCEAS instrument (2.3 s), it observed greater conversion of NO to NO 2 when the high NO concentration standard was sampled in an ambient air matrix.On the other hand, when the lower concentration NO x standard was sampled in O 3 -free air, the two instruments were in better agreement (Fig. 6a).
Shown in Fig. 7a is a scatter plot of all of the NO x calibration data.This plot shows a linear relationship with a slope of 1.11 ± 0.01, an intercept of 0.2 ± 0.1 ppbv, and r 2 = 0.975.The slope of the line in Fig. 7a reveals an 11 % systematic difference in the measurement of the calibration standard.
The median NO 2 mixing ratio during ORCA was 1.24 ppbv.Shown in Fig. 6b is a time series of a subset of ambient air NO 2 measurements.The IBBCEAS NO 2 data clearly exhibit more scatter than the CRDS NO 2 data and occasionally fall below zero.For ambient air data collected on 18-19 July 2015, the scatter plot of IBBCEAS and CRDS NO 2 data (Fig. 7b) has the same slope (1.12 ± 0.02) as that shown in Fig. 7a, while the scatter plot of the entire cam-paign has a slope of 1.02 ± 0.01 (Fig. 7c).When data are filtered by removing all IBBCEAS points below a factor of 3 times the standard deviation of a blank measurement for 60 s (∼ 1 ppbv; see Sect.4.3), the slope of the scatter plot, when forced through an intercept of zero, is 1.18 ± 0.01 (r 2 = 0.70), suggesting that the IBBCEAS data are biased high.
4.2.2Iodine species (I 2 , OIO) Figure 5b shows an example retrieval for I 2 .The extinction spectrum was recorded during the ORCA campaign on 18 July 2015, when I 2 from a permeation source had been added to the inlet and represents a 60 s average.A mixing ratio of (10.1 ± 0.2) ppbv was retrieved in this example.
Figure S5 shows spectra of laboratory-generated I 2 samples and their respective absorption spectra, fits, and residuals.Fit results near 20 pptv showed large residuals (±27 %), indicating that the IBBCEAS instrument was near its detection limit.
A sample time series of laboratory I 2 measurements (30 s averaged data) is shown in Fig. S6.Gas streams containing I 2 were generated using four permeation tubes of different wall thickness, which were exchanged while the output of the permeation was bypassed and the IBBCEAS sampled zero air.The IBBCEAS instrument responded rapidly to concentration changes: for example, at 22:05, the response 30 s after "zeroing" was 10 % of the preceding value, indicating that memory effects (from slow adsorption and desorption kinetics of iodine on the inner walls of the inlet) were negligible.
Figure S7 shows a sample time series of I 2 mixing ratios during the ORCA campaign.Concentrations of I 2 (and of OIO, not shown) in ambient air were below the instrument's detection limits.At 22:30 UTC, I 2 from a diffusion source was added to the inlet.While the transmission of iodine through inlets was not systematically investigated in this work, the square-wave response and quick rise and fall times suggest the absence of inlet transmission losses.
The limit of detection (LOD) values were assessed through Allan deviation analyses (Werle et al., 1993).For NO 2 sam-  pled at a flow of 1.5 slpm and cell pressure of 890 hPa, the Allan deviation was ±360 pptv for 10 s data, ±135 pptv for 60 s data, and ±63 pptv for 5 min data averages (Fig. 8a).At the higher sample flow of 5 slpm and reduced cell pressure of 466 hPa used during the ORCA campaign, the Allan deviation was ±164 pptv for 334 s data (Fig. S8).
For I 2 in laboratory-generated samples, typical 1σ fit errors for 60 s averaged data were ±7 pptv when 21 pptv I 2 was sampled and ±9 pptv when 544 pptv I 2 was sampled (Fig. S6).The accuracy of I 2 data is, in principle, of similar magnitude to that of the NO 2 data, except that it also depends on knowledge of inlet transmission efficiency, which was not assessed in this work.The Allan deviation plot in Fig. 8b demonstrates 1σ measurement precisions for I 2 of ±49 pptv for 60 s and ±22 pptv for 5 min averaged data.During ORCA, the 1σ precisions were ±120 pptv for 60 s and ±60 pptv for 5 min data, respectively.
For OIO, the Allan deviation analysis gives 1σ measurement precisions of ±5.7 pptv for 60 s and ±2.3 pptv for 5 min averaged data (Fig. S9b) in the laboratory.

Discussion
The IBBCEAS instrument described in this work adds to a growing number of instruments designed for the measurement of atmospheric trace gases (Table 2).It differentiates   itself foremost through its wavelength region, where several trace gases of atmospheric interest (e.g., NO 2 , I 2 , IO, and OIO) absorb (Fig. 1).Measurements of optical absorption in the cyan region enable simultaneous quantification of I 2 , and, potentially, IO and OIO in a single channel.Currently, such measurements require multiple detection channels, for example the quantification of I 2 and OIO mixing ratios using the 525-555 nm range and those of IO using the 420-460 nm window in a separate channel (Vaughan et al., 2008).
A reduction of channels and LEDs corresponds to savings in space and power requirements, important considerations in the field.As such, this paper has laid some of the ground work for future measurements in the cyan spectral region, in that we measured relevant scattering and absorption cross sections of pure gases (see Supplement).Such information is needed to determine mirror reflectivity and is useful to probe the consistency of scattering cross sections reported for this wave-length region (of which there have been relatively few) and, hopefully, improve upon their accuracy, to better describe transmission of radiation through the atmosphere.
Scattering cross section measurements for non-absorbing gases agree with recent literature values (Fig. 4 and Table 1).For example, the IBBCEAS scattering cross section measurements agree, within the combined experimental uncertainties (1σ ), with the CRDS data by Sneep and Ubachs (2005) at 532.2 nm for N 2 , Ar, CO 2 , and CH 4 , and with the nephelometer measurements of Shardanand and Rao (1977) at 488.0 and 514.5 nm for N 2 , Ar, O 2 , and CH 4 .For CO 2 , the IBBCEAS data closely match the data by Bideau-Mehu et al. (1973) and recent data by He et al. (2018); the older data by Shardanand and Rao (1977), in contrast, appear high (Fig. 4e).
In the case of CH 4 , to a lesser extent, O 2 and air, the analysis is complicated by absorption lines.For O 2 , our data and derived collisional-induced absorption cross sections agree with literature (Table S1).For regions that appear to be free of absorption bands, the IBBCEAS data for CH 4 (Fig. 4f) agree with cross section measurements of Shardanand and Rao (1977) but not with the more recent work of Hohm (1993) whose data appear high.The CRDS data point by Sneep and Ubachs (2005) is on a shoulder of a large absorption band and hence not a valid measurement of scattering but rather of the extinction cross section; their data agree with the extinction cross section observed in this work.Accurate knowledge of the scattering cross section of CH 4 is important in the study of planetary bodies with atmospheric content of CH 4 such as Saturn's moon Titan.Future studies should re-examine the scattering cross section measurements of CH 4 to resolve the reported differences.
The detection limit for NO 2 achieved under laboratory conditions in this work (63 pptv for 5 min data) is of similar magnitude as those by instruments operated in other wavelength regions (Table 2).However, the IBBCEAS measurement precision in this work was surpassed by the more mature blue diode CRDS, though future upgrades (see below) may improve the IBBCEAS precision.On the other hand, the IBBCEAS may ultimately be more accurate.Unlike IB-BCEAS, for which all absorbing molecules are incorporated in the fit, the measurement of NO 2 by blue diode laser CRDS is prone to potential interference from molecules that absorb at 405 nm, such as glyoxal and methyl glyoxal (Fuchs et al., 2009).For example, Fuchs et al. (2009) estimated that the presence of glyoxal could introduce an interference of ∼ 1 % in polluted and up to 10 % in forested environments.The low precision of the data in the scatter plot shown in Fig. 7c prevents us from drawing a definitive conclusion as to the magnitude of such interferences.Still, future comparisons of IBBCEAS and blue diode CRDS measurements of ambient NO 2 should be conducted.The data presented in this paper show that accurate measurements of NO 2 in ambient air by a cyan IBBCEAS are possible.
An important parameter in any IBBCEAS measurement is the effective cavity length.Duan et al. (2018) recently summarized past practices in its determination and found that these practices vary considerably between groups.In this work, the necessary correction (R L = 1.28 at 2.0 slpm sample flow rate) differed substantially from the physical dimension of the cell (1.39), underlining its importance for accurate measurements.However, it is unclear to what extent R L varies between molecules, which will need to be investigated.To improve the accuracy of future measurements, we will deploy more frequent zeroing and verify R L periodically, for example by periodic sampling of molecular oxygen and measuring O 4 concentration.
The sensitivity for I 2 and OIO (22 and 2.3 pptv for 5 min averaged data under laboratory conditions) suffices for their quantification in environmental chamber studies (Dixneuf et al., 2009).Further, these LOD values are below maximum I 2 and OIO abundances reported at Mace Head, Ireland, of 94 and 13 pptv (Bitter et al., 2005) and 61 and 9.2 pptv (Peters et al., 2005), respectively, but above the maximum I 2 level of 4 pptv reported in California (Finley and Saltzman, 2008).This implies that iodine species on the west coast of British Columbia, Canada, might have been detected if the instrument had been operated optimally.Even then, the LOD for I 2 of ∼ 100 pptv (2σ , 60 s average) is larger than the LOD of 26 pptv (2σ , 60 s average) reported by Vaughan et al. (2008).This suggests that improvements are possible.
One area for improvement is the thermal management of the LED.Its intensity and wavelength drift with temperature, necessitating temperature stabilization to achieve a constant emission profile.In this work, the cyan LED was stabilized by mounting a single Peltier thermoelectric module and thermocouple between the heat sink and the LED.When evaluated in the climate-controlled laboratory, this yielded a cavity emission profile that varied along the entire wavelength range (446.9-563.2nm) of the spectrometer by ±15 counts (Fig. S10).At 500 nm, this corresponded to an absorption coefficient of ±1 × 10 −9 cm −1 .During the field deployment, however, the LED was subjected to greater output fluctuations, since the trailer temperature was not as well controlled as in the laboratory, and thermal gradients between the LED and thermocouple may have come into play.
When deployed in the field, the performance of the IB-BCEAS instrument was also compromised by variable cell pressure which added noise to the optical extinction, caused in part by the large pressure drop from ∼ 1013 to ∼ 467 hPa.Furthermore, we discovered post-campaign that the alignment of the round-to-linear fibre bundle is very sensitive to vibrations, which would have added additional noise during the field campaign.In future deployments, improved pressure, temperature, and vibrational stabilization of the sample cell (as is common in aircraft deployments, for example) and instrument will be paramount.

Figure 1 .
Figure 1.Absorption coefficients for atmospheric absorbers in the cyan region at typical tropospheric mixing ratios (stated in brackets after each species).Absorption cross sections were obtained from Keller-Rudek et al. (2013) and are based on the following: H 2 O (Coheur et al., 2002), O 3(Burkholder and Talukdar, 1994), NO 2(Voigt et al., 2002), I 2(Spietz et al., 2006), IO and OIO(Spietz et al., 2005), and OBrO(Knight et al., 2000).The literature cross sections for H 2 O (shown in blue) were convolved with the resolution of the IBBCEAS (shown in red).The absorption cross section of O 4 is shown in Fig.S4in the Supplement.

Figure 2 .
Figure 2. (a) Set-up of the cyan IBBCEAS (not to scale).(b) Set-up of NO 2 calibration gas delivery for instrument characterization experiments.Abbreviations: LED -light emitting diode, BS -quartz beam splitter, CCD -charge-coupled device, ZA -zero air, USB -universal serial bus, and MFC -mass flow controller.

Figure 3 .
Figure 3. (a) Transmission spectra (averaged over 15 min) observed when the sample cell was filled with He, Ar, or N 2 at ambient pressure (890 hPa).The LED output spectrum (in arbitrarily scaled units) is superimposed.(b) Mirror reflectivity and effective path length (based on Ar/He) calculated from the data shown in panel (a).

Figure 5 .
Figure 5. Spectral fits to 60 s IBBCEAS signal, including a thirdorder polynomial and fit residual, for (a) NO 2 (ambient) and (b) I 2 (calibration) during the ORCA field campaign.

Figure 6 .
Figure 6.(a) Sample time series of continuous IBBCEAS and CRDS data collected during the ORCA campaign showing measurements of ambient air, NO 2 -free zero air (every 30 min), and hourly standard additions of NO 2 (∼ 10 ppbv) in zero air and NO (∼ 130 ppbv) to ambient air.The areas shaded in grey indicate times when both instruments sampled zero air.When high concentrations of NO were added, more NO 2 was observed by CRDS due to this instrument's longer inlet residence time.(b) Sample time series of ambient air NO 2 mixing ratios observed by IBBCEAS and CRDS during ORCA.

Figure 7 .
Figure 7. Scatter plots and straight-line fits of IBBCEAS and CRDS NO 2 mixing ratios for (a) standard additions and calibrations, (b) ambient air measurements on 18-19 July 2015, and (c) ambient air measurements of the entire campaign.For the fit line shown in red, IBBCEAS retrievals with NO 2 mixing ratios < 1 ppbv were excluded.All data were averaged to 60 s.

Figure 8 .
Figure8.Allan deviation plots of data collected while the IB-BCEAS was sampling zero air to determine the optimum integration time of (a) the calculated NO 2 mixing ratios at 890 hPa (1.5 slpm sample rate) and 298 K and (b) the calculated I 2 mixing ratios under the same conditions.

Table 2 .
Selected CEAS detectors for quantification of NO 2 in the near-UV and visible region of the electromagnetic spectrum.NA stands for not available.