Abstract We test two methods of 3-D acoustic source localization on volcanic explosions and small-scale jetting events at Karymsky Volcano, Kamchatka, Russia. Recent infrasound studies have provided evidence that volcanic jets produce low-frequency aerodynamic sound (jet noise) similar to that from man-made jet engines. For man-made jet noise, noise sources localize along the turbulent jet flow downstream of the nozzle. Discrimination of jet noise sources along the axis of a volcanic jet requires high resolution in the vertical dimension, which is very difficult to achieve with typical volcano-acoustic network geometries. At Karymsky Volcano, an eroded edifice (Dvor Caldera) adjacent to the active cone provided a platform for the deployment of five infrasound sensors in July 2012 with intra-network relief of ~ 600 m. The network was designed to target large-scale jetting, but unfortunately only small-scale jetting and explosions were recorded during the 12-day experiment. A novel 3-D inverse localization method, srcLoc, is tested and compared against a more common grid-search semblance technique. Simulations using synthetic signals show that srcLoc is capable of determining vertical solutions to within ± 150 m or better (for signal-to-noise ratios ≥ 1) for this network configuration. However, srcLoc locations for explosions and small-scale jetting at Karymsky Volcano show a persistent overestimation of source elevation and underestimation of sound speed. The semblance method provides more realistic source locations, likely because it uses a fixed, realistic sound speed of ~ 340 m/s. Explosion waveforms exhibit amplitude relationships and waveform distortion strikingly similar to those theorized by modeling studies of wave diffraction around the crater rim. We suggest that the delay of acoustic signals and apparent elevated source locations are due to raypaths altered by topography and/or crater diffraction effects, implying that topography in the vent region must be accounted for when attempting 3-D volcano acoustic source localization. Though the data presented here are insufficient to resolve small-scale jet noise sources, similar techniques may be successfully applied to large volcanic jets in the future.

Abstract Data on the geology, petrography, mineralogy, and geochemistry of rocks from Kamen Volcano (Central Kamchatka Depression) are presented and compared with rocks from the neighbouring active volcanoes. The rocks from Kamen and Ploskie Sopky volcanoes differ systematically in major elemental and mineral compositions and could not have been produced from the same primary melts. The compositional trends of Kamen stratovolcano lavas and dikes are clearly distinct from those of Klyuchevskoy lavas in all major and trace element diagrams as well as in mineral composition. However, lavas of the monogenetic cones on the southwestern slope of Kamen Volcano are similar to the moderately high-Mg basalts from Klyuchevskoy and may have been derived from the same primary melts. This means that the monogenetic cones of Kamen Volcano represent the feeding magma for Klyuchevskoy Volcano. Rocks from Kamen stratovolcano and Bezymianny form a common trend on all major element diagrams, indicating their genetic proximity. This suggests that Bezymianny Volcano inherited the feeding magma system of extinct Kamen Volcano. The observed geochemical diversity of rocks from the Klyuchevskaya group of volcanoes can be explained as the result of both gradual depletion over time of the mantle N-MORB-type source due to the intense previous magmatic events in this area, and the addition of distinct fluids to this mantle source.

Abstract For the past few decades, Bezymianny volcano has erupted once to twice per year. Here, I examine eight eruptive events between 2006 and 2010. This is the first time period for which proximal or broadband seismic data have been recorded at Bezymianny. Several recurring patterns are demonstrated in advance of eruptions. Eruptions are generally preceded by 12–36 h of tremor energy elevated by 2 to 3 orders of magnitude. Locatable earthquake activity is quite erratic in the days before eruptions. For eruptions of juvenile magma, however, the cumulative moment magnitude increases with the repose time since the previous eruption. Though tenuous, this relationship is statistically significant and could improve forecasts of Bezymianny eruptions. The most energetic eruptions demonstrate increasing multiplet activity in the run-up, followed by a rapid cessation at the time of eruption. When present, this behavior marks increasing pressure in the conduit system as degassing eclipses the capacity for venting. Very long period seismicity (> 20 s periods) accompanies some eruptions. These tend to be the same short-lived high-energy eruptions that exhibit multiplet precursors. Four eruptions are examined in detail to illustrate the variety in eruption mechanisms. Lava dome collapses, sustained eruptions, singular paroxysmal explosions and post-explosion lava flows occur in different combinations demonstrating that more than one eruption trigger is regulating Bezymianny. Compared to Bezymianny's fifty-year modern history, recent eruptions have been shorter-lived and separated by longer repose times. Some evidence suggests that these eruptions may be increasingly explosive—a speculation that carries significant hazard implications. If true, however, this threat is tempered by solid evidence that the most explosive eruptions are preceded by the clearest precursors, suggesting an ability to improve the already excellent eruption forecasts available for Bezymianny.

