Hybrid seismic modeling based on discrete-wave number and finite-difference methods

Any calculation of seismic wave propagation comprising the seismic source, the travel path, and the receiver site in a single finite-difference (FD) model requires a considerable amount of computer time and memory. Moreover, the methods currently available for including point sources in the 2D FD calculations are far-field approximations only. Therefore we have developed a new hybrid method for treating the seismic wave fields at localized 2D near-surface structures embedded in a 1D background medium, and excited by a point source. The source radiation and propagation in the background model is solved by the discrete-wave number (DW) method, while the propagation in the local 2D structure is calculated by the FD method. The coupling between the two sets of calculations is performed on a rectangular excitation box surrounding the local structure. We show the usefulness of the method in ground-motion studies where both near-field source effects and local site effects are important. Technical problems connected with the inconsistency between the 3D source radiation and the 2D FD calculation are minor for the relatively distant in-plane point explosive sources, but are more serious for the in-plane dislocation sources.     

研究“源兆”的若干问题及途径分析

从位场理论的观点讨论了场和源的概念,建议将“场兆”改称为构造前兆,将“源兆”称为震源前兆。文中还分析了震源体与各种前兆观测方法探测的深度,并根据地磁与重力观测方法不同的物理基础和地质基础,指出震磁效应对于研究震源前兆的特殊意义,文中还介绍了１９７６唐山地震观测到震源前兆的部分实际地磁资料。     

The Epi－continental arc of Southeast China and relevant earthquakes

ＴｈｅＥｐｉ┐ｃｏｎｔｉｎｅｎｔａｌａｒｃｏｆＳｏｕｔｈｅａｓｔＣｈｉｎａａｎｄｒｅｌｅｖａｎｔｅａｒｔｈｑｕａｋｅｓＪＩＡ－ＷＥＩＸＵ（徐嘉炜）ＤｅｐａｒｔｍｅｎｔｏｆＲｅｓｏｕｒｃｅａｎｄＥｎｖｉｒｏｎｍｅｎｔａｌＳｃｉｅｎｃｅｓ,ＨｅｆｅｉＵ...     

Blind prediction of the site effects at Ashigara Valley,Japan, and its comparison with reality

Weak and strong ground motions were numerically predicted for three stations of the Ashigara Valley test site. The prediction was based on the records from a rock-outcrop station, one weak-motion record from a surface-sediments station, and the standard geotechnical model. The data were provided by the Japanese Working Group on the Effects of Surface Geology as a part of an international experiment. The finite-difference method for SH waves in a 2-D linear viscoelastic medium (a causalQ model) was employed.Comparison with the real records shows that at two stations the predictions fit better than at the third one. Strangely, the two better predictions were for stations situated at larger distances from the reference rock station (one station was on the surface, the other in a borehole). The strong ground motion (the peak acceleration of about 200 cm s^–2) was not predicted qualitatively worse than the weak motion (8 cm s^–2). A less sophisticated second prediction (not submitted during the experiment), in which we did not attempt to fit the available weak-motion record at the sedimentary station, agrees with the reality significantly better.     

The emplacement of intermediate volume ignimbrites: A case study from Roccamonfina Volcano,Southern Italy

A model is presented for the emplacement of intermediate volume ignimbrites based on a study of two 6 km³ volume ignimbrites on Roccamonfina Volcano, Italy. The model considers that the flows were slow moving, and quickly deflated from turbulent to non-turbulent conditions. Yield strength and density increased whereas fluidisation decreased with time and runout of the pyroclastic flows. In proximal locations, on the caldera rim, heterogeneous exposures including discontinuous lithic breccias, stratified and cross-stratified units interbedded with massive ignimbrite suggest deposition from turbulent flows. In medial locations thick, massive ignimbrite occurs associated with three types of co-ignimbrite lithic breccia which we interpret as being emplaced by non-turbulent flows. Multiple grading of different breccia/lithic concentration types within single flow units indicates that internal shear occurred producing overriding or overlapping of the rear of the flow onto the slower-moving front part. This overriding of different parts of non-turbulent pyroclastic flows could be caused by at least two different mechanisms: (1) changes in flow regime, such as hydraulic jumps that may occur at breaks in slope; and (2) periods of increased discharge rate, possibly associated with caldera collapse, producing fresh pulses of lithic-rich material that sheared onto the slower-moving part of the flow in front.We propose that ground surge deposits enriched in pumice compared with their associated ignimbrite probably formed by a flow separation mechanism from the top and front of the pyroclastic flow. These turbulent clouds moved ahead of the non-turbulent lower part of the flow to form stratified pumice-rich deposits. In distal regions well-developed coarse, often clast-supported, pumice concentrations zones and coarse intra-flow-unit lithic concentrations occur within the massive ignimbrite. We suggest that the flows were non-turbulent, possessed a relatively high yield strength and may have moved by plug flow prior to emplacement.     

