EarthRef.org Reference Database (ERR)
Development and Maintenance by the EarthRef.org Database Team

Detailed Reference Information
Spera 1984
Spera, F.J. (1984). Some numerical experiments on the withdrawal of magma from crustal reservoirs. Journal of Geophysical Research 89: doi: 10.1029/JB080i010p08222. issn: 0148-0227.

Compositional heterogeneity in the form of continuous or discontinuous chemical and thermal gradients in lava and/or pyroclastic flows is very common. An understanding of the dynamics of magma withdrawal is essential to palinspastic reconstruction of intensive variable gradients in magma reservoirs. Important parameters governing the extent of subterranean magma mixing triggered by an eruption include the vertical structure of density and viscosity within the chamber, the discharge the size and shape of the chamber, and whether eruption takes place along a sublinear fissure, a ring fracture, or a central vent. A numerical model has been set up to study isoviscous magma withdrawal from a central vent conduit as a function of the Reynolds number, the reservoir to conduit width ratio reservoir aspect ratio (width/depth) and differing kinematic boundary conditions. Both open (magma recharge) and closed (caldera collapse) system behavior are considered. Finite difference solutions to the vorticity transport and Poisson equations enable determination of vorticity, stream function and velocity fields as a function of time. The most petrologically significant output is the stream function (particle trajectories) and evacuation isochron diagrams. An evacuation isochron represents the locus of points within the chamber such that magma parcels along a given isochron arrive at the bottom of the volcanic conduit concurrently. Open systems evolve toward a time invariant state (fully developed flow). Spin-up times depend on chamber aspect ratio (Br/Dr), reservois/conduit width ratio (Br/Bc), and Reynolds number (Re), Br, Dr, Bc represent chamber half-width, depth, and conduit half-width, respectively. Spin-up times are relatively small (1/10 to 1/5) fractions of typical eruption duraitons. The shape and orientation of evacuation isochrons (FI) depend on Re (increasing Re decreases withdrawal depth) and geometric factors (increasing Br/Bc) at constant Re and Br/Dr or decreasing Br/Dr or decreasing Br/Dr at constant Br/Bc and Re increases withdrawal depth). A significant amount of roofward magma can remain untapped in a chamber even for long duration eruptions. Systems driven by caldera collapse also involve juxtaposition of roofward and deep-seated magma during the course of an eruption. Relative to the magma recharge (open system) situation. EI's are laterally elongated. The extent of vertical mixing is thus smaller although small significant in this case. Maximum withdrawal depths very monotonically in both cases. There is excellent qualitative agreement between predictions based on the numerical experiments and Fe-Ti oxide temperatures for a thermally zoned ash flow deposit south of Mono Lake in eastern California (Bishop Tuff).

BACKGROUND DATA FILES

Abstract

Journal
Journal of Geophysical Research
http://www.agu.org/journals/jb/
Publisher
American Geophysical Union
2000 Florida Avenue N.W.
Washington, D.C. 20009-1277
USA
1-202-462-6900
1-202-328-0566
service@agu.org
Click to clear formClick to return to previous pageClick to submit