We present the results of time-dependent local earthquake tomography for the Kluchevskoy group of volcanoes in Kamchatka, Russia. We consider the time period from 1999 to 2009, which covers several stages of activity of Kluchevskoy and Bezymianny volcanoes. The results are supported by synthetic tests that recover a common 3D model based on data corresponding to different time windows. Throughout the period, we observe a robust feature below 25 km depth with anomalously high Vp/Vs values (up to 2.2). We interpret this feature as a channel bringing deep mantle materials with high fluid and melt content to the bottom of the crust. This mantle channel directly or indirectly determines the activity of all volcanoes of the Kluchevskoy group. In the crust, we model complex structure that varies over time. During the pre-eruptive period, we detected two levels of potential magma storage: one in the middle crust at 10–12 km depth and one close to the surface just below Kluchevskoy volcano. In 2005, a year of powerful eruptions of Kluchevskoy and Besymiyanny volcanoes, we observe a general increase in Vp/Vs throughout the crust. In the relaxation period following the eruption, the Vp/Vs values are generally low, and no strong anomalous zones in the crust are observed. We propose that very rapid variations in Vp/Vs are most likely due to abrupt changes in the stress and deformation states, which cause fracturing and the active transport of fluids. These fluids drive more fracturing in a positive feedback system that ultimately leads to eruption. We envision the magma reservoirs beneath the Kluchevskoy group as sponge-structured volumes that may quickly change the content of the molten phases as fluids pulse rapidly through the system.

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.

We present a new dataset for whole-rock major, trace, isotopic, and phenocryst compositions indicating a genetic link between andesites of the Holocene eruptions of the Bezymianny stratovolcano (the Bezymianny stage), the andesitic to dacitic Late Pleistocene lava dome complex (the pre-Bezymianny stage), and the magnesian to high-alumina basalts of the adjacent Kliuchevskoi Volcano. We demonstrate that volcanic products from the Bezymianny stage of volcano evolution are most likely the products of magma mixing between silicic products of the earliest stages of magma fractionation and the less evolved basaltic andesite parental melts periodically injected into the magma reservoir. In contrast, the intermediate and silicic magmas of the pre-Bezymianny stage together with basalts from Kliuchevskoi much more closely resemble the liquid line of descent and may represent a unique prolonged and continuous calc-alkaline trend of magma evolution from high-magnesian basalt to dacite. As a result of the geothermobarometry, we recognize variable conditions of magma fractionation and magma storage beneath Bezymianny for different magma types during its evolution since the Late Pleistocene: (1) 1100–1150 °C, 500–640 MPa, 1–2.5 wt. H2O for parental basaltic andesite; (2) 1130–1050 °C, 700–600 MPa, 2.5–5 wt. H2O for two-pyroxene andesites; (3) 1040–990 °C, 560–470 MPa, 5–6.5 wt. H2O for orthopyroxene-bearing andesites; (4) 950–1000 °C, 450–150 MPa, 3.5–5.5 wt. H2O for hornblende-bearing andesites; and (5) 950–900 °C, 410–250 MPa, 6–7 wt. H2O for dacites. Repeated basalt injections and magma fractionation combined with internal mixing in the magma chamber are the main processes responsible for both the complex petrography and the geochemical trends observed in the lavas of Bezymianny Volcano.

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.

