1 Introduction
A large number of inorganic and organic species are emitted into the atmosphere from anthropogenic and biogenic sources, making it difficult to detect all of them simultaneously in field experiments and in air quality monitoring stations . These compounds are chemically transformed in the atmosphere by oxidation reactions and thereby form secondary pollutants such as ozone and particles. Most of them react with the primary oxidant in the atmosphere, the hydroxyl radical (), which is formed primarily by the photolysis of ozone and the subsequent reaction of the excited oxygen atom () with water vapour. Therefore, atmospheric measurements of the reactivity, the inverse lifetime of the radical, can be used as a proxy for the total amount of chemically active compounds. The reactivity () is defined as
1 where is the reaction rate coefficient of the compound at a concentration of . As the reactant concentrations are weighted by the reaction rate coefficient, the reactivity describes the total chemical turnover of both the radical and the reactive trace gases and therefore gives the potential for the formation of secondary pollutants from oxidation .
reactivity has been measured in field campaigns for more than 20 years , providing valuable complementary information to individual trace gas measurements. Measurements of single compounds could explain the measured reactivity in some campaigns
Two main methods have been developed to measure reactivity. The comparative reactivity method (CRM) compares the consumption of radicals when either an artificially introduced reactant that is not typically present in the atmosphere (most commonly pyrrole) or reactants in the sampled air react with artificially produced radicals in a reaction volume . The higher the concentration of co-reactants in the sampled air, the less of the artificial reactant is consumed. The artificial reactant is most commonly measured by proton-transfer-reaction mass spectrometry (PTR-MS) , but gas chromatography has also been used . A challenge of the CRM technique is the need for large corrections to account for secondary chemistry in the reaction volume
In the flash photolysis and laser-induced fluorescence method, the loss of radicals is directly measured in a flow tube through which air containing the reactants is sampled . Some instruments use a movable injector to inject artificially produced radicals, allowing the reaction time to be varied . Most instruments, however, produce by flash photolysis of ozone using a short laser pulse at a wavelength of from a quadrupled laser . The following decay is observed with a high time resolution by laser-induced fluorescence after excitation by a pulsed, high-frequency dye laser system providing radiation at a wavelength of 308 . Chemical ionisation mass spectrometry has also been used to detect the radicals .
Instruments for reactivity measurements used in field campaigns were compared in chamber experiments in the large outdoor chamber SAPHIR in 2015 and 2016, which allowed for a systematic investigation of the performances of the instruments under controlled conditions . The results showed that all instruments gave accurate results. Instruments using flash photolysis and detection by laser-induced fluorescence showed the highest precision and accuracy as a high repetition rate of measurements is possible and no corrections are required due to secondary chemistry in the instruments for typical atmospheric conditions. However, as there is no commercial instrument available to detect radicals, this method is currently only used by groups that also have instruments to measure atmospheric concentrations. The complexity of this method has certainly prevented wider use of this type of instrument.
Figure 1
Schematic of the flow tube, through which air containing reactants flows, and the detection cell for measuring radicals by fluorescence. The photolysis and the detection laser beams are directed by turning mirrors into the flow tube and the detection cell, respectively. The laser beams are expanded by lenses to diameters of 30 (photolysis laser) and 8 (detection laser). All parts of the instrument including all utilities, parts, and pumps were built in a 19 in. double rack for the flights on the NASA DC-8 aircraft during the AEROMMA campaign in 2023.
[Figure omitted. See PDF]
In this work, it is shown that it is possible to reduce the technical complexity for the flash photolysis and laser-induced fluorescence method such that it can run autonomously. As a proof of concept, the improved instrument was deployed on the NASA DC-8 aircraft during the AEROMMA (Atmospheric Emissions and Reactions Observed from Megacities to Marine Areas) campaign led by the National Oceanic and Atmospheric Administration (NOAA) in summer 2023.
2 Measurement of reactivity by flash photolysis and laser-induced fluorescenceThe instrument described in this work is based on an instrument that was first used in a field campaign in the Pearl River Delta, China , and in experiments in the SAPHIR atmospheric simulation chamber at Forschungszentrum Jülich, Germany . Previous versions shared several parts, such as the laser system for the detection, with an instrument measuring contemporary radical concentrations . However, the new instrument used on the aircraft was designed to be used as a stand-alone instrument.
The instrument consists of two main parts: (1) a flow tube, through which air containing reactants flows continuously and in which a high concentration of radicals is generated by ozone flash photolysis applying a short laser-pulse at a wavelength of 266 at a low repetition rate (around 1 ), and (2) the detection of radicals by laser-induced fluorescence using a laser at a wavelength of 308 operated at a high repetition rate (13 ) (Fig. ).
Air is sampled in a flow tube through a stainless-steel inlet tube (inner diameter: 8 , variable length of up to several metres) coated with SilcoNert® to minimise losses of reactive species in the inlet. The flow tube made of anodised aluminium has an inner diameter of 40 and a length of 50 . The flow is controlled by a calibrated mass flow controller (Bronkhorst, Low P Series) downstream of the flow tube that is backed up by a scroll pump (Agilent, IDP-3). The flow rate is chosen such that the residence time of air in the flow tube is approximately 2 s. For ambient conditions, the typical flow rate is between 13 and 20 .
radicals are produced in the flow tube by the photolysis of ozone at 266 forming excited oxygen atoms , which subsequently react with water vapour forming 2 radicals on a timescale of nanoseconds for conditions in the flow tube:
In field experiments, ozone and water vapour concentrations in the sampled ambient air are usually high enough to produce a sufficiently large concentration, but in laboratory experiments, in which synthetic air is used, ozone and humidity have to be added. For this purpose, oxygen is photolysed by 185 radiation from a low-pressure mercury lamp in a custom-built ozoniser. Water vapour is added using either a water bubbler or a controlled evaporator mixing (CEM) system (Bronkhorst), in which Milli-Q® water is evaporated. The CEM system allows precise control of the water vapour mixing ratio. The addition of ozone and water vapour requires the availability of bottled synthetic air. If water vapour and/or ozone is added to the sampled ambient air, the dilution of the ambient air needs to be considered in the evaluation. Sensors measure the pressure (Honeywell, precision pressure transducer (PPT)) and relative humidity together with temperature (Vaisala, HUMICAP) at the outlet of the flow tube.
The 266 radiation is generated by a compact quadrupled laser (LUMIBIRD, Ultra 100) that delivers short laser pulses (10 ) with a pulse energy of 20 . The laser operates at a low repetition rate (0.93 to 1 ). The exact frequency is set to minimise the overhead time between two consecutive decay measurements, taking into account the duration over which the decay is observed and the time required to transfer the data from the photon-counting electronics to the computer. The laser beam is expanded by a lens telescope to a diameter of 30 . For typical atmospheric ozone (20 to 50 ) and water vapour (0.2 % to 1.8 %) mixing ratios, the initial concentration is on the order of a few billion per cubic centimetre. The photolysis laser is aligned to illuminate almost the full cross section of the flow tube, so the radicals are approximately homogeneously distributed in the flow tube.
Near the end of the flow tube, a small fraction of the air (1 , litres per minute at standard conditions) is sampled into a low-pressure detection cell through a conical nozzle (pinhole diameter: 0.4 ) that sticks into the centre of the flow tube (Fig. ). The pressure of the detection cell is typically 2.5 to 4 for atmospheric pressure in the flow tube. In the detection cell, radicals are excited at their rotational absorption line of the band transition by a short laser pulse (20 ) at a wavelength of 308 . The 308 laser radiation is generated by a custom-built dye laser (dye: Rhodamine 519) that is pumped by a commercial frequency-doubled laser (Spectra-Physics, Talon). The 616 light produced by the dye laser is frequency doubled to the UV by a barium borate (BBO) crystal inside the laser cavity. The average laser power at a frequency of 13 is up to 200 . Mirrors guide the laser light to the detection cell. The beam size is expanded to a diameter of approximately 8 before entering the detection cell.
Perpendicular to the axis of the air flow and the axis of the laser beam, fluorescence photons are collected by a set of condenser lenses which direct the photons to a microchannel plate photomultiplier (Photek, MCP 325). Opposite of the lens system, a spherical mirror reflects the photons towards the detection system, almost doubling the solid angle from which the photons are detected. As the fluorescence wavelength is resonant to the laser excitation wavelength, the detector is electronically gated, while the laser pulse is applied and photon counting starts shortly after with a delay of approximately 100 . Single photons are counted by photon-counting electronics (Becker and Hickl, PMS-400) for 500 .
After application of the photolysis laser pulse, the (initial concentration ) reacts away in the flow tube. As the reactants have much higher concentrations than the radicals, the time behaviour of the concentration can be described by a pseudo-first-order loss process:
2
The corresponding measured photon counts () include a constant background signal (), which is mainly caused by scattered laser light and detector noise: 3
After the application of the photolysis laser, the decay of the photon counts is recorded for 1 with a time resolution of 1 . The reactivity is calculated from a single exponential fit to the measured photon counts using a Levenberg–Marquardt minimisation procedure. Depending on the detection sensitivity of the measurement and the concentration produced, several measurements are summed or averaged before the fit is applied. A minimum amplitude of the decay curve around 40 counts is sufficient to obtain a reliable fit result.
radicals are also lost in wall reactions on the surface of the flow tube. This is mainly diffusion-limited, and the corresponding decay can be described by a single exponential function. The resulting instrumental zero decay rate () is typically between 1 and 3 . Its value needs to be regularly determined by measuring the loss rate in pure synthetic air containing only ozone and water vapour. The reactivity from the ozone added to the sampled air is negligible ( 0.1 ).
As discussed in , photolysis of reactants does not significantly affect the measurements for typical atmospheric conditions. Deviations from a single exponential decay of the concentration are possible if is produced from secondary chemistry in the flow tube on the timescale of the loss rate . For example, radicals are produced by the reaction of hydroperoxy radicals () with nitric oxide (). However, this can only become relevant under exceptional conditions with high concentrations (e.g. 20 , ) and rapid production of , for example, in the reaction of with carbon monoxide (). In most cases, however, regeneration of does not play a role, even at high concentrations, because the overall reactivity is typically high at these conditions, so the radical lifetime is much shorter than the timescale of the production, and the regeneration is too slow to affect the results. can also be regenerated if the reaction with the reactant forms an adduct which then can decompose and eliminate an radical. An example is the reaction of with isoprene hydroxy hydroperoxide (), which is photochemically produced in the oxidation of isoprene mainly emitted by vegetation . The contribution of such species to the total atmospheric reactivity is typically small, so the underestimation of the reactivity due to regeneration in the flow tube is usually negligible.
Figure 2
Schematic of the simplified dye laser system (wavelength 616 ) pumped by a commercial laser (wavelength 532 ) used for the detection. The dye laser cavity consists of the two end mirrors and two additional mirrors that fold the beam path of the cavity. The dye cuvette is positioned at the Brewster angle. Prisms are used to select the wavelength of the dye laser, which can be tuned by the horizontal position of the end mirror, which is mounted in a motorised mirror mount. The dye laser radiation is frequency converted to a wavelength of 308 by a BBO crystal inside the laser cavity. The UV light is directed to the detection cell in a lens tube system via deflection mirrors (Fig. ).
[Figure omitted. See PDF]
3 Improvements and characterisation of reactivity measurements for (airborne) field campaigns3.1 Improvements of the laser stability and sensitivity
In order to measure reactivity specifically on an aircraft, the instrument needs to be robust against vibration, pressure, and temperature changes. In addition, an autonomous operation is advantageous and a high time resolution in the range of a few seconds is desired, as rapid changes in the ambient reactivity are expected due to the high speed of the aircraft, e.g. when flying through a pollution plume. An instrument meeting these requirements will also be suitable for easy deployment in ground-based field campaigns and may also be used in air quality monitoring stations.
The major challenge of reactivity measurements using the direct detection of radicals by fluorescence is the detection laser, which is typically a dye laser system (Sect. ). In previous versions, the reactivity instrument was part of a measurement system that also included concentration measurements and therefore shared several parts, e.g. the laser system . The aim of developing an reactivity instrument for aircraft applications was to have a stand-alone instrument. Therefore, the previously developed laser system was further developed and optimised for the measurement of reactivity only. A reduced complexity of the system is possible because no absolute concentration needs to be measured and the initial concentration is high.
For the relative time-resolved measurements required in the reactivity instrument, the laser wavelength does not need to be tuned on and off the absorption line, since background signals, e.g. from detector noise or the fluorescence of species other than , do not need to be subtracted from the total signal. They only appear as an offset in the decay curve (Eq. ) as long as they do not change over the time of the decay and/or are small compared to the fluorescence counts. Due to the high concentration, also a lower detection sensitivity due to a lower excitation efficiency than for concentration measurements is acceptable.
In the instrument used for concentration measurements, the tuning of the laser wavelength and the narrow spectral width of the laser in the order of the Doppler-broadened absorption (approximately 3 ) are achieved by a movable etalon in the dye laser resonator . The alignment is sensitive to temperature variations and vibrations, which can especially occur during flights. As the laser wavelength tuning is not required for reactivity measurements, the etalon is removed from the optical design in the stand-alone instrument. In this laser design, the laser wavelength is determined by the optical path through the prisms in the dye laser cavity (Fig. ). This results in a broad spectral width of approximately 0.03 , much wider than the absorption line, so the radical is effectively excited, even if the central wavelength drifts slightly, making the setup robust to small changes in the laser alignment.
The peak wavelength of the dye laser can be tuned by the horizontal position of the cavity end mirror of the dye laser, which is mounted in a motorised mirror mount (Newport, Picomotor). Changes are expected for example in the aircraft, as the cabin pressure is reduced after take-off, changing the refractive index of the air and therefore the centre wavelength of the dye laser output. The laser wavelength is monitored by a high-resolution spectrometer (Ocean Insight, HR4000, resolution: 0.03 ), which allows automatic tuning to the absorption line by software.
In addition, autonomous and stable operation of the dye laser is achieved by the following steps:
-
heating the plate, on which the dye laser cavity is mounted to a slightly higher temperature than the ambient to avoid temperature drifts in the alignment;
-
mounting the mirror that directs the pump laser beam and the BBO crystal in motorised mirror mounts (Newport, Picomotor) to remotely and automatically tune and optimise their position for maximum laser output;
-
mounting all other mirrors of the laser cavity in ultra-stable mirror mounts (Thorlabs, Polaris).
These improvements make the detection robust to small changes in the cavity alignment, allowing the laser to operate at a high performance without operator intervention. For example, the dye laser power achieved during the deployment on the NASA DC-8 aircraft was at least 150 .
The total number of fluorescence photons in the new reactivity instrument is maximised by a high laser repetition rate of 13 , 50 % higher than in previous versions of the instrument. This is possible because the pump laser used in this system (Newport, Talon) delivers a nearly constant pulse energy up to this repetition rate, so the dye laser power scales with the repetition rate. In addition, potential interferences from artificial production in the detection cell, which increase with the laser repetition rate such as photolysis processes of, for example, ozone, are not important . They only slightly increase the background signal in the measured decay, since it can be assumed that the concentration of a species causing the interference does not change over the time of a decay curve of 1 .
The fluorescence yield is further enhanced by back-reflecting a large fraction of the 308 laser at the exit of the detection cell using a mirror with a low transmission of 10 %, so the photon density in the cell is almost doubled. The transmitted laser light is used to monitor the laser power. Again, this is only possible because small artificial photolytic sources of inside the detection cell and an increased background signal from laser scattering do not affect the measured decay. An example for a decay curve from measurements on the NASA DC-8 aircraft is shown in Fig. , demonstrating the high precision of measurements at a high time resolution (here: average of 10 decay curves resulting in a time resolution of 10 ) that could be achieved with the optimised instrument design.
Figure 3
Example of a measured decay during the AEROMMA campaign on 2 August 2023 at an altitude around 500 on the NASA DC-8 aircraft. Nine traces were summed before the fit was applied, resulting in a time resolution of 10 due to some overhead time. A fit of the fluorescence counts to a single exponential function gives . In the measurement shown, the laser power of the 308 detection laser was 147 .
[Figure omitted. See PDF]
3.2 Precision of reactivity measurements using the improved designFigure shows an Allan deviation plot derived from reactivity measurements in humidified synthetic air (1.5 % water vapour mixing ratio) with added ozone resulting in a mixing ratio of 60 ozone. The initial concentration is approximately 8 109 in the flow tube and corresponds to an amplitude of the fluorescence signal of 24 counts for one photolysis laser shot. This can be converted to a sensitivity of the detection of 0.002 counts per 106 radicals per laser power of the 308 detection laser. This number is approximately a factor of 10 lower than the sensitivities achieved in instruments for the measurement of radical concentrations due to the much broader spectral width of the laser used in the new reactivity instrument.
Figure 4
Allan deviation plot (b) of reactivity measurements (a) from measurements after subtracting the zero reactivity value in a mixture of humidified synthetic air (water vapour mixing ratio: 1.5 %) with ozone (60 ). The dashed line gives the Allan deviation expected from Gaussian noise. The distribution (c) of zero measurements for 1 data (red) shows deviations from a Gaussian distribution (lines) that are reduced when three traces are summed before the single exponential fit is applied (blue).
[Figure omitted. See PDF]
The Allan deviation demonstrates a high precision of approximately 0.3 of the reactivity measurement at a time resolution of 1 . An even higher precision of, for example, 0.07 is obtained for an integration time of 10 . The distribution of reactivity measurements (Fig. ) shows deviations from a Gaussian distribution when individual decay curves (1 integration time) are evaluated as seen by the fraction of values that deviate from zero by more than 0.5 . The number of outliers can be significantly reduced if at least three decay curves are summed before applying the exponential fit because small systematic deviations from a single exponential behaviour are smoothed out. This demonstrates that an integration time of at least 3 is recommended to ensure a statistical error of the reactivity measurements.
The produced initial concentration in ambient air may be lower than in the laboratory measurements, as the mixing ratios of ozone and water vapour are highly variable. However, the Allan deviation of measurements in synthetic air shows that a high time resolution in the range of seconds can still be easily achieved by summing several traces. The detection limit is much lower than the reactivity in ambient air, which typically has minimum values around 1 even in clean environments (Sect. ).
For the evaluation of the data collected during the AEROMMA campaign on board the NASA DC-8 aircraft, typically 5 to 10 decay curves were summed before calculating the reactivity from the decay curve, in order to achieve a sufficiently high precision of the data (e.g. minimum amplitude of 40 counts). As the repetition rate of the photolysis laser beam was 0.93 (Sect. ), this results in a time resolution of measurements of 5.5 to 11 . For conditions of very low water vapour mixing ratios of less than 0.1 % encountered at high altitudes or in dry areas, up to 30 to 40 decay curves had to be added up. At even lower water vapour mixing ratios of less than 0.05 %, experienced in the free troposphere at altitudes above 8 , the instrument cannot operate without the addition of water vapour, which was not foreseen in this campaign.
3.3 Considerations for the use on an aircraftWhen the instrument is operated on board an aircraft, the air is drawn into the flow tube from outside the aircraft; therefore the pressure in the flow tube is similar to the ambient pressure, which decreases with altitude. On the NASA DC-8 aircraft, a standard inlet system provided by NASA was used. A restrictor that was part of the inlet system resulted in an approximately 100 higher pressure than ambient pressure in the flow tube during the flight. A reduced air pressure in the flow tube compared to ground conditions has several consequences, some of which require adjustments to the operation of the instrument during the flight:
-
In order to keep a similar residence time of the air in the flow tube, the sampling flow rate controlled by a mass flow controller is automatically adjusted, so the volumetric flow rate becomes similar to the changing pressure. If this is not done, the residence time of the air in the flow tube can become shorter than the time between two successive photolysis laser shots, so the decay curve would be affected by the flushing out of the radicals.
-
The diffusion rate of radicals increases, which may affect the zero decay rate.
-
The initial concentration in the flow tube is reduced because the concentrations of ozone and humidity are reduced due to the reduced number density of molecules, leading to a lower amplitude of the fluorescence counts per shot of the photolysis laser.
-
A critical nozzle is used to sample air from the flow tube into the detection cell, which inherently results in a constant volumetric flow rate into the cell. Since the volume flow rate through the flow tube is also constant, the fraction of air sampled into the detection cell remains independent of the flight altitude.
-
As the mass flow sampled into the detection cell decreases, the pressure in the detection cell decreases if the power of the pump downstream of the cell is kept constant, leading to a decrease in the detection sensitivity, if this is optimised for ambient pressure on the ground (typical detection cell pressure: 1.5 to 4 ). The pressure could be increased by adjusting the flow restriction using a butterfly valve between the detection cell and the pump (Fig. ). However, only a manual valve was installed in the AEROMMA campaign.
Overall, there is a significant reduction in the detection sensitivity at high altitude of the aircraft, which can be compensated for by a longer integration time. This can typically be accepted as the variability of the reactivity is expected to be small at high altitude.
3.4 Characterisation of the instrument zero decay rate
The instrument's zero is the loss of radicals in the absence of an reactant due to the wall loss of (Sect. ). It can be described as a pseudo-first-order loss process and must be characterised thoroughly for an accurate determination of reactivity measurements.
In the laboratory characterisation experiments high-purity synthetic air mixed from evaporated liquid nitrogen and oxygen (purity 99.9999 %, Linde) and ultra-pure Milli-Q® water was used. This makes it unlikely that the zero decay is caused by the introduction of reactants in the laboratory measurements. Tests in the field during the AEROMMA campaign were performed with ultra-pure bottled synthetic air (Linde) and Milli-Q® water but resulted partly in up to 1 higher zero values than measured in the laboratory. The exact value depended on the specific synthetic air bottle, so the higher values were likely due to impurities. Therefore, zero measurements from the laboratory were used for the evaluation of measurements of the AEROMMA campaign.
The dependence of the zero decay rate on pressure and humidity was characterised in laboratory measurements. The pressure in the flow tube was varied by inserting a valve in the inlet line acting as a variable flow restrictor. The volume flow rate in the flow tube was kept constant during these tests by automatically adjusting the setpoint of the mass flow controller downstream of the flow tube, as done when operating the instrument on an aircraft. The humidity of the synthetic air was varied by changing the amount of water that was evaporated in the humidification system (Bronkhorst, CEM).
The overall dependence of the zero decay rate on humidity and pressure is small for pressures higher than 300 and water vapour mixing ratios higher than 0.5 % (Fig. ). The values are between 1.6 and 1.9 . Only at very low pressures of 200 , lower than the pressure in the flow tube experienced during the AEROMMA campaign ( 350 ), does the zero decay time increase to values higher than 2.5 . The small dependence of the zero decay time on pressure ( 300 ) and the water vapour mixing ratio can be expressed as
4
Fitting the values obtained in the laboratory experiments gives values of the parameters of 0.09 , 2.7 10−4 , 5.4 10−4 , and 2.1 . The differences between the parametrisation and the measured values are less than 0.3 , which is within the reproducibility of the zero decay measurement, giving a lower limit of the accuracy of the reactivity measurements.
Figure 5
Zero decay values depending on the pressure (). The colours give the water vapour mixing ratio during the zero measurement. The dependence can be described by Eq. () for pressure values higher than 300 . The measurements of the higher zero decay values below 300 were reproducible and indicate that other effects like the increased diffusion of radicals led to a higher wall loss. During the AEROMMA campaign on the NASA DC-8 aircraft, the pressure in the flow tube was always higher than 350 .
[Figure omitted. See PDF]
The zero decay rate is most likely caused by the loss of radicals on the wall of the flow tube. Due to its high reactivity, it can be assumed that the probability of wall loss on metal surfaces, such as in the reactivity instrument, is very high for radicals, so the total loss rate is mainly limited by diffusion . However, this description only holds if the initial concentration is homogeneously distributed in the flow tube. Deviations from this can occur if either the photon density of the 266 photolysis laser is not homogeneous or the laser beam is not well aligned.
With decreasing pressure, the diffusion of radicals towards the wall of the flow tube increases as the diffusion coefficient is inversely proportional to the pressure. Therefore, an increase in the zero decay rate is only clearly visible at lower pressures for low water vapour mixing ratios. The increase in the zero decay rate with increasing water vapour mixing ratio could be caused by a higher probability of the loss at the wall. However, it cannot be fully excluded that there was a small contamination of the water.
Only one other reactivity instrument has been deployed on an aircraft . This instrument uses a movable injector inside a flow tube, through which a small amount of humidified air containing and radicals produced by water photolysis at a wavelength of 185 from a mercury lamp is injected. In this instrument, the wall loss is assumed to be independent of pressure, but an loss due to impurities in the injected air was observed . Since the mass flow rate of the injected air was the same at all altitudes but the volume flow rate of sampled air decreased with height, a pressure dependence of the zero decay rate appeared due to the change in the dilution of contaminants . The effect changed between the campaigns as the contaminant concentrations varied. As there is no air injection in the instrument in this work, the type of pressure-dependent zero decay rate such as observed by and does not apply.
3.5 Validation of the reactivity measurements using a propane gas standardThe accuracy of the measurements was tested by providing well-defined concentrations of propane in humidified synthetic air. The gas mixture was prepared in a canister by mixing propane (purity 99.5 %, Linde) with nitrogen. The resulting mixing ratio of (2166 22) was measured using the total organic carbon method described in detail in , in which all carbon is converted to on a heated palladium catalyst. The is measured with a high accuracy using cavity ring-down spectroscopy (Picarro, G1301) allowing the propane concentration to be back-calculated.
For the tests with the reactivity instrument, a small flow of a few millilitres per minute of this propane gas standard was mixed with a large flow of humidified synthetic air that was prepared in the same way as done for the zero measurements. The inlet of the instrument was overflowed with this well-defined mixture. Contaminations from other reactants are therefore not expected to affect the measurements. All flows were provided by calibrated mass flow controllers (Bronkhorst, El-Flow Series).
Figure shows the measured reactivity, when the propane concentration was varied between 1 1012 and 30 1012 , resulting in reactivity values between 1.5 and 34 . Two sets of measurements were performed at different pressures (997 and 600 ) that were in the range of typical pressure values experienced during the flights of the AEROMMA campaign. Due to the dilution at lower pressure, propane concentrations and reactivity values were also lower in this set of measurements. The slope of a regression analysis gives the reaction rate constant of the reaction with propane. The accuracy of the value includes the accuracy of the propane concentration in the flow tube of 1.4 %, which takes into account the uncertainties in the flows and the propane concentration in the canister, and the statistical error of the slope of 2.7 %. The total accuracy is therefore 4.1 %.
Figure 6
Determination of the bimolecular rate coefficient of the reaction with propane in air from reactivity measurements when the inlet of the instrument is overflowed with a mixture of propane in synthetic air. The rate coefficient can be described by an Arrhenius expression and is independent of pressure . The rate coefficient is determined at two different values of pressure to test the accuracy of measurements for the operational conditions on the aircraft.
[Figure omitted. See PDF]
Table 1Rate constant of the reaction with propane determined at two pressure values () by the reactivity instrument () compared to recommendations by IUPAC (, ) and NASA-JPL (, ) at the respective temperature () of the measurement. The rate coefficient can be described by an Arrhenius expression and is independent of pressure. The rate coefficient is determined at two different values of pressure () to test the accuracy of measurements for the operational conditions on the aircraft.
/ | / | /10−12 | /10−12 | /10−12 |
---|---|---|---|---|
997 | 299 | 1.11 0.05 | 1.08 0.2 | 1.12 0.04 |
602 | 296 | 1.13 0.05 | 1.05 0.2 | 1.09 0.04 |
The reaction rate coefficient of the reaction with propane has been investigated in several studies. Calculations using the Arrhenius expressions in the recommendations show excellent agreement within 8 % (IUPAC, ) and 4 % (NASA-JPL, ) with the values obtained from the reactivity measurements for both pressure values (Table ). This demonstrates the high accuracy of the new instrument's reactivity measurements, as shown for previous versions of the instrument . It also shows that this high accuracy is also achieved for measurements on board an aircraft when the sampled ambient air is at low pressure.
4 reactivity measurements on board the DC-8 aircraft during the AEROMMA campaignThe improved reactivity instrument was deployed on the NASA DC-8 aircraft during the AEROMMA campaign in summer 2023. The aircraft carried about 30 different instruments to measure gas-phase species and aerosol properties. Flights were conducted over the Pacific Ocean and over major urban areas in North America, including Chicago, Toronto, New York City, and Los Angeles. The detailed analysis of the measurements will be presented in separate papers. Here, two flights over the Pacific Ocean are shown to demonstrate the performance of the instrument when low reactivity values are observed, and only a few reactants are expected to contribute significantly to the reactivity; therefore these measurements can demonstrate the high precision and accuracy of the measurements.
The measured reactivity differs from the value outside the aircraft. One reason is the slightly different ( 5 %) number densities of the reactants in the flow tube, where the pressure is close to ambient pressure and the temperature is at cabin temperature (295 to 305 ). The other reason is the possible influence of temperature and pressure on the rate constants.
Figure 7
Time series and map of reactivity (1 average) measured during the flights on 21 and 23 June 2023 over the Pacific Ocean starting from Palmdale, CA, USA (red circle). Coloured areas indicate contributions from measured reactants. Biogenic volatile organic compounds (BVOCs) include dimethyl sulfide (DMS), isoprene, and monoterpenes, and oxygenated volatile organic compounds (OVOCs) include formaldehyde, acetaldehyde, ethanol, methyl vinyl ketone, methacrolein, nonanal, octanal, dimethyl sulfoxide (DMSO), and hydroperoxymethyl thioformate (HPMTF). The reactivity is given for conditions inside the flow tube of the instrument, where the pressure is close to ambient pressure, but the temperature is at cabin temperature (295 to 305 ). The zero decay value is subtracted. Error bars are the 1 precision of measurements, and the grey lines indicate the total uncertainty of the reactivity values of 0.5 . The reactivity measurements started after the transit at high altitude (8 ), where water vapour mixing ratios were too low to produce sufficiently high concentrations in the flow tube.
[Figure omitted. See PDF]
Table 2Instruments used for measurements shown in Fig. or used to analyse reactivity data from the AEROMMA campaign.
Method | Species | Data | Instrument/reference |
---|---|---|---|
version | |||
Laser flash photolysis/laser-induced fluorescence | reactivity | R0∗ | this work |
Diode laser hygrometer | RA∗ | ||
Off-axis integrated cavity output spectroscopy | R0∗ | LGR F-N2O/CO-23r, | |
Cavity ring-down spectroscopy | Picarro 2401-m, | ||
Chemiluminescence | R0∗ | ||
Laser-induced fluorescence | R0∗ | ||
Laser-induced fluorescence | , | R0∗ | |
Laser-induced fluorescence | R0∗ | ||
Chemical ionisation mass spectrometry | peroxyacetyl nitrate (PAN), | R0∗ | , |
hydroperoxymethyl thioformate (HPMTF) | R1∗ | ||
Chemical ionisation mass spectrometry | dimethyl sulfoxide (DMSO) | R0∗ | |
Proton-transfer-reaction mass spectrometry | dimethyl sulfide (DMS), | R1∗ | |
monoterpenes, isoprene, | |||
methyl vinyl ketone + methacrolein, | |||
ethanol, nonanal, octanal | |||
acetaldehyde |
∗ .
The zero decay rate is calculated from the parametrisation of Eq. (), using the ambient humidity measurements from a diode laser hygrometer on board the aircraft and the pressure measurements in the flow tube of the reactivity instrument, resulting in values between 1.65 and 1.85 . Figure shows a map and the time series of the measured reactivity after subtracting the zero decay value. The measured total reactivity is compared to calculations of the reactivity (Fig. ) using reactant measurements from several instruments listed in Table . Measurements included inorganic compounds such as carbon monoxide (), sulfur dioxide (), nitrogen dioxide (), nitric oxide (), and ozone () and organic compounds such as methane and formaldehyde ().
As would be expected in the clean, marine air, the reactivity is low, with values between 1 and 1.5 over the Pacific Ocean and over land when the aircraft flew at high altitude. The measured reactivity can largely be explained by the presence of carbon monoxide and methane, having mixing ratios of around 100 and 2 , respectively, which gives a reactivity of approximately 1 . Individual contributions from other measured reactants are less than 0.2 , of which formaldehyde (mixing ratios around 200 ) and dimethyl sulfide (mixing ratios around 100 ) are the largest.
The measurements demonstrate that the reactivity can be measured with high precision and accuracy on an aircraft using the new instrument. Due to the large number of reactants in the atmosphere, the total reactivity measurements are expected to be rather higher than the calculations using the limited number of reactant measurements. The difference between measured total reactivity and calculations using individual reactants is in most cases less than 0.4 , which is less than the accuracy of the measurement due to the uncertainty of the zero decay rate. Since the reactivity of measured non-methane hydrocarbons was very small, the contribution of unmeasured oxygenated VOCs (OVOCs) that could be expected from their oxidation was likely to be small. Differences are higher with values of up to 0.6 after the start of the flight and show a decreasing trend. This could be due to a slight drift in the value of the zero decay rate, which could be caused by contamination if dirty air was sampled immediately after take-off, leading to an increased wall loss.
At high altitudes (approximately 4 ), the water vapour mixing ratio drops below 0.5 %, so the initial concentration produced per photolysis laser shot is only around 1 108 . At the same time, the expected reactivity becomes very low. Although up to 40 traces are summed before applying the fit procedure, resulting in an amplitude similar to traces acquired at higher water vapour mixing ratios, the scatter of the data is significantly increased (e.g. 17:30 to 18:00 PST local time, 21 June 2023, Fig. ). At altitudes of more than 8 , the produced concentrations are too low to evaluate the decay.
provided the first reactivity measurements from an aircraft over marine environments during the INTEX-B measurement campaign in April 2006. An instrument was used, in which the decay was also measured directly by laser-induced fluorescence and in which a movable injector for radicals was used to vary the reaction time. The same group measured reactivity again with this instrument during the four ATom campaigns between 2016 and 2018 . The largest uncertainty in their measurements was due to the uncertainty in the zero decay value, as in the measurements in this work. In order to reduce the uncertainty, and adjusted the zero decay rate so that the measured reactivity agreed with calculations of reactants for certain parts of the flight at high altitudes over the oceans, assuming that there were no relevant unmeasured reactants in this clean air.
The reactivity observed during the INTEX-B and ATom campaigns over the oceans was maximally 2 and dropped to low values around 0.2 at high altitudes in the free troposphere. As during the AEROMMA campaign over the Pacific Ocean, the reactivity was mainly due to and methane, with small contributions from oxygenated organic compounds. Differences between the measured reactivity and calculations using reactant measurements were also similar and on the order of the uncertainty of the zero decay rate. attempted to estimate the contribution of unmeasured species to the reactivity by a statistical approach and found that an reactivity between 0.4 and 0.7 cannot be explained by measured reactants in the marine boundary layer.
5 ConclusionsAn instrument for measuring reactivity using laser flash photolysis and the direct detection of the decay by laser-induced fluorescence has been further developed in this work for use in field experiments in challenging environments such as on board an aircraft. This instrument can operate largely autonomously and with a high sensitivity, providing a high precision of less than 0.3 with a time resolution on the order of seconds. The accuracy of the measurements is mainly limited by the uncertainty in the zero decay rate of 0.5 . Validation with a well-defined mixture of propane in synthetic air at two different pressures demonstrates that the measured reactivity values give reaction rate coefficients in excellent agreement with values calculated from the Arrhenius expression recommended by IUPAC within 8 % and NASA-JPL within 4 % .
The effort required to operate the instrument has been greatly reduced compared to previous versions. By simplifying the dye laser system used to detect the radicals, the instrument is robust against vibrations and changes in the temperature and the pressure. If necessary, motorised mounts for the optical elements can automatically compensate for small changes in the laser alignment.
The good performance is demonstrated during the AEROMMA campaign, where measurements were conducted on board the NASA DC-8 aircraft. Measurements in a clean environments above the Pacific Ocean gave low reactivity values in the range 1.5 to 2.0 . These low values are well explained by measured reactants with major contributions from carbon monoxide and methane within the uncertainty of the zero decay value that needs to be subtracted from the measurements.
Overall, the complexity of the new reactivity instrument could be significantly reduced compared to previous versions of the instrument and therefore has the potential for a wide application in laboratory and field experiments. Widespread use of reactivity measurements would provide valuable information on the load of pollutants in the atmosphere and the potential for the formation of secondary pollutants from their chemical transformation .
Data availability
Data from the 21 June 2023 flight are available at
Author contributions
HF wrote the manuscript. AS, FB, AN, and MF performed reactivity measurements, and FB and RB performed TOC measurements. FH and AW contributed to the application of the instrument in the AEROMMA campaign. KHB, SSB, MMC, GSD, GIG, CMJ, JP, MAR, AWR, AW, NBS, RHS, CES, PRV, CW, EMW, LX, and KZ contributed to measurements used for the analysis of the measurements during the AEROMMA campaign. All co-authors discussed the content of the paper and contributed to the writing.
Competing interests
At least one of the (co-)authors is a member of the editorial board of Atmospheric Measurement Techniques. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
Acknowledgements
The authors thank the Department of Engineering and workshops at Forschungszentrum Jülich for the great support to develop the instrument. The authors also thank the entire team of the AEROMMA campaign for the support and opportunity to perform measurements during the campaign, specifically for the provision of additional data for comparing measured reactivity with calculations using single reactant concentrations (formaldehyde: Glenn M. Wolfe, NASA Goddard Space Flight Center; PAN: Gordon Novak, NOAA).
Financial support
This research has been supported by the Klaus Tschira Stiftung (Klaus Tschira Boost Fund) and the National Oceanic and Atmospheric Administration (grant no. NA22OAR43200151). The article processing charges for this open-access publication were covered by the Forschungszentrum Jülich.
Review statement
This paper was edited by Hartwig Harder and reviewed by two anonymous referees.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2025. This work is published under https://creativecommons.org/licenses/by/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Hydroxyl radical (
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details




















