New estimates of sulfate diffusion rates in the EPICA Dome C ice core

Abstract. To extract climatically relevant chemical signals from the deepest, oldest Antarctic ice, we must first investigate the degree to which chemical ions diffuse within solid ice. Volcanic sulfate peaks are an ideal target for such an investigation because they are high-amplitude, short-duration (∼3 years) events with a quasi-uniform structure. Here we present an analysis of the EPICA Dome C sulfate record over the last 450 kyr. We identify volcanic peaks and isolate them from the non-sea-salt sulfate background to reveal the effects of diffusion: amplitude damping and broadening of peaks in the time domain with increasing depth and age. Sulfate peak shape is also altered by the thinning of ice layers with depth that results from ice flow. Both processes must be simulated to derive effective diffusion rates. This is achieved by running a forward model to diffuse idealised sulfate peaks at different rates while also accounting for ice thinning. Our simulations suggest a median effective diffusion rate of sulfate ions of 2.4±1.7×10-7 m2 yr−1 in Holocene ice, slightly faster than suggested by previous work. The effective diffusion rate observed in deeper ice is significantly lower, and Holocene ice shows the highest rate of the last 450 kyr. Beyond the Holocene, there is no systematic difference between the effective diffusion rates of glacial and interglacial periods despite variations in soluble ion concentrations, dust loading, and ice grain radii. Effective diffusion rates for 40 to 200 ka are relatively constant and of the order 1×10-8 m2 yr−1. Our results suggest that the diffusion of sulfate ions within volcanic peaks is relatively fast initially, perhaps through an interconnected vein network, but slows significantly after 40 kyr. In the absence of clear evidence for a controlling influence of temperature on sulfate diffusivity with depth and age, we hypothesise that the rapid decrease in effective diffusion rate from the time of deposition to ice of 50 ka age may be due to a switch in the mechanism of diffusion resulting from the changing location of sulfate ions within the ice microstructure and/or interconnectedness of veins and grain boundaries.