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.

Fossil diatom assemblages in a sediment core from a small lake in Central Kamchatka (Russia) were used to reconstruct palaeoenvironmental conditions of the late Holocene. The waterbody may be a kettle lake that formed on a moraine of the Two-Yurts Lake Valley, located on the eastern slope of the Central Kamchatka Mountain Chain. At present, it is a seepage lake with no surficial outflow. Fossil diatom assemblages show an almost constant ratio between planktonic and periphytic forms throughout the record. Downcore variations in the relative abundances of diatom species enabled division of the core into four diatom assemblage zones, mainly related to changes in abundances of Aulacoseira subarctica, Stephanodiscus minutulus, and Discostella pseudostelligera and several benthic species. Associated variations in the composition and content of organic matter are consistent with the diatom stratigraphy. The oldest recovered sediments date to about 3220 BC. They lie below a sedimentation hiatus and likely include reworked deposits from nearby Two-Yurts Lake. The initial lake stage between 870 and 400 BC was characterized by acidic shallow-water conditions. Between 400 BC and AD 1400, lacustrine conditions were established, with highest contributions from planktonic diatoms. The interval between AD 1400 and 1900 might reflect summer cooling during the Little Ice Age, indicated by diatoms that prefer strong turbulence, nutrient recycling and cooler summer conditions. The timing of palaeolimnological changes generally fits the pattern of neoglacial cooling during the late Holocene on Kamchatka and in the neighbouring Sea of Okhotsk, mainly driven by the prevailing modes of regional atmospheric circulation.

A geometric technique is introduced to estimate the height of volcanic eruption columns using the generally discarded near-limb portion of geostationary imagery. Such oblique observations facilitate a height-by-angle estimation method by offering close-to-orthogonal side views of eruption columns protruding from the Earth ellipsoid. Coverage is restricted to daytime point estimates in the immediate vicinity of the vent, which nevertheless can provide complementary constraints on source conditions for the modeling of near-field plume evolution. The technique is best suited to strong eruption columns with minimal tilting in the radial direction. For weak eruptions with severely bent plumes or eruptions with expanded umbrella clouds the radial tilt/expansion has to be corrected for either visually or using ancillary wind profiles. Validation on a large set of mountain peaks indicates a typical height uncertainty of ±500 m for near-vertical eruption columns, which compares favorably with the accuracy of the common temperature method.

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.

We determine a 3-D P-wave velocity structure of the mantle down to 700 km depth under the Kamchatka peninsula using 678 P-wave arrival times collected from digital seismograms of 75 teleseismic events recorded by 15 portable seismic stations and 1 permanent station in Kamchatka. The subducting Pacific slab is imaged clearly that is visible in the upper mantle and extends below the 660-km discontinuity under southern Kamchatka, while it shortens toward the north and terminates near the Aleutian–Kamchatka junction. Low-velocity anomalies are visible beneath northern Kamchatka and under the junction, which are interpreted as asthenospheric flow. A gap model without remnant slab fragment is proposed to interpret the main feature of high-V anomalies. Combining our tomographic results with other geological and geophysical evidences, we consider that the slab loss may be induced by the friction with surrounding asthenosphere as the Pacific plate rotated clockwise at about 30 Ma ago, and then it was enlarged by the slab-edge pinch-off by the asthenospheric flow and the presence of Meiji seamounts. As a result, the slab loss and the subducted Meiji seamounts have jointly caused the Pacific plate to subduct under Kamchatka with a lower dip angle near the junction, which made the Sheveluch and Klyuchevskoy volcanoes shift westward.