The fate of NOx emissions due to nocturnal oxidation at high latitudes: 1-D simulations and sensitivity experiments

The fate of nitrogen oxide pollution during high-latitude winter is controlled by reactions of dinitrogen pentoxide (N2O5) and is highly affected by the competition between heterogeneous atmospheric reactions and deposition to the snowpack. MISTRA (MIcrophysical STRAtus), a 1-D photochemical model, simulated an urban pollution plume from Fairbanks, Alaska to investigate this competition of N2O5 reactions and explore sensitivity to model parameters. It was found that dry deposition of N2O5 made up a significant fraction of N2O5 loss near the snowpack, but reactions on aerosol particles dominated loss of N2O5 over the integrated atmospheric column. Sensitivity experiments found the fate of NOx emissions were most sensitive to NO emission flux, photolysis rates, and ambient temperature. The results indicate a strong sensitivity to urban area density, season and clouds, and temperature, implying a strong sensitivity of the results to urban planning and climate change. Results suggest that secondary formation of particulate (PM2.5) nitrate in the Fairbanks downtown area does not contribute significant mass to the total PM2.5 concentration, but appreciable amounts are formed downwind of downtown due to nocturnal NOx oxidation and subsequent reaction with ammonia on aerosol particles

Model study of multiphase DMS oxidation with a focus on halogens

We studied the oxidation of dimethylsulfide (DMS) in the marine boundary layer (MBL) with a one-dimensional numerical model and focused on the influence of halogens. Our model runs show that there is still significant uncertainty about the end products of the DMS addition pathway, which is especially caused by uncertainty in the product yield of the reaction of the intermediate product methyl sulfinic acid (MSIA) with OH. BrO strongly increases the importance of the addition branch in the oxidation of DMS even when present at mixing ratios smaller than 0.5pmol mol^-1. The inclusion of halogen chemistry leads to higher DMS oxidation rates and smaller DMS to SO₂ conversion efficiencies. The DMS to SO₂ conversion efficiency is also drastically reduced under cloudy conditions. In cloud-free model runs between 5 and 15% of the oxidized DMS reacts further to particulate sulfur, in cloudy runs this fraction is almost 100%. Sulfate production by HOCl_aq and HOBr_aq is important in cloud droplets even for small Br^- deficits and related small gas phase halogen concentrations. In general, more particulate sulfur is formed when halogen chemistry is included. A possible enrichment of HCO₃^- in fresh sea salt aerosol would increase pH values enough to make the reaction of S(IV)^* (=SO_2,aq+HSO₃^-+SO₃^2-) with O₃ dominant for sulfate production. It leads to a shift from methyl sulfonic acid (MSA) to non-sea salt sulfate (nss-SO₄^2-) production but increases the total nss-SO₄^2- only somewhat because almost all available sulfur is already oxidized to particulate sulfur in the base scenario. We discuss how realistic this is for the MBL. We found the reaction MSA_aq+OH to contribute about 10% to the production of nss-SO₄^2- in clouds. It is unimportant for cloud-free model runs. Overall we find that the presence of halogens leads to processes that decrease the albedo of stratiform clouds in the MBL

Modelling iodide – iodate speciation in atmospheric aerosol: Contributions of inorganic and organic iodine chemistry

The speciation of iodine in atmospheric aerosol is currently poorly understood. Models predict negligible iodide concentrations but accumulation of iodate in aerosol, both of which is not confirmed by recent measurements. We present an updated aqueous phase iodine chemistry scheme for use in atmospheric chemistry models and discuss sensitivity studies with the marine boundary layer model MISTRA. These studies show that iodate can be reduced in acidic aerosol by inorganic reactions, i.e., iodate does not necessarily accumulate in particles. Furthermore, the transformation of particulate iodide to volatile iodine species likely has been overestimated in previous model studies due to negligence of collision-induced upper limits for the reaction rates. However, inorganic reaction cycles still do not seem to be sufficient to reproduce the observed range of iodide – iodate speciation in atmospheric aerosol. Therefore, we also investigate the effects of the recently suggested reaction of HOI with dissolved organic matter to produce iodide. If this reaction is fast enough to compete with the inorganic mechanism, it would not only directly lead to enhanced iodide concentrations but, indirectly via speed-up of the inorganic iodate reduction cycles, also to a decrease in iodate concentrations. Hence, according to our model studies, organic iodine chemistry, combined with inorganic reaction cycles, is able to reproduce observations. The presented chemistry cycles are highly dependent on pH and thus offer an explanation for the large observed variability of the iodide – iodate speciation in atmospheric aerosol

