Предложены математическая модель и метод расчета течений магмы с вязкостью, зависящей от температуры, по питающему цилиндрическому каналу вулкана. Приведен пример расчета: радиус канала 10 м, глубина магматического очага 30 км, избыточное давление в очаге 20 бар, температура магмы в очаге 1300°С, вязкость магмы в нем 104, 105, 106 Пас. Рассмотрен начальный этап деятельности канала длительностью 10 лет. Выяснены условия замерзания канала (останавливающаяся экструзия), возникновения квазистационарного режима (устойчивое истечение), течения с интенсивным прогревом стенок.

A mathematical model and a method of computation are developed for the flow of magma with temperature-dependent viscosity in the conduit of a volcano. An example involves the following parameters: conduit radius is 10 m, depth to the magma chamber is 30 km, the overpressure, magma temperature and viscosity in the chamber are 20 bars, 1300° C, and 104, Ю5, 106Pa s, respectively. The initial phase of conduit operation lasting less than 10 years is considered. Conditions are determined under which the conduit freezes (stopping the extrusion), a quasistationary behavior sets in (steady flow), and the flow heats the conduit wall intensively.

Статья представляет результаты исследования Трещинного Толбачинского извержения им. 50-летия ИВиС ДВО РАН (ТТИ-50) 2012-2013 гг. по данным дистанционного зондирования.
Нами оценены количественные характеристики ТТИ-50: величины вертикальных смещений, площадь лавовых полей, их мощность и объем. Значения вертикальных смещений оценивались по серии радиоинтерферометрических пар для зоны извержения. Пары снимков соответствуют заключительной фазе извержения, когда величины смещения были небольшими. Вертикальные смещения рассчитаны для участков лавовых полей, значение когерентности которых превышает 0,4. Полученная серия значений вертикальных смещений отражает преимущественно процесс остывания лавы, для которого характерны просадки поверхности. Максимальные величины смещений составили 27 см за 24 дня.
Вычисление мощности лавовых полей выполнялось на основе анализа разновременных ЦМР. Высотные профили, измеренные геодезическими приемниками GPS в ходе полевых работ в августе 2013 года, были использованы для оценки точности ЦМР: общедоступных SRTM, SRTM-X, ASTER GDEM и ЦМР, построенной ИТЦ СКАНЭКС по двум оптическим стереопарам SPOT 6 (от 18.07.2013 и 11.10.2013). Среднеквадратическая погрешность определения абсолютных высот по ЦМР SRTM-X и SPOT6, по сравнению с данными наземных съемок, не превышает 5 м. Это делает возможным оценку мощности лавовых потоков по разности высот SRTM-X и SPOT6. ЦМР SPOT6 за две даты использовались совместно для исключения ошибок, связанных с облачностью и свежевыпавшим снегом. Максимальные значения мощности превышают 80 м. Вычисленный объем извержения - 0,521±0.25 км3.

This paper presents results of study of the 2012-2013 Tolbachik fissure eruption on the basis of remote sensing
techniques.
We have calculated values of vertical displacements, lava thickness and the volume of the erupted lava. Values of
vertical displacements were estimated using a series of
radar interferometric pairs for the Tolbachik eruption zone.
These pairs correspond to the concluding phase of the erupti
on, when vertical displacements were relatively small.
Vertical displacements were calculated for parts of lava fields with coherence value over 0,4. The obtained values of
vertical displacement are typical for subsidence caused by cooling lava flows. The maximum value of subsidence is
27 cm for 24 days. The calculation of lava thickness was based on comparison of multitemporal DEMs. Height profiles measured by geodetic GPS receivers during fieldwork in August 2013 were used to estimate the quality of DEMs, derived from satellite imagery: freely available SRTM, SRTM-X, ASTER GDEM and the DEMs calculated at RDC ScanEx from two stereopairs of SPOT6 images (of 18.07.2013 and 11.10.2013). The RMS error for heights of SRTM-X and
SPOT 6 in relation to GPS data is within ±5 m. This enables to estimate the total thickness of new lava fields on the
basis of height differences between SRTM-X and SPOT 6 DEMs. Both SPOT 6 DEMs were used together to eliminate errors caused by clouds and snow. The maximum lava thickness is over 80 m. The volume of the erupted lava is 0,521±0,25 km3.

Приведен алгоритм восстановления параметров пепловых выбросов – суммарной массы и ее распределения по высотам – при эксплозивных извержениях. Решение обратной задачи строится на основе метода множественной регрессии, при минимальной априорной информации о характере эксплозивного процесса. В качестве примера, рассмотрено сильное эксплозивное событие на вулкане Безымянный 24.12.2006 г., для которого распределение массы пеплового выброса по высотам, согласно расчетам, частично контролировалось выносом пеплового материала в облаках пирокластических потоков. Данная особенность проявилась в характерном двухмодальном распределении массы выброса с максимумами на высотах средней тропосферы и нижней стратосферы.

