Since Bezymianny Volcano resumed its activity in 1956, eruptions have been frequent; recently with up to 1–2 explosive events per year. To investigate deformation related to this activity we installed a GPS network of 8 continuous and 6 campaign stations around Bezymianny. The two striking observations for 2005–2010 are (1) rapid and continuous network-wide subsidence between 8 and 12 mm/yr, which appears to affect KAMNET stations more than 40 km away where we observe 4–5 mm/yr of subsidence, and (2) only the summit station BZ09 shows slight deviations from the average motion in the north component at times of eruptions.
The network-wide subsidence cannot be explained by tectonic deformation related to the build-up of interseismic strain due to subduction of the Pacific plate. A first order model of surface loading by eruptive products of the Kluchevskoy Group of Volcanoes also explains only a fraction of the subsidence. However, a deep sill at about 30 km under Kluchevskoy that constantly discharges material fits our observations well. The sill is constrained by deep seismicity which suggests 9.5 km width, 12.7 km length, and a 13° dip-angle to the south-east. We infer a closing rate of 0.22 m/yr, which implies a volume loss of 0.027 km3/yr (0.16 m/yr and 0.019 km3/yr considering surface loading). Additional stations in the near and far field are required to uniquely resolve the spatial extent and likely partitioning of this source.
We explain the eruption related deformation at BZ09 with a very shallow reservoir, likely within Bezymianny's edifice at a depth between 0.25 km and 1.5 km with a volume change of 1–4 × 10− 4 km3. Much of the material erupted at Bezymianny may be sourced from deeper mid-crustal reservoirs with co-eruptive volume changes at or below the detection limit of the GPS network. Installation of more sensitive instruments such as tiltmeters would allow resolving of subtle co-eruptive motion.

Temporal irregularity of the output of volcanic material is studied for the sequence of large (V ≥ 0.5 km3, N = 29) explosive eruptions on Kamchatka during the last 10,000 years. Informally, volcanic productivity looks episodic, and dates of eruptions cluster. To investigate the probable self-similar clustering behavior of eruption times, we determine correlation dimension Dc. For intervals between events 800 and 10,000 years, Dc ≈ 1 (no self-similar clustering). However, for shorter delays, Dc = 0.71, and the significance level for the hypothesis Dc < 1 is 2.5%. For the temporal structure of the output of volcanic products (i.e., for the sequence of variable-weight points), a self-similar “episodic” behavior holds over the entire range of delays 100–10,000 years, with Dc = 0.67 (Dc < 1 at 3.4% significance). This behavior is produced partly by the mentioned common clustering of event dates, and partly by another specific property of the event sequence, that we call “order clustering”. This kind of clustering is a property of a time-ordered list of eruptions, and is manifested as the tendency of the largest eruptions (as opposed to smaller ones) to be close neighbors in this list. Another statistical technique, of “rescaled range” (R/S), confirms these results. Similar but weaker-expressed behavior was also found for two other data sets: historical Kamchatka eruptions and acid layers in Greenland ice column. The episodic multiscaled mode of the output of volcanic material may be a characteristic property of a sequence of eruptions in an island arc, with important consequences for climate forcing by volcanic aerosol, and volcanic hazard.

In a companion paper (Horváth et al., 2021), we introduced a new technique to estimate volcanic eruption column height from extremely oblique near-limb geostationary views. The current paper demonstrates and validates the technique in a number of recent eruptions, ranging from ones with weak columnar plumes to subplinian events with massive umbrella clouds and overshooting tops that penetrate the stratosphere. Due to its purely geometric nature, the new method is shown to be unaffected by the limitations of the traditional brightness temperature method, such as height underestimation in subpixel and semitransparent plumes, ambiguous solutions near the tropopause temperature inversion, or the lack of solutions in undercooled plumes. The side view height estimates were in good agreement with plume heights derived from ground-based video and satellite stereo observations, suggesting they can be a useful complement to established techniques.

Abstract We provide results of a detailed study of the first peridotite xenoliths of proven mantle origin reported from Bezymyanny volcano in the Klyuchevskoy Group, northern Kamchatka arc. The xenoliths are coarse spinel harzburgites made up mainly of Mg-rich olivine as well as subhedral orthopyroxene (opx) and Cr-rich spinel, and also contain fine-grained interstitial pyroxenes, amphibole and feldspar. The samples are unique in preserving the evidence for both initial arc mantle substrate produced by high-degree melt extraction and subsequent enrichment events. We show that the textures, modal and major oxide compositions of the Bezymyanny xenoliths are generally similar to those of spinel harzburgite xenoliths from Avacha volcano in southern Kamchatka. However, coarse opx from the Bezymyanny harzburgites has higher abundances of light and medium rare earth elements and other highly incompatible elements than coarse opx from the Avacha harzburgites. We infer that (1) the sub-arc lithospheric mantle beneath both Avacha and Bezymyanny (and possibly between these volcanoes) consists predominantly of harzburgitic melting residues, which experienced metasomatism by slab-related fluids or low-fraction, fluid-rich melts and (2) the degrees of metasomatism are higher beneath Bezymyanny. By contrast, xenolith suites from Shiveluch and Kharchinsky volcanoes 50–100 km north of the Klyuchevskoy Group include abundant cumulates and products of reaction of mantle rocks with silicate melts at high melt/rock ratios. The high melt flux through the lithospheric mantle beneath Shiveluch and Kharchinsky may be related to the asthenospheric flow around the northern edge of the sinking Pacific plate; lateral propagation of fluids in the mantle wedge south of the plate edge may contribute to metasomatism in the mantle lithosphere beneath the Klyuchevskoy Group volcanoes.