Air–snowpack exchange of bromine, ozone and mercury in the springtime Arctic simulated by the 1-D model PHANTAS – Part 2: Mercury and its speciation

Atmospheric mercury depletion events (AMDEs) refer to a recurring depletion of mercury occurring in the springtime Arctic (and Antarctic) boundary layer, in general, concurrently with ozone depletion events (ODEs). To close some of the knowledge gaps in the physical and chemical mechanisms of AMDEs and ODEs, we have developed a one-dimensional model that simulates multiphase chemistry and transport of trace constituents throughout porous snowpack and in the overlying atmospheric boundary layer (ABL). This paper constitutes Part 2 of the study, describing the mercury component of the model and its application to the simulation of AMDEs. Building on model components reported in Part 1 ("In-snow bromine activation and its impact on ozone"), we have developed a chemical mechanism for the redox reactions of mercury in the gas and aqueous phases with temperature dependent reaction rates and equilibrium constants accounted for wherever possible. Thus the model allows us to study the chemical and physical processes taking place during ODEs and AMDEs within a single framework where two-way interactions between the snowpack and the atmosphere are simulated in a detailed, process-oriented manner. Model runs are conducted for meteorological and chemical conditions that represent the springtime Arctic ABL characterized by the presence of "haze" (sulfate aerosols) and the saline snowpack on sea ice. The oxidation of gaseous elemental mercury (GEM) is initiated via reaction with Br-atom to form HgBr, followed by competitions between its thermal decomposition and further reactions to give thermally stable Hg(II) products. To shed light on uncertain kinetics and mechanisms of this multi-step oxidation process, we have tested different combinations of their rate constants based on published laboratory and quantum mechanical studies. For some combinations of the rate constants, the model simulates roughly linear relationships between the gaseous mercury and ozone concentrations as observed during AMDEs/ODEs by including the reaction HgBr + BrO and assuming its rate constant to be the same as for the reaction HgBr + Br, while for other combinations the results are more realistic by neglecting the reaction HgBr + BrO. Speciation of gaseous oxidized mercury (GOM) changes significantly depending on whether or not BrO is assumed to react with HgBr to form Hg(OBr)Br. Similarly to ozone (reported in Part 1), GEM is depleted via bromine radical chemistry more vigorously in the snowp