2-D crustal Poisson's ratio from seismic travel time inversion in Changbaishan Tianchi volcanic region

Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L1 are determined simultaneously with deep seismic sounding data in Changbaishan Tianchi volcanic region, and then its Poisson's ratio is obtained. Calculated results show that this technique overcomes some defects of traditional forward calculation method, and it is also very effective to determine Poisson's ratio distribution of deep seismic sounding profile, especially useful for study on volcanic magma and crustal fault zone. Study result indicates that there is an abnormally high Poisson's ratio body that is about 30 km wide and 12 km high in the low velocity region under Tianchi crater. Its value of Poisson's ratio is 8% higher than that of surrounding medium and it should be the magma chamber formed from melted rock with high temperature. There is a high Poisson's ratio zone ranging from magma chamber to the top of crust, which may be the uprise passage of hot substance. The lower part with high Poisson's ratio, which stretches downward to Moho, is possibly the extrusion way of hot substance from the uppermost mantle. The conclusions above are consistent with the study results of both tomographic determination of 3D crustal structure and magnetotelluric survey in this region.     

再论各向异性介质转换波地震学,第一部分:理论

我们业已研发了计算各向异性、非均质介质中P- SV转换波(C-波)的转换点和旅行时的新理论。据此 可以利用诸如相似性分析、迪克斯模型建模、克契 霍夫求和等常规方法来完成各向异性的处理和各向 异性处理,并使各向异性的处理成为可能。这里将 我们的新发展分作两部分来介绍。第一部分为理 论,第二部分为对速度分析和参数计算的应用。第 一部分理论包括转换点的计算和动校正的分析。     

Converted-wave seismology in anisotropic media revisited,Part I: Basic theory

We have developed new basic theories for calculating the conversion point and the travel time of the P-SV converted wave (C-wave) in anisotropic, inhomogeneous media. This enables the use of conventional procedures such as semblance analysis, Dix-type model building and Kirchhoff summation, to implement anisotropic processing, and makes anisotropic processing affordable. Here we present these new developments in two parts: basic theory and application to velocity analysis and parameter estimation. This part deals with the basic theory, including both conversion-point calculation and moveout analysis. Existing equations for calculating the PS-wave (C-wave) conversion point in layered media with vertical transverse isotropy (VTI) are strictly limited to offsets about half the reflector depth (an offset-depth ratio, xlz, of 0.5), and those for calculating the C-wave traveltimes are limited to offsets equal to the reflector depth (x/z=l.0). In contrast, the new equations for calculating the conversion-point extend into offsets about three-times the reflector depth (x/z=3.0), those for calculating the C-wave traveltimes extend into offsets twice the reflector depth (x/z=2.0). With the improved accuracy, the equations can help in C-wave data processing and parameter estimation in anisotropic, inhomogeneous media. This work is funded by the Edinburgh Anisotropy Project (EAP) of the British Geological Survey. First author: Xiangyang Li, Mr. Li is currently a professorial research seismologist (Grade 6) and technical director of the Edinburgh Anisotropy Project in the British Geological Survey. He also holds a honorary professorship in multicomponent seismology at the School of Geosciences, University of Edinburgh. He received his BSc(1982) in Geophysics from Changchun Geological Institute, China, an MSc (1984) in applied geophysics from East China Petroleum Institute (now known as the China University of Petroleum), and a PhD (1992) in seismology from the University of Edinburgh. During 1984–1987, he worked as a lecturer with the East China Petroleum Institute. Since 1991, he has been employed by the British Geological Survey. His research interests include seismic anisotropy and multicomponent seismology.     

Ground deformation of Nisyros Volcano (Greece) for the period 1995–2002: Results from DInSAR and DGPS observations

Differential GPS (DGPS) and Differential Interferometric Synthetic Aperture Radar (DInSAR) analyses were applied to the Kos-Yali-Nisyros Volcanic Field (SE Hellenic Volcanic Arc) to quantify the ground deformation of Nisyros Volcano. After intense seismic activity in 1996, a GPS network was installed in June 1997 and re-occupied annually up to 2002. A general uplift ranging from 14 to 140 mm was determined at all stations of the network. The corresponding horizontal displacements ranged from 13 to 53 mm. The displacement vectors indicate that the island is undergoing extension towards the East, West and South. A two-source “Mogi” model combined with assumed motion along the Mandraki Fault was constructed to fit the observed deformation. The best-fit model assumes sources at a depth of 5500 m NW of the centre of the island and at 6500 m offshore ESE of Yali Island. DInSAR analysis using four pairs of images taken between May 1995 and September 2000 suggests that deformation was occurring during 1995 before the start of the seismic crisis. An amplitude of at least 56 mm along the slant range appeared for the period 1996 through 1999. This deformation is consistent with the two-source model invoked in DGPS modelling. Surface evidence of ground deformation is expressed in the contemporaneous reactivation of the Mandraki Fault. In addition, a 600 m long N-S trending irregular rupture in the caldera floor was formed between 2001 and 2002. This rupture is interpreted as the release of surface stress in the consolidated epiclastic and hydrothermal sediments of the caldera floor.     

