Multifractal and Chaotic Analysis of Vrancea (Romania) Intermediate-depth Earthquakes: Investigation of the Temporal Distribution of Events

The Vrancea seismic region contains an isolated cluster of events beneath the Carpathian Arc Bend in Romania, dipping to about 200 km depth. Seismic activity mainly occurs at intermediate depths (h > 60 km). The main goal of the paper is to perform an in-depth, complex analysis of the occurrence times of these intermediate-depth events. We also try to show the versatility of the methods used to characterize different aspects of the seismicity evolution and to offer a user-friendly software toolbox to do most of the related computations. The earthquake catalog used in this study spans from 1974 to 2002 and includes only the intermediate-depth events. In the first part of the paper, we analyze the multifractal characteristics of the temporal distribution of earthquakes. The study reveals two distinct scaling regimes. At small scales we found a clear nonhomogeneous, multifractal pattern, while at large scales the temporal distribution of events shows a monofractal, and close to Poissonian (random), behavior. The multifractal behavior at small scales (minutes-hours) is shown to be clearly an effect of the short aftershock sequences that occurred after some major Vrancea earthquakes. In the second part of the paper we analyze whether our temporal series shows a persistent (or anti-persistent) long-term behavior, by using the Detrended Fluctuation Analysis (DFA) method. The results suggest that the analyzed temporal series of Vrancea earthquakes is a non-correlated process. In part three of the paper we seek to determine whether the dynamics of our earthquake system (described by the occurrence time of Vrancea earthquakes) is deterministically chaotic, deriving from a rather simple evolution law, or whether it is stochastic and is generated by a system that possesses many degrees of freedom. The results suggest that our signal is stochastic (probably does not possess an attractor). The limited time-span of the catalog and the analysis performed in this paper cannot rule out the emergence of an interesting, quasi-deterministic and low-dimensional structure in the case of major Vrancea earthquakes.Acknowledgement One of the authors (BE) is grateful to the Japanese Ministry of Education for providing him a Monbusho scholarship for studying in DPRI, Kyoto University. We thank Z.R. Struzik, M Holschneider, J. Mori and D. Kaplan for their useful comments, and acknowledge the support of the staff of DPRI, Kyoto University and the National Inst. for Earth Physics, Bucharest. We thank the two reviewers, M.B. Geilikman and M. Anghel, for their useful suggestions which improved the quality of this work.     

地球磁极倒转的分形混沌研究

地球磁场多次发生南北(正负)磁极位置的变换,这已为大家所公认。但造成这种异常现象的原因,则是迄今未能很好解答的一个难题。应用分形混沌理论对地球磁极倒转进行了分析,认为发生地球磁极倒转时间段在时间轴上具有分形分布的性质;地球磁场系统通过不断与外界交换物质和能量,维持一种空间或时间的有序结构。地球磁场极性的随机倒转具有混沌运动的自逆转特性,混沌理论给地磁极性倒转提出了一个简明的动力机制解释。     

Finite element modeling of heating phenomena of cracks excited by high-intensity ultrasonic pulses

A three-dimensional thermo-mechanical coupled finite element model is built up to simulate the phenomena of dynamical contact and frictional heating of crack faces when the plate containing the crack is excited by high-intensity ultrasonic pulses. In the finite element model, the high-power ultrasonic transducer is modeled by using a piezoelectric thermal-analogy method, and the dynamical interaction between both crack faces is modeled using a contact-impact theory. In the simulations, the frictional heating taking place at the crack faces is quantitatively calculated by using finite element thermal-structural coupling analysis, especially, the influences of acoustic chaos to plate vibration and crack heating are calculated and analysed in detail. Meanwhile, the related ultrasonic infrared images are also obtained experimentally, and the theoretical simulation results are in agreement with that of the experiments. The results show that, by using the theoretical method, a good simulation of dynamic interaction and friction heating process of the crack faces under non-chaotic or chaotic sound excitation can be obtained.     

双混沌映射的土家织锦图案数字加密算法

结合混沌序列易生成性、对初始条件敏感性、可完全重现性等特点,给出了一种基于混沌映射的土家织锦图案数字加密算法。首先对原始图案进行离散小波变换(DWT),变换的矩阵奇行与Logistic映射产生的二值序列异或,偶行与Hybrid映射产生的二值序列异或,然后对产生的新矩阵进行DWT逆变换,即得到加密的织锦图案。实验结果表明,本文算法具有较高的安全性和较好的加密效果。     

Dynamics of Two Planets in the 2/1 Mean-Motion Resonance

The dynamics of two planets near a first-order mean-motion resonance is modeled in the domain of the general three-body planar problem. The system studied is the pair Uranus-Neptune (2/1 resonance). The phase space of the resonance and near-resonance regions is studied by means of surfaces of section and spectral analysis techniques. After a thorough investigation of the topology of the phase space, we find that several regimes of motion are possible for the Uranus-Neptune system, and the regions of transition between the regimes of motion are the seats of chaotic motion. This revised version was published online in July 2006 with corrections to the Cover Date.     

