Radiocarbon Age of Soil Organic Matter Fractions Buried by Tephra in Alaska

Radiocarbon ages were determined on different fractions extracted from buried paleosols in south-central Alaska as an experiment to establish best practices for analysis of low-organic-matter paleosols. Seven samples were collected from directly beneath tephra deposits to determine the eruption frequency of Mount Spurr Volcano, Alaska. Soil development near the volcano is poor due to the high-latitude climate and frequent burial of soil surfaces by tephra. Contamination of soils by local wind-blown material is a concern. The humic acid 14C ages are consistently younger than both the bulk soil and residue after extraction ages. The difference in ages between the humic acid extract and bulk soil range from 60–1130 14C yr BP and 180–4110 14C yr BP, respectively, for residue. Previous observations from dating different soil fractions show that residue ages are typically younger than humic acid extracts presumably because they contain a fraction of younger plant material including roots. We attribute the older ages to contamination by old carbon from eolian charcoal particles. This study supports the use of accelerator mass spectrometry (AMS) 14C dating of the humic acid fraction in order to estimate the age of soil that presumably marks the age of burial and avoids suspected contamination by old carbon.

Adjusting particle-size distributions to account for aggregation in tephra-deposit model forecasts

Volcanic ash transport and dispersion (VATD) models are used to forecast tephra deposition during volcanic eruptions. Model accuracy is limited by the fact that fine-ash aggregates (clumps into clusters), thus altering patterns of deposition. In most models this is accounted for by ad hoc changes to model input, representing fine ash as aggregates with density ρagg, and a log-normal size distribution with median μagg and standard deviation σagg. Optimal values may vary between eruptions. To test the variance, we used the Ash3d tephra model to simulate four deposits: 18 May 1980 Mount St. Helens; 16–17 September 1992 Crater Peak (Mount Spurr); 17 June 1996 Ruapehu; and 23 March 2009 Mount Redoubt. In 192 simulations, we systematically varied μagg and σagg, holding ρagg constant at 600 kg m−3. We evaluated the fit using three indices that compare modeled versus measured (1) mass load at sample locations; (2) mass load versus distance along the dispersal axis; and (3) isomass area. For all deposits, under these inputs, the best-fit value of μagg ranged narrowly between ∼ 2.3 and 2.7φ (0.20–0.15 mm), despite large variations in erupted mass (0.25–50 Tg), plume height (8.5–25 km), mass fraction of fine ( < 0.063 mm) ash (3–59 %), atmospheric temperature, and water content between these eruptions. This close agreement suggests that aggregation may be treated as a discrete process that is insensitive to eruptive style or magnitude. This result offers the potential for a simple, computationally efficient parameterization scheme for use in operational model forecasts. Further research may indicate whether this narrow range also reflects physical constraints on processes in the evolving cloud.

Adjusting particle-size distributions to account for aggregation in tephra-deposit model forecasts

Volcanic ash transport and dispersion models (VATDs) are used to forecast tephra deposition during volcanic eruptions. Model accuracy is limited by the fact that fine ash aggregates, altering patterns of deposition. In most models this is accounted for by ad hoc changes to model input, representing fine ash as aggregates with density ρagg, and a log-normal size distribution with median μagg and standard deviation σagg. Optimal values may vary between eruptions. To test the variance, we used the Ash3d tephra model to simulate four deposits: 18 May 1980 Mount St. Helens; 16–17 September 1992 Crater Peak (Mount Spurr); 17 June 1996 Ruapehu; and 23 March 2009 Mount Redoubt. In 158 simulations, we systematically varied μagg and σagg, holding ρagg constant at 600 kg m−3. We evaluated the fit using three indices that compare modeled versus measured (1) mass load at sample locations; (2) mass load versus distance along the dispersal axis; and (3) isomass area. For all deposits, under these inputs, the best-fit value of μagg ranged narrowly between ~ 2.1–2.5 φ (0.23–0.18 mm), despite large variations in erupted mass (0.25–50 Tg), plume height (8.5–25 km), mass fraction of fine (< 0.063 mm) ash (3–59 %), atmospheric temperature, and water content between these eruptions. This close agreement suggests that aggregation may be treated as a discrete process that is insensitive to eruptive style or magnitude. This result offers the potential for a simple, computationally-efficient parameterization scheme for use in operational model forecasts. Further research may indicate whether this narrow range also reflects physical constraints on processes in the evolving cloud.