Physical Characterization of Frozen Saltwater ... - ACS Publications

May 23, 2018 - cryosphere; heterogeneous chemistry; model; photochemistry; polar boundary layer; quasi-liquid layer; snowpack. The Supporting Informat...
0 downloads 0 Views 1001KB Size
Subscriber access provided by Kaohsiung Medical University

Article

Physical Characterization of Frozen Saltwater Solutions using Raman Microscopy Philip Patrick Anthony Malley, Subha Chakraborty, and Tara F Kahan ACS Earth Space Chem., Just Accepted Manuscript • DOI: 10.1021/ acsearthspacechem.8b00045 • Publication Date (Web): 23 May 2018 Downloaded from http://pubs.acs.org on June 3, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

Physical Characterization of Frozen Saltwater Solutions using Raman Microscopy Philip P. A. Malley, Subha Chakraborty, and Tara F. Kahan*

1-014 Center for Science and Technology, Syracuse University 111 College Place Syracuse, NY 13244 *Author to whom correspondence should be addressed: [email protected], (315) 443-3285. Fax: (315) 443-4070. 1-014 Center for Science and Technology, Syracuse University, Syracuse, NY 13244.

1 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 26

Abstract: Ice is an important but poorly understood atmospheric reaction medium. Reactions in ice and at air-ice interfaces are often modelled using rate constants measured in liquid aqueous solution, despite evidence that reactivity in these two media can be very different. This approach may be valid at high ionic strengths (e.g. in sea ice) due to the formation of liquid brine. However, recent experiments indicate uneven solute distribution at ice surfaces, suggesting that liquid water does not completely wet ice surfaces at environmentally-relevant solute concentrations. We have investigated the distribution of liquid solution, solid ice, and solid salt (NaCl.2H2O, “hydrohalite”) at the surface of frozen aqueous sodium chloride (NaCl) solutions and frozen sea water using Raman microscopy. At temperatures above the eutectic (–21.1 °C), the ice surfaces were incompletely wetted except at the highest temperatures (approximately –5 °C). Liquid water at the surface took the form of either isolated patches or channels, depending on salt concentration and sample temperature; liquid fractions ranged from approximately 11% to 85%. 3-Dimensional (“volume”) maps showed similar liquid fractions and channel widths at all depths investigated (up to 100 µm) as at the surface for each sample composition. Below –21.1 °C, no liquid was observed in any sample. Instead, hydrohalite was observed with surface coverages ranging from 13% to 100% depending on salt concentration; surface coverage was independent of temperature between –30 °C and –22 °C. Accounting for the presence of two distinct reaction environments at the surface of salty ice might improve predictions of physical and chemical processes in snow-covered regions.

Key words: heterogeneous chemistry, polar boundary layer, cryosphere, snow pack, Sargasso Sea, photochemistry, quasi-liquid layer, model

2 ACS Paragon Plus Environment

Page 3 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

Introduction Ice has long been recognized as an important reaction medium in diverse fields including food preservation, astrochemistry, and stratospheric ozone depletion. In the mid-1990s, researchers discovered that reactions within terrestrial snowpacks could greatly alter the composition of the atmosphere.1-3 Atmospheric models have consequently incorporated reactions in snowpacks and at air-ice interfaces to improve predictions of atmospheric composition in snow-covered regions. In these models, ice is often treated as an inert substance coated by a layer of liquid water; all snowpack chemistry occurs within this liquid phase (e.g. References 4-5). This simplification allows researchers to use rate constants measured in aqueous solutions to model reactions in snow and ice. While this approach significantly reduces the complexity of modelling snowpack chemistry, its validity has been questioned. Specifically, the assumption that reactivity in snow (or at air-ice interfaces) is always well-described by kinetics measured in aqueous solution is not strongly supported.6-9 Many physical and chemical processes occur very differently at liquid and frozen water surfaces. For example, a number of aromatic pollutants have been reported to photolyze more rapidly at ice surfaces than at liquid water surfaces, and some aromatic pollutants that do not photolyze in aqueous solution can photolyze at ice surfaces.10-13 Heterogeneous reaction kinetics have also been reported to differ at frozen and liquid water surfaces: Ozonation of several aromatic species has been reported to be much faster at ice surfaces than at liquid water surfaces, while hydroxyl radicals (OH), which react rapidly with aromatic species at liquid water surfaces, have been reported to be unreactive toward a range of aromatic species at ice surfaces.10, 14-17 If solutes cause the surface to be coated in a liquid “brine,” then these physical and chemical processes will occur as though the surface is a liquid rather than a solid, and can be modelled as

3 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 26