The article provides an algorithm for recovery of parameters of ash emissions (total volume and its height distribution) during explosive eruptions. The solution for the corresponding inverse task uses a multiple regression approach with minimal a prior information on the eruption dynamics. As an example, we consider a strong explosive event at Bezymianny Volcano, Kamchatka, on 24.12.2006. The estimations showed that the mass distribution for ash emission with heights was partially controlled by the emission of ash material inside the clouds from pyroclastic flows. This peculiarity was revealed as a bimodal distribution of the emission mass with maximums at the mid tropospheric and low stratospheric heights.

Приводятся результаты гравиметрических, сейсмических и электромагнитных исследований. Выявлены основные особенности глубинного строения района. Создана комплексная геолого-геофизическая модель земной коры под Авачинским вулканом. Она включает коровый магматический очаг на глубине ~15-25 км, перекрывающую его интрузию и периферический очаг под конусом вулкана на глубине ~0-2 км, а также зону, насыщенную жидкими флюидами в Авачинском грабене. Рассмотрены возможные геодинамические процессы, протекающие в земной коре в настоящее время. Важная роль отводится коровой проницаемой зоне, содержащей флюиды. Даны рекомендации для бурения глубокой скважины в районе Авачинского грабена с целью поисков геотермального месторождения.

Results are presented from gravity, seismic and electromagnetic studies. Main features of the deep structure of the area have been identified. A multidisciplinary geologic-geophysical model has been developed for the crust beneath Avacha Volcano. The model involves a crustal magma chamber at a depth of about 15-25 km, an intrusion that overlies it, and a peripheral chamber under the volcanic cone at 0-2 km depth, as well as a fluid-saturated zone in the Avacha Graben. We discuss possible geodynamic processes that are going on in the crust at present. Importance is attached to the fluid-containing crustal permeable zone. Recommendations are provided for drilling a deep well in the Avacha Graben area to search for a geothermal field.

Рассмотрены методика и результаты магнитотеллурического зондирования в модификациях АМТЗ и МТЗ. Аудиомагнитотеллурическое зондирование (АМТЗ) было проведено впервые в районе современного извержения Толбачинского вулкана. Результаты анализа магнитотеллурических параметров свидетельствуют, что геоэлектрическую среду, в связи с региональным разломом, можно аппроксимировать в виде двумерно-неоднородной модели. В качестве основных для интерпретации приняты продольная и поперечная кривые зондирований. Совместный анализ этих кривых и псевдоразрезов фаз импеданса свидетельствуют о геоэлектрической неоднородности среды в районе прорыва магматических расплавов им. С.И. Набоко. По данным бимодальной инверсии кривых АМТЗ получен геоэлектрический разрез, содержащий проводящую неоднородность, связываемую с разломом, по которому флюиды поступают к дневной поверхности. Наряду с АМТЗ для изучения глубинной электропроводности использованы МТЗ в расширенном диапазоне, по которым выделяется коровая проводящая аномалия на глубинах 15–35 км. По данным АМТЗ, МТЗ и другой геолого-геофизической информации составлена концептуальная модель района, характеризующая возможную природу аномальных зон. Даны приближенные оценки пористости пород в разломной зоне, по которой магматические расплавы поступали в вышележащие толщи в районе прорыва им. С.И. Набоко.

Рассмотрены методика и результаты магнитотеллурических зондирований в районе Авачинско-Корякской группы вулканов. Геоэлектрический разрез изучен в диапазоне периодов 0.0001 – 1000 с и более. Выполнено численное двумерное моделирование. Предварительно с помощью пробных моделей изучены возможные искажения кривых зондирований. По данным качественного анализа магнитотеллурических параметров определен характер геоэлектрических неоднородностей. В качестве основных приняты кривые по простиранию и вкрест простирания структур Камчатки. Продольные кривые, в меньшей степени подверженные влиянию берегового эффекта, использованы совместно с поперечными кривыми для создания геоэлектрической модели с помощью численного двумерного моделирования магнитотеллурического поля. Полученная геоэлектрическая модель содержит в верхней части разреза проводящий слой, связанный с осадочно-вулканогенным чехлом. Глубинная часть модели включает субвертикальные проводящие зоны, отражающие зону глубинных разломов. Рассматривается возможная природа выявленных аномалий и приближенная оценка пористости пород в проводящих зонах.

The article presents the methods and results of the magnetotelluric sounding within the Avacha-Koryaksky group of volcanoes. Geoelectrical section was studied within the period range from 0.0001 to 1000 seconds and above. The authors performed a numerical two-dimensional modeling. Initially, we used test models for possible distortions of curves. The analysis of the magnetotelluric parametres allowed us to characterize the geoelectrical inhomogenuities. Curves along the strike and across the strike were used as main curves. Since longitudinal curves are less prone to coast effect, they were used with transverse curves in order to create a geoelectrical model based on a 2D magnetotelluric field numerical modeling. The created geolectrical model has a conductive bed in the upper part of the section that is connected with an igneous-sedimentary cover. The deep part of the model includes near-vertical conductive zones, which denote a zone with deep faults. The paper describes possible nature of the revealed anomalies and provides rough estimation of rock porosity in the conductive zones.