Geodynamics as wave dynamics of the medium composed of rotating blocks Текст научной статьи по специальности «Науки о Земле и смежные экологические науки»

GEODYNAMICS & TECTONOPHYSICS

2015 VOLUME 6 ISSUE 3 PAGES 345-364

ISSN 2078-502X

http://dx.doi.org/10.5800/GT-2015-6-3-0185

Geodynamics as wave dynamics of the medium composed of rotating blocks

А. V. Vikulin

Institute of Volcanology and Seismology, Far East Branch of RAS, Petropavlovsk-Kamchatsky, Russia

Abstract: The geomedium block concept envisages that stresses in the medium composed of rotating blocks have torque and thus predetermine the medium's energy capacity (in terms of [Ponomarev, 2008]). The present paper describes the wave nature of the global geodynamic process taking place in the medium characterized by the existence of slow and fast rotation strain waves that are classified as a new type of waves. Movements may also occur as rheid, superplastic and/or superfluid motions and facilitate the formation of vortex geological structures in the geomedium.

Our analysis of data on almost 800 strong volcanic eruptions shows that the magma chamber’s thickness is generally small, about 0.5 km, and this value is constant, independent of the volcanic process and predetermined by properties of the crust. A new magma chamber model is based on the idea of 'thermal explosion’/'self-acceleration' manifested by intensive plastic movements along boundaries between the blocks in conditions of the low thermal conductivity of the geomedium. It is shown that if the solid rock in the magma chamber is overheated above its melting point, high stresses may occur in the surrounding area, and their elastic energy may amount to 1015 joules per 1 km3 of the overheated solid rock. In view of such stresses, it is possible to consider the interaction between volcano’s magma chambers as the migration of volcanic activity along the volcanic arc and provide an explanation of the interaction between volcanic activity and seismicity within the adjacent parallel arcs.

The thin overheated interlayer/magma chamber concept may be valid for the entire Earth's crust. In our hypothesis, properties of the Moho are determined by the phase transition from the block structure of the crust to the nonblock structure of the upper mantle.

Key words: geodynamics, force moment, rotational waves, rheid flow, gravitational waves.

For citation: Vikulin A.V. 2015. Geodynamics as wave dynamics of the medium composed of rotating blocks. Geodynamics & Tectonophysics 6 (3), 345-364. doi:10.5800/GT-2015-6-3-0185.

Институт вулканологии и сейсмологии ДВО РАН, Петропавловск-Камчатский, Россия

Аннотация: В работе развивается концепция блоковой геосреды применительно к геодинамическому (сейсмическому и вулканическому) процессу. Работу, в целом, можно представить в виде четырех частей: введения, шести разделов, кратко описывающих разработанную ранее автором ротационную модель сейсмотектонического процесса, трех разделов, посвященных развитию блоковой концепции применительно к вулканическому процессу, и обсуждения результатов.

Во введении на примере анализа исследований, проводимых в течение последних десятилетий сотрудниками Института земной коры СО РАН, показывается, что концепция тектонофизического процесса Байкальской рифтовой зоны, разрабатываемая в рамках разломной тектоники на региональном уровне, не позволяет видеть всю картину в целом. Региональные концепции существенным образом опираются на представления

Геодинамика как волновая динамика блоковой

ВРАЩАЮЩЕЙСЯ СРЕДЫ

А. В. Викулин

345

Recent Geodynamics

A.V. Vikulin: Geodynamics as wave dynamics...

о разломах, представления о блоковой среде используются чисто формально. И, несмотря на большие успехи [Sherman, 2014], уже на первых этапах построения региональной модели ученые вынуждены вводить взаимосвязи между ее параметрами, тем самым резко ограничивая возможности интерпретации модели на заключительных этапах исследования. Отмечается, что принцип Сен-Венана в применении к задачам сейсмологии и геодинамики не может рассматриваться как фундаментальный.

В разделе «Напряжения с моментом силы» проводится построение механически очевидной модели движения блока, являющегося частью вращающейся среды - геосреды. Показывается (рис. 1), что в блоковой вращающейся и передвигающейся вдоль поверхности Земли геосреде генерируется упругое поле с моментом силы, которое действует на блоки через их поверхности. Такие свойства упругого поля являются следствием закона сохранения момента количества движения. Движение блока во вращающейся системе координат механически эквивалентно движению блока в невращающейся (инерциальной) системе координат под действием собственного момента силы (спина), который в окружающем блок пространстве создает упругое поле с моментом силы. Такие напряжения в геосреде накапливаются, что и может объяснить ее известное свойство - энергонасыщенность [Ponomarev, 2008].

В разделе «Близкодействие и дальнодействие ротационного упругого поля» показывается, что поле упругих напряжений с моментом в геосреде описывается симметричным тензором напряжений. Оно характеризуется близкодействием - «моментным» взаимодействием рядом расположенных блоков, и дальнодействием -«энергетическим» взаимодействием всех блоков сейсмического пояса, протягивающегося на десятки тысяч километров, что может являться отражением общего физического принципа - корпускулярно-волнового дуализма. В ротационной концепции не требуется привлекать модель Коссера, которая является математической, не физической.

В разделах «О ротационных волнах в блоковых вращающихся геосредах» и «Новый тип геодинамических колебаний» описывается разработанная ранее автором в рамках блоковой концепции геосреды модель ротационного сейсмотектонического процесса [Vikulin, 2011]. Характерными скоростями волновой модели являются с0 = y^VRVs, с0 ~ 1 — 10 см/с (8) и (9) и 1^1-10 км/с (10). Значение первой с0 - «солитонной» предельной скорости - определяется только угловой скоростью вращения Земли вокруг своей оси, откуда и название модели - ротационная. Вторая Vex - «экситонная» предельная скорость - является скоростью упругих сейсмических волн. Отмечается, что такие же, по сути, деформационные волны [Khachai О.А., Khachai О.Yu., 2012; Khachai et ai., 2013] инструментально зарегистрированы в шахтах. Проводится сопоставление теоретических модельных решений с данными о скоростях миграции очагов землетрясений. Отмечается и качественное и количественное совпадение с такими же и решениями и значениями характерных скоростей, соответственно. Формулируется вывод о том, что ротационные волны - это новый тип волн, являющихся для блоковых вращающихся сред такими же характерными, как и продольные и поперечные сейсмические волны.

В разделе «Реидные свойства геосреды» отмечено, что медленные ротационные волновые движения могут быть ответственными за сверхпластичные течения [Leonov, 2008] геосреды - ее реидные [Carey, 1954] свойства [Vikulin, Ivanchin, 2013а]. Отмечается, что физическим аналогом такого движения может являться сверхтекучесть.

В разделе «Ротационные и маятниковые волны» показывается, что описание волнового геодинамического процесса в рамках ротационной модели близко такому же, по сути, блоковому описанию с позиции концепции «композиции и декомпозиции вещества Земли» [Oparin et ai., 2010]. Это позволяет ротационные и маятниковые, а также и деформационные [Khachai О.А., Khachai ОУи, 2012; Khachai et ai., 2013] волны отнести к одному классу явлений - взаимодействию блоков геосреды между собой посредством упругого поля с моментом силы.

В разделе «Блоковая геосреда и вулканизм» отмечается, что отражением блокового строения среды может являться магматический очаг, питающий извержения действующего вулкана.

В разделе «О параметрах магматических очагов» с использованием большого объема данных об извержениях вулканов планеты показывается, что толщина магматического очага является инвариантом наиболее общих статистических распределений, характеризующих вулканический процесс (рис. 3 и 4), и, как следствие, мало изменяющейся величиной 4h=0.5±0.1 км (19). Это позволяет толщину магматического очага рассматривать как постоянную величину, не зависящую от вулканического процесса, которая определяется геодинамическими параметрами земной коры, ее блоковым строением. При размерах кальдер до 0=10-15 км и более форма очага близка блинообразной: D>>Ah.

В разделе «Магматический очаг как состояние земной коры» для блоковой геосреды, вдоль границ которой возможны интенсивные пластические деформации [Magnitskii et ai., 1998; Turcotte, Schubert, 1985], с использованием разработанной в материаловедении концепции «теплового самоускорения» [Ivanchin, 1982] проводится построение модели «тонкого» (блинообразного) магматического очага. В пределах земной коры в результате малой теплопроводности горных пород и интенсивной пластической деформации вдоль границы залегания кристаллического фундамента на глубине 6 км (Камчатка) может реализоваться состояние перегретого выше точки плавления вещества, находящегося в твердой фазе. При локальных плавлениях в очаге и увеличении его объема в окружающей очаг земной коре создаются упругие напряжения, величина которых может достигать 1015 Дж на 1 км3 перегретой породы. Такие напряжения сопоставимы с напряжениями в очагах сильнейших тектонических землетрясений с магнитудами около 8 и более. «Энергетическая» близость магматических и сейсмических рядом расположенных очагов в рамках модели блоковой гесореды позволяет объяснить и взаимодействие вулканов между собой (миграцию вулканической активности [Vikulin et ai., 2010, 2012, 2013]), и взаимодействие вулканизма, сейсмичности и тектоники [Vikulin, 2011; Vikulin et ai., 2013]. Развитые представления распространяются на границу Мохо. Формулируется гипотеза, согласно которой подошва земной коры может представлять собой фазовую поверхность, ниже которой геосреда не является блоковой, способной к сдвиговому течению.

346

Geodynamics & Tectonophysics 2015 Volume 6 Issue 3 Pages 345-364

В разделе «Обсуждение результатов» рассматривается ряд возможностей ротационной модели. Проблема взаимосвязи энергонасыщенности, реидности и способности двигаться вихревым способом, с одной стороны, и сильной нелинейности геосреды - с другой, сформулирована как фундаментальная задача геодинамики и тектонофизики. Свойство энергонасыщенности геосреды, ее нахождение в напряженном состоянии [Rykunov et al., 1979], способном ее разрушить [Bogdanovich, 1909; Ponomarev, 2008], показывают, что землетрясение, скорее всего, происходит не в соответствии с теорией Ф. Рейда, т.е. не в результате создания локальных напряжений в очаге будущего землетрясения и преодоления предела прочности горных пород. Землетрясение в рамках развиваемой автором ротационной концепции является результатом дальнодействующего взаимодействия всех блоков земной коры и создания в области очага будущего землетрясения и прилегающих к ней блоков условий для их близкодействующего взаимодействия, которое может сопровождаться образованием свободной поверхности разрыва и излучением сейсмических волн. В рамках ротационной концепции механизм «зацепления» блоков и плит друг за друга и «выделения» тепла за счет трения их границ становится маловероятным.

Ключевые слова: геодинамика, момент силы, ротационные волны, реидное течение, гравитационные волны.

1. Introduction

The Baikal region and its neighbouring territories are the best studied among many areas subject to traditional geodynamic studies, including comprehensive geological and geophysical observations with regards to tectonophysical concepts. The research team under the leadership by Prof. S.I. Sherman from the Institute of the Earth's Crust, SB RAS has conducted detailed studies in this region for several decades. They have proposed and developed a tectonophysical model of the deep structure of faults in Central Asia [Sherman et al., 1992, 1994]. Results of their pioneering studies include the following:

- Large faults in lithosphere extension zones are physically modelled, and quantitative parameters of

There is no “relativity of rotation.” A rotating system is not an inertial frame, and the laws of physics

are different.

R. Feynman [Feynman et al., 1964]

Относительности вращения не существует. Вращательная система - не инерциальная система, и законы физики в ней другие.

Р. Фейнман [Feynman et al., 1964]

Geological time scale is close to the scale of the Universe. Geologists own chronicle, which recorded the events of history of the Earth, as well as the Universe.

D.V. Nalivkin [Nalivkin, 1969]