such. There is some evidence that this might be the case: The photolysis rate constant of the aromatic dye harmine was larger at frozen freshwater surfaces than in aqueous solution, but no enhancement compared to in aqueous solution was observed at the surface of frozen NaCl solutions when initial NaCl concentrations exceeded 0.2 M). These results were interpreted as indicating that the reaction environment experienced by harmine transitioned from that of an ice surface to a liquid brine as the NaCl concentration increased.18 While the fractional volume of liquid water in ice can be predicted based on solute concentration and temperature, this does not provide enough information to determine the composition of salty ice surfaces.19-20 It has been suggested that ice surfaces may not be completely wetted even at high solute concentrations.7 Incomplete wetting would result in an ice surface consisting of discrete regions of bare ice and liquid brine. In support of this, several experimental studies have reported inhomogeneous solute distribution at salty ice surfaces.21-24 Scanning electron microscopy (SEM) studies of sintered frozen NaCl solutions showed isolated liquid patches at ice surfaces and at grain boundaries above the eutectic temperature, and “webs” of solid salt below the eutectic temperature, while environmental SEM was recently used to image evaporating laboratory-prepared frost flowers.22, 25 In the latter work, liquid brine was inferred to coat thin (~2 µm) fingers of ice at high temperatures, and salt (likely NaCl or NaCl•2H2O) crystals were observed even at temperatures as high as –6 °C, which is above the eutectic temperature of the NaCl-water system (–21.1 °C).26 Another environmental SEM study showed that a uranyl salt was distributed in channels at ice surfaces, likely in liquid regions at grain boundaries.21 Atomic force microscopy (AFM) and optical images indicate the presence of isolated “clumps” or islands of solid salt surrounded by bare ice below the eutectic temperature.24, 27 X-ray fluorescence (XRF) has recently been used to investigate the spatial

4 ACS Paragon Plus Environment

Page 5 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

distribution of solutes in frozen samples. Significant spatial heterogeneity of several elements including Br and Cl was reported in sea ice cores at depths ranging from the top centimeter to almost 2 m.28 Another study reported that transition metals were distributed in channels when they were frozen from aqueous solutions containing NaCl.23 If surface wetting of frozen NaCl solutions is incomplete, as these studies suggest, at least two distinct reaction environments would be present at ice surfaces – pure ice and a liquid brine. As discussed above, reactivity can differ greatly at ice and liquid water surfaces, and atmospheric models might be improved by accounting for multiple distinct reaction environments at salty ice surfaces. Raman microscopy has become a popular technique to study the spatially-resolved composition of atmospheric condensed phases. Unlike most traditional surface science techniques, it can operate at atmospheric pressure, and can be used with liquid and solid samples, including high-vapour pressure compounds. It is less likely to cause beam damage than other techniques such as electron microscopy and x-ray microscopy. Raman microscopy has been used to study the composition of laboratory-generated and ambient particulate matter, and recently it has been used to investigate the composition of micro-inclusions in sea ice and to investigate the speciation of nitrate in and at the surface of frozen aqueous solutions.29-31 In addition to detecting and identifying solutes in ice and at ice surfaces, Raman spectroscopy can be used to distinguish between liquid water and ice, as discussed in Reference 32. In this work we use Raman microscopy to detect liquid water (or brine) in frozen aqueous NaCl solutions at a range of concentrations and temperatures relevant to Earth’s surface. We determine the fraction of liquid at the ice surface and in the bulk sample, as well as the spatial distribution of the liquid. These results will help inform future atmospheric models that include ice as a reaction medium.

5 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 26

Methods Materials: Sodium chloride (Sigma Aldrich ≥99.5%) solutions were prepared by dissolving a measured quantity of NaCl in 18.2 MΩ cm deionized water to a concentration of either 0.6 M or 0.02 M. Frozen sea ice was prepared from Sargasso Sea water. The seawater was collected in 2012, was filtered to remove particulate matter, and was refrigerated until use.33 Sample Preparation and Raman Imaging: A custom-built cold chamber was used for all experiments. The chamber consisted of a stainless steel plate and a tightly-fitting stainless steel cover with a quartz window in the ceiling. A Peltier device attached to the bottom of the chamber controlled the temperature. A drop of solution (approximately 100 µL) was placed on the chamber floor, and a thermocouple was inserted into the sample through a port on the side of the chamber. All samples were initially frozen at a temperature between –25 and –30 °C, then slowly (0.5 °C min-1) warmed to the desired temperature. This was done to ensure consistent results, as different sample preparation methods (rapid sample warming or freezing the sample directly to the desired temperature) resulted in inconsistent liquid fractional surface coverage. Frozen samples were allowed to equilibrate at the set temperature for 10 minutes prior to acquiring spectra. Raman spectra were acquired with a Renishaw InVia Raman microscope using a 532 nm continuous wave laser with a 50× long working distance objective. The microscope was optically focused at the sample surface prior to each experiment, and an optical image was acquired. Scattered light was collected through the microscope objective, then passed through a 1800 lines/mm visible range grating. The spectrum was imaged onto a 1015 pixel CCD camera. Spectra were acquired between 2502 and 3824 cm-1. The depth of focus in our experiments is ~4 µm and the measurement diameter is ~2 µm. Several experiments were performed to ensure that

6 ACS Paragon Plus Environment

Page 7 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

exposure to the laser beam did not alter the physical properties of the ice surface. These are discussed in the Supplementary Information (SI). 150 µm × 150 µm maps were acquired with 3 µm steps along the x and y axes with an acquisition time of 0.25 s per spectrum. Collection of a single map took ~41 minutes. Threedimensional (“volume”) maps were obtained using the same parameters as the surface area maps but were repeated in 10 µm steps along the z-axis from the ice surface (z = 0 µm) to a depth of 100 µm. Since these step sizes are larger than the depth of focus and measurement diameter of our system, we do not expect adjacent pixels to contribute to spectra acquired at a given location. Figure 1a shows spectra of deionized water at room temperature, at –17 °C, and at –26 °C. Distinct differences are observed in the spectra of liquid and solid samples. Stronger hydrogen bonding between water molecules in ice shifts the Raman peak to lower energies (~3135 cm-1, compared to ~3426 cm-1 in liquid water). There is little difference in the spectra acquired at the two temperatures in the frozen samples. Figure 1b shows spectra of a 0.6 M aqueous NaCl solution acquired at similar temperatures. While the room temperature spectrum looks very similar to that of deionized water, the spectra acquired in ice have very different appearances. At –14.7 °C, the spectrum of the NaCl solution resembles that of liquid water much more so than ice, with peak intensity at 3425 cm-1, likely due to the formation of a liquid brine. We note that spectra of frozen aqueous NaCl solutions never perfectly matched those of liquid water; a small “ice-like” peak at 3135 cm-1 was always observed. This feature is likely due to contributions from nearby ice within the observational volume described above. The spectrum in Figure 1b acquired at –25 °C resembles that of ice, but two additional peaks are observed at 3421 cm-1 and 3539 cm-1. These peaks are due to the formation of NaCl•2H2O (“hydrohalite”), which

