MODELING AND FORECASTING TEPHRA ACCUMULATION WITH TEPHRA2

CONNOR, Charles, Department of Geology, University of South Florida, SCA 528, 4202 E. Fowler Ave, Tampa, FL 33620, CONNOR, Laura, Department of Geology, University of South Florida, 4202 E. Fowler Ave, SCA 528, Tampa, FL 33620, BONADONNA, Costanza, Section des Sciences de la Terre, Université de Genève, 13, rue des Maraîchers, Geneva, 1205, Switzerland and VOLENTIK, Alain C.M., Department of Geology, Univsersity of South Florida, SCA 528, 4202 E. Fowler Ave, Tampa, FL 33620, cconnor@cas.usf.edu

TEPHRA2 is a parallel code written in C and MPI and is used to forecast tephra accumulation following explosive volcanic eruptions. The code uses a closed-form solution of the advection-diffusion equation, particle fall velocities that depend on local Reynold’s number, and stratified wind field to forecast tephra accumulation, usually expressed as kilogram per cubic meter, particle size distribution at specific locations from the vent, and maximum clast size expected as a function of distance from the vent. In practice, deposits of specific eruptions can be modeled if sufficient field data are available using TEPHRA2 and a nonlinear inversion method, the downhill simplex algorithm, to estimate best fit eruption parameters. We demonstrate wide applicability of TEPHRA2 with examples from Cerro Negro, Nicaragua (small volume basaltic eruption), Colima, Mexio (short-lived basaltic-andesite vulcanian eruption), and Chaiten, Chile (comparatively small volume Plinian eruption of rhyolite composition). In each case models of the eruption based on inversion of the deposit thickness and/or granulometry compares well with independent observations. An advantage of this inversion procedure is that it is possible to assign numerical uncertainty to estimates of deposit volume, eruption column height, eruption duration, and similar eruption input parameters. TEPHRA2 is used as a forecasting tool with Monte Carlo simulation. Input parameters are specified using probability density functions for eruption mass, column height and grainsize distribution. Wind data are specified using historical records or reanalysis data. Output is cast as a probability density function of tephra accumulation (or maximum grainsize) at a specific point, or the probability that tephra accumulation will exceed some threshold value. One challenge in the continued development of codes like TEPHRA2 involves model parsimony. How complex need a code become in order to forecast activity from a hazard perspective?

Session No. 82

T64. Physics of Volcanic Eruptions: Implications for Hazards

Monday, 19 October 2009: 8:00 AM-12:00 PM

B116 (Oregon Convention Center)

Geological Society of America Abstracts with Programs. Vol. 41, No. 7, p.232

© Copyright 2009 The Geological Society of America (GSA), all rights reserved. Permission is hereby granted to the author(s) of this abstract to reproduce and distribute it freely, for noncommercial purposes. Permission is hereby granted to any individual scientist to download a single copy of this electronic file and reproduce up to 20 paper copies for noncommercial purposes advancing science and education, including classroom use, providing all reproductions include the complete content shown here, including the author information. All other forms of reproduction and/or transmittal are prohibited without written permission from GSA Copyright Permissions.

Back to: T64. Physics of Volcanic Eruptions: Implications for Hazards

<< Previous Abstract | Next Abstract

2009 Portland GSA Annual Meeting (18-21 October 2009)

MODELING AND FORECASTING TEPHRA ACCUMULATION WITH TEPHRA2