Modeling iodide ? iodate speciation in atmospheric aerosol

International audienceThe speciation of iodine in atmospheric aerosol is currently poorly understood. Models predict negligible iodide concentrations and accumulation of iodate in aerosol, both of which is not confirmed by recent measurements. We present an updated aqueous phase iodine chemistry scheme for use in atmospheric chemistry models and discuss sensitivity studies with the marine boundary layer model MISTRA. These studies show that iodate can be reduced in acidic aerosol by inorganic reactions, i.e., iodate does not necessarily accumulate in particles. Furthermore, the transformation of particulate iodide to volatile iodine species likely has been overestimated in previous model studies due to negligence of collision-induced upper limits for the reaction rates. However, inorganic reaction cycles still do not seem to be sufficient to reproduce the observed range of iodide ? iodate speciation in atmospheric aerosol. Therefore, we also investigate the effects of the recently suggested reaction of HOI with dissolved organic matter to produce iodide. If this reaction is fast enough to compete with the inorganic mechanism, it would not only directly lead to enhanced iodide concentrations but, indirectly via speed-up of the inorganic iodate reduction cycles, also to a decrease in iodate concentrations. Hence, according to our model studies, organic iodine chemistry, combined with inorganic reaction cycles, is able to reproduce observations. The presented chemistry cycles are highly dependent on pH and thus offer an explanation for the large observed variability of the iodide ? iodate speciation in atmospheric aerosol

Accommodation coefficient of HOBr on deliquescent sodium bromide aerosol particles

Uptake of HOBr on sea salt aerosol, sea salt brine or ice is believed to be a key process providing a source of photolabile bromine (Br₂) and sustaining ozone depletion cycles in the Arctic troposphere. In the present study, uptake of HOBr on sodium bromide (NaBr) aerosol particles was investigated at an extremely low HOBr concentration of 300 cm^-3 using the short-lived radioactive isotopes ^83-86Br. Under these conditions, at maximum one HOBr molecule was taken up per particle. The rate of uptake was clearly limited by the mass accommodation coefficient, which was calculated to be 0.6 ± 0.2. This value is a factor of 10 larger than estimates used in earlier models. The atmospheric implications are discussed using the box model "MOCCA'', showing that the increase of the accommodation coefficient of HOBr by a factor of 10 only slightly affects net ozone loss, but significantly increases chlorine release

Modeling chemistry in and above snow at Summit, Greenland – Part 1: Model description and results

Sun-lit snow is increasingly recognized as a chemical reactor that plays an active role in uptake, transformation, and release of atmospheric trace gases. Snow is known to influence boundary layer air on a local scale, and given the large global surface coverage of snow may also be significant on regional and global scales. We present a new detailed one-dimensional snow chemistry module that has been coupled to the 1-D atmospheric boundary layer model MISTRA. The new 1-D snow module, which is dynamically coupled to the overlaying atmospheric model, includes heat transport in the snowpack, molecular diffusion, and wind pumping of gases in the interstitial air. The model includes gas phase chemical reactions both in the interstitial air and the atmosphere. Heterogeneous and multiphase chemistry on atmospheric aerosol is considered explicitly. The chemical interaction of interstitial air with snow grains is simulated assuming chemistry in a liquid-like layer (LLL) on the grain surface. The coupled model, referred to as MISTRA-SNOW, was used to investigate snow as the source of nitrogen oxides (NOx) and gas phase reactive bromine in the atmospheric boundary layer in the remote snow covered Arctic (over the Greenland ice sheet) as well as to investigate the link between halogen cycling and ozone depletion that has been observed in interstitial air. The model is validated using data taken 10 June–13 June, 2008 as part of the Greenland Summit Halogen-HOx experiment (GSHOX). The model predicts that reactions involving bromide and nitrate impurities in the surface snow can sustain atmospheric NO and BrO mixing ratios measured at Summit, Greenland during this period