7 ACS Paragon Plus Environment

ACS Earth and Space Chemistry

is thermodynamically stable at temperatures below the eutectic point of the NaCl-H2O system (–

Normalized intensity

21.1 °C).26

Normalized intensity

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 26

(a)

(b)

Figure 1. Raman spectra of (a) deionized water acquired at 25 °C, −16.7 °C, and –25.6 °C; and (b) an aqueous 0.6 M NaCl solution acquired at 24.6 °C, −14.7 °C, and −24.6 °C. The intensity of each spectrum is normalized to that at 3136 cm−1.

Spectra were analyzed in two ways. In the first (“Method 1”), we examined covariance of each spectrum with a reference spectrum of pure ice and normalized the covariance by the product of the total counts in the entire sample spectrum and the reference spectrum. This normalized covariance, ρ, represents the spectral deviation from pure ice and hence provides qualitative information about the presence of any species (including liquid water) other than pure ice at a particular coordinate. We plot surface maps of (1 – ρ) over the entire scan-area after scaling between 0 and 1 such that (1 – ρ) = 0 and (1 – ρ) = 1 on the scan-area would correspond to minimum and maximum spectral deviation from pure ice. A description of the calculations is provided in the SI. In the second method (“Method 2”), each spectrum was assigned as liquid 8 ACS Paragon Plus Environment

Page 9 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

water (or a liquid brine), ice, or hydrohalite based on the relative intensities at energies corresponding to peaks in the spectra of liquid water, ice, and hydrohalite. The criteria used for this analysis are shown in the SI. This method does not provide information about the relative contribution of each potential component to the spectra; it classifies each spectrum based on the dominant component. Fractional liquid surface coverage was determined as the ratio of “water-like” pixels to the total number of pixels in a map. For maps generated using Method 2, this entailed dividing the number of spectra classified as water (i.e. teal pixels) by the total number of pixels in a map. For maps generated using Method 1, a threshold value is required to distinguish between “waterlike” and “ice-like” spectra. The threshold value chosen can greatly affect the calculated liquid fraction. To choose the threshold value, we calculated the expected liquid volume fraction of frozen aqueous 0.6 M NaCl solutions at -15 ºC and -8 ºC using Equation 1.19 Φ௟ (ܶ) =

[ே௔஼௟]എ

[1]

[ே௔஼௟]೅

where Φl(T) is the fractional liquid coverage at a given temperature, [NaCl]o is the concentration of NaCl in liquid solution, and [NaCl]T is the NaCl concentration in the liquid fraction of the ice at a given temperature. [NaCl]T, which is greater than [NaCl]o due to freeze exclusion, is determined from Equation 2.19 [ܰܽ‫= ் ]݈ܥ‬

ଵ଴଴଴୼ு೑എ

ெಹమ ೀ ோ்೑



்೑ ି் ்



[2]

Here, ΔHfo is the enthalpy of fusion, MH2O is the molar mass of water, R is the gas constant, and Tf is the temperature of fusion. We calculated Φl in bulk ice at -15 ºC and -8 ºC at a range of threshold values, and chose the threshold value that yielded liquid fractions closest to those predicted by Equation 1 at those temperatures. The threshold value used in this work was 0.85,

9 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 26

meaning that spectra that deviated by at least 15% from a spectrum of pure ice were considered to be liquid water. Liquid fractions determined from maps created using Method 2 consistently yielded greater values of Φl than those determined from maps created using Method 1, both at the air-ice interface and in bulk ice. Method 1 was therefore used for most analysis. Method 2 was used for experiments in which maps were acquired below the eutectic temperature, as Method 2 can distinguish between liquid water, ice, and hydrohalite. A scaling factor of 0.85 was applied to the values of Φl determined from these experiments to bring them into agreement with the values determined using Method 1. The widths of liquid channels at air-ice interfaces and within bulk ice were estimated by measuring the width approximately perpendicular to the channel direction at several locations along the length of the channel.

Results and Discussion Surface Maps of Salt Solutions Above the Eutectic Temperature Figure 2 shows Raman maps acquired at –15 °C at the surface of frozen 0.02 M and 0.6 M NaCl solutions, which correspond to NaCl concentrations expected in snow contaminated with road salt and in sea ice, respectively.34 Figure 2 also shows a map of a frozen seawater sample at the same temperature. Distinct regions of solid ice and liquid solution are apparent at the surface of each sample. No liquid was observed at the surface of frozen deionized water in the absence of salt, as shown in the SI. The deviations from ice observed in Figure 2 are due to the formation of a liquid brine, and not to the so-called “disordered surface region” or “quasiliquid layer” that exists at the top few layers of ice.6 Unlike techniques such as sum-frequency

