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.

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.

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.