Atmospheric distribution of HCN from satellite observations and 3-D model simulations

Hydrogen cyanide (HCN) is an important tracer of biomass burning, but there are significant uncertainties in its atmospheric budget, especially its photochemical and ocean sinks. Here we use a tracer version of the TOMCAT global 3-D chemical transport model to investigate the physical and chemical processes driving the abundance of HCN in the troposphere and stratosphere over the period 2004-2016. The modelled HCN distribution is compared with version 4.1 of the Atmospheric Chemistry Experiment Fourier transform spectrometer (ACE-FTS) HCN satellite data, which provide profiles up to around 42 km, and with ground-based column measurements from the Network for the Detection of Atmospheric Composition Change (NDACC). The long-term ACE-FTS measurements reveal the strong enhancements in upper-tropospheric HCN due to large wildfire events in Indonesia in 2006 and 2015. Our 3-D model simulations confirm previous lower-altitude balloon comparisons that the currently recommended NASA Jet Propulsion Laboratory (JPL) reaction rate coefficient of HCN with OH greatly overestimates the HCN loss. The use of the rate coefficient proposed by in combination with the HCN oxidation by O(D-1) gives good agreement between ACE-FTS observations and the model. Furthermore, the model photochemical loss terms show that the reduction in the HCN mixing ratio with height in the middle stratosphere is mainly driven by the O(1D) sink with only a small contribution from a reaction with OH. From comparisons of the model tracers with ground-based HCN observations we test the magnitude of the ocean sink in two different published schemes . We find that in our 3-D model the two schemes produce HCN abundances which are very different to the NDACC observations but in different directions. A model HCN tracer using the scheme overestimates the HCN concentration by almost a factor of 2, while a HCN tracer using the scheme underestimates the observations by about one-third. To obtain good agreement between the model and observations, we need to scale the magnitudes of the global ocean sinks by factors of 0.25 and 2 for the schemes of and , respectively. This work shows that the atmospheric photochemical sinks of HCN now appear well constrained but improvements are needed in parameterizing the major ocean uptake sink.