10 ACS Paragon Plus Environment

Page 11 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

vibrational spectroscopy or glancing-angle Raman spectroscopy, our method is not sensitive or selective to this region.32, 35

(a)

(b)

(a) (b) (c)

(c)

11 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 26

Figure 2. Correlation map of spectra acquired by raster-scan over the surface of a (a) 0.02 M NaCl solution at −13.7 °C; (b) 0.6 M NaCl solution at −13. 5°C; and (c) seawater at −13.5 °C. The scan-size is 150 µm × 150 µm, with a 3 µm step size. The normalized covariance ρ of the spectra at each coordinate of the samples and that of DI water at the coordinate (0, 0) is evaluated and (1 – ρ) is plotted as a surface-plot as shown in the color-bar. The value of (1 – ρ) indicates the dissimilarity between the acquired spectra and the reference spectrum corresponding to DI water at – 18 °C; increasing brightness corresponds to increasing liquid-like nature of the surface.

The deviation from ice is generally greatest in the middle of the liquid regions, whether they are channels or “islands.” This may indicate that liquid regions are deepest in the centre of the feature, as reported by a recent XRF study.23 The deviation from ice is much greater in the 0.6 M solution and the seawater than in the 0.02 M solution, perhaps suggesting deeper liquid channels at the higher salt concentrations. The maps shown in Figure 2 are also shown in the SI with binary analysis (i.e. classifying each spectrum as either liquid water or ice). This analysis was used to calculate liquid surface fraction for each experimental condition. Less liquid is observed at the surface of 0.02 M NaCl solutions than at the surface of 0.6 M solutions; average fractional liquid surface coverages at 15 °C were 13 ± 3% and 19 ± 6% respectively. The fractional liquid surface coverage of frozen seawater was 22 ± 4%, which is similar to that measured at the surface of frozen 0.6 M NaCl solutions at the same temperature. This suggests that NaCl is responsible for the majority of surface melting observed in sea ice, despite the presence of other solutes. These results are in agreement with previous studies that suggest incomplete wetting of ice surfaces by NaCl, but this is the first study to show this at seawater NaCl concentrations.22-24, 27-28 Figure 3 shows the fraction of the ice surface covered by liquid as a function of temperature for 0.02 and 0.6 M NaCl solutions. Maps acquired at each temperature for an individual sample are also shown. As the temperature increases, isolated patches of liquid extend and merge to form channel-like structures in the 0.6 M samples. Liquid at the surface of frozen 0.02 M NaCl solutions generally took the form of patches or islands rather than channels at all

12 ACS Paragon Plus Environment

Page 13 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

temperatures. Our results are consistent with suggestions that solute-induced melting occurs at grain boundaries between ice crystals.21, 23 A positive deviation from the liquid volume fractions predicted by Equations 1 and 2 is observed for the 0.02 M solutions and for the warmest 0.6 M solutions. We hesitate to ascribe this to a greater liquid fraction at the air-ice interface than predicted by thermodynamics, however. First, liquid fractions measured within the bulk agreed with those at the surface for a given sample (see below). Second, significant uncertainty is associated with choosing a threshold value to distinguish between liquid water and ice. We therefore focus on the relative amounts of liquid at the air-ice interface under different conditions rather than on the exact liquid fractions calculated. -20.3 oC 0.6 M

-15.5 oC 0.6 M

-10.2 oC 0.6 M

-15.1 oC 0.02 M

-5.1 oC 0.6 M

-4.9 oC 0.02 M

13 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 26

Figure 3. Temperature dependence of the percent coverage of liquid at air-ice interfaces. Raman maps acquired at temperatures corresponding to the data points in the Figure are also shown. The 4 upper maps are from a 0.6 M NaCl solution and the 2 lower maps are from a 0.02 M NaCl solution. The maps are aligned along the x-axis with the corresponding data point on the graph. Maps were analyzed using Method 1. Black indicates ice and teal indicates liquid water.

Figure 4 shows the average width of liquid channels observed at the surface of frozen aqueous 0.6 M NaCl solutions. Over the span of 10 ºC (from -20.3 to -10.2 ºC), the average width increased by approximately a factor of 3. This is similar to the 2.5× increase in the fractional liquid surface coverage observed over this temperature range.

30 20 10 0 -22

-20

-18

-16

-14

-12

-10

-8

-6

-4

-2

Temperature (oC) Figure 4. Average channel width at the surface of frozen aqueous 0.6 M NaCl solutions. Error bars indicate the standard deviation about the mean of at least 3 trials.

Volume Maps We acquired volume (3-dimensional) maps of frozen 0.02 M and 0.6 M aqueous NaCl solutions and of seawater. Figure 5 shows maps of each sample type at temperatures of approximately -15.5 °C. Liquid regions were observed at all depths. We calculated the liquid fractions of frozen 0.6 M NaCl solutions at depths of –40 and –80 µm; they were 16 ± 6% and 23 ± 8%, which are in good agreement with the fractional coverage of 19 ± 6% measured at the surface. Channel widths at the two depths were 10 ± 9 µm and 11 ± 6 µm, which are also in excellent agreement with the width of 8 ± 4 µm measured at the surface at the same temperature. 14 ACS Paragon Plus Environment

Page 15 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