Масштаб геологического времени близок к масштабу Вселенной. Геологи владеют летописью, в которой записаны события истории Земли, а также и Вселенной.

Д.В. Наливкин [Nalivkin, 1969]

deformation in zones of large faults are established [Sherman et al., 2001];

- An original geodynamic model is developed to show how rift basins in Pribaikalie and Transbaikalie develop in space and time [Lunina et al., 2009];

- A seismic zone is modelled in terms of tectonophy-sics [Sherman, 2009]; the model shows that sources of rare strong earthquakes are clustered in linear or arcshaped systems [Sherman, 2013] and migration velocities of earthquakes К>12 (M>4-5) range from 1 to 100 km per year [Sherman, 2013; Sherman, Gorbunova, 2008]; according to [Vikulin et al., 2012a, 2012b], these values do not contradict with the relevant global data;

- According to above-mentioned model, faults are activated by slow strain waves of excitation which are generated by interplate and interblock movements of

347

A.V. Vikulin: Geodynamics as wave dynamics...

the lithosphere [Sherman, Gorbunova, 2008] and can also occur in zones wherein seismicity is slowly migrating [Sherman, 2009, 2013; Sherman et al., 2011]; a general assumption is that L>l, wherein L is strain wave length and I is length of a fault activated by strain waves [Sherman et aI., 2008];

- It is revealed that two major fault zones in the Baikal rift zone and the Amur River region are related; active faults are identified in both regions; it is established that the fault activation is manifested by seismicity and caused by trigger mechanisms with the leading role of slow strain waves [Sherman, 2013; Sherman et al., 2011].

The above-mentioned tectonophysical model of the wide geodynamic zone, including Central Asia and the Amur River region, is based on the concept of the fault-block geomedium and the assumption that faults are activated in real time by slow strain waves, and the fault activation is accompanied by seismic events that sequentially take place along the faults. It is believed that developing the tectonophysical model of the seismic process can “give an insight to regularities of earthquake location in space and time and facilitate seismic forecasting” [Sherman, 2009]. The current model is based on commonly accepted geological, geophysical and tectonophysical concepts of the geodynamic process. In general, it seems promising and compliant with the major concepts that are internationally accepted in modern geoscience. Undoubtedly, this model is an important contribution to the Earth sciences.

The above-mentioned research results are consolidated in [Sherman, 2014]. S.I. Sherman and his colleagues have proposed and develop the tectonophysi-cal concept of the Baikal rift zone development. In his book, S.I. Sherman states that their concept does not depend on physical concepts of seismicity (p. 6) and is thus universally applicable.

However, a review of the tectonophysical model of Central Asia with respect to the planetary scale [Vikulin et al., 2012b] reveals a seemingly insignificant inconsistency and raises a number of important questions.

As mentioned above, the model estimates of velocities of earthquake sources migration along the faults are generally consistent with the global data. Nonetheless, the detailed analysis of migration velocities in geodynamic settings of all active zones of the Earth [Vikulin et al., 2012a] shows that this is not exactly the case. The model estimates of velocities of earthquake sources migration along the faults in Central Asia do not contradict with similar plots constructed for the Pacific margin and the Alpine-Hymalayan belt (i.e. compression zones), but are in conflict with the plots constructed for the Middle Atlantic ridge (i.e. extension zone). Thus, the earthquake migration estimates for Central Asia at the regional level are contradicting with

the estimates at the planetary level. Otherwise, it has to be accepted that the studied region of Central Asia is not a rift, or it should be clarified that the earthquake migration estimates for the Baikal region refer to one of its flanges rather than to the entire rift zone. This poses the questions: - Which flange in particular? -What is the difference between the flanges? The regional model fails to provide answers to these questions.

In the regional tectonophysical models, including the model of Central Asia, the concept envisaging the block structure of the geomedium is applied formally or 'rhetorically' - such models are typically constructed with reference to boundaries between the blocks and block length values rather than blocks themselves and their volumes, and the models thus consider waves propagating along the boundaries between the blocks (i.e. along the faults), but not the waves propagating inside the block medium. This terminological ‘swapping’ - ‘speaking about blocks, while thinking about faults’- is not just a habitual use of ‘fault' as a commonly accepted term with a ‘standard’ reference to the local stress accumulation mechanism (according to

H. Rheid) and its modifications. In the regional tectonic models, the accumulated energy is not released - it is redistributed in the form of slow strain waves, while in seismic models, the energy is released by slow strain waves. The notions are substituted ‘automatically’ or ‘by analogy’ as a result of ignoring the major consequence of the geomedium block concept which states that stresses are redistributed in the geomedium's volume (p. 193-197, 332-334 in [Sadovsky, 2004]). The idea of stress redistribution is supported by detailed seismic monitoring data obtained recently in mining projects in Russia [Mulev et al., 2013]: “While a combined machine is moving inside the mine, the seismic energy release plot is changing; once the machine is stopped, the corresponding isolines in the plot are immediately ‘frozen’ in real time”. In other words, even when an insignificant rock volume is ‘extracted’ from the massif, stresses are significantly redistributed1. The 'volume' mechanism of stress redistribution in the geomedium is in conflict with the tectonophysical (actually, 'fault-based') interpretation of the data obtained for Central Asia. In fact, an obscure issue is how a nonradiating fault, being only activated [Sherman, Gor-

It follows from the instrumental seismological data published in [Mulev et al., 2013] that the Saint-Venant's principle (being local, according to [Sedov, 1973], p. 364-365) is not applicable to solving the geophysical problems. Seismologists and geophysicists consider the Saint-Venant's principle as fundamental (page 11 in [Rizni-chenko, 1985], though a sufficient justification of this approach is lacking. This problem is discussed in detail in [Vikulin, 2014a]. Below in the present paper, a geodynamic justification of the Saint-Venant's principle is provided under the rotation concept, and it is suggested that long-range interaction takes place between the geomedium blocks.

348

Geodynamics & Tectonophysics 2015 Volume 6 Issue 3 Pages 345-364

bunova, 2008], can 'be aware' of the length of the wave passing through this fault.

The geomedium block structure is a challenging subject for both the Earth sciences and practical applications, and the key problem is to specify what part of the geomedium should be regarded as a block and what part is a boundary of this block. At the conference held in Irkutsk in August 2014, the report on this problem [Tveritinova, Vikulin, 2014] sparked off a dispute that was lively, although not completed by any specific conclusion, and the participants agreed that the definitions need to be sorted out.

The seismic and geodynamic setting of the entire Baikal-Amur River zone as a global interplate boundary [Sherman et al., 2011] has not been fully clarified yet, and the knowledge is mainly based on the regional and hypothetical conclusions that are 'cross-linked' only with regard to some separate ideas. With account of the current uncertainties, ‘paving a pathway to regional earthquake forecasting’ [Sherman, 2009,

2014; Sherman, Gorbunova, 2008; Sherman et al., 1992, 1994, 2001, 2011] may appear not that straight forward.

The 'regional' approach fails to provide a comprehensive view, and relationships between model parameters have to be introduced to design a model, which means setting up strong limitations on potential interpretations of consequences following the model at the final stage of modelling which is most important for geodynamics. A more general approach can facilitate finding original pathways to reviewing and solving geodynamic problems in the Earth sciences [Vikulin et al., 2012a, 2012b].

2. Stresses with force momentum in the geomedium

COMPOSED OF BLOCKS

An important scientific achievement of the past decades is the justification of the hypothesis of the block structure of the geological medium [Peive, 1961] / geomedium (p. 5-20 in [Nikolaev, 1987]) / geophysical medium (p. 332-334 in [Sadovsky, 2004]).

The Earth crust is subject to motions and changes as the crustal blocks and plates migrate on the Earth surface (Fig. 1). Properties of the crust are predetermined by the Earth rotation and its structure composed of the blocks that are subject to translational motions at the Earth surface [Vikulin, 2010, 2011; Vikulin, Ivanchin, 2013а].

The coordinate system, that is fixed for the body (i.e. the Earth), rotates with angular velocity fl that does not depend on the coordinate system as it is - in the given moment of time, all systems of the kind are rotating around axes (that are parallel to each other) with absolute velocity fl [Landau, Lifshitz, 1973]. Thus, re-

gardless of its size, any block/plate is characterized by impulse momentum М that is independent of the block/plate size and origin and directed parallel to the body's (i.e. Earth's) rotation axis: М = Ifl, where I is inertia momentum of the block/plate. As shown in Fig. 1, a, crustal movements from position М1 to position М 2 lead to changes in the direction of momentum, М^М 2 because of the rigid relationship between the block and the Earth that rotates with angular velocity fl. Under the impulse momentum preservation law, when the direction and, consequently, value of М is changed, the block is subject to force momentum K.

The value and direction of force momentum K can be estimated theoretically as follows. It is assumed that the block is a uniform ball-shaped body that is stopped in position М 2 and subject to elastic stresses with force momentum Р2 (see Fig. 1, b). Then elastic stresses with force momentum Р1 are applied to rotate the block to its initial state with momentum М1. With the assumption that kinetic energy of the block's rotation is converted to elastic stresses and vice versa without any loss of energy (| Р11 = | Р21 = |Р|), the cosine law is applied to force momentum K:

\K\= 2|P|sin£/2. (1)

It should be noted that elastic stresses with force momentum K are applied to the block across its surface from the side of its surrounding crust and lithosphere [Vikulin, Ivanchin, 2013а].

Therefore, in the rotation model [Vikulin, 2009, 2010, 2011; Vikulin, Ivanchin, 2013а; Vikulin et al., 2010,

2013], the block's movement with angular velocity fl is mechanically equivalent to its movement in the non-rotating/inert system of coordinates under the impact of its own force momentum K that generates an elastic field with force momentum (equation 1) in the medium surrounding the block (as a consequence of the impulse momentum conservation law).

The 'internal' and/or own [Peive, 1961] momentum М, that is actually a macrospin (p. 146-148 in [Sedov, 1973]), has a specific feature - under the impulse momentum conservation law, the geomedium cannot be 'deprived' of it by any means, including plastic deformation. Thus, as a result of the translational motion of the block along the Earth surface (due to increasing block rotation angle в and constancy of the angle a, Fig. 1), rotation stresses with force momentum (equation 1) occur in the crust. The accumulation of such stresses provides an explanation [Vikulin, Ivanchin, 2013а] of such a property of the crust as energy capacity (in terms of [Ponomarev, 2008]). In fact, a similar conclusion was drawn by geologists at the end of the 19th century and in the early 20th century on the basis of observations of auto-destruction of rock samples

349

A.V. Vikulin: Geodynamics as wave dynamics...

Fig. 1. While the crustal block moves from position Mi to position М2 (i.e. rotates by angle в) (о), stresses with force momentum K are 'generated' in the crust and the lithosphere, force momentum K is applied to the block from its surrounding medium (b). a - the angle between the momentum M and the vertical to the Earth surface, which remains constant during the block movement. See explanations in the text.

Рис. 1. Движение блока земной коры из положения с моментом импульса Mi в положение М2 (поворот блока на угол в) (о), сопровождающееся «генерацией» в земной коре и литосфере напряжений с моментом силы К, прикладываемым к блоку со стороны окружающей его среды (b). Пояснения в тексте.

extracted to the ground surface [Bogdonovich, 1909]. Besides, the fact that the Earth's crust is permanently subject to stresses is confirmed by the available geophysical data obtained by instrumental measurements [Rykunov et ol., 1979; Nikoloev, 2003; Chebotareva, 2011; Khachai О.А., Khachai O.Yu., 2012].

3. Short- and long-range impacts of the

ROTATIONAL ELASTIC FIELD