1 Institute of Climate and Energy Systems, ICE-3: Troposphere, Forschungszentrum Jülich GmbH, Jülich, Germany; Department of Physics, University of Cologne, Cologne, Germany
2 Institute of Climate and Energy Systems, ICE-3: Troposphere, Forschungszentrum Jülich GmbH, Jülich, Germany
3 NOAA Chemical Sciences Laboratory, Boulder, Colorado, USA; Cooperative Institute for Research in Environmental Sciences, University of Colorado Boulder, Boulder, Colorado, USA; now at: Department of Mechanical Engineering, University of Colorado, Boulder, Colorado, USA
4 NOAA Chemical Sciences Laboratory, Boulder, Colorado, USA; Department of Chemistry, University of Colorado, Boulder, Colorado, USA
5 NOAA Chemical Sciences Laboratory, Boulder, Colorado, USA
6 NASA Langley Research Center, Hampton, Virginia, USA
7 NOAA Chemical Sciences Laboratory, Boulder, Colorado, USA; Cooperative Institute for Research in Environmental Sciences, University of Colorado Boulder, Boulder, Colorado, USA
8 NOAA Chemical Sciences Laboratory, Boulder, Colorado, USA; Cooperative Institute for Research in Environmental Sciences, University of Colorado Boulder, Boulder, Colorado, USA; now at: NOAA Global Monitoring Laboratory, Boulder, Colorado, USA
9 NOAA Chemical Sciences Laboratory, Boulder, Colorado, USA; now at: National Center for Atmospheric Research, Boulder, Colorado, USA
10 NOAA Chemical Sciences Laboratory, Boulder, Colorado, USA; Cooperative Institute for Research in Environmental Sciences, University of Colorado Boulder, Boulder, Colorado, USA; now at: Department of Energy, Environmental and Chemical Engineering, Washington University in St. Louis, St. Louis, Missouri, USA