This demonstrates that not all solutes are excluded to the ice surface, in agreement with recent reports that the concentrations of solutes such as Cl- and NO3- measured at ice surfaces are lower than expected if all of the solute was excluded to the ice surface.30-31, 36-37 Recent x-ray tomography studies have reported liquid volume fractions below 5% in sea ice samples at –20 °C.38 We did not acquire volume maps at –20 °C, but liquid fractional surface coverage (which appears to be similar to bulk liquid volume fraction) was 12 ± 7% for 0.6 M NaCl solutions at – 20.3 °C. This discrepancy could be due to different histories and formation processes of the ice samples, which have been reported (and which we have observed) to affect liquid volume fraction.39 It could also be due to the choice of threshold used to distinguish between water and ice in the two studies. The thermodynamically-predicted liquid volume fraction at this temperature (based on Equation 1) is 12%, in good agreement with our experimental value.

15 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 26

Figure 5. Volume maps of (a) a 0.02 M aqueous NaCl solution at –15.5 °C, (b) a 0.6 M aqueous NaCl solution at – 15.4 °C, and (c) seawater at –15.6 °C. Maps were analyzed using Method 1. Teal indicates liquid water.

X-ray tomography studies of frozen sea ice show vertical liquid channels extending to depths of several mm.38, 40 Conversely, a different x-ray tomography study showed that cesium chloride in frozen aqueous solutions was concentrated in distinct unconnected pockets rather than in channels.39 Our results are consistent with both studies. We observed liquid channels in frozen 0.6 M aqueous NaCl solutions and in frozen seawater (some of which extended at least 100 µm in depth), while liquid water in frozen 0.02 M NaCl solutions was primarily in isolated 16 ACS Paragon Plus Environment

Page 17 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

patches. The study investigating cesium chloride distributions used 1.0 mM concentrations.39 At this low concentration, there is likely not enough liquid present to form the channels that are observed at seawater NaCl concentrations. Finally, we note that under our experimental conditions less than 1% of the volume of frozen deionized water was in liquid form (as shown in the SI), which indicates that the liquid observed in the volume maps presented in Figure 5 is due to solute-induced melting.

Surface Maps of Frozen Salt Solutions Below the Eutectic Temperature At temperatures below the eutectic (–21.1 °C), NaCl–H2O mixtures at equilibrium will consist of a combination of solid ice and solid salt (including the hydrohalite, NaCl•2H2O, which we are able to detect using Raman spectroscopy as shown in Figure 1b).26 Liquid water (or brine) is not expected to be present, although an NMR study reported the presence of liquid water at temperatures as low as –45 °C in frozen 0.5 M NaCl solutions.19 Figure 6 shows surface maps of frozen NaCl solutions and seawater at approximately – 31°C. At this temperature, we do not see any evidence of liquid (which would be represented as teal pixels on the maps if it was detected). Instead, we see distinct regions of pristine ice and hydrohalite. The surface coverage of hydrohalite varied significantly between samples and ranged from 40 to 100% for the 0.6 M NaCl solutions and from 45 to 80% for the seawater. Coverage was much lower for the 0.02 M NaCl solutions, ranging from 13 to 20%. At temperatures lower than the eutectic, the ice is opaque due to the presence of hydrohalite, and volume maps could not be obtained.

17 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 26

Figure 6. Surface maps of (a) 0.6 M NaCl (–32.0 °C), (b) Sargasso Sea Water (–29.5 °C), and (c) 0.02 M NaCl (– 30.9 °C). Maps were analyzed using Method 2. Black indicates ice, teal indicates liquid water (not observed), and magenta indicates hydrohalite.

Surface Maps of Frozen Salt Solutions Near the Eutectic Temperature Figure 7 shows a surface map of a 0.6 M NaCl solution acquired just below the eutectic temperature, as well as a map of the same sample region acquired after warming the sample to just above the eutectic temperature. The transition from solid hydrohalite to liquid brine when the temperature rises above the eutectic temperature is evident. At all temperatures investigated, hydrohalite and liquid were never observed simultaneously in the same sample.

18 ACS Paragon Plus Environment

Page 19 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

Figure 7. Raman maps of a 0.6 M NaCl solution at –21.8 °C (top) and –18.7 °C (bottom). Maps were analyzed using Method 2. Black indicates ice, teal indicates liquid water, and magenta indicates hydrohalite.

We performed a series of experiments to investigate the transition between liquid brine and hydrohalite near the eutectic temperature using 0.6 M NaCl solutions. For these experiments, a Raman map was acquired at a temperature below the eutectic, and then the temperature was slowly raised by a few degrees. After equilibrating the sample at the new temperature, another map was acquired, and this process was repeated until several maps were acquired at temperatures ranging from approximately –32 °C to –15 °C. Then maps were acquired as the temperature was slowly decreased. Figure 8a shows the fraction of the ice surface containing hydrohalite at each temperature, and Figure 8b shows the fractional liquid surface coverage. The hydrohalite coverage remained constant at temperatures between –32 and –21.7 °C. As the temperature increased above –21.1 °C, the extent of hydrohalite coverage decreased rapidly. Negligible hydrohalite was detected at temperatures above the eutectic. As we cooled the sample, we observed a reproducible hysteresis: hydrohalite was not observed until the temperature decreased below –22 °C – a full degree lower than the eutectic temperature. The 19 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 26

same hysteresis reported for hydrohalite formation is observed for fractional liquid surface coverage (Fig. 8b). Whereas the surface coverage of hydrohalite was constant at temperatures below the eutectic, at temperatures greater than the eutectic the liquid fraction increased with increasing temperature.