The rotation concept (p. 124-136 in [Vikulin, 2011]) is applied to solve the following two-point boundary-value problem: a ball-shaped block has radius R0; its own momentum is given by equation 1 for the inert non-rotating spherical system of coordinates (r, V 6 with the starting point (r=0) in the block's centre; angle Ф is calculated in the plane located perpendicular to the block's rotation axis (0=tc/2); displacements at the infinite are zero; the non-zero-value force momentum does not depend on the block's size. According to theoretical estimations for the geomedium surrounding the block (i.e. r > R0) [Vikulin, 2010], the stress field is estimated as follows:

= <Г(рГ = 4QR0r 3J|“ sin0sin^/2 (2)

(values of other components are zero). The stress field is symmetrical2. Its energy is W0 = 16 npQ2 Rq sin2 ^/2, and the force momentum/ seismic momentum is

K0 = —8n3/2QRqsin^/2, where p=3 g/cm3 and

G=10n n/m2 are density and shear modulus of the crust, respectively, П=7.3-105 rad/sec is angular velocity of the Earth rotation around its axis, R0 is typical size of the block/earthquake source with М ~ 8, and suffix number (0) refers to the block's 'own' energy and torque.

In [Vikulin, 2010, 2011], the model is applied to study the interaction between ony two blocks with sizes R01 and R02, which are located at distance l from each other. The interaction energy and force momentum are

A fundamental role of 'stress tensor symmetry' for mechanics of continua, including the geomedium, is discussed in [Vikulin, Ivanchin, 2013a]. They also draw attention to the fact that the nonphysical character of the model proposed by Cosserat E. and F. [Cosserat E. et F., 1909] was noted immediately after its publication (p. 26 in [Hirth, Lothe, 1968]).

350

Geodynamics & Tectonophysics 2015 Volume 6 Issue 3 Pages 345-364

estimated as follows: Wint = 2npQ.2R0j1R0j2l 3cos$ and Kint = — 2npQ2RQ1Rg2l-3sin^>, respectively. In case of

equal-size blocks (Roi =Ro2), interaction energy Wint and force momentum Кш become similar to their own energy W0 and own force momentum K0 (Wint=W0, and Kint=K0 [Vikulin, 2010]) at distances W and I0K estimated in [Vikulin, Ivanchin, 2013а] as follows:

10Wf »102R 0, (3)

10Kя R0. (4)

Equation (3) shows that the interaction between the blocks, which controls their energy exchange, is manifested across large distances exceeding the block sizes and can be thus termed the long-range3 interaction. In seismology, such interaction is well-known and evidenced by the migration of earthquake sources along seismic belts for long distances amounting to dozens of thousand kilometres [Vikulin et al., 2012a, 2012b], as well as by remote fore- and aftershocks [Prozorov, 1978] and coupled earthquakes (p. 119-123 in [Vikulin,

2011]).

The instantaneous interaction between the blocks controls their momentum exchange and, according to equation (4), propagates to small distances not exceeding the block sizes. It can thus be called the short-range interaction. Such interaction is also well-known and evidenced by strong earthquakes, including coupled and multiple earthquakes with long-length sources (1000 km and more). Their sub-sources (with lengths ranging from 100 to 300 km) are involved in the instantaneous interaction and generate intensive own oscillations of the planet (p. 244-258 in [Vikulin, 2011]).

The above-mentioned theoretical solutions obtained for stress field o, elastic energies W0 and Wint and mo-mentums К0 and Кш are applicable to any block and any couple of the geomedium blocks. The summary of the stress fields may show the planet's elastic field that is 'self-congruent' and predetermines a rotation component in motions of each geomedium block. The stress fields may vary by scale and have variable complex shapes in various geodynamic polygons due to the irregular structure of the geomedium and different modes of the geomedium block motions. In view of the above, the block rotation pattern in Central Asia is chaotic and multi-directional (Fig. 2). In other regions studied in other scales with insignificant local inconsistencies, stress field patterns may seem regular and manifested in geophysical fields by vortex, ring-shaped

and other non-linear geological structures (see Figures 1, 2, 7 and 8 in [Vikulin, 2010]).

In physics, short- and long-range actions are often related to corpuscular and wave interactions taking place across boundaries of particles and across the medium wherein the particles are located, respectively. Under the geomedium block concept, the blocks can be viewed as 'elementary' particles. Therefore, in terms of physics, the geodynamics of block interactions in the rotation model can be viewed as a manifestation of the general physical principle of corpuscular-wave dualism envisaging that movements of geophysical blocks, tectonic plates and geological structures have both corpuscular and wave features, as discussed below.

4. Rotation waves in the geomedium composed of rotating blocks

For the block that generates an elastic field with force momentum (1) and interacts with elastic fields generated by other blocks of similar sizes in a mass chain, the movement law is established as the sine-Gordon equation (p. 85-95 in [Vikulin, 2011]). In this case, a seismic belt is modelled by a unidimentional chain of the crustal blocks/earthquake sources interacting with each other. Each block is characterized by inertia momentum I and volume V = 4/3(nR3).

The block motion can be given by the following д2 В

equation: I-^ = K0 + K1, where K0 is force momentum

corresponding to the elastic stress field generated by an individual block; K1 is force momentum controlling the interaction of the block with other blocks included in the chain. It is generally assumed that K1 is proportional to both the elastic energy accumulated due to

д2 В

motions of the corresponding block (V r^, where z is

distance along the chain of mass/blocks) and the elastic energy of other blocks included in the mechanical chain. The block motion can be thus given by the following dimensionless equation:

d2e

d2e .

П =sin0-

(5)

where В = fi/2, % = k0z and r) = v0k0t are dimensionless coordinates, and t is time. Taking into account that the wave length is close to the block size (A,«R0), and wave number is k0 = 2n/R0, the following equation is obtained for process development velocity v0 typical of the motion (equation 5):

3 As noted above, with account of the long-range interaction between the geomedium blocks, the local Saint-Venant principle becomes inapplicable to solving the problems of geodynamics. Flexural-torsional strain problems cannot be solved on the basis of this principle either (p. 33 in [Bezukhov, 1953]).

Vo

15

8п245л

QRo I-

^2vrvs = 0.27^, (6)

351

A.V. Vikulin: Geodynamics as wave dynamics...

Fig. 2. GPS data obtained at 323 observation points in Central Asian polygon in the period from 1995 to 2005 [Kuzikov, Mukhamediev, 2010].

Dots show the observation points. a - displacement vectors of the entire polygon; b - differentiated rotation of blocks (viewed as rigid bodies), according to data from clusters of observation points characterised by zero horizontal displacements against each other.

Рис. 2. Данные 323 пунктов GPS-наблюдений на Центрально-Азиатском полигоне в 1995-2005 гг., обозначенных точками [Kuzikov, Mukhamediev, 2010].

a - вектора смещений всего полигона, в целом; b - дифференцированные вращения блоков как жестких образований, в каждом блоке выявленные по таким совокупностям наблюдательных пунктов, имеющих нулевые горизонтальные перемещения друг относительно друга.

where VR = QR0, and Vs ~ 4 km/sec is velocity of a transversal seismic wave.

The form of equation (5) is predetermined by equation (1) used to estimate the force momentum of the elastic field generated by an individual block. It is thus assumed that sine-Gordon equation (5) and its right part (sinff) can be viewed, in terms of mathematics, as the law of impulse momentum preservation in wave processes. This is a principal assumption - when applied to the rotation problem of interrelated blocks arranged in a chain, it allows us not to consider any interaction between blocks due to friction along their boundaries, while the interaction needs to be considered in the elasticity theory (as shown, for example, in [Nikolaevsky, 1996]). Logically, this approach facilitates a physically 'transparent' interpretation of the typical velocity (equation 6) of the geophysical process described by sine-Gordon equation (5), if the rotation problem of stress fields around the block, that is rotating due to its own momentum, is solved under the classical theory of elasticity [Landau, Lifshitz, 2003] with a symmetric stress tensor (equation 2).

Equation (6) shows that if physical parameters G, p and Rо are fixed, velocity v0 depends only on angular velocity Q (via inertia momentum I and volume V of the

Ы.-ШИ

block). This means that the actual cause of such strain is Earth rotation (p. 237-243 in [Vikulin, 2011]). This is why the model introduced in [Vikulin, 2011] is called "rotation model". Given the above-mentioned parameters of the crust, the typical velocity is v0 = 10 — 102 m/sec.

In this case, we analyse a chain of irregularly rotating blocks characterised by deviations of force momentums that diverge from those in equilibrium positions (u). In general, such a chain corresponds to the seismic process in nature. Friction af along the block boundaries is taken into account, though it is viewed not as a mechanism of interaction between the blocks due to their 'hitch-up' to each other (as envisaged in the elasticity theory), but as a dissipative factor that hinders the rotation interaction between the blocks. Based on the above, the law of block motions in the chain can be thus shown by a modified sine-Gordon equation (p. 237243 in [Vikulin, 2011]):

0 — 0 = sin0 + af % + p8 ()sine. (7)

A numerical solution of equation (7) is found by the McLaughlin-Scott method. In this case, 5(<f) is the Dirac

Geodynamics & Tectonophysics 2015 Volume 6 Issue 3 Pages 345-364

function. In the model, initial conditions are given with regard to an average strain velocity in seismically active regions, and coefficients of friction ay- and inhomogeneity p correspond to those of real faults. If blocks/earth-quake sources are interacting mainly due to creeping (i.e. in case of the slowly developing seismic process), an asymptotic value of rotation strain transfer velocity is c0 — 1 — 10 cm/sec (p. 237-243 in [Vikulin, 2011]).

Under the model envisaging that the non-linear geomedium is composed of blocks, it is acceptable that velocity {v0, c0}, that is typical of the transfer of soli-tonic rotation strain (i.e. stress with torque), can be estimated as follows [Vikulin, Ivanchin, 2013а]:

Co = Y^VrVs, Co - 1 — Ю cm/sec, (8)

where у = к-1 — 10~4 is non-linear parameter characterising a chain of blocks/cluster of earthquake sources that make up a seismic belt; the block/cluster is variable in size; its rotation is irregular due to friction; к — 104(103 — 105) is geomedium non-linearity coefficient that is equal to the ratio between elastic constants of the third order and elastic constants of the second order, i.e. linear moduli of elasticity [Nikolaev, 1987].

It is noteworthy that instrumental records in mines have also revealed 'slow' and 'fast' strain waves [Kha-chai et al., 2013] and oscillation waves (u-waves) [Oparin et al., 2010].

5. A NEW TYPE OF GEODYNAMIC OSCILLATIONS

In our study, we compare experimental data on the migration of earthquake sources [Vikulin et al., 2012a, 2012b] and theoretical solutions of the sine-Gordon equation (p. 124-136 in [Vikulin, 2011]) [Davydov, 1982] and find two solutions that correspond to limiting energies [Vikulin, 2010], i.e. slow solitons (sol) and fast exitons (ex) (in terms of [Davydov, 1982])4. The slow solitons control the global migration of earthquake sources within an entire seismic belt (i.e. ‘long-range action of energy’) with the following limiting velocity:

Vsol ^ c0 « 1 - 10 cm/sec. (9)

The fast exitons control the local migration of fore-and aftershocks in sources of strong earthquakes (i.e. ‘instantaneous’ short-range interaction) with a limiting seismic (s) velocity as follows:

cо < Vex ^ Vs« 1 - 10 km/sec. (10)

4 Definitions of exitons differ in the literature [Davydov, 1976; Zim-an, 1974, and others]. In this paper, solitons and exitons are viewed as excitations according to [Davydov, 1982].

The slow solitons (9) and fast exitons (10), as a reflection of the general physical principle of wave-particle duality, are, in fact, a new type of geodynamic fluctuations. They are characteristic of the block rotating medium /geomedium in the same way as longitudinal and transverse seismic waves for a ‘normal’ body.

Based on plate dimensions and movement velocities of plate boundaries, two relationships are established between plate boundary length L and its movement velocity V, L12(V) (p. 313-316 and 317-322 in [Vikulin,

2011] ). Relationship LgL1 — (0.7 ± 0.3)LgV1 is close to LgL — LgV that characterizes spreading and subduc-tion [Zharkov, 1983; New Global Tectonics, 1974] and actually corresponds to the ‘long-range’ soliton action (9). The second relationship LgL2 — (0.4 ± 0.2)LgV2 is close to LgL — 0.3LgV that determines fault activation in Central Asia [Sherman, Gorbunova, 2008]. In our opinion [Vikulin, 2010], it corresponds to ‘instantaneous action’ (10). Tectonic solitons and exitons, as excitations corresponding to relationships L12(V), control changes of the tectonic activity at the Earth surface, while seismic disturbances (9) and (10) predetermine the seismic activity migration along the belts.

The conclusion about the existence of “a new type of solitary waves” with “limiting velocity values” is also stated for non-local elastic solids [Pamyatnykh, Ursulov,

2012] . The idea of a ‘slow’ mode “propagating with a velocity that is much lower than the velocity of sound in liquids, solid granules and gas” is theoretically and experimentally supported by results of theoretical and experimental studies [Rudenko et al., 2012]. “The slow dynamics” and its impact on “the elastic properties of materials” are specified in [Korobov et al., 2013]. The idea of the new type of rotation waves [Vikulin, 2010, 2011], as the consequence of the concept envisaging that the geomedium is composed of rotating blocks, is consistent with results of theoretical and instrumental acoustic studies of ‘normal’ solids [Korobov et al., 2013; Pamyatnykh, Ursulov, 2012; Rudenko et al., 2012].

6. Rheid properties of the geomedium

There is much evidence of movements at the Earth's surface from earthquake focal areas which are mani-festted by ‘humps’ and often accompanied by rotation [Vikulin, 2008]. Such a case is described in [Ionina, Ku-beev, 2013] - during the catastrophic earthquake in China on 16 December 1920, “a missionary had to put his legs wide apart like a drinker trying to remain steady on his feet as he experienced strong rotational movements underneath” (p. 124). In 1930s, geophysical and geological data on ‘slow’ movements of the geomedium (typically lasting for 10-1013 sec) were

353

A.V. Vikulin: Geodynamics as wave dynamics...

consolidated, and terms of 'rheid deformation’ [Geological Dictionary, 1978; Carey, 1953] and 'superplastic deformation’ of the Earth were introduced to denote “a flow of the material in the solid state” [Leonov, 2008]. It is shown below that the above-described state of the geomedium is a direct consequence of its motions due to rotation.

A review of definitions of the rheological properties of the Earth [Vikulin, 2009] shows that the Debye temperature 9d of the geomedium can be given as follows [Zharkov, 1983]:

0d * 10~3V(H)3Vp(H), (11)

where V is average velocity of excitations in the geomedium, cm/sec; p is density of the geomedium, g/cm3; Н is depth, km. If an average velocity (determined by longitudinal and transverse seismic velocities) for the lithosphere and upper mantle ranges from 1 to 10 km/sec, the Debye temperature is high enough. If Н=100 km, Qd * 660 °C * 1000 °K, and these values correlate well with the common model of the Earth physics [Zharkov, 1983].

The situation radically changes in case of transition to rotation mode c о (equation 8) that is determined by motions of the geomedium blocks, tectonic plates and geological structures in the aggregate5. According to equation 8, limiting velocity value c0 is at least five orders lower than the longitudinal and transverse seismic velocities, and the Debye temperature (equation 11) is thus very low:

0d * 10~2K.

This provides for potential quantum, frictionless superfluid motions of the geomedium [Vikulin, 2013],

i.e. its rheidity [Vikulin, Ivanchin, 2013а] and/or superplastic flow in the solid state [Leonov, 2008; Carey, 1953].

The Debye temperature is proportional to the maximum possible frequency of oscillations of the complete set of particles composing the medium [Ziman, 1974] or a set of mesovolumes of a solid body or all the geomedium blocks, tectonic plates and other geological structures of the Earth. According to [Vikulin, 2011] (p. 244-258), in case of the Earth, the frequency is given by Chandler frequency typical of the aggregate

5 The Earth crust is a solid body. Therefore, applying the Debye theory to geodynamic problems is reasonable in terms of physics, to prove a possibility of transition under the Debye theory to rotation mode C0, a transition from elastic phonons to tectonic solitons (equation 9). An elastic phonon is a quantum of the interaction between 'point-size' atoms (i.e. atoms of zero size, without elastic own momentum) in the crystal lattice. A tectonic soliton is a carrier of the geodynamic interaction between (rotating) blocks of the geomedium. It has a considerable size and, respectively, its own torque.

oscillation of all the blocks in the seismotectonic belt. Actually, the oscillation of the entire belt is determined by the energy of 'zero' oscillations, i.e. energy E о (see. Fig. 6 in [Vikulin, 2010]).

7. Rotation waves and oscillation (ц) waves

Under the geomedium block approach, considering fracturing and “composition and decomposition of the Earth material” as the major processes [Oparin, Vostrikov, 2010], oscillation (p) waves in the geomedium are viewed as waves that predetermine geodynamic processes [Oparin et al., 2010]. Velocities of both oscillation (p) and rotation waves are lower than velocities of longitudinal waves. To experimentally substantiate the existence of p-waves, oscillation is analysed in chains of rigid massive blocks that are analogous to the chain of blocks in the rotation model. Two types of waves are revealed in both the rotation model and the chains [Oparin et al., 2007]. According to estimations based on laboratory experiments [Oparin, Vostrikov, 2010], velocities of oscillation (p) waves amount to 102-103 m/sec and are close to seismic wave velocities. Oscillation (p) wave velocities measured on site range from 1 to 10 cm/sec [Oparin, Vostrikov, 2010] and are close to the typical rotation wave velocity, eg (equation 8).

As shown above, in the framework of the wave geodynamics, both the approach based on the rotation concept [Vikulin, 2009, 2010, 2011; Vikulin, Ivanchin, 2013а] and the approach considering “composition and decomposition of the Earth material” [Oparin, Vostrikov, 2010; Oparin et al., 2007, 2010] yield similar results, while being developed independently of each other. Therefore, the rotation and oscillation (p) waves can be considered as a phenomenon of one and the same class - interaction between the geome-dium/rotating medium blocks in the field of elastic stresses with a force momentum. The methods used to study the oscillation (p) waves [Oparin et al., 2007] and strain waves [Khachai et al., 2013] may prove useful for developing a technology to ensure proper recording of the rotation waves, which are currently detected (similar to tectonic waves [Bykov, 2005]) by indirect methods, and instrumental methods for this purpose are still lacking.

8. THE BLOCK GEOMEDIUM AND VOLCANISM

The above concepts of waves that occur during the seismotectonic process in the block geomedium seem to be applicable to volcanic processes. In fact, volcanic belts of the Earth, as well as seismically active zones are the largest linear structures of the planet. They occurred in the Early Cretaceous and developed in quite a

354

Geodynamics & Tectonophysics 2015 Volume 6 Issue 3 Pages 345-364

synchronous pattern (p. 460 in [Krasny, 2004]). In studies of the three most active volcanic belts of the Earth -Pacific, Alpine-Himalayan and Mid-Atlantic belts, it is revealed that the volcanic activity migration, as well as the migration of seismic/tectonic activity, is a manifestation of the wave geodynamic process [Vikulin et al., 2010,2012].

A reflection of the 'block character’ of the volcanic process and a ‘quantum’ characteristic of a single volcanic eruption is a magma chamber feeding the volcanic eruption which is an analogue of an earthquake in the seismic process. Several definitions of a magma chamber /reservoir/ are given in [Vlodavets, 1984]. According to the most general definition that is not contradicting with the others, a magma chamber is an isolated volume feeding the volcano [Geological Dictionary, 1978]. Processes taking place in the magma chamber have not been fully clarified yet [Vlodavets, 1984; Luchitsky, 1971; Macdonald, 1975]. Nevertheless, in case of a sufficiently intense volcanic eruption when a large volume of volcanic material is conveyed to the ground surface, the size of a magma chamber can be estimated from the size of a caldera, i.e. a large circularshaped area (with a diameter of 10 to 15 km and more) of the volcanic depression with sloping walls (that can be several hundred metres high) which formed due to eruption (p. 123 in [Vlodavets, 1984]).

Many notions in volcanology have not been clearly defined yet, and a precise definition of 'caldera’ is also lacking (p. 123-129 in [Vlodavets, 1984]). The available data on calderas are reviewed in many publications, including [Luchitsky, 1971; McDonald, 1975; Leonov, Grib, 2004; Laverov, 2005; Vikulin, Akmanova, 2014; Vikulin, Ivanchin, 2015]. In our study, without going into detail, we refer to the above-mentioned general definition of ‘caldera’.

Our attempt to ‘incorporate’ the volcanic process in the geomedium block concept under the rotational approach is described below.

9. The parameters of magma chambers

Plotting the numbers and sizes of volcanic eruptions of the recent Kuril-Kamchatka volcanoes. In the region under study, the most complete and regionally consistent database, including numbers of recent (past 200 thousand years) eruptions and sizes of volcanic forms, is available for the Kuril-Kamchatka arc [Vikulin, Akmanova, 2014]. Based on the data on 70 volcanoes that erupted 676 times (N) in the past 9.5 thousand years, the eruption repeatability is given in the range of volcanic explosivity indices 2<W<7 [Siebert et al., 2010] as follows [Vikulin et al., 2012a, 2012b]:

LgN=3.60-(0.48±0.06)W. (12)

With account of the relationship between the ‘energy’ characteristics of eruption and the volume of erupted material (W and V, respectively): W=LgV+5, [V]=km3 [Siebert et al., 2010], the eruption repeatability (equation 12) can be given in terms of ‘volume’:

LgN=1.15-(0.48±0.06)LgV, [V]=km3. (13)

By applying the mid-square method to analyse the data from [Laverov, 2005] on recent volcanic forms, including calderas and large craters on cones and flat tops of underwater volcanoes (N=287), the following distribution of volcanic form numbers N by square areas S is yielded:

LgN=(2.32±0.16)-(0.47±0.14)LgS, [S]=km2. (14)

The above equations show that the ‘inclination angle’ of the ‘energy’ and/or ‘volume’ distribution (W and V in equations 12 and 13, respectively) amounts to 0.48±0.06, which is practically coincident with the ‘inclination angle’ of the distribution of the volcanic form numbers by their square areas S (equation 14): 0.47±0.14. It should be noted that the volcanic eruptions number (N=676) and the volcanic forms number (N=278) are statistically significant. The energy range and the sizes of calderas in equations (12), (13) and (14) cover the whole spectrum of volcanic eruptions, including those with maximum values, Wmax=7, Vmax~100 km3, and Smax=20x25 km2.

Essentially, equations (12), (13) and (14) show fairly general statistical distributions that are typical of the volcanic process taking place in a wide region within a large time period. The similar inclination angles of the energy (12) and/or volume (13) and square area (14) distributions suggest the following hypothesis: the ratio of erupted material volume V to area S of the resultant volcanic form (that is large enough) is constant:

V/S=Ah=const. (15)

Being actually an invariant of more general statistical distributions, hypothetical equation (15) has a fundamental importance for volcanology [Vikulin,

Akmanova, 2014].

Eruptions of volcanoes of the planet. The inclination angle of the volcanic eruption recurrence plot [Vikulin et al., 2012a] is given as follows:

LgN«-(0.52±0.05) W. (16)

It is equal to the inclination angle of a similar plot constructed for the Kuril-Kamchatka volcanoes (equation 12). With reference to the data on volcanoes of the planet [Vikulin, Ivanchin, 2015], the volcanic eruption numbers are plotted against the areas of resultant calderas and volumes of erupted material.

355

A.V. Vikulin: Geodynamics as wave dynamics...

I Fig. 3. Numbers of calderas (points) (N=373) and large craters on cones (triangles) (N=175) versus areas S (km2).

Рис. 3. Распределение чисел N=373 кальдер (точки) и крупных воронок N=175 на конусах (треугольники) по величинам их площадей S (км2), образованных на вулканах планеты.

Figure 3 shows a plot of numbers of calderas versus their areas: LgN=f(LgS), N=373. In Figure 4, numbers of eruptions are shown versus erupted volumes: LgN=f(LgV), N=125. The mid-square method is applied to analyse relationships shown in Figures 3 and 4 and given by equations (17) and (18), respectively (see below); correlation coefficients are 0.9 and 0.7, respectively. Intervals for averaging the initial values of S and V are increased with increasing values of S and V in order to provide a sufficiently uniform spacing of averaged values in the logarithmic scale.

LgN=(1.86±0.11)-(0.49±0.05)LgS, [S]=km2, (17)

I Fig. 4. Numbers of eruptions versus material volumes V (km3) erupted by volcanoes of the planet (N=125). I

I Рис. 4. Распределение чисел извержений как функция объемов выброшенного материала V (км3), N=125, при извержениях вулканов планеты.

LgN=(1.41±0.19)-(0.42±0.09)LgV, [V]=km3. (18)

Parameters of calderas - diameter D, area S, and erupted material volume V - are variable in wide ranges: D=2^150 km, S=2^4.648 km2, V=0.3-3000 km3. Nonetheless, it is established that the inclination angles of the plots showing the distribution of areas S (equation 17) and erupted material volumes V (equation 18) are similar, and the statistically significant reliability is high, minimum 0.7. This finding is actually similar to the conclusion that the ratio is constant in a wide range of scales of V and S.

Thus, the available data on the volcanic eruptions in the past 33 Ma support the hypothesis that Ah=const (equation 15), i.e. does not depend on parameters of the volcanic process.

Thickness of a magma chamber. Based on the above, it is possible to estimate a magma chamber thickness from parameters of a volcanic eruption. For this purpose, equation (15) is used with data on volcanic eruptions with known values of S and V (N=125). The volumes of volcanic eruptions of the planet are recalculated by the method proposed by I.V. Melekestsev, taking into account the fact that densities of magma at different focal depths and densities of eruption products on the ground surface are different [Laverov, 2005]. It is shown that equation (15) is valid [Vikulin, Akmanova, 2014]:

Ah=0.5±0.1 km. (19)

Minimum and maximum values of magma chamber thickness estimated from data on 125 volcanic eruptions (each with known caldera area S and erupted material volume V) range from a few meters to a few kilometres. These values do not contradict with thicknesses of molten granitoid layers formed while the rocks were irreversibly deformed during folding of the crust in the Pamirs [Magnitskii et al., 1998].

For a caldera with an average diameter D ranging from 10 to 15 km, a magma chamber which thickness is estimated from (19) is represented by a thin 'pancakeshaped' layer: D>>Ah.

Noteworthy are the following 'coinciding' features. Firstly, thicknesses of magma chamber Ah are close to heights of caldera sides h, and the minimum and maximum thickness values are generally close (Ah«h), which gives grounds to consider this fact in support of the hypothesis of a thin magmatic chamber with a permanent thickness. Secondly, as shown in Fig. 3, analytical equation (17) cannot be significantly changed by data on large craters on cinder cones (marked by triangles) (N=175, i.e. almost a half of the number of calderas, N=373). As such, distinctions between large craters on cinder cones and calderas are eliminated under the thin magma chamber concept, as

356

Geodynamics & Tectonophysics 2015 Volume 6 Issue 3 Pages 345-364

shown in equations (15) to (19). In other words, the very concept of a caldera-generating eruption viewed as a 'specific' class of eruptions becomes meaningless.

Equation (19) is based on analyses of different data sets including (1) data on the entire planet and individual regions, specifically the Kuril-Kamchatka arc, Kamchatka, Kuril Islands and individual volcanoes located in the Kuril Islands and Kamchatka; (2) data covering different time periods, specifically data on volcanic eruptions of the Earth during the period of 33 Ma, and recent volcanic eruptions in the Kuril-Kamchatka region; (3) data from regions with different geotectonic settings; and (4) data on different types of volcanoes and different eruptions. Thus, there are grounds to consider that the small thickness of the magma pocket in equation (19) is a constant that does not depend on the volcanic process [Vikulin, Akmanova,

2014].

10. A MAGMATIC CHAMBER AS THE STATE OF THE CRUST

The mode! of a 'pancake-shaped' thin magma chamber in the crust. Recent instrumental geophysical studies conducted at volcanoes in the Kuril-Kamchatka region show that their magma chambers are generally located at depths from 5 to 30 km [Anosov et al., 1990; Balesta, 1981; Fedotov, 1984; Ermakov, Shteinberg, 1999; Fedotov, Masurenkov, 1991], i.e. within the crust composed of blocks. Therefore, the magma chamber thickness, which is independent of the volcanic process, may depend only on geodynamic movements of the blocks.

The following data can provide prerequisites for a statement of the problem. Heat generation in the crust due to mechanical movements of the crustal blocks was described in general terms by P.N. Kropotkin [1948]. Potential melting of the material in fault zones due to dissipative heating during shifting of the crustal blocks was noted in many publications on geodynamic modelling (for instance, p. 308-310 in [Turcotte, Schubert, 1985]). The first quantitative model of granitoid melting in the crust as a result of its intensive deformation during folding (a case of the Pamirs) was proposed in [Magnitskii et a!., 1998]. It should be emphasized that under this approach it is not required to refer to the fluid concept in discussion of melting. The potential existence of polingenic magma chambers in the crust was substantiated in [Ermakov, 1977] and confirmed by results of instrumental geophysical observations [Balesta, 1981].

The thin ‘pancake-shaped’ magma chamber model is based on principles established in the material science of solids with intensive plastic deformation; in case of low heat transfer, the plastic deformation

zone can become considerably heated under the impact of such deformations and even destroyed [Ivanchin, 1982].

The main provisions of the model [Vikulin, Ivanchin,

2015] are as follows. Local plastic deformation s may be high (s«1) in sliding bands and low outside the bands (s«0). As the strain rate is exponentially dependent on stress and temperature, the temperature in the area subject to plastic deformation may increase due to heat emission. In case of a small thermal conductivity factor, when heating is not compensated, the temperature in the zone of intensive plastic deformation may continuously grow until transition to the thermal self-acceleration mode, which typically results in partial destruction of the deforming body, and in many cases such destruction results from melting.

A heated band may form under the following physically transparent conditions: on the one hand, deformation should be fast enough not to give time for the generated heat to be removed by thermal conductivity; on the other hand, it should be slow enough to provide a sufficient time for the plastic deformation rate and the specific heat capacity to decrease due to stress relaxation.

With account of the low thermal conductivity of the crust, thermodynamic calculations are obtained for the following conditions: a magma chamber is a solid that is heated to the temperature above its melting point; it is located at the depth of the crystalline basement boundary, Н=5-6 km (Kamchatka) and composed of aluminium which thermal properties are similar to those of magma. (All the thermodynamic parameters of aluminium and their dependence on temperature and pressure are known). Based on the thermodynamic calculations, it is shown that when local melting takes place and the magma chamber's volume is increasing, an additional pressure is generated around the magma chamber and the field of elastic stresses occurs; its energy is about 1015 joules per 1 km3 of the overheated rock. Considering the elastic energy accumulated around magma pockets, overheated magma chambers of relatively strong (W>4; V>0.1 km3) eruptions appear to be similar to sources of the strongest earthquakes (M=8 and above) which lengths range from 100 to 200 km and above. In the block geomedium model, the similarity of the neighbouring magma chambers and seismic sources in terms of energy can be viewed as an explanation of both the interaction between volcanoes (i.e. migration of volcanic activity according to [Vikulin et al., 2010,2012a, 2012b]) and the interaction between volcanism, seismicity and tectonics [Vikulin, 2011; Vikulin et al., 2013].

The above general conclusion is supported by data on the three adjacent large cones of the northern eruption of the Tolbachik volcano in Kamchatka, which

357

A.V. Vikulin: Geodynamics as wave dynamics...

formed one after another during the eruption in 1975 and 'migrated' along the entire Tolbachik zone of cinder cones [Fedotov, 1984]. It is also supported by data on calderas (D=10 km and above) that migrated along the 500-km long Snake River valley in the USA; the calderas were formed during eight strongest eruptions (W>6-7; V>10-100 km3) in the past 16 Ma [Koro-novsky, 2012], data on volcanoes that migrated along the central Pacific Ocean for distances from 800 to 1000 km and formed the chain of Hawaii Islands [Mcdonald, 1985], and data on volcanic activity migration at the Earth surface [Luchitsky, 1971].

The Mohorovicic discontinuity (Moho). It is stated in [Peive, 1961] that the “most important conclusion of modern tectonics, which forces us to review our concepts, is that the crust ... is partitioned into blocks not only by a system of steeply dipping and vertical tectonic surfaces but also by gently dipping and horizontal surfaces. ... Deep tangential tectonic zones located in areas of high pressures and temperatures are zones with ‘plastic flow’ ... and primary magma chambers. ... Mechanical movements of the crustal material are the main source of this energy”. According to [Peive, 1961; Magnitskii et al., 1998], mechanisms of plastic flow in thin layers along surfaces of the blocks are quite real and can be applied to other tectonic boundaries in the crust.

Among all the layers within the upper mantle, only the Moho is quite reliably established in all the regions of the Earth. Therefore, the concept of an overheated transition layer, which depth is constant, sounds reasonable with reference to the Moho, in the first instance.

It is not known whether the upper mantle located below the base of the crust is composed of blocks. We can only be certain that the crust has the block structure. Under the р-Т conditions at the Moho, it can be expected that the geomedium below the crust may not be composed of blocks. In this case, another reason against the block structure is the assumption that a volume flow [Leonov, 2008] is more typical of the substance below the Moho than a shear flow. Under the concepts of the rotational mechanics of the block geomedium [Vikulin, 2011; Vikulin et al., 2013] and the model of 'thermal explosion' and 'thermal self-acceleration' [Ivanchin, 1982; Vikulin, Akmanova, 2014; Vikulin, Ivanchin, 2015], it can be assumed that the base of the Earth's crust is represented by a phase surface, below which the geomedium (including the lithosphere and upper mantle) is not composed of blocks, and the shear flow cannot occur in it or may occur but to a lower extent than in the crust [Vikulin, Ivanchin, 2013b]. Our approach described above is close to the Moho model described in [Pavlenkova, 2013].

Thus, the conclusion, which is fundamental for both volcanology and geodynamics, is that a magma pocket

is a thin layer, and its thickness depends on properties of the crust. Our conclusion makes it possible within the framework of the rotation model [Vikulin, 2011; Vikulin, Ivanchin, 2013b; Vikulin et al., 2013] to consistently ‘mate’ the views on the volcanic process [Vikulin, Akmanova, 2014; Vikulin, Ivanchin, 2015], the block geomedium concept and the idea of the wave geodynamic process.

11. Discussion of results

Our review shows that the regional tectonophysical studies of the Baikal fault zone [Sherman, 2014] have some major shortcomings. In the reviewed regional studies and analyses, the fault structure of the geomedium was mainly referred to, while the block geomedium concept was not actually applied. The migration of earthquake sources was taken into account in the tectonophysical modelling of the regional processes, and an important conclusion was stated: “Deforma-tion/strain waves should be introduced in the scope of the geophysical medium concept and viewed as its dynamic parameter” (p. 283 in [Sherman, 2014]). However, in the cited publication, the strain wave process was considered mainly on the basis of the fault tectonics concept without any reference to the mechanism of stress redistribution in the geomedium's volume.

The present paper describes the rotation concept that can be useful for eliminating the difficulties associated with the need to take into account both the block structure of the geomedium and the wave process taking place in the geomedium. It is shown that the geomedium is a body composed of rotating blocks, and stresses with torque are generated in the geomedium. As a consequence of the momentum conservation law, such stresses predetermine the energy capacity (in terms of [Ponomarev, 2008]) of the geomedium. In our study, the new type of rotation waves is reviewed, and it is shown that slow and fast strain waves (and/or tectonic waves [Bykov, 2005] and/or oscillation (u) waves [Oparin, Vostrikov, 2010] and/or strain waves [Khachai et al., 2013]) can occur in the geomedium composed of rotating blocks as a consequence of the wave-particle duality. The fast and slow strain waves are as typical for the geomedium composed of rotating blocks as longitudinal and transverse seismic waves for the ‘normal’ solid.

In the geomedium, movements can also take place as rheid and/or super-fluid flow [Vikulin, 2013] and/or super-plastic flow [Leonov, 2008; Carey, 1954] which is most probably associated with the energy capacity of the geomedium and its ability to generate vortex, ring-shaped and other non-linear geological structures [Vikulin, 2014b]. In terms of physics (p. 5-20 in

358

Geodynamics & Tectonophysics 2015 Volume 6 Issue 3 Pages 345-364

[Nikolaev, 1987], this state of the geomedium corresponds to its geometric rotational-structural nonlinearity.

In geodynamics and tectonophysics, a fundamental challenge and experiment crucis, i.e. crucial/critical experiment [Vikulin, 2013, 2014a, 2014b], is to reveal the relationship between the geomedium's energy capacity, rheidity and vortex motion ability, on the one hand, and it's clear nonlinearity [Nikolaev, 1987], on the other hand.

Considering the geomedium's energy capacity [Ponomarev, 2008] and its state of stresses [Nikolaev, 2003; Rykunov et al., 1979], which can destroy the geomedium [Bogdanovich, 1909; Ponomarev, 2008], it can be suggested that an earthquake is most likely to occur out of compliance with the Raid theory, i.e. not as a result of local stresses in the future earthquake source in excess of the tensile strength of rocks. Under the rotation concept proposed by the author, an earthquake in the energy-saturated medium (ready for destruction) can result from the long-range interaction between all blocks of the crust through a long/geo-logical time span. Besides, it can occur as a consequence of specific conditions in the future earthquake source and its adjacent blocks which provide for the short-range interaction between the blocks and may be accompanied by the formation of a free surface of a fault and emission of seismic waves.

Magma chambers feeding volcanic eruptions are a reflection of the 'block character’ of the volcanic process. Our analysis of the data on almost 800 strong volcanic eruptions, that were followed by the formation of calderas and the transition of large volumes of volcanic material to the ground surface, shows that the magma chamber's thickness is generally small, about 0.5 km, and this value is constant, independent of the volcanic process and predetermined by properties of the crust.

We propose the new magmatic chamber model based on the ideas of ‘thermal explosion’ and ‘thermal self-acceleration’. It can supplement and develop the well-known concepts of potential heat generation in the crust [Kropotkin, 1948; Ermakov, 1977; Turcotte, Schubert, 1985; Magnitskii et al., 1998]. The model envisages intensive plastic movements along the boundaries of the blocks in conditions of the low thermal conductivity of the geomedium. If the solid material in the magma chamber is heated to the temperature above the melting point, elastic stresses with energy up to 1015 joules per 1 km3 of the overheated rock can be accumulated in the area around the magma chamber. Such stresses are comparable to stresses in the sources of the strongest earthquakes. This can provide an explanation of the interaction between volcanoes (manifested by the migration of volcanic activity along the volcanic arc) and the interaction between the volcanic

and seismic processes taking place in the adjacent and parallel belts.

Under the above concepts, the magma chambers are viewed as ‘normal’ blocks of the crust which have the same rotational properties as the blocks/earthquake sources.

This idea of thin overheated interlayers/volcanic chambers is valid for the entire crust of the Earth. In our hypothesis, it is stated that properties of the Moho are determined by the phase transition from the block structure of the crust to the non-block structure of the upper mantle.

12. Conclusions

Our conclusions [Vikulin et al., 2012b], that are based on the physical and mathematical models developed by the author and his colleagues, are of fundamental importance for geodynamics. A brief review of the conclusions is given below.

1. In geodynamics, a currently popular idea envisages the mechanism of 'hitching' of blocks and plates with each other and heat 'emission' due to friction at their boundaries. However, it becomes 'useless' under the rotation concept - this mechanism is unlikely to occur. It is suggested by models of ‘thermal explo-sion’/'self-acceleration’ and intensive deformation due to geological processes [Magnitskii et al., 1998] that overheated areas may occur in the crust, and phase transitions (from solids to liquids with free gas separation) may take place in such areas.

2. The concept of the geomedium composed of rotating blocks can facilitate solving the problems of geodynamics without any need to employ the currently popular models of magma ascent from the mantle (possibly, core) depths.

3. The problem of the Earth thermals and ‘hot spots’ can be considered from other positions. Within such zones, the kinetic energy generated by the rotation of the geomedium blocks and plates is partially converted into elastic stresses with torque (equation 1). It can be released not only by earthquakes, volcanic eruptions and tectonic plate movements, but also by the generation of heat that is redistributed inside the Earth and brought to the ground surface by various mechanism, including rotation waves with specific wave velocity eg (equation 8).

13. Acknowledgement

The author is grateful to Sh.A. Mukhamediev who thoroughly analysed the initial paper and thus encouraged the author to review and significantly improve it.

359

A.V. Vikulin: Geodynamics as wave dynamics...

14. References

Anosov G.I., Argentov V.V., Abdurakhmanov A.I. et al., 1990. The deep structure of the Uratman volcano. Vulkanologiya i Seismologiya (4), 85-91 (in Russian) [Аносов Г.И., Аргентов В.В. Абдурахманов А.И. и др. Глубинное строение вулкана Уратман // Вулканология и сейсмология. 1990. № 4. С. 85-91].

Balesta S.T., 1981. The Earth Crust and Magma Pockets in Recent Volcanism Areas. Nauka, Moscow, 135 p. (in Russian) [Балеста С.Т. Земная кора и магматические очаги областей современного вулканизма. М.: Наука, 1981. 135 с.].

Bezukhov N.I., 1953. The Theory of Elasticity and Plasticity. State Publishing House of Techn.-Theor. Literature, Moscow, 420 p. (in Russian) [Безухов Н.И. Теория упругости и пластичности. М.: Госуд. изд-во техн.-теор. лит-ры, 1953. 420 с.].

Bogdanovich R.I., 1909. Earthquakes in Messina and San Francisco. St. Petersburg, 165 p. (in Russian) [Богданович К.И. Землетрясения в Мессине и Сан Франциско. СПб., 1909. 165 с.].

Bykov V.G., 2005. Strain waves in the Earth: Theory field data, and models. Russian Geology and Geophysics 46 (11), 1176-1190.

Carey S.W., 1953. The rheid concept in geotectonics. Journal of the Geological Society of Australia 1 (1-2), 67-117. http://dx.doi.org/10.1080/14400955308527848.

Chebotareva I.Y., 2011. Methods for passive study of the geological environment using seismic noise. Acoustical Physics 57 (6), 857-865. http://dx.doi.org/10.1134/S1063771011060042.

CosseratE. etF, 1909. Theorie des corps deformables. Hermann, Paris. 226 p.

DavydovА.S., 1976. Physics of Solid Bodies. Nauka, Moscow, 640 p. (in Russian) [ДавыдовА.С. Физика твердого тела. М.: Наука, 1976. 640 с.].

Davydov A.S., 1982. Solitons in quasi-one-dimensional molecular structures. Soviet Physics Uspekhi 25 (12), 898-918. http://dx.doi.org/10.1070/PU1982v025n12ABEH005012.

Ermakov V.А., 1977. Formational Segmentation of Quaternary Igneous Rocks. Nedra, Moscow, 223 p. (in Russian) [Ермаков В.А. Формационное расчленение четвертичных вулканических пород. М.: Недра, 1977. 223 с.].

Ermakov V.А., Shteinberg G.S., 1999. The Kudryavy volcano and the evolution of the Medvezhaya caldera, Iturup Island, Kuril Islands. Vulkanologiya i Seismologiya (3), 19-40 (in Russian) [Ермаков ВА, Штейнберг Г.С. Вулкан Кудрявый и эволюция кальдеры Медвежья (о-в Итуруп, Курильские о-ва) // Вулканология и сейсмология. 1999. № 3. С. 19-40].

Fedotov S.A. (Ed.), 1984. The Major Fracture-Type Eruption of Tolbachik Volcano. Kamchatka. 1975-1976, Nauka, Moscow, 638 p. (in Russian) [Большое трещинное Толбачинское извержение. Камчатка. 1975-1976 / Ред. С.А. Федотов. М.: Наука, 1984. 638 с.].

Fedotov S.A., Masurenkov Yu.P. (Eds.), 1991. Active Volcanoes of Kamchatka. Two volumes. Nauka, Moscow, vol. 1, 320 p.; vol. 2, 415 p. (in Russian) [Действующие вулканы Камчатки. В 2-х т. / Ред. С.А. Федотов, Ю.П. Масуренков. М.: Наука, 1991. Т. 1. 320 с. Т. 2. 415 с.].

Feynman R.P., Leighton R.B., Sands M., 1964. The Feynman lectures on physics. Vol. 2. Addison-Wesley Publishing Company, Inc., Massachusetts, Palo Alto, London, 296 p.

Geological Dictionary, 1978. Nedra, Moscow, vol. 2, 456 p. (in Russian) [Геологический словарь. М.: Недра, 1978. Т. 2. 456 с.].

Hirth J.P., LotheJ. 1968. Theory of Dislocations. McGraw-Hill, New York, 600 p.

Ionina NA, Kubeev M.N., 2013. The 100 Major Catastrophes. Veche, Moscow, 432 p. (in Russian) [Ионина Н.А., Кубеев М.Н. 100 великих катастроф. М.: Вече, 2013. 432 с.].

Ivanchin A.G., 1982. The Role of Cooperative Effects in the Energy Dissipation During Movements of Dislocations. Candidate of Physics and Mathematics Thesis. Institute of Atmospheric Optics, SB RAS, Tomsk, 131 p. (in Russian) [Иванчин А.Г. Роль кооперативных эффектов при движении дислокаций в диссипации энергии: Дис. ... канд. физ.-мат. наук. Томск: ИОА СО АН СССР, 1982. 131 с.].

Khachai О.А., Khachai O.Yu, 2012. About the non-linear dynamic modelling of responses of the heterogenic medium to strong explosive impacts. In: Tectonophysics and Top Issues of Earth Sciences. Proceedings of the 3rd Tectono-physical Conference, 8-12 October 2012. IPE RAS, Moscow, vol. 2, p. 415-418 (in Russian) [Хачай О.А., Хачай О.Ю. О построении нелинейной динамической модели отклика гетерогенной среды на сильные взрывные воздействия // Тектонофизика и актуальные вопросы наук о Земле: Материалы третьей Тектонофизической конференции. Т. 2. М.: ИФЗ РАН, 2012. С. 415-418].

Khachai ОА, Khachai O.Yu., Klimko V.K. et al., 2013. Kinematic and dynamic characteristics of slow strain waves in the natural massif as a response to blast impacts. In: Geodynamics and the state of stresses of the Earth interior. Proceedings of the 20th All-Russia Scientific Conference. Mining Institute SB RAS, Novosibirsk, p. 38-42 (in Russian) [Хачай ОА., Хачай О.Ю., Климко В.К. и др. Кинематические и динамические характеристики медленных де-

360

Geodynamics & Tectonophysics 2015 Volume 6 Issue 3 Pages 345-364

формационных волн в породном массиве как отклик на взрывные воздействия // Геодинамика и напряженное строение недр Земли: Труды ХХ Всероссийской научной конференции. Новосибирск: ИГД СО РАН, 2013. С. 38-42].

Korobov A.I., Odina N.I., Mekhedov D.M., 2013. Effect of slow dynamics on elastic properties of materials with residual and shear strains. Acoustical Physics 59 (4), 387-392. http://dx.doi.org/10.1134/S106377101304009X.

Koronovsky N., 2012. The Yellowstone supervolcano. Nauka i Zhizn (Science and Life) (5), 110-113 (in Russian) [Коро-новский Н. Йеллоустонский супервулкан // Наука и жизнь. 2012. № 5. С. 110-113].

Krasny L.I. (Ed.), 2004. Earth Planet. Encyclopaedia. Tectonic and Geodynamics. VSEGEI, St. Petersburg, 652 p. (in Russian) [Планета Земля. Энциклопедический справочник. Том «Тектоника и геодинамика» / Ред. Л.И. Красный. СПб.: ВСЕГЕИ, 2004. 652 с.].

Kropotkin P.N., 1948. Main problems of energy of tectonic processes. Izvestiya AN SSSR (5), 89-104 (in Russian) [Кропоткин П.Н. Основные проблемы энергетики тектонических процессов // Известия АН СССР. 1948. № 5. С. 89-104].

Kuzikov S.I., Mukhamediev S.A., 2010. Structure of the present-day velocity field of the crust in the area of the Central-Asian GPS network. Izvestiya, Physics of the Solid Earth 46 (7), 584-601. http://dx.doi.org/10.1134/S106935131 0070037.

Landau L.D., Lifshits E.M., 1973. Mechanics. Nauka, Moscow, 208 p. (in Russian) [Ландау Л.Д., Лифшиц Е.М. Механика. М.: Наука, 1973. 208 с.].

Landau L.D., Lifshits E.M., 2003. Theory of Elasticity. Nauka, Moscow, 264 p. (in Russian) [Ландау Л.Д., Лифшиц Е.М. Теория упругости. М.: Наука, 2003. 264 с.].

Laverov N.P. (Ed.), 2005. Recent Volcanism in the Territory of Russia. Nauka, Moscow, 604 p. (in Russian) [Новейший и современный вулканизм на территории России / Ред. Н.П. Лаверов. М.: Наука, 2005. 604 с.].

Leonov M.G., 2008. Tectonics of Consolidated Crust. Nauka, Moscow, 457 p. (in Russian) [Леонов М.Г. Тектоника консолидированной коры. М.: Наука, 2008. 457 с.].

Leonov M.G., Grib E.N., 2004. Structural Positions and Volcanism of Quaternary Calderas of Kamchatka. Dalnauka, Vladivostok, 189 p. (in Russian) [Леонов В.Л., Гриб Е.Н. Структурные позиции и вулканизм четвертичных кальдер Камчатки. Владивосток: Дальнаука, 2004. 189 с.].

Luchitsky I.V., 1971. Fundamentals of Paleovolcanology. V. 1. Nauka, Moscow, 480 p. (in Russian) [Лучицкий И.В. Основы полеовулканологии. Т. 1. М.: Наука, 1971. 480 с.].

Lunina O.V., GladkovA.S., Nevedrova N.N., 2009. Rift Basins of Pribaikalie: Tectonic Structure and Development History. GEO, Novosibirsk, 316 p. (in Russian) [Лунина О.В., Гладков А.С., Неведрова Н.Н. Рифтовые впадины Прибайкалья: тектоническое строение и история развития. Новосибирск: ГЕО, 2009. 316 с.].

Magnitskii V.A., Mukhamediev S.A., Khasanov R., 1998. The possibility of rock melting in the Earth's crust during intense folding: evidence from the Pamirs. Doklady Earth Sciences 363A (9), 1292-1296.

McDonald G, 1975. Volcanoes. Mir, Moscow, 432 p. (in Russian) [Макдональд Г. Вулканы. М.: Мир, 1975. 432 с.].

Mulev S.N., Bondarev A.V., Panin S.F., 2013. The experience of the implementation of seismic monitoring systems in mines in Russia. In: Geodynamics and the state of stresses of the Earth interior. Proceedings of the 20th All-Russia Conference with participation of researchers from other countries (7-11 October 2013, Novosibirsk). Mining Institute SB RAS, Novosibirsk, p. 416-421 (in Russian) [Мулев С.Н., Бондарев А.В., Панин С.Ф. Опыт внедрения систем сейсмического мониторинга на шахтах и рудниках России // Геодинамика и напряженное состояние недр Земли: Труды ХХ Всероссийской конференции с участием иностранных ученых, Новосибирск, 711 октября 2013. Новосибирск: ИГД СО РАН, 2013. С. 416-421].

Nalivkin D.V, 1969. Hurricanes, Storms and Tornadoes. Nauka, Leningrad, 487 p. (in Russian) [Наливкин Д.В. Ураганы, бури и смерчи. Л.: Наука, 1969. 487 с.].

New Global Tectonics (Plate Tectonics), 1974. Mir, Moscow. 472 p. (in Russian) [Новая глобальная тектоника (тектоника плит). М.: Мир, 1974. 472 с.].

Nikolaev A.V. (Ed.), 1987. Problems of Non-Linear Seismics. Nauka, Moscow, 228 p. (in Russian) [Проблемы нелинейной сейсмики / Ред. А.В. Николаев. М.: Наука, 1987. 228 с.].

Nikolaev A.V. (Ed.), 2003. Problems of Geophysics in the 21st Century. Two volumes. Nauka, Moscow, (in Russian) [Проблемы геофизики XXI века / Ред. А.В. Николаев. В 2-х кн. М.: Наука, 2003].

Nikolaevsky V.N., 1996. Geodynamics and Fluid Dynamics. Nedra, Moscow, 447 p. (in Russian) [Николаевский В.Н. Геодинамика и флюидодинамика. М.: Недра, 1996. 447 с.].

Oparin V.N., Annin B.D., Chugui Yu.V. et al., 2007. Methods and Measurement Equipment for Modelling and Field Studies of Non-Linear Strain and Wave Processes in Rock Blocks. SB RAS Publishing House, Novosibirsk, 320 p. (in Russian) [Опарин В.Н., Аннин Б.Д., Чугуй Ю.В. и др. Методы и измерительные приборы для моделирования и натурных исследований нелинейных деформационно-волновых процессов в блочных массивах горных пород. Новосибирск: Изд-во СО РАН, 2007. 320 с.].

361

A.V. Vikulin: Geodynamics as wave dynamics...

Oparin V.N., Simonov B.F., Yushkin V.F. et al., 2010. Geomechanical and Technical Bases for Enhancement of Oil Recovery from Formations with Application of Vibration Wave Technology. Nauka, Novosibirsk, 404 p. (in Russian) [Опарин В.Н., Симонов Б.Ф., Юшкин В.Ф. и др. Геомеханические и технические основы увеличения нефтеотдачи пластов в виброволновых технологиях. Новосибирск: Наука, 2010. 404 с.].

Oparin V.N., Vostrikov VI, 2010. The energy criterion of the volume destruction of focal zones, and oscillation waves. In: Methods and Systems for Monitoring of Seismic Deformation Caused by Technogenic Earthquakes and Rock Shocks. V. 2: Integration Projects of SB RAS. Issue 25. SB RAS Publishing House, Novosibirsk, p. 93-128 (in Russian) [Опарин

В.Н., Востриков В.И. Энергетический критерий объемного разрушения очаговых зон и волны маятникового типа // Методы и системы сейсмодеформационного мониторинга техногенных землетрясений и горных ударов. Т. 2. Интеграционные проекты СО РАН. Вып. 25. Новосибирск: Изд-во СО РАН, 2010. C. 93-128].

Pamyatnykh E.A., Ursulov A.V, 2012. Nonlinear solitary waves in nonlocal elastic solids. Acoustical Physics 58 (2), 160-165. http://dx.doi.org/10.1134/S1063771012010162.

Pavlenkova N.l, 2013. The origin of the Moho according to geophysical data. In: Proceedings of the XLV Tectonic Meeting. GEOS, Moscow, p. 138-141 (in Russian) [Павленкова Н.И. Природа границы М по геофизическим данным // Материалы XLV Тектонического совещания. М.: ГЕОС, 2013. С. 138-141].

Peive A.V., 1961. Tectonics and magmatism. Izvestia AN SSSR. Geological Series (3), 36-54 (in Russian) [Пейве А.В. Тектоника и магматизм // Известия АН СССР. Серия геологическая. 1961. № 3. С. 36-54].

Ponomarev V.S., 2008. Energy Capacity of the Geological Medium. Nauka, Moscow, 379 p. (in Russian) [Пономарев В.С. Энергонасыщенность геологической среды. М.: Наука, 2008. 379 с.].

Prozorov A.G., 1978. On a lower probability of strong shocks in the space-and-time vicinity of strong global earthquakes. In: Issues of earthquake forecasting and the Earth structure. Computational Seismology. Issue 11. Nauka, Moscow, p. 35-41 (in Russian) [Прозоров А.Г. О пониженной вероятности сильных толчков в некоторой пространственно-временной окрестности сильных землетрясений мира // Вопросы прогноза землетрясений и строение Земли. Вычислительная сейсмология. Вып. 11. М.: Наука, 1978. С. 35-41].

Riznichenko Yu.V., 1985. Problems of Seismology. Selected Works. Nauka, Moscow. 408 p. (in Russian) [Ризниченко Ю.В. Проблемы сейсмологии. Избранные труды. М.: Наука, 1985. 408 с.].

Rudenko O.V., Sobisevich A.L., Sobisevich L.E., Hedberg C.M., Shamaev N.V., 2012. Nonlinear model of a granulated medium containing viscous fluid layers and gas cavities. Acoustical Physics 58 (1), 99-106. http://dx.doi.org/ 10.1134/S1063771012010198.

Rykunov L.N., Khavroshkin O.B., Tsiplakov V.V, 1979. Temporary variations of high-frequency seismic noise. Izvestiya AN SSSR, Fizika Zemli (11), 72-77 (in Russian) [Рыкунов Л.Н., Хаврошкин О.Б., Циплаков В.В. Временные вариации высокочастотных сейсмических шумов // Известия АН СССР. Физика Земли. 1979. № 11. С. 72-77].

Sadovsky M.A., 2004. Selected Works. Geophysics and Blast Physics. Nauka, Moscow, 440 p. (in Russian) [Садовский М.А. Избранные труды: Геофизика и физика взрыва. М.: Наука, 2004. 440 с.].

Sedov L.l., 1973. Solid Medium Mechanics. V. 1. Nauka, Moscow. 536 p. (in Russian) [Седов Л.И. Механика сплошной среды. Т. 1. М.: Наука, 1973. 536 с.].

Sherman S.l., 2009. A tectonophysical model of a seismic zone: experience of development based on the example of the Baikal rift system. Izvestiya, Physics of the Solid Earth 45 (11), 938-941. http://dx.doi.org/10.1134/S1069351309 110020.

Sherman S.l., 2013. Deformation waves as a trigger mechanism of seismic activity in seismic zones of the continental lithosphere. Geodynamics & Tectonophysics 4 (2), 83-117 (in Russian) [Шерман С.И. Деформационные волны как триггерный механизм сейсмической активности в сейсмических зонах континентальной литосферы // Геодинамика и тектонофизика. 2013. Т. 4. № 2. С. 83-117]. http://dx.doi.org/10.5800/GT-2013-4-2-0093.

Sherman S.l., 2014. Seismic Process and the Forecast of Earthquakes: Tectonophysical Conception. Academic Publishing House "Geo", Novosibirsk, 359 p. (in Russian) [Шерман С.И. Сейсмический процесс и прогноз землетрясений: тектонофизическая концепция. Новосибирск: Академическое издательство «Гео», 2014. 359 с.].

Sherman S.l., Cheremnykh A.V., Bornyakov S.A., Shishkina L.P., 2001. Modeling of large faults in zones of lithospheric extension and numerical constraints on deformation. Geologiya i Geofizika (Russian geology and geophysics) 42 (7), 1052-1057.

Sherman S.l., Gorbunova E.A, 2008. Wave origin of fault activation in Central Asia on the basis of seismic monitoring. Fizicheskaya Mezomekhanika 11 (1), 115-122 (in Russian) [Шерман С.И., Горбунова Е.А. Волновая природа активизации разломов Центральной Азии на базе сейсмического мониторинга // Физическая мезомеханика. 2008. Т 11. № 1. С. 115-122].

Sherman S.l., Seminsky K.Zh., Bornyakov S.A., Buddo V.Yu., Lobatskaya R.M., Adamovich A.N., Truskov V.A., BabichevA.A., 1992. Faulting in the Lithosphere. Tensile Stress Zones. Nauka, Siberian Branch, Novosibirsk, 227 p. (in Russian) [Шерман С.И., Семинский К.Ж., Борняков С.А., Буддо В.Ю., Лобацкая Р.М., Адамович А.Н., Трусков В.А., Бабичев А.А. Разломообразование в литосфере. Зоны растяжения. Новосибирск: Наука. Сибирское отделение, 1992. 227 с.].

Sherman S.l., Seminsky K.Zh., Bornyakov S.A., Buddo V.Yu., Lobatskaya R.M., Adamovich A.N., Truskov V.A., BabichevA.A., 1994. Faulting in the Lithosphere. Compression Zones. Nauka, Siberian Branch, Novosibirsk, 262 p. (in Russian)

362

Geodynamics & Tectonophysics 2015 Volume 6 Issue 3 Pages 345-364

[Шерман С.И., Семинский К.Ж., Борняков С.А., Буддо В.Ю., Лобацкая Р.М., Адамович А.Н., Трусков В.А., Бабичев А.А. Разломообразование в литосфере. Зоны сжатия. Новосибирск: Наука. Сибирское отделение, 1994. 262 с.].

Sherman S.I., Sorokin A.P., Sorokina A.T., Gorbunova E.A., Bormotov V.A., 2011. New data on the active faults and zones of modern lithosphere destruction in the Amur region. Doklady Earth Sciences 439 (2), 1146-1151. http://dx.doi.org/10.1134/S1028334X11080186.

Siebert L., Simkin T., Kimberly P, 2010. Volcanoes of the World. Third edition. Smithsonian Institution. Washington DC, University of California Press, 551 p.

Turcotte D., Schubert G., 1985. Geodynamics: Applications of Continuum Physics to Geological Problems. Part 1. Mir, Moscow, 376 p. (in Russian) [Теркот Д., Шуберт Дж. Геодинамика: Геологические приложения физики сплошных сред. Ч. 1. М.: Мир, 1985. 376 с.].

Tveritinova T.Yu., Vikulin A.V., 2014. The geomedium block concept and the problem of block boundaries. In: Faulting in the lithosphere and related processes: tectonophysical analysis. Abstracts of the All-Russia Conference with visiting researchers from other countries (11-16 August 2014, Irkutsk). IEC SB RAS, Irkutsk. P. 35-36 (in Russian) [Тверитинова Т.Ю., Викулин А.В. Концепция блоковой геосреды и проблема границ блоков // Разломобра-зование в литосфере и сопутствующие процессы: тектонофизический анализ: Тезисы докладов Всероссийского совещания с участием приглашенных исследователей из других стран. Иркутск, 11-16 августа, 2014. Иркутск: ИЗК СО РАН, 2014. С. 35-36].

Vikulin A.V., 2008. Vortex Motion World. KSTU, Petropavlovsk-Kamchatsky, 230 p. (in Russian) [Викулин А.В. Мир вихревых движений. Петропавловск-Камчатский: КГТУ, 2008. 230 с.].

Vikulin A.V., 2009. Physics of the Earth and Geodynamics. Textbook. Kamchatka State University, Petropavlovsk-Kamchatsky, 463 p. (in Russian) [Викулин А.В. Физика Земли и геодинамика: Учебное пособие. Петропав-ловск-Камчатский: КамГУ, 2009. 463 с.].

Vikulin A.V., 2010. New type of elastic rotational waves in geo-medium and vortex geodynamics. Geodynamics & Tectonophysics 1 (2), 119-141 (in Russian) [Викулин А.В. Новый тип упругих ротационных волн в геосреде и вихревая геодинамика // Геодинамика и тектонофизика. 2010. Т. 1. № 2. С. 119-141]. http://dx.doi.org/10. 5800/GT-2010-1-2-0010.

Vikulin A.V., 2011. Seismicity. Volcanism. Geodynamics. Collection of Works. Kamchatka State University, Petropavlovsk-Kamchatsky, 407 p. (in Russian) [Викулин А.В. Сейсмичность. Вулканизм. Геодинамика: Сборник трудов. Петропавловск-Камчатский: КамГУ, 2011. 407 с.].

VikulinA.V., 2013. Non-linearity-fractality or rheidity-power capacity: what categories are closer to geology? (Review of the article by N.V. Koronovsky and A.A. Naimark "Methods of dynamical geology at the threshold of applicability"). Vestnik KRAUNC (Kamchatka Regional Association Education-Science Centre). Earth Sciences (1), 152-157 (in Russian) [Викулин А.В. Нелинейность-фрактальность или реидность-энергонасыщенность: какие категории ближе геологии? (Отзыв на статью Н.В. Короновского, А.А. Наймарка «Методы динамической геологии на критическом рубеже применимости») // Вестник КРАУНЦ. Науки о Земле. 2013. № 1 (21).

С. 152-157].

Vikulin A.V., 2014a. The Earth Acoustics. In: The system of Earth Planet. The 20th anniversary of the system of Earth Planet workshop. LENARD, Moscow, p. 202-242 (in Russian) [Викулин А.В. Акустика Земли // Система «Планета Земля»: ХХ лет семинару «Система "Планета Земля"». М.: ЛЕНАРД 2014. С. 202-242].

Vikulin A.V., 2014b. The problem of vortex motions in geology (Review of the article by E.G. Mirlin and Yu.V. Mironov "The role of vortex motions in geodynamics of the Aegean sea (based on comparative analysis with the geodynamics of the Woodlark basin)". Vestnik KRAUNC (Kamchatka Regional Association Education-Science Centre). Earth Sciences (2), 186-190 (in Russian) [Викулин А.В. Проблема вихревых движений в геологии (Отзыв на статью Е.Г. Мирлина, Ю.В. Миронова «Роль вихревого движения в геодинамике Эгейского моря (на основе сравнительного анализа с геодинамикой котловины Вудларк)») // Вестник КРАУНЦ. Науки о Земле. 2014. № 2 (24). С. 186-190].

Vikulin A.V., Akmanova D.R., 2014. A magma focus as a property of the crust. Vestnik KRAUNC (Kamchatka Regional Association Education-Science Centre). Earth Sciences (1), 213-230 (in Russian) [Викулин А.В., Акманова Д.Р. Магматический очаг как свойство земной коры // Вестник КРАУНЦ. Науки о Земле. 2014. № 1 (23). С. 213230].

Vikulin A.V., Akmanova D.R., Osipova N.A., 2010. Volcanism as an indicator of geodynamic processes. Litosfera (3), 5-11 (in Russian) [Викулин А.В., Акманова Д.Р., Осипова Н.А. Вулканизм как индикатор геодинамических процессов // Литосфера. 2010. № 3. С. 5-11].

Vikulin A.V., Akmanova D.P., Vikulina S.A., Dolgaya A.A, 2012a. Migration of seismic and volcanic activity as display of wave geodynamic process. Geodynamics and Tectonophysics 3 (1), 1-18. http://dx.doi.org/10.5800/GT-2012-3-1-0058.

Vikulin A.V., Ivanchin A.G., 2013a. Modern concept of block hierarchy in the structure of geomedium and its implications in geosciences. Journal of Mining Science 49 (3), 395-408. http://dx.doi.org/10.1134/S1062739149030076.

Vikulin A.V., Ivanchin A.G., 2013b. On the Moho origin as a phase boundary between the crust composed of blocks and the non-block upper mantle. In: Proceedings of the XLV Tectonic meeting. GEOS, Moscow, p. 38-42 (in Russian)

363

A.V. Vikulin: Geodynamics as wave dynamics...

[Викулин А.В., Ивнчин А.Г. О природе Мохо как фазовой границе раздела между блоковой земной корой и неблоковой верхней мантией // Материалы XLV Тектонического совещания. М.: ГЕОС, 2013. С. 38-42].

Vikulin A.V., Ivanchin A.G., 2015. A new insight to the origin of a magma pocket. In: The system of Earth Planet. The 200th anniversary of the Holy Alliance (1815-2015). LENARD, Moscow, p. 293-321 (in Russian) [Викулин А.В., Иванчин А.Г. Новый взгляд на природу магматического очага // Система «Планета Земля»: 200 лет Священному союзу (1815-2015). М.: ЛЕНАРД, 2015. С. 293-321].

Vikulin A.V., Melekestsev I.V., Akmanova D.R., Ivanchin A.G., Vodinchar G.M., Dolgaya A.A., Gusyakov V.K., 2012b. The IT system for modelling of seismic and volcanic processes as a basis for studies of wave geodynamic phenomena. Computational Technologies 17 (3), 34-54 (in Russian) [Викулин А.В., Мелекесцев И.В., Акманова Д.Р. и др. Информационно-вычислительная система моделирования сейсмического и вулканического процессов как основа изучения волновых геодинамических явлений // Вычислительные технологии. 2012. Т. 17. № 3. С. 34-54].

Vikulin A.V., Tveritinova T.Yu., Ivanchin A.G., 2013. Wave moment geodynamics. Acta Geophysica 61 (2), 245-263. http://dx.doi.org/10.2478/s11600-012-0079-8.

Vlodavets V.I., 1984. Volcanology Guidebook. Nauka, Moscow, 340 p. (in Russian) [Влодавец В.И. Справочник по вулканологии. М.: Наука, 1984. 340 с.].

Zharkov V.N., 1983. The Interior Structure of the Earth and Other Planets. Nauka, Moscow, 416 p. (in Russian) [Жарков В.Н. Внутреннее строение Земли и планет. М.: Наука, 1983. 416 с.].

Ziman J.M., 1974. Principles of the Theory of Solids. Mir, Moscow, 472 p. (in Russian) [Займан Дж. Принципы твердого тела. М.: Мир, 1974. 472 с.].

Vikulin, Alexander V., Doctor of Physics and Mathematics, Lead Researcher

Institute of Volcanology and Seismology, Far East Branch of the Russian Academy of Sciences

9 Boulevard Piip, Petropavlovsk-Kamchatsky 683006, Russia

И e-mail: vik@kscnet.ru

Викулин Александр Васильевич, д.ф.-м.н, в.н.с.

Институт вулканологии и сейсмологии ДВО РАН 683006, Петропавловск-Камчатский, б-р Пийпа, 9, Россия И e-mail: vik@kscnet.ru