Abstract. Snow in the environment acts as a host to rich chemistry and provides a matrix for physical exchange of contaminants within the ecosystem. The goal of this review is to summarise the current state of knowledge of physical processes and chemical reactivity in surface snow with relevance to polar regions. It focuses on a description of impurities in distinct compartments present in surface snow, such as snow crystals, grain boundaries, crystal surfaces, and liquid parts. It emphasises the microscopic description of the ice surface and its link with the environment. Distinct differences between the disordered air–ice interface, often termed quasi-liquid layer, and a liquid phase are highlighted. The reactivity in these different compartments of surface snow is discussed using many experimental studies, simulations, and selected snow models from the molecular to the macro-scale. Although new experimental techniques have extended our knowledge of the surface properties of ice and their impact on some single reactions and processes, others occurring on, at or within snow grains remain unquantified. The presence of liquid or liquid-like compartments either due to the formation of brine or disorder at surfaces of snow crystals below the freezing point may strongly modify reaction rates. Therefore, future experiments should include a detailed characterisation of the surface properties of the ice matrices. A further point that remains largely unresolved is the distribution of impurities between the different domains of the condensed phase inside the snowpack, i.e. in the bulk solid, in liquid at the surface or trapped in confined pockets within or between grains, or at the surface. While surface-sensitive laboratory techniques may in the future help to resolve this point for equilibrium conditions, additional uncertainty for the environmental snowpack may be caused by the highly dynamic nature of the snowpack due to the fast metamorphism occurring under certain environmental conditions. Due to these gaps in knowledge the first snow chemistry models have attempted to reproduce certain processes like the long-term incorporation of volatile compounds in snow and firn or the release of reactive species from the snowpack. Although so far none of the models offers a coupled approach of physical and chemical processes or a detailed representation of the different compartments, they have successfully been used to reproduce some field experiments. A fully coupled snow chemistry and physics model remains to be developed.

A comparative study of OH, O3, and H2O equilibrium aqueous solvation and gas-phase accommodation on liquid water at 300 K is performed using a combination of ab initio calculations and molecular dynamics simulations. Polarizable force fields are developed for the interaction potential of OH and O3 with water. The free energy profiles for transfer of OH and O3 from the gas phase to the bulk liquid exhibit a pronounced minimum at the surface, but no barrier to solvation in the bulk liquid. The calculated surface excess of each oxidant is comparable to calculated and experimental values for short chain, aliphatic alcohols. Driving forces for the surface activity are discussed in terms of the radial distribution functions and dipole orientation distributions for each molecule in the bulk liquid and at the surface. Simulations of OH, O3, and H2O impinging on liquid water with a thermal impact velocity are used to calculate thermal accommodation (S) and mass accommodation (alpha) coefficients. The values of S for OH, O3, and H2O are 0.95, 0.90, and 0.99, respectively. The approaching molecules are accelerated toward the liquid surface when they are approximately 5 angstroms above it. The molecules that reach thermal equilibrium with the surface do so within 2 ps of striking the surface, while those that do not scatter into the gas phase with excess translational kinetic energy in the direction perpendicular to the surface. The time constants for absorption and desorption range from approximately 35 to 140 ps, and the values of alpha for OH, O3, and H2O are 0.83, 0.047, and 0.99, respectively. The results are consistent with previous formulations of gas-phase accommodation from simulations, in which the process occurs by rapid thermal and structural equilibration followed by diffusion on the free energy profile. The implications of these results with respect to atmospheric chemistry are discussed.

The reaction of ozone with aqueous sodium bromide particles is investigated with a combination of aerosol chamber experiments, kinetics modeling, and molecular dynamics simulations. The molecular bromine production in the chamber experiments is approximately an order of magnitude greater than that predicted by known chemistry in the gas and bulk aqueous phases with use of a comprehensive computer kinetics model. Molecular dynamics simulations indicate that ozone has significant residence time at the air−solution interface, while making frequent contacts with bromide ions for as long as 50 ps in the surface layer of a 6.1 M NaBr solution. The formation of a complex between ozone and bromide ion, [O3···Br-], which can lead to production of Br2 by reaction at the air−water interface, is therefore feasible. Experimentally observed Br2 is well predicted by including an interface process with a reaction probability of [1.9 ± 0.8] × 10-6 (1 s) as the first step in a surface mechanism to produce additional gas-phase Br2. An estimate of the impact of this interface reaction on bromine formation in the marine boundary layer shows that several ppt of bromine could potentially be produced during the night from this proposed surface chemistry.

Interaction of the hydroxyl radical with the liquid water surface was studied using classical molecular dynamics computer simulations. From a series of scattering trajectories, the thermal and mass accommodation coefficients of OH on liquid water at 300 K were determined to be 0.95 and 0.83, respectively. The calculated free energy profile for transfer of OH across the air-water interface at 300 K exhibits a minimum in the interfacial region, with the free energy of adsorbtion (DeltaGa) being about 1 kcal/mol more negative than the hydration free energy (DeltaGs). The propensity of the hydroxyl radical for the air-water interface manifests itself in partitioning of OH radicals between the bulk water and the surface. The enhancement of the surface concentration of OH relative to its concentration in the aqueous phase suggests that important OH chemistry may be occurring in the interfacial layer of water droplets, aqueous aerosol particles, and thin water films adsorbed on solid surfaces. This has profound consequences for modeling heterogeneous atmospheric chemical processes.