Modeling the Sulfate Deposition to the Greenland Ice Sheet From the Laki Eruption

Using the state of the art Goddard Institute for Space Studies (GISS) modelE general circulation model, simulations were conducted of the chemistry and transport of aerosols resulting from the 1783-84 Laki (64°N) flood lava eruption. A set of 3 ensemble simulations from different initial conditions were conducted by injecting our estimate of the SO2 gas into the atmosphere by the 10 episodes of the eruption and allowing the sulfur chemistry model to convert this gas into sulfate aerosol. The SO2 gas and sulfate aerosol is transported by the model and wet and dry deposition is calculated over each grid box during the simulation. We compare the resulting sulfate deposition to the Greenland Ice Sheet in the model to 23 ice core measurements and find very good agreement. The model simulation deposits a range of 169 to over 300 kg/km2 over interior Greenland with much higher values along the coastal areas. This compares to a range of 62 to 324 kg/km2 for the 23 ice core measurements with an average value of 158 kg/km2. This comparison is one important model validation tool. Modeling and observations show fairly large spatial variations in the deposition of sulfate across the Greenland Ice Sheet for the Laki eruption, but the patterns are similar to those we modeled for the 1912 Katmai and 1991 Pinatubo eruptions. Estimates of sulfate loading based on single ice cores can show significant differences, so ideally several ice cores should be combined in reconstructing the sulfate loading of past volcanic eruptions, taking into account the characteristic spatial variations in the deposition pattern.