Figure 8. Average fraction of the surface covered with (a) hydrohalite and (b) liquid water as a function of sample temperature. Red triangles indicate increasing temperature, and blue circles indicate decreasing temperature. Error bars represent the standard deviation about the mean of at least 3 trials. The red and blue arrows indicate the direction of the temperature ramp (warming or cooling). Fractional surface coverages were determined using Method 2.

The observed hysteresis is likely due to the fact that hydrohalite formation requires the proper conformation of salt and water molecules. Similar hystereses are reported for several atmospheric phase changes such as the efflorescence of aqueous aerosols and the homogeneous

20 ACS Paragon Plus Environment

Page 21 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

freezing of water.41-43 This hysteresis was observed whether the solution was cooled slowly between map acquisition temperatures, or at a rapid rate (0.5 °C min-1 or 0.5 °C s-1).

Atmospheric Implications The results of this work have two major implications for the way air-ice interactions are treated in atmospheric models. The first is that we observe an even distribution of liquid water throughout the entire ice sample, whereas models generally position all of the water (and solutes) at the air-ice interface. Only solutes in liquid regions at the air-ice interface (or connected to it) can participate in heterogeneous reactions. The volume of the liquid brine at the air-ice interface, and the number of solute molecules in this brine, could be adjusted in atmospheric models. The second implication of our results is that ice surfaces will not be completely coated in liquid water, even at high solute concentrations and high temperatures. Some fraction of the surface will always consist of bare ice. This spatial heterogeneity may affect gas-surface interactions, as two distinct environments – pure ice and liquid brine (or pure ice and solid hydrohalite at low temperatures) – will be present. Physical and reactive uptake coefficients of gases to these two environments, as well as reaction kinetics of adsorbed / absorbed species, may be very different. Altering the volume of liquid accessible to the air in atmospheric models is relatively straightforward. Accounting for the effects this will have on reactivity is not. For one thing, the distribution of solutes between different liquid regions may not be uniform. For example, nitrate was recently shown to reside in liquid compartments within frozen aqueous solutions containing seawater levels of NaCl.30 Nitrate’s surface concentration was ~20% greater than that in the bulk in frozen solutions prepared from 0.5 M NaCl, but was 80% greater in solutions prepared from an aquarium sea salt mixture. This suggests that the concentration of reactive species at air-ice

21 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 26

interfaces may depend on the identity of ionic species present, and not solely on the extent of brine formation. It is also possible that some solutes will not reside within liquid regions in the ice. X-ray fluorescence mapping of ice core samples from Antarctica suggested that Cl, K, and Ca were located primarily within grain boundaries, but Br and Sr had different lateral distributions at all ice core depths investigated.28 Finally, understanding that two distinct reaction environments (bare ice and liquid brine) exist at air-ice interfaces under most environmental conditions may be helpful for interpreting and predicting heterogeneous reaction kinetics in snow-covered regions, but more information is necessary to incorporate this into models in a meaningful way. Specifically, the distribution of reactive solutes between these two media should be known, as should the physical and reactive uptake coefficients of gas-phase species. Further work is required to better understand reactivity on salty ice surfaces and to improve predictions of chemistry in snow-covered regions.

Acknowledgments: This work was funded by NSF award 1454959. The authors thank Dr. David J. Kieber for providing Sargasso Sea Water for the experiments.

Supporting Information: Quality control experiments, description of spectral analysis methods, surface and volume maps of deionized water.

22 ACS Paragon Plus Environment

Page 23 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