The general solution of the HENON–HEILES problem

The general solution of the Henon–Heiles system is approximated inside a domain of the (x, C) of initial conditions (C is the energy constant). The method applied is that described by Poincaré as ‘the only “crack” permitting penetration into the non-integrable problems’ and involves calculation of a dense set of families of periodic solutions that covers the solution space of the problem. In the case of the Henon–Heiles potential we calculated the families of periodic solutions that re-enter after 1–108 oscillations. The density of the set of such families is defined by a pre-assigned parameter ε (Poincaré parameter), which ascertains that at least one periodic solution is computed and available within a distance ε from any point of the domain (x, C) for which the approximate general solution computed. The approximate general solution presented here corresponds to ε = 0.07. The same solution is further improved by “zooming” into four square sub-domain of (x, C), i.e. by computing sufficient number of families that reduce the density parameter to ε = 0.003. Further zooming to reduce the density parameter, say to ε = 10⁻⁶, or even smaller, although easily performable in both areas occupied by stable as well as unstable solutions, was found unnecessary. The stability of all members of each and all families computed was calculated and presented in this paper for both the large solution domain and for the sub-domains. The correspondence between areas of the approximate general solution occupied by stable periodic solutions and Poincaré sections with well-aligned section points and also correspondence between areas occupied by unstable solutions and Poincaré sections with randomly scattered section points is shown by calculating such sections. All calculations were performed using the Runge-Kutta (R-K) 8th order direct integration method and the large output received, consisting of many thousands of families is saved as “Atlas of the General Solution of the Henon–Heiles Problem,” including their stability and is available at request. It is concluded that approximation of the general solution of this system is straightforward and that the chaotic character of its Poincaré sections imposes no limitations or difficulties.     

Collective resonant phenomena on small bodies in the solar system

For both asteroids and meteor streams, and also for comets, resonances play a major role for their orbital evolutions but on different time scales. For asteroids both mean motion resonances and secular resonances not only structure the phase space of regular orbits but are mainly at the origin for the inherent chaos of planet crosser objects.For comets and their chaotic routes temporary trapping into orbital resonances is a well known phenomenon. In addition for slow diffusion through the Kuiper belt resonances are the only candidates for originating a slow chaos.Like for asteroids, resonances with Jupiter play a major role for the orbital evolution of meteor streams. Crossing of separatrix like zones appears to be crucial for the formation of arcs and for the dissolution of streams. In particular the orbital inclination of a meteor stream appears to be a critical parameter for arc formation. Numerical results obtained in an other context show that the competition between the Poynting-Robertson drag and the gravitational interaction of grains near the 2/1 resonance might be very important in the long run for the structure of meteor streams.     

地质系统的复杂性--地质科学的基本问题(Ⅱ)

笔者根据其所提出的“地质作用与时空结构是一切地质现象的本质与核心”的自然哲学理念将地质科学的学科体系划分为地球物质的成分与结构、地质作用、地质学场与地质系统的演化等四大基本领域,它与9个重大基础地质问题以及与其相关的基础理论问题相对应（表2）,笔者将其已发表的地质科学的复杂性理论“地质作用的临界过程动力学－地质系统在混沌边缘分形生长”应用于研究9个重大基础地质问题,着重对于其中所包含的主要基础理论问题的实质及其理论与研究方法作较深入而简要的剖析。目的是通过将地质系统的复杂性作为地质科学的基本问题的新视角对古老而常新的地质科学进行再认识,将重大基础地质问题的研究提高到非线性科学和复杂性理论的层次,并实现地质科学向精确科学的跨越,取得突破性进展。     

东亚地区天气吸引子特征的研究

利用东亚地区１１个有代表性测站（１９７９－１９８８年）连续１０年的５００ｈＰａ逐日位势高度时间序列，全面计算了它们的浑沌特征量。结果表明，各测站相应的关联维数都是分数；完整的Ｌｙａｐｕｎｏｖ指数谱中，不仅都有正值、零和负值，而且Ｌｙａｐｕｎｏｖ指数之和都小于零。因此，从整体来看，相应于东亚地区天气尺度的吸引子是由一个有限窨支撑起来的。东亚地区的天气吸引子所显示的系统平均Ｋｏｌｍｏｇｏｌｏｖ熵约为０     

储集层测井信号的非线性混沌特性

对测井信号的特征提取采用非线性分析方法。利用非平稳信号分析方法——游程检验法,分析松辽盆地某井储集层测井数据的非平稳特性;在此基础上,采用非线性混沌方法对测井数据进行最大Lyapunov指数提取。结果表明:储集层测井序列具有非平稳特性;最大的Lyapunov指数大于零,表明了储集层测井序列具有混沌特征。为此,储集层测井信息的提取可采用非线性混沌理论方法,以便更有效地描述测井信息特征。