Tracer Measurements in Growing Sea Ice Support Convective Gravity Drainage Parameterizations

Gravity drainage is the dominant process redistributing solutes in growing sea ice. Modeling gravity drainage is therefore necessary to predict physical and biogeochemical variables in sea ice. We evaluate seven gravity drainage parameterizations, spanning the range of approaches in the literature, using tracer measurements in a sea ice growth experiment. Artificial sea ice is grown to around 17 cm thickness in a new experimental facility, the Roland von Glasow air‐sea‐ice chamber. We use NaCl (present in the water initially) and rhodamine (injected into the water after 10 cm of sea ice growth) as independent tracers of brine dynamics. We measure vertical profiles of bulk salinity in situ, as well as bulk salinity and rhodamine in discrete samples taken at the end of the experiment. Convective parameterizations that diagnose gravity drainage using Rayleigh numbers outperform a simpler convective parameterization and diffusive parameterizations when compared to observations. This study is the first to numerically model solutes decoupled from salinity using convective gravity drainage parameterizations. Our results show that (1) convective, Rayleigh number‐based parameterizations are our most accurate and precise tool for predicting sea ice bulk salinity; and (2) these parameterizations can be generalized to brine dynamics parameterizations, and hence can predict the dynamics of any solute in growing sea ic

Modeling chemistry in and above snow at Summit, Greenland – Part 2: Impact of snowpack chemistry on the oxidation capacity of the boundary layer

The chemical composition of the boundary layer in snow covered regions is impacted by chemistry in the snowpack via uptake, processing, and emission of atmospheric trace gases. We use the coupled one-dimensional (1-D) snow chemistry and atmospheric boundary layer model MISTRA-SNOW to study the impact of snowpack chemistry on the oxidation capacity of the boundary layer. The model includes gas phase photochemistry and chemical reactions both in the interstitial air and the atmosphere. While it is acknowledged that the chemistry occurring at ice surfaces may consist of a true quasi-liquid layer and/or a concentrated brine layer, lack of additional knowledge requires that this chemistry be modeled as primarily aqueous chemistry occurring in a liquid-like layer (LLL) on snow grains. The model has been recently compared with BrO and NO data taken on 10 June–13 June 2008 as part of the Greenland Summit Halogen-HOx experiment (GSHOX). In the present study, we use the same focus period to investigate the influence of snowpack derived chemistry on OH and HOx + RO2 in the boundary layer. We compare model results with chemical ionization mass spectrometry (CIMS) measurements of the hydroxyl radical (OH) and of the hydroperoxyl radical (HO2) plus the sum of all organic peroxy radicals (RO2) taken at Summit during summer 2008. Using sensitivity runs we show that snowpack influenced nitrogen cycling and bromine chemistry both increase the oxidation capacity of the boundary layer and that together they increase the midday OH concentrations. Bromine chemistry increases the OH concentration by 10–18 % (10 % at noon LT), while snow sourced NOx increases OH concentrations by 20–50 % (27 % at noon LT). We show for the first time, using a coupled one dimensional snowpack-boundary layer model, that air-snow interactions impact the oxidation capacity of the boundary layer and that it is not possible to match measured OH levels without snowpack NOx and halogen emissions. Model predicted HONO compared with mistchamber measurements suggests there may be an unknown HONO source at Summit. Other model predicted HOx precursors, H2O2 and HCHO, compare well with measurements taken in summer 2000, which had lower levels than other years. Over 3 days, snow sourced NOx contributes an additional 2 ppb to boundary layer ozone production, while snow sourced bromine has the opposite effect and contributes 1 ppb to boundary layer ozone loss