References 1. Domine, F.; Shepson, P. B., Air-Snow Interactions and Atmospheric Chemistry. Science 2002, 297, 1506-1510. 2. Grannas, A. M.; Jones, A. E.; Dibb, J.; Ammann, M.; Anastasio, C.; Beine, H. J.; Bergin, M.; Bottenheim, J.; Boxe, C. S.; Carver, G.; Chen, G.; Crawford, J. H.; Domine, F.; Frey, M. M.; Guzman, M. I.; Heard, D. E.; Helmig, D.; Hoffmann, M. R.; Honrath, R. E.; Huey, L. G.; Hutterli, M.; Jacobi, H. W.; Klan, P.; Lefer, B.; McConnell, J.; Plane, J.; Sander, R.; Savarino, J.; Shepson, P. B.; Simpson, W. R.; Sodeau, J. R.; von Glasow, R.; Weller, R.; Wolff, E. W.; Zhu, T., An Overview of Snow Photochemistry: Evidence, Mechanisms and Impacts. Atmos. Chem. Phys. 2007, 7, 4329-4373. 3. Simpson, W. R.; von Glasow, R.; Riedel, K.; Anderson, P.; Ariya, P.; Bottenheim, J.; Burrows, J.; Carpenter, L.; Friess, U.; Goodsite, M. E.; Heard, D.; Hutterli, M.; Jacobi, H.-W.; Kaleschke, L.; Neff, B.; Plane, J.; Platt, U.; Richter, A.; Roscoe, H.; Sander, R.; Shepson, P.; Sodeau, J.; Steffen, A.; Wagner, T.; Wolff, E., Halogens and Their Role in Polar BoundaryLayer Ozone Depletion. Atmos. Chem. Phys. Discuss. 2007, 7, 4285-4403. 4. Toyota, K.; McConnell, J. C.; Staebler, R. M.; Dastoor, A. P., Air-Snowpack Exchange of Bromine, Ozone and Mercury in the Springtime Arctic Simulated by the 1-D Model Phantas Part 1: In-Snow Bromine Activation and Its Impact on Ozone. Atmos. Chem. Phys. 2014, 14, 4101-4133. 5. Thomas, J. L.; Stutz, J.; Lefer, B.; Huey, L. G.; Toyota, K.; Dibb, J. E.; von Glasow, R., Modeling Chemistry in and above Snow at Summit, Greenland - Part 1: Model Description and Results. Atmos. Chem. Phys. 2011, 11, 4899-4914. 6. Bartels-Rausch, T.; Jacobi, H. W.; Kahan, T. F.; Thomas, J. L.; Thomson, E. S.; Abbatt, J. P. D.; Ammann, M.; Blackford, J. R.; Bluhm, H.; Boxe, C.; Domine, F.; Frey, M. M.; Gladich, I.; Guzman, M. I.; Heger, D.; Huthwelker, T.; Klan, P.; Kuhs, W. F.; Kuo, M. H.; Maus, S.; Moussa, S. G.; McNeill, V. F.; Newberg, J. T.; Pettersson, J. B. C.; Roeselova, M.; Sodeau, J. R., A Review of Air-Ice Chemical and Physical Interactions (Aici): Liquids, Quasi-Liquids, and Solids in Snow. Atmos. Chem. Phys. 2014, 14, 1587-1633. 7. Domine, F.; Bock, J.; Voisin, D.; Donaldson, D. J., Can We Model Snow Photochemistry? Problems with the Current Approaches. J. Phys. Chem. A 2013, 117, 47334749. 8. Kahan, T. F.; Malley, P. P. A.; Grossman, J. N.; Stathis, A. A., Photochemistry in Model Aqueous-Organic Condensed Phases. In Multiphase Atmospheric Chemistry in the Environment, Hunt, S.; Nizkorodov, S. A.; Laskin, A., Eds. ACS Symposium Series, accepted. 9. Kahan, T. F.; Wren, S. N.; Donaldson, D. J., A Pinch of Salt Is All It Takes: Chemistry at the Frozen Water Surface. Accounts Chem. Res. 2014, 47, 1587-1594. 10. Kahan, T. F.; Donaldson, D. J., Photolysis of Polycyclic Aromatic Hydrocarbons on Water and Ice Surfaces. J. Phys. Chem. A 2007, 111, 1277-1285. 11. Kahan, T. F.; Donaldson, D. J., Benzene Photolysis on Ice: Implications for the Fate of Organic Contaminants in the Winter. Environ. Sci. Technol. 2010, 44, 3819-3824. 12. Malley, P. P. A.; Kahan, T. F., Nonchromophoric Organic Matter Suppresses Polycyclic Aromatic Hydrocarbon Photolysis in Ice and at Ice Surfaces. J. Phys. Chem. A 2014, 118, 16381643. 13. Stathis, A. A.; Hendrickson-Stives, A. K.; Kahan, T. F., Photolysis Kinetics of Toluene, Ethylbenzene, and Xylenes at Ice Surfaces. J. Phys. Chem. A 2016, 120, 6693-6697. 23 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 26

14. Kahan, T. F.; Donaldson, D. J., Heterogeneous Ozonation Kinetics of Phenanthrene at the Air-Ice Interface. Environ. Res. Lett. 2008, 3, 045006. 15. Kahan, T. F.; Zhao, R.; Donaldson, D. J., Hydroxyl Radical Reactivity at the Air-Ice Interface. Atmos. Chem. Phys. 2010, 10, 843-854. 16. Ray, D.; Liskova, H.; Klan, P., Kinetics of Heterogeneous Reactions of Ozone with Representative Pahs and an Alkene at the Air-Ice Interface at 258 and 188 K. Environ. Sci. Proc. Imp. 2014, 16, 770-776. 17. Ray, D.; Malongwe, J. K.; Klan, P., Rate Acceleration of the Heterogeneous Reaction of Ozone with a Model Alkene at the Air-Ice Interface at Low Temperatures. Environ. Sci. Technol. 2013, 47, 6773-6780. 18. Kahan, T. F.; Kwamena, N.-O. A.; Donaldson, D. J., Different Photolysis Kinetics at the Surface of Frozen Freshwater Vs. Frozen Salt Solutions. Atmos. Chem. Phys. 2010, 10, 1091710922. 19. Cho, H.; Shepson, P. B.; Barrie, L. A.; Cowin, J. P.; Zaveri, R., Nmr Investigation of the Quasi-Brine Layer in Ice/Brine Mixtures. J. Phys. Chem. B 2002, 106, 11226-11232. 20. Kuo, M. H.; Moussa, S. G.; McNeill, V. F., Modeling Interfacial Liquid Layers on Environmental Ices. Atmos. Chem. Phys. 2011, 11, 9971-9982. 21. Krausko, J.; Runstuk, J.; Nedela, V.; Klan, P.; Heger, D., Observation of a Brine Layer on an Ice Surface with an Environmental Scanning Electron Microscope at Higher Pressures and Temperatures. Langmuir 2014, 30, 5441-5447. 22. Blackford, J. R.; Jeffree, C. E.; Noake, D. F. J.; Marmo, B. A., Microstructural Evolution in Sintered Ice Particles Containing Nacl Observed by Low-Temperature Scanning Electron Microscope. Proc. Inst. Mech. Eng. L J. Mater. Des. Appl. 2007, 221, 151-156. 23. Tokumasu, K.; Harada, M.; Okada, T., X-Ray Fluorescence Imaging of Frozen Aqueous Nacl Solutions. Langmuir 2016, 32, 527-533. 24. Doppenschmidt, A.; Butt, H. J., Measuring the Thickness of the Liquid-Like Layer on Ice Surfaces with Atomic Force Microscopy. Langmuir 2000, 16, 6709-6714. 25. Yang, X.; Nedela, V.; Runstuk, J.; Ondruskova, G.; Krausko, J.; Vetrakova, L.; Heger, D., Evaporating Brine from Frost Flowers with Electron Microscopy and Implications for Atmospheric Chemistry and Sea-Salt Aerosol Formation. Atmos. Chem. Phys. 2017, 17, 62916303. 26. Lange, N. A.; Forker, G. M., Lange's Handbook of Chemistry. 9 ed.; Handbook Publishers Inc: 1956. 27. Krepelova, A.; Huthwelker, T.; Bluhm, H.; Ammann, M., Surface Chemical Properties of Eutectic and Frozen Nacl Solutions Probed by Xps and Nexafs. Chemphyschem 2010, 11, 38593866. 28. Obbard, R. W.; Lieb-Lappen, R. M.; Nordick, K. V.; Golden, E. J.; Leonard, J. R.; Lanzirotti, A.; Newville, M. G., Synchrotron X-Ray Fluorescence Spectroscopy of Salts in Natural Sea Ice. Earth and Space Science 2016, 3, 463-479. 29. Eichler, J.; Kleitz, I.; Bayer-Giraldi, M.; Jansen, D.; Kipfstuhl, S.; Shigeyama, W.; Weikusat, C.; Weikusat, I., Location and Distribution of Micro-Inclusions in the Edml and Neem Ice Cores Using Optical Microscopy and in Situ Raman Spectroscopy. Cryosphere 2017, 11, 1075-1090. 30. Morenz, K. J.; Donaldson, D. J., Chemical Morphology of Frozen Mixed Nitrate-Salt Solutions. J. Phys. Chem. A 2017, 121, 2166-2171.

