Gordeychik Boris, Churikova Tatiana, Kronz Andreas, Simakin Alexander, Wörner Gerhard First data on magma ascent and residence times retrieved from Fe-Mg and trace element zonation in olivine phenocrysts from Kamchatka basalts // Geophysical Research Abstracts. 2016. Vol. 18. P. EGU2016-12839.
Compositional zonation in olivine phenocrysts and diffusion modelling have been used in the last ten years to estimate magma residence times and the duration of magma ascent. The fundamental assumption is that mixing with newly injected magma into a reservoir triggers diffusional exchange between mafic olivine crystals and more evolved magma and that this magma mixing eventually triggers eruption. If depth of mixing is known, this translates to ascent rates of magmas to the surface. We applied this approach to a series of different arc basalt lavas from Kamchatka to constrain the rates of magma ascent and magma resident in what is one of the most active subduction zones in the world that is also dominated by an abundance of unusually mafic magmas. Our sample collection cover the principal modes of arc magmatism in Kamchatka: from different volcanic complexes (stratovolcano, dikes, summit eruptions, monogenetic cones), of different age (from Late-Pleistocene to Holocene and recent eruptions), from different magmatic regimes (long-lived volcanoes vs. monogenetic eruptions) and different major element composition (from basalt to basaltic andesite of different geochemical character including LILE enrichments). We analyzed and modelled zonation profiles for a range of elements with different diffusivities (e.g. Mg-Fe, Ca, Ni, Mn, Cr) to assess the role of variable diffusivities as a function of major and trace elements in the olivines from different P-T conditions. First data were obtained on samples from the Klyuchevskoy, Shiveluch and Tolbachik, including recent most eruption in 2012/2013. These data show that for some samples the zonation patterns are much more complex than is usually observed: high-Mg olivines at different volcanoes have very different zonation patterns, including normally, reversely zoned grains or even show highly complex repetitive zonation that indicate large compositional changes in the surrounding magma at very short time scales (years). Thus in some Kamchatka basalts, we observe unusual Mg-Fe zonations that are linked to complex mixing, possibly resorption and subsequent crystal growth processes that are generally not preserved due to fast diffusion of Mg-Fe. Based on a first assessment of our measured profiles, the values for diffusion times in Fo-rich olivines (88 to 92% Fo) vary from only a few months to years and thus magma ascent from deep magma sources must have been fast.
Gordeychik Boris, Churikova Tatiana, Kronz Andreas, Simakin Alexander, Wörner Gerhard Olivine zoning in high-Mg basalts of the Shiveluch volcano (Kamchatka) // Geophysical Research Abstracts. 2017. Vol. 19. P. EGU2017-10473.
Shiveluch volcano located in northern Kamchatka erupted mainly high-Mg andesites during Holocene times. However, tephrochronologists found two Holocene tephra layers that are unusual for this volcano: a high-Mg middle-K basalts with an age of 7600 yr BP and high-Mg high-K basalt with an age of 3600 yr BP [Volynets et al, 1997]. The proximal outcrops for these two tephra deposits were discovered just recently [Churikova et al., 2010; Gorbach & Portnyagin, 2011]. Our study of olivines from the high-Mg basalts documents unusual Mg-Fe zonation [Gordeychik et al., 2016]: Inner cores of olivines from both eruptions show Fo87-92, falling to the rim to Fo75-85. In the outer cores of both basalt tephra, forsterite decreases linearly abruptly changing to a steeper gradient towards the rim. Electron microprobe element maps reveal the complex and highly unusual zoning features of these olivines.
The inner cores of the olivines of 7600 yr BP tephra have bell-shaped distributions for forsterite and nickel. The maximum forsterite in their core can be up to Fo92, decreasing outward to the outer core to Fo86. At the same time, the trace elements in the inner core remain constant. Such element distribution is consistent with diffusion of Fe, Mg, and Ni in the initially uniform high Mg cores after the phenocrysts were changed to non-equilibrium in a less mafic melt. The shape of the inner cores suggests partial dissolution after magma mixing. The interfaces between the inner and outer cores are marked by abundant melt/fluid inclusions. The inner cores were overgrown by olivine with Fo90 when the crystals moved to the high-Mg melt. As result some olivine grains have the maximum forsterite values in the outer core. The specific feature of the olivine outer cores from basalt of the 7600 yr BP tephra eruption are concentric zones with higher values of Ca, Cr, Al, P. One of the crystals has five distinct growth zones with high Cr concentrations. The width of these zones can be only a few microns and thus such zones are often missed in typical quantitative point measurements in microprobe profiles.
Inner cores of olivines from the 3600 yr BP tephra are uniform in forsterite and nickel. However, Al and Ca element distribution maps show in inner cores higher concentrations with rather smooth contours. This suggests that initially the olivines were formed from high-Al and high-Ca melt, then were dissolved and the overgrowth zonation has been smoothed out due to faster Mg-Fe diffusion. Only Ca and Al with low diffusivity were conserved. The concentric zones with higher element concentrations are not so well expressed in olivines from the 3600 yr BP tephra, but some distinct growth zones are also shown in Ca, Cr, and P.
Information extraction and decoding of the elemental maps allow seeing highly complex growth-dissolutiondiffusion history of magma mixing processes prior to eruption. This research was supported by RFBR-DFG grant # 16-55-12040.
Gorelchik V.I., Levina V.I. Space-Time and Dinamic Characteristics of Earthquakes Associated with the 1974 Eruption of Klyuchevskoi // Volcanology and Seismology. 1989. Vol. 7. № 6. P. 943-974.
Gorelchik V.I., Shirokov V.A., Firstov P.P., Chubarova O.S. Shiveluch volcano: seismicity, deep structure and forecasting eruptions (Kamchatka) // Journal of Volcanology and Geothermal Research. 1997. Vol. 78. № 1–2. P. 121 - 137. doi: 10.1016/S0377-0273(96)00108-4.
The deep structure, Wadati-Benioff zone (focal zone) geometry and the magma feeding system of Shiveluch volcano are investigated based on 1962–1994 detailed seismic surveillance. A focal zone beneath Shiveluch is dipping at an angle of 70° at depths of 100–200 km. Based on the revealed interrelations between seismicity at depths of 105–120 km and an extrusive phase of its eruptions in 1980 through 1994, it is inferred that primary magmas, periodically feeding the crustal chamber, are melted at depths of at least 100 km. An upsurge of extrusive-explosive activity at the volcano is preceded and accompanied by the increasing number and energy of both volcanic earthquakes beneath the dome and tectonic or volcano-tectonic earthquakes in the zones of NW-striking crustal faults near the volcano.The eruption of April 1993 has been the most powerful since 1964. It was successfully predicted based on interactive use of all seismic data. At the same time the influence of seismicity at depths of 105–120 km under the volcano on the style (and consequently on prediction) of its activity is decisive.
Gorokhova N.V., Melnik O.E., Plechov P.Yu., Shcherbakov V.D. Numerical simulation of plagioclase rim growth during magma ascent at Bezymianny Volcano, Kamchatka // Journal of Volcanology and Geothermal Research. 2013. Vol. 263. P. 172 - 181. doi: 10.1016/j.jvolgeores.2013.03.020.
Slow CaAl-NaSi interdiffusion in plagioclase crystals preserves chemical zoning of plagioclase in detail, which, along with strong dependence of anorthite content in plagioclase on melt composition, pressure, and temperature, make this mineral an important source of information on magma processes. A numerical model of zoned crystal growth is developed in the paper. The model is based on equations of multicomponent diffusion with diagonal cross-component diffusion terms and accounts for mass conservation on the melt–crystal interface and growth rate controlled by undercooling. The model is applied to the data of plagioclase rim zoning from several recent Bezymianny Volcano (Kamchatka) eruptions. We show that an equilibrium growth model cannot explain crystallization of naturally observed plagioclase during magma ascent. The developed non-equilibrium model reproduced natural plagioclase zoning and allowed magma ascent rates to be constrained. Matching of natural and simulated zoning suggests ascent from 100 to 50 MPa during 15–20 days. Magma ascent rate from 50 MPa to the surface varies from eruption to eruption: plagioclase zoning from the December 2006 eruption suggests ascent to the surface in less than 1 day, whereas plagioclase zoning from March 2000 and May 2007 eruptions are better explained by magma ascent over periods of more than 30 days). Based on comparison of diffusion coefficients for individual elements a mechanism of atomic diffusion during plagioclase crystallization is proposed.
Gorshkov G.S. Catalogue of the active volcanoes of the world including solfatara fields. Part VII - Kurile Islands / Ed. Signore Francesco Napoli, Italy: the International Volcanological Association. 1958. 100 p.
Gorshkov G.S. Determination of the explosion energy in some volcanoes according to barograms // Bulletin Volcanologique. 1960. Vol. 23. № 2. P. 141-144.