Crystallization and Melt Removal at Arenal Volcano, Polytopic Vector Analysis

Tephra sequences ET3 and ET4 from Arenal volcano in Costa Rica have recently been interpreted to be a product of crystal fractionation by Bolge and coworkers in a series of papers (2004, 2006). The two tephra units are part of a sequence of 22 tephra units that represent a 7000 year span of the Arenal volcano activity. The tephro- stratigraphy has been described extensively by Melson (1982; 1994). The ET3 and ET4 tephras were interpreted (based on major- and trace-element, isotopic analyses of whole rocks and microchemical analyses of individual phases) as clear evidence of crystal separation by gravity settling (Bolge et al., 2004, 2006). The lower ET4 tephra sequence (andesitic and crystal poor) and the upper ET3 tephra (basaltic and crystal rich) represent an inverted snapshot of the magma chamber with contrasting geochemical properties. The ET3 sequence (deeper part of the magma chamber) has nearly constant composition with only a few elements varying stratigraphically (best represented by CaO). This is consistent with gradually decreasing amounts of melt in the upper part of ET3. The lower ET4 tephra (upper part of the magma chamber) contains large chemical gradients in both incompatible and compatible elements. In the present study we use whole-rock geochemical data from the recent tephra sequences ET3 and ET4 as inputs to Polytopic Vector Analysis (PVA) (for a review of this method see Vogel and coworkers, in press). With this method we produce a three end member solution that is consistent with crystallization of Olivine, plagioclase and pyroxene from the most mafic end member (EM1) resulting in a crystal rich mush zone. As crystallization progresses the compositions of the liquids are driven towards an intermediate end member (EM3), which has an intermediate composition liquid. At EM3 composition, rapid depletion of FeO, MgO and TiO2 by crystallization of Fe-Ti oxides, rapidly drives the liquid composition towards the silicic EM1 (incompatible element enriched end member). Using PVA we refine the interpretations of Bolge and coworkers and show that melt from the crystalline rich ET3 tephra was removed and ponded in the magma chamber above the crystalline mush (top part of ET4 unit). Thus when the eruption occurred the most evolved tephra (ET4) were deposited first followed by the least evolved tephra (ET3), which resulted in sampling of a chemically zoned magma chamber. Using PVA on stratigraphically controlled whole-rock analyses of tephra samples, we can unambiguously identify processes and end members that are involved in crystal accumulation and liquid separation processes. Thus PVA is a rigorous analytical tool that uses only whole-rock chemical data to produce robust results that can be used with other analytical techniques to test petrological models.