Automatic equipment to monitor groundwater radon

Prototype instrumentation, able to automatically measure groundwater radon content variations, is presented. The equipment is made of stainless steel and has spherical valves with automatic and pneumatic control. The deemanation of the gases from the water is obtained by evacuating a suitable expansion chamber. The instrumentation can make discrete sampling ranging from 1 per hour to 1 per 99 hours. The equipment was tested in the laboratory: the efficiency was measured by means of a²⁶⁶Ra solution. A mean value of (0.65±0.07) count/s/Bq was obtained. A calibration test was carried out by comparing countings from the automatic equipment with those obtained by the standard laboratory cell. Results of an operational check over a period of approximately one year indicate that variations in radon at the calibration site are attributable more to meteorological than to tectonic causes.     

1983年菏泽5.9级地震前后沂水泉氡的变化

沂水泉位于山东省沂水一汤头断裂带上。多年来的观测分析表明,影响沂水泉氡值变化的最大干扰是降雨。由于这种干扰存在着“记忆”滞后影响,所以用一般的相关分析难以将这种影响排除干净。本文对最小二乘法为建立模型的判据,用动态灰箱分析法拟合并推估观测系统的变化。在分析中用了四个状态变量：长趋势成分Ｌ（ｔ）、年周期变化成分Ｐ（ｔ）、降雨影响部分Ｒ（ｔ）和随机变化部分Ｓ（ｔ）。经过分析表明,沂水泉氡的长趋势部分１     

The arrangement in mean elements space of the periodic orbits close to that of the Moon

In a simplified model of the Earth-Moon-Sun system based on the restricted circular 3-dimensional 3-body problem, it is possible to find numerically a set of 8 periodic orbits whose time evolutions closely resemble that of the Moon's orbit. These orbits have a period of 223 synodic months (i.e. the period of the Saros cycle known for more than two millennia as a means of predicting eclipses), and are characterized by a secular rotation of the argument of perigee . Periodic orbits of longer durations exhibiting this last feature are very abundant in Earth-Moon-Sun dynamical models. Their arrangement in the space of the mean orbital elements- for various values of the lunar mean motion is presented.     

Generalized Lyapunov exponents indicators in Hamiltonian dynamics: An application to a double star system

The Lyapunov characteristic numbers (LCNs) which are defined as the mean value of the distribution of the local variations of the tangent vectors to the flow (=ln _k ⁱ ) (see Froeschlé, 1984) have been found to be sensitive indicators of stochasticity. So we computed the distribution of these local variations and determined the moments of higher order for the integrable and stochastic regions in a binary star system with =0.5.     

Asteroid motion near the 3 : 1 resonance

A mapping model is constructed to describe asteroid motion near the 3 : 1 mean motion resonance with Jupiter, in the plane. The topology of the phase space of this mapping coincides with that of the real system, which is considered to be the elliptic restricted three body problem with the Sun and Jupiter as primaries. This model is valid for all values of the eccentricity. This is achieved by the introduction of a correcting term to the averaged Hamiltonian which is valid for small values of the ecentricity.We start with a two dimensional mapping which represents the circular restricted three body problem. This provides the basic framework for the complete model, but cannot explain the generation of a gap in the distribution of the asteroids at this resonance. The next approximation is a four dimensional mapping, corresponding to the elliptic restricted problem. It is found that chaotic regions exist near the 3 : 1 resonance, due to the interaction between the two degrees of freedom, for initial conditions close to a critical curve of the circular model. As a consequence of the chaotic motion, the eccentricity of the asteroid jumps to high values and close encounters with Mars and even Earth may occur, thus generating a gap. It is found that the generation of chaos depends also on the phase (i.e. the angles andv) and as a consequence, there exist islands of ordered motion inside the sea of chaotic motion near the 3 : 1 resonance. Thus, the model of the elliptic restricted three body problem cannot explain completely the generation of a gap, although the density in the distribution of the asteroids will be much less than far from the resonance. Finally, we take into account the effect of the gravitational attraction of Saturn on Jupiter's orbit, and in particular the variation of the eccentricity and the argument of perihelion. This generates a mixing of the phases and as a consequence the whole phase space near the 3 : 1 resonance becomes chaotic. This chaotic zone is in good agreement with the observations.