Abstract Interferometric synthetic aperture radar (InSAR) images reveal a pre-eruption deformation signal at Kizimen volcano, Kamchatka, Russia, where an ongoing eruption began in mid-November, 2010. The previous eruption of this basaltic andesite-to-dacite stratovolcano occurred in 1927–1928. InSAR images from both ascending and descending orbital passes of Envisat and ALOS PALSAR satellites show as much as 6 cm of line-of-sight shortening from September 2008 to September 2010 in a broad area centered at Kizimen. About 20 cm of opening of a nearly vertical dike provides an adequate fit to the surface deformation pattern. The model dike is approximately 14 km long, 10 km high, centered 13 km beneath Kizimen, and strikes NE–SW. Time-series analysis of multi-temporal interferograms indicates that (1) intrusion started sometime between late 2008 and July 2009, (2) continued at a nearly constant rate, and (3) resulted in a volume expansion of 3.2 × 107 m3 by September 2010, i.e., about two months before the onset of the 2010 eruption. Earthquakes located above the tip of the dike accompanied the intrusion. Eventually, magma pressure in the dike exceeded the confining strength of the host rock, triggering the 2010 eruption. Our results provide insight into the intrusion process that preceded an explosive eruption at a Pacific Rim stratovolcano following nearly a century of quiescence, and therefore have implications for monitoring and hazards assessment at similar volcanoes elsewhere.

Volcanic eruptions from Kamchatka have deposited many unique tephra layers over a large region within the North Pacific, providing important isochrons between key sites such as marine ODP core 883 (Pacific Ocean, Detroit Seamount) and Elgygytgyn Lake (Chukotka, eastern Siberia). Here we present a compilation of C14 dates on major Holocene tephras from the volcanically highly active region, based on decades of detailed stratigraphical fieldwork on Shiveluch, Kliuchevskoy, and other volcanoes.The 12-m thick tephra sequence at the Kliuchevskoy slope has been continuously accumulating during the last ∼11 ka. It contains over 200 visible individual tephra layers and no datable organic material. The section is dominated by dark-gray mafic cinders related to Kliuchevskoy activity. In addition, it contains 30 light-colored thin layers of silicic tephra from distant volcanoes including 11 layers from Shiveluch volcano located only 65 km to the north. We have used EMPA glass analysis to correlate most of the marker tephra layers to their source eruptions dated earlier by C14 (Braitseva et al., 1997; Ponomareva et al., 2007), and in this way linked Kliuchevskoy tephra sequence to sequences at other volcanoes including Shiveluch. The C14 dates and tephras from the northern Kamchatka are then combined into a single Bayesian framework taking into account stratigraphical ordering within and between the sites. This approach has allowed us to enhance the reliability and precision of the estimated ages for the eruptions. Age-depth models are constructed to analyse changes in deposition rates and volcanic activity throughout the Holocene. This detailed chronology of the eruptions serves as a basis for understanding temporal patterns in the eruption sequence and geochemical variations of magmas. This research could prove important for the long-term forecast of eruptions and volcanic hazards.

This study presents data on the geochemical composition of boiling mud pools at the Mutnovsky volcano. The physicochemical characteristics of the pools and the concentrations of major, minor and trace elements in pool solutions vary widely. A comparison of the geochemical compositions of host rocks and solutions indicates that leaching from rocks is not the only source of chemicals in thermal solutions. Geophysical studies reveal the inner structure of thermal fields, which reflect the shapes of the underground reservoirs and feed channels. Using geophysical methods (electrical resistivity tomography and frequency domain investigations), it was shown that the vertical structure and complex geochemical zonation of the feed channels leads to a high contrast in the compositions of the mud solutions. These findings answer questions about the origin and composition of surface manifestations. To elucidate the mechanisms of solution formation, an attempt was made to describe the magmatic fluid evolution and the resulting mixing of waters by physical and mathematical models. The model illustrates fluid migration from a magma chamber to the surface. It is shown that the formation of brines corresponding to the mud pool composition is possible during secondary boiling.