2-D crustal Poisson’s ratio from seismic travel time inversion in Changbaishan Tianchi volcanic region

Based on the inversion method of 2D velocity structure and interface, the crustal velocity structures of P-wave and S-wave along the profile L ₁ are determined simultaneously with deep seismic sounding data in Changbaishan Tianchi volcanic region, and then its Poisson’s ratio is obtained. Calculated results show that this technique overcomes some defects of traditional forward calculation method, and it is also very effective to determine Poisson’s ratio distribution of deep seismic sounding profile, especially useful for study on volcanic magma and crustal fault zone. Study result indicates that there is an abnormally high Poisson’s ratio body that is about 30 km wide and 12 km high in the low velocity region under Tianchi crater. Its value of Poisson’s ratio is 8% higher than that of surrounding medium and it should be the magma chamber formed from melted rock with high temperature. There is a high Poisson’s ratio zone ranging from magma chamber to the top of crust, which may be the uprise passage of hot substance. The lower part with high Poisson’s ratio, which stretches downward to Moho, is possibly the extrusion way of hot substance from the uppermost mantle. The conclusions above are consistent with the study results of both tomographic determination of 3D crustal structure and magnetotelluric survey in this region. Foundation item: Key Project from China Earthquake Administration and the Project (95-11-02-01) from Ministry of Science and Technology (2001DIA10003). Contribution No. RCEG200401, Geophysical Exploration Center, China Earthquake Administration.     

Seismic signature of a phreatic explosion: hydrofracturing damage at Karthala volcano, Grande Comore Island, Indian Ocean

Karthala volcano is a basaltic shield volcano with an active hydrothermal system that forms the southern two-thirds of the Grande Comore Island, off the east coat of Africa, northwest of Madagascar. Since the start of volcano monitoring by the local volcano observatory in 1988, the July 11th, 1991 phreatic eruption was the first volcanic event seismically recorded on this volcano, and a rare example of a monitored basaltic shield. From 1991 to 1995 the VT locations, 0.5<M_l<4.3, show a crack shaped pattern (3 km long, 1 km wide) within the summit caldera extending at depth from –2 km to +2 km relative to sea level. This N-S elongated pattern coincides with the direction of the regional maximum horizontal stress as deduced from regional focal mechanism solutions. This brittle signature of the damage associated with the 1991 phreatic eruption is a typical pattern of the seismicity induced by controlled fluid injections such as those applied at geothermal fields, in oil and gas recovery, or for stress measurements. It suggests the 1991 phreatic eruption was driven by hydraulic fracturing induced by forced fluid flow. We propose that the extremely high LP and VT seismicity rates, relative to other effusive volcanoes, during the climax of the 1991 phreatic explosion, are due to the activation of the whole hydrothermal system, as roughly sized by the distribution of VT hypocenters. The seismicity rate in 1995 was still higher than the pre-eruption seismicity rate, and disagrees with the time pattern of thermo-elastic stress readjustment induced by single magma intrusions at basaltic volcanoes. We propose that it corresponds to the still ongoing relaxation of pressure heterogeneity within the hydrothermal system as suggested by the few LP events that still occurred in 1995.Editorial responsibility: H Shinohara     

Caldera formation by magma withdrawal from a reservoir beneath a volcanic edifice

Caldera formation has been explained by magma withdrawal from a crustal reservoir, but little is known about the conditions that lead to the critical reservoir pressure for collapse. During an eruption, the reservoir pressure is constrained to lie within a finite range: it cannot exceed the threshold value for eruption, and cannot decrease below another threshold value such that feeder dykes get shut by the confining pressure, which stops the eruption. For caldera collapse to occur, the critical reservoir pressure for roof failure must therefore be within this operating range. We use an analytical elastic model to evaluate the changes of reservoir pressure that are required for failure of roof rocks above the reservoir with and without a volcanic edifice at Earth's surface. With no edifice at Earth's surface, faulting in the roof region can only occur in the initial phase of reservoir inflation and affects a very small part of the focal area. Such conditions do not allow caldera collapse. With a volcanic edifice, large tensile stresses develop in the roof region, whose magnitude increase as the reservoir deflates during an eruption. The edifice size must exceed a threshold value for failure of the roof region before the end of eruption. The largest tensile stresses are reached at Earth's surface, indicating that faulting starts there. Failure affects an area whose horizontal dimensions depend on edifice and chamber dimensions. For small and deep reservoirs, failure conditions cannot be achieved even if the edifice is very large. Quantitative predictions are consistent with observations on a number of volcanoes.     

Interevent Time Distribution in Seismicity: A Theoretical Approach

This paper presents an analysis of the distribution of the time τ between two consecutive events in a stationary point process. The study is motivated by the discovery of unified scaling laws for τ for the case of seismic events. We demonstrate that these laws cannot exist simultaneously in a seismogenic area. Under very natural assumptions we show that if, after rescaling to ensure Eτ =1, the interevent time has a universal distribution F, then F must be exponential. In other words, Corral’s unified scaling law cannot exist in the whole range of time. In the framework of a general cluster model we discuss the parameterization of an empirical unified law and the physical meaning of the parameters involved. An erratum to this article is available at .