24 ACS Paragon Plus Environment

Page 25 of 26 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Earth and Space Chemistry

31. Wren, S. N.; Donaldson, D. J., Exclusion of Nitrate to the Air-Ice Interface During Freezing. J. Phys. Chem. Lett. 2011, 2, 1967-1971. 32. Kahan, T. F.; Reid, J. P.; Donaldson, D. J., Spectroscopic Probes of the Quasi-Liquid Layer on Ice. J. Phys. Chem. A 2007, 111, 11006-11012. 33. Kieber, D. J.; Keene, W. C.; Frossard, A. A.; Long, M. S.; Maben, J. R.; Russell, L. M.; Kinsey, J. D.; Tyssebotn, I. M. B.; Quinn, P. K.; Bates, T. S., Coupled Ocean-Atmosphere Loss of Marine Refractory Dissolved Organic Carbon. Geophys. Res. Lett. 2016, 43, 2765-2772. 34. Labadia, C. F.; Buttle, J. M., Road Salt Accumulation in Highway Snow Banks and Transport through the Unsaturated Zone of the Oak Ridges Moraine, Southern Ontario. Hydrol. Process. 1996, 10, 1575-1589. 35. Wei, X.; Shen, Y. R., Vibrational Spectroscopy of Ice Interfaces. Applied Phys. B Lasers and Optics 2002, 74, 617-620. 36. Walker, R. L.; Searles, K.; Willard, J. A.; Michelsen, R. R. H., Total Reflection Infrared Spectroscopy of Water-Ice and Frozen Aqueous Nacl Solutions. J. Chem. Phys. 2013, 139, 8. 37. Kahan, T. F.; Wren, S. N.; Donaldson, D. J., A Pinch of Salt Is All It Takes: Chemistry at the Frozen Water Surface. Accounts Chem. Res. 2014, 47, 1587-1594. 38. Lieb-Lappen, R. M.; Golden, E. J.; Obbard, R. W., Metrics for Interpreting the Microstructure of Sea Ice Using X-Ray Micro-Computed Tomography. Cold Reg. Sci. Tech. 2017, 138, 24-35. 39. Hullar, T.; Anastasio, C., Direct Visualization of Solute Locations in Laboratory Ice Samples. Cryosphere 2016, 10, 2057-2068. 40. Lieblappen, R. M.; Kumar, D. D.; Pauls, S. D.; Obbard, R. W., A Network Model for Characterizing Brine Channels in Sea Ice. Cryosphere 2018, 12, 1013-1026. 41. Wang, J.; Hoffmann, A. A.; Park, R. J.; Jacob, D. J.; Martin, S. T., Global Distribution of Solid and Aqueous Sulfate Aerosols: Effect of the Hysteresis of Particle Phase Transitions. J. Geophys. Res. Atmos. 2008, 113, 11. 42. Colberg, C. A.; Krieger, U. K.; Peter, T., Morphological Investigations of Single Levitated H2so4/Nh3/H2o Aerosol Particles During Deliquescence/Efflorescence Experiments. J. Phys. Chem. A 2004, 108, 2700-2709. 43. Morishige, K.; Kawano, K., Freezing and Melting of Water in a Single Cylindrical Pore: The Pore-Size Dependence of Freezing and Melting Behavior. J. Chem. Phys 1999, 110, 48674872.

25 ACS Paragon Plus Environment

ACS Earth and Space Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 26

For TOC only

26 ACS Paragon Plus Environment