Abstract. Two methods were used to quantify the flux of volcanic sulphur (as the equivalent mass of SO2) to the stratosphere over different timescalses during the Holocene. A combination of satellite-based measurements of sulphur yields from recent explosive volcanic eruptions with an appropriate rate of explosive volcanism for the past 200 years constrains the medium-term (~102 years) flux of volcanic sulphur to the stratosphere to be ~1 Mt a-1, with lower and upper bounds of 0.3 and 3 Mt a-1. The short-term (10- to 20-year) flux due to small magnitude (1010-1012 kg) eruptions is of the order of 0.4 Mta-1. At any time the instantaneous levels of sulphur in the stratosphere are dominated by the most recent (0-3 years) volcanic events. The flux calculations do not attempt to address this very short timescalses variability. Although there are significant errors associated with the raw sulphur emission data on which this analysis is based, the approach present is general and may be readily modified as the quantity and quality of the data improve. Data from a Greenland ice core support these conclusions. Integration of the sulphate signals from presumed volcanic sources recorded in the GISP2 core provides a minimum estimate of the 103-year volcanic SO2 flux to the stratosphere of 0.5-1 Mt a-1 over the past 9000 years. The short-term flux calculations do not account for the impact of rare, large events. The ice-core record does not fully account for the contribution from small, frequent events. |
Table 3 Volcanic SO2 injection to the stratosphere
Eruption magnitude (k) |
Time (years) for >95% probablity of an eruptiona |
Number of erption reaching the stratosphere per yearb |
SO2 injection (Mt a-1; bounds in parentheses)c |
Eruptive mass flux of tephra (Mt a-1)d |
3 |
1.0 |
0.72 |
0.08 (0.03-0.3) |
20-30 |
4 |
6.4 |
0.47 |
0.29 (0.09-0.9) |
130-200 |
5 |
40 |
0.08 |
0.26 (0.08-0.8) |
210-320 |
6 |
250 |
0.01 |
0.24 (0.08-0.8) |
340-510 |
Total |
0.9 (.03-3) |
700-1060 |
a From Poisson statistics, with l from Eq. (1)
b From Eq. (1) or Fig. 2
c From Eq. (4), allowing ±0.5 log unit range in the intercept
value; see text
d Calcuted mass flux from Pyle (1995), assuming bulk tephra with
density of 1000-1500 kg m-3