The tropospheric processing of acidic gases and hydrogen sulphide in volcanic gas plumes as inferred from field and model investigations

Improving the constraints on the atmospheric fate and depletion rates of acidic compounds persistently emitted by non-erupting (quiescent) volcanoes is important for quantitatively predicting the environmental impact of volcanic gas plumes. Here, we present new experimental data coupled with modelling studies to investigate the chemical processing of acidic volcanogenic species during tropospheric dispersion. Diffusive tube samplers were deployed at Mount Etna, a very active open-conduit basaltic volcano in eastern Sicily, and Vulcano Island, a closed-conduit quiescent volcano in the Aeolian Islands (northern Sicily). Sulphur dioxide (SO₂), hydrogen sulphide (H₂S), hydrogen chloride (HCl) and hydrogen fluoride (HF) concentrations in the volcanic plumes (typically several minutes to a few hours old) were repeatedly determined at distances from the summit vents ranging from 0.1 to ~10 km, and under different environmental conditions. At both volcanoes, acidic gas concentrations were found to decrease exponentially with distance from the summit vents (e.g., SO₂ decreases from ~10 000 μg/m³at 0.1 km from Etna's vents down to ~7 μg/m³ at ~10 km distance), reflecting the atmospheric dilution of the plume within the acid gas-free background troposphere. Conversely, SO₂/HCl, SO₂/HF, and SO₂/H₂S ratios in the plume showed no systematic changes with plume aging, and fit source compositions within analytical error. Assuming that SO₂ losses by reaction are small during short-range atmospheric transport within quiescent (ash-free) volcanic plumes, our observations suggest that, for these short transport distances, atmospheric reactions for H₂S and halogens are also negligible. The one-dimensional model MISTRA was used to simulate quantitatively the evolution of halogen and sulphur compounds in the plume of Mt. Etna. Model predictions support the hypothesis of minor HCl chemical processing during plume transport, at least in cloud-free conditions. Larger variations in the modelled SO₂/HCl ratios were predicted under cloudy conditions, due to heterogeneous chlorine cycling in the aerosol phase. The modelled evolution of the SO₂/H₂S ratios is found to be substantially dependent on whether or not the interactions of H₂S with halogens are included in the model. In the former case, H₂S is assumed to be oxidized in the atmosphere mainly by OH, which results in minor chemical loss for H₂S during plume aging and produces a fair match between modelled and measured SO₂/H₂S ratios. In the latter case, fast oxidation of H₂S by Cl leads to H₂S chemical lifetimes in the early plume of a few seconds, and thus SO₂ to H₂S ratios that increase sharply during plume transport. This disagreement between modelled and observed plume compositions suggests that more in-detail kinetic investigations are required for a proper evaluation of H₂S chemical processing in volcanic plumes

Quantification of the depletion of ozone in the plume of Mount Etna

Volcanoes are an important source of inorganic halogen species into the atmosphere. Chemical processing of these species generates oxidised, highly reactive, halogen species which catalyse considerable O3 destruction within volcanic plumes. A campaign of ground-based in situ O3, SO2 and meteorology measurements was undertaken at the summit of Mount Etna volcano in July/August 2012. At the same time, spectroscopic measurements were made of BrO and SO2 columns in the plume downwind. Depletions of ozone were seen at all in-plume measurement locations, with average O3 depletions ranging from 11–35 nmol mol 1 (15–45 %). Atmospheric processing times of the plume were estimated to be between 1 and 4 min. A 1-D numerical model of early plume evolution was also used. It was found that in the early plume O3 was destroyed at an approximately constant rate relative to an inert plume tracer. This is ascribed to reactive halogen chemistry, and the data suggests the majority of the reactive halogen that destroys O3 in the early plume is generated within the crater, including a substantial proportion generated in a high-temperature “effective source region” immediately after emission. The model could approximately reproduce the main measured features of the ozone chemistry. Model results show a strong dependence of the near-vent bromine chemistry on the presence or absence of volcanic NOx emissions and suggest that near-vent ozone measurements can be used as a qualitative indicator of NOx emission