Analytical models for bubble growth during decompression of high viscosity magmas

Analytical models for decompressional bubble growth in a viscous magma are developed to establish the influence of high magma viscosity on vesiculation and to assess the time-scales on which bubbles respond to decompression. Instantaneous decompression of individual bubbles, analogous to a sudden release of pressure (e.g. sector collapse), is considered for two end-member cases. The infinite melt model considers the growth of an isolated bubble before significant bubble interaction occurs. The shell model considers the growth of a bubble surrounded by a thin shell and is analogous to bubble growth in a highly vesicular magmatic foam. Results from the shell model show that magmas less viscous than 10⁹ Pa s can freely expand without developing strong overpressures. The timescales for pressure re-equilibration are shortened by increased ratios of bubble radius to shell thickness and by larger decompression. Time-scales for isolated bubbles in rhyolitic melts (infinite melt model) are significantly longer, implying that such bubbles could experience internal pressures greater than the ambient pressure for at least a few hours following a sudden release of pressure. The shell model is developed to assess bubble growth during the linear decompression of a magma body of constant viscosity. For the range of decompression rates and viscosities associated with actual volcanic eruptions, bubble growth continues at approximately the equilibrium rate, with no attendant excess of internal pressure. The results imply that viscosity does not have any significant role in preventing the explosive expansion of high viscosity foams. However, for viscosities of >10⁹ Pa s there is the potential for a viscosity quench under the extreme decompression rates of an explosive eruption. It is proposed that the typical vesicularities of pumice of 0.7–0.8 are a consequence of the viscosity of the degassing magmas becoming sufficiently high to inhibit bubble expansion over the characteristic time-scale of eruption. For fully degassed silicic lavas with viscosities in the range 10¹⁰ to 10¹² Pa s time-scales for decompression of isolated bubbles can be hours to many months.     

邯郸地形变与地下水位动态相互关系研究

依据邯郸地区地质构造及20多年的观测资料,系统研究了地下水位、地壳形变的正常与异常动态变化及其相互关系.本区地形变的主控因素是西部太行山隆起和东部华北平原大幅度沉降的继承性运动,造成地形变向东南方向倾斜;在大区域的周期性强降水后由于地质构造因素地形变向西北向倾斜,形成水位和形变的正常动态变化.采用有效的定量方法提取异常,显示一部分反向西北倾斜的异常与华北的中强震有较好的对应关系.故强降水引起的地倾斜变化不应看作干扰,而是一种引起地壳变形地震前兆的外动力因素,应予高度关注.     

Progress in studies on the structure and physical properties of the mantle in China

The studies on the structure and physical properties of the mantle by Chinese geophysicists from 2003 to 2007 are reviewed in this report. It mainly contains studies on the seismic velocity structure of the mantle,anisotropy of the mantle,mantle discontinuities,mantle convection and the physical properties of mantle. The review concerns mainly the contents,the methods used and the results of the studies. It can be found that new progress in the study on the structure and physical properties of mantle has been made in the last four years in China. In some preexis-tent areas much progress has been made,advanced methods have been adopted,extensive international co-operation has been conducted in many ways,and the scope of the co-operation has gradually expanded. More-over,some new fields appear as well.     

利用地球简正模耦合研究上地幔过渡区方位各向异性

地震方位各向异性广泛存在于地球上地幔中,目前利用地震体波或面波分析研究上地幔各向异性的地球物理方法有很多种,但是由于各自的局限性均难以分析上地幔过渡区中的各向异性特征.方位各向异性可导致球形简正模和环形简正模之间发生耦合.地球长周期自由振荡的简正模可深入到上地幔过渡区.本文利用各向异性地球模型计算各向异性简正模耦合深度敏感核,表明长周期(250~400 s)简正模各向异性耦合(如₀S₂₀-₀T₂₁ 和₀S₂₅-₀T₂₅)的敏感度峰值在400~600 km之间.在不受地球自转影响的台站,如位于南极极点的QSPA站,仍然可以观测到强烈的简正模耦合现象.本文的研究表明:只有在地震观测台站靠近长周期球形振荡的节点时,才能在其观测数据中观测到各向异性耦合现象,许多各向异性耦合在震后18~24 h期间最强,并可导致垂直方向的环形简正模的振幅大于球形耦合简正模的振幅.这些特征是在地震观测数据中寻找各向异性耦合的重要线索.长周期简正模的方位各向异性耦合为我们提供了一个新的认识上地幔过渡区各向异性的窗口.     

六安地震台TJ-Ⅱ井下体应变观测与资料初析

介绍了六安地震台TJ-Ⅱ体应变仪安装井的地质构造特点与井孔岩芯物理特性,以及TJ-Ⅱ体应变仪在六安地震台的运行情况。通过对六安台TJ-Ⅱ一年多来观测资料映震效应与区域前兆异常分析,初步分析了区域地下形变对大别山区未来中强地震预测研究的实际意义。     

中国地幔结构及物性研究的进展

对2004——2007年中国地球物理学家在地幔内部结构和物质性质方面的研究工作进行了总结. 对地幔结构的地震波速度成像,地幔介质的各向异性,地幔间断面,地幔对流,以及地幔介质物性进行了综述,指出了各方面的主要内容,使用的主要方法和主要结果. 从这4年的研究可以发现,一些原有的研究领域工作更加深入,方法更加先进,而且进行了广泛的国际合作,合作的范围也逐渐扩大,方式多样,并出现了一些新的研究方向.      

震后GPS观测数据揭示的日本M_W9.0地震周边地区地幔黏滞性结构垂向变化

本文利用大范围的震后GPS数据和黏弹性球形地球位错理论,定量研究了日本M_W9.0地震周边地区地幔黏滞性结构的垂向变化.首先结合陆地和海底的GPS观测数据,以及基于球形地球位错理论格林函数和贝叶斯反演方法,反演了该地震的同震滑动分布,发现其最大错动量高达59m.然后在均一地幔黏滞性结构的假设前提下,确定了震源周边地区地幔黏滞因子的最优解,发现依据该地幔黏滞因子获得的理论远场震后位移和GPS观测结果之间的均方根误差高达0.81cm,不能解释远场观测结果.为解决上述问题,本文对震中周边地区地幔黏滞性结构沿垂向方向进行分层,建立了一个随深度变化的地幔黏滞性构造模型,然后综合利用远近场的GPS数据对该地区地幔黏滞因子进行反演研究,结果表明,震源周边地区岩石圈弹性层厚度最优解为40km,40～220km深度的地幔黏滞因子最优解为6×10~(18)Pa·s,220～670km深度之间的地幔黏滞因子最优解为1.5×10~(19)Pa·s.上述地幔黏滞性构造使远场的均方根误差降为0.12cm,仅为利用均一地幔黏滞性构造所得均方根误差值的15%,大大提高了远场模拟结果的准确性.最后,观测值和模拟值之间的均方根误差分析表明,近场震后形变数据主要约束浅层的地幔黏滞性结构,而远场震后形变数据主要约束深部的地幔黏滞性结构.     

Influence of the Load Wavelength on the Permeability of a Viscosity Interface in the Mantle

Assuming a radially stratified Newtonian mantle in a steady-state approximation, we demonstrate that the permeability of a viscosity interface at 660-km depth strongly depends on the wavelength of buoyancy forces driving the flow. The flow induced by long-wavelength loads penetrates through the boundary freely even if the viscosity increases by two orders. In contrast, the boundary is practically impermeable for short-wavelength loads located in the upper mantle. Thus, a stepwise increase of viscosity is a significant obstacle for small descending features in the upper mantle, but huge upper mantle downwellings, or upwellings formed in the-lower mantle can overcome it easily. This indicates that certain care is necessary in interpreting the seismic structure of the mantle by means of flow models. The global tomographic image includes only the first few degrees of the harmonic series and, consequently, its interpretation in terms of a present-day flow field results in a predominantly whole-mantle circulation even for extreme viscosity contrasts.     

The Ms = 6.2, June 15, 1995 Aigion earthquake (Greece): evidence for low angle normal faulting in the Corinth rift

We present the results of a multidisciplinary study of the M_s = 6.2, 1995, June 15, Aigion earthquake (Gulf of Corinth, Greece). In order to constrain the rupture geometry, we used all available data from seismology (local, regional and teleseismic records of the mainshock and of aftershocks), geodesy (GPS and SAR interferometry), and tectonics. Part of these data were obtained during a postseismic field study consisting of the surveying of 24 GPS points, the temporary installation of 20 digital seismometers, and a detailed field investigation for surface fault break. The Aigion fault was the only fault onland which showed detectable breaks (< 4 cm). We relocated the mainshock hypocenter at 10 km in depth, 38 ° 21.7 N, 22 ° 12.0 E, about 15 km NNE to the damaged city of Aigion. The modeling of teleseismic P and SH waves provides a seismic moment M_o = 3.4 10¹⁸ N.m, a well constrained focal mechanism (strike 277 °, dip 33 °, rake – 77°), at a centroidal depth of 7.2 km, consistent with the NEIC and the revised Harvard determinations. It thus involved almost pure normal faulting in agreement with the tectonics of the Gulf. The horizontal GPS displacements corrected for the opening of the gulf (1.5 cm/year) show a well-resolved 7 cm northward motion above the hypocenter, which eliminates the possibility of a steep, south-dipping fault plane. Fitting the S-wave polarization at SERG, 10 km from the epicenter, with a 33° northward dipping plane implies a hypocentral depth greater than 10 km. The north dipping fault plane provides a poor fit to the GPS data at the southern points when a homogeneous elastic half-space is considered: the best fit geodetic model is obtained for a fault shallower by 2 km, assuming the same dip. We show with a two-dimensional model that this depth difference is probably due to the distorting effect of the shallow, low-rigidity sediments of the gulf and of its edges. The best-fit fault model, with dimensions 9 km E–W and 15 km along dip, and a 0.87 m uniform slip, fits InSAR data covering the time of the earthquake. The fault is located about 10 km east-northeast to the Aigion fault, whose surface breaks thus appears as secondary features. The rupture lasted 4 to 5 s, propagating southward and upward on a fault probably outcropping offshore, near the southern edge of the gulf. In the shallowest 4 km, the slip – if any – has not exceeded about 30 cm. This geometry implies a large directivity effect in Aigion, in agreement with the accelerogram aig which shows a short duration (2 s) and a large amplitude (0.5 g) of the direct S acceleration. This unusual low-angle normal faulting may have been favoured by a low-friction, high pore pressure fault zone, or by a rotation of the stress directions due to the possible dip towards the south of the brittle-ductile transition zone. This fault cannot be responsible for the long term topography of the rift, which is controlled by larger normal faults with larger dip angles, implying either a seldom, or a more recently started activity of such low angle faults in the central part of the rift.     

Mass heterogeneities and convection in the earth's mantle inferred from gravity and core-mantle boundary irregularities

Mass heterogeneities in the earth's mantle are retrieved from the gravity data and the topography of the core-mantle boundary as well as the topography of the earth's surface. A mantle circulation induced by the heterogeneities is modelled by solving the Stokes problem for incompressible Newtonian fluid. The derived models of mantle motions correlate well with the plate tectonics and point at a close relation between the surface tectonic activity and the processes in the vicinity of the core-mantle boundary.     

Constraints on the age of basement and crustal growth in Tianshan Orogen by Nd isotopic composition

Based on study of Nd isotopic composition for 101 rocks of various types from Tianshan Orogen, the age and character of basement and continental crustal evolution of the Tianshan Orogen were proposed. It is deduced that the continental crustal basement of the Tianshan Orogen was formed 1. 8 Ga ago. The protolith of its metamorphic rocks was derived from long-term depleted mantle source in the ancient are tectonic setting probably. The Tianshan Orogen is obviously different from the North Tarim Block in age of basement and post-evolution history. It was also shown that Paleozoic continental crustal growth happened extensively in the Tianshan Orogen, which is distinguished from Yangtse Block and Cathaysia Block in eastern China. Project supported by the National Natural Science Foundation of China (Grant No. 49633250). It belongs to the National “305” Project in Xinjiang, which is one of the National Key Projects in the Ninth Five-Year Plan (96-915-07-05A).     

First-principles investigation of the concentration effect on equilibrium fractionation of Ca isotopes in forsterite

Previous theoretical studies have found that the concentration variations within a certain range have a prominent effect on inter-mineral equilibrium isotope fractionation(10^3 lna).Based on the density functional theory,we investigated how the average Ca–O bond length and the reduced partition function ratios(10^3 lnb)and103lna of 44 Ca/40 Ca in forsterite(Fo)are affected by its Ca concentration.Our results show that Ca–O bond length in forsterite ranges from 2.327 to 2.267 A with the Ca/(Ca+Mg)varying between a narrow range limited by an upper limit of 1/8 and a lower limit of 1/64.However,outside this narrow range,i.e.,Ca/(Ca+Mg)is lower than1/64 or higher than 1/8,Ca–O bond length becomes insensitive to Ca concentration and maintains to be a constant.Because the 10^3 lnb is negatively correlated with Ca–O bond length,the 10^3lnb significantly increases with decreasing Ca/(Ca+Mg)when 1/64相似文献   

Upper mantle structure of the Apulian plate from Rayleigh waves

Phase velocities of Rayleigh waves for the Adriatic Sea area are obtained in the period range 25–190 sec along the path (l'Aquila-Trieste) AQU-TRI and 20–167 sec along the path (Trieste-Bari) TRI-BAI.The phase velocities are systematically higher than the known values for the surrounding regions. The data inversion indicates the presence of a lithosphere typical of stable continental areas with clear high-velocity lid (V _s 4.6 km/sec) overlying a well developed low velocity zone (V _s 4.2 km/sec).P. F. Geodinamica C.N.R., Roma Pubbl. N. 189.     

On the stable growth of segregated ice in freezing soil under negligible overburden pressure

The stable growth condition of a segregated ice layer is studied by the use of the principle of mass and heat conservation. This condition evidently depends upon the properties of a thin transitional zone, which is believed to exist between the boundary of an ice layer and a 0[°C] isotherm. All probable models of the transitional zone are classified and the conditions for each model is derived. The effect of the small amount of soil minerals contained in an ice layer is also studied.     

Mining soil databases for landscape‐scale patterns in the abundance and size distribution of hillslope rock fragments

Landscape‐scale variation in rock fragments on soil‐mantled hillslopes is poorly understood, despite the potential importance of rock fragments in soil weathering and coarse sediment supply to river networks. We explored the utility of soil survey databases for data mining, with the goals of identifying landscape‐scale patterns in the abundance and size distribution of rock fragments (diameter D > 2 mm) and potential controls on grain size production. We focus on data from three regions: the Hawaiian Islands, and the Sierra Nevada and Cascade Mountains, where elevation transects span a range of environmental conditions. We selected pedons from pits dug on hillslopes with active soil production and transport. For the 27 pedons selected, we constructed depth‐averaged grain size distributions and calculated the mass fraction of rock fragments (F_RF) and the median rock fragment grain size (D_50RF). We also categorized as bimodal, size distributions with a clear ‘breakpoint’ between fine and coarse modes. Several strong patterns emerge from the data. We find rock fragments in 85% of the pedons, primarily in distinct coarse modes within bimodal size distributions. Values of F_RF and D_50RF are strongly correlated, although the best‐fit power law scaling between F_RF and D_50RF differs between the warmer Hawaiian, and colder Sierra Nevada and Cascade Mountain sites. We also find a regional contrast in the variation in F_RF with elevation; F_RF declines with elevation in Hawaii, but increases in the mainland sites. Although this contrast could be an artifact of variable lithology, precipitation may influence many patterns in the data. Lower mean‐annual precipitation correlates with higher F_RF, dominantly bimodal distributions and surface enrichment in the vertical distribution of rock fragments. These observations may be useful in refining models of coarse sediment supply to rivers, and suggest opportunities for future work to test mechanistic hypotheses for rock fragment production on soil‐mantled hillslopes. Copyright © 2011 John Wiley & Sons, Ltd.     

基于TRT系统的地质构造三维成像技术及其改进方法

TRT技术是利用地震波反射原理来实现对灾害体的辨识,利用震源点与接收点共同组成的随时间变化的共焦椭球体来实现对灾害体的定位.根据TRT成果图像异常区域与隧道实际灾害的对应情况,将TRT图像异常区域划分为集群异常区、单体异常区、错断异常区、互层异常区四种类型.针对TRT技术自身存在的直达波速计算不准确、背景波速取定不合理等问题,提出了相应的改进方法,一方面通过计算出直达波在隧道洞身的实际传播的最小轨迹,求得震源点到接收点的实际直达波速;另一方面根据实际直达波速的分布特征,求出正态分布条件下,实际直达波速在置信水平为1-a的置信区间,并将置信区间的上下限值作为背景波速录入TRT系统,得出两组TRT成果图像,这样就改变了人工取定背景波速所存在的单一性与盲目性,使得改进后的灾害预测有了一定的区间考量.通过实例验证,发现改进后的预测效果明显优于改进前采用人工取定背景波速时的预测效果,并在一定程度上避免了灾害漏报的可能.     

南北地震带北段的远震P波层析成像研究

本文利用“中国地震科学台阵”探测项目在南北地震带北段布设的678个流动地震台站在2013年10月至2015年4月期间记录到的远震波形数据,经过波形互相关拾取到473个远震事件共130309条P波走时残差数据,通过远震层析成像研究获得了该区（30°N-44°N,96°E-110°E）下方0.5°×0.5°的P波速度扰动图像.结果显示,研究区下方P波速度结构显示强烈的不均一性和显著的分区、分块特征.岩石圈速度结构具有显著的东西差异：祁连、西秦岭和松潘甘孜地块组成的青藏东北缘地区显示明显的低速异常,而属于克拉通性质的鄂尔多斯地块和四川盆地则显示高速异常,表明东部克拉通块体对青藏高原物质的东向挤出起到了强烈的阻挡作用.阿拉善地块显示出弱高速和局部弱低速的异常并存的特征.阿拉善地块西部显示低速异常,而东部与鄂尔多斯相邻的地区显示高速异常,可能表明该地区的岩石圈的变形主要受到青藏高原东北缘的挤压作用.在鄂尔多斯和四川盆地之间的秦岭下方100~250 km深度上表现为明显的低速异常,表明该处可能存在软流圈物质的运移通道.鄂尔多斯北部的河套裂陷盆地下方在100~500 km深度内低速异常表现明显,说明该区有深部热物质上涌且至少来源于地幔过渡带.青藏东北缘上地幔显示低速异常且地幔过渡带中出现明显的高速异常,这种结构模式暗示了在青藏高原东北缘可能发生了岩石圈拆沉作用,而高速异常体可能是拆沉的岩石圈地幔.

     

Redox state of amphibole-bearing mantle peridotite from Nüshan, Anhui Province in eastern China and its implications

Using secondary spinel standard method, we have measured precisely the compositions of spinels of amphibole-bearing mantle peridotite xenoliths from Nüshan in eastern China, and calculated the mantle oxygen fugacities recorded by the xenoliths. Results indicate that the mantle metasomatism for forming amphiboles in Nüshan region of Anhui has resulted in the decrease of mantle redox, which is in contrast with theoretical estimation and previous research results from other areas around the world. Combining with related studies on the mantle of eastern China, we give a reasonable explanation to the ‘new finding’ and further elucidate the compositions and nature of mantle fluids in eastern China.     

Equation of state of gold and its application to the phase boundaries near 660 km depth in Earth’s mantle

The pressure-volume-temperature equation of state (EOS) of gold is fundamental to high-pressure science because of its widespread use as an internal pressure standard. In particular, the EOS of gold has been used in recent in situ multi-anvil press studies for determination of phase boundaries related to the 660-km seismic discontinuity. These studies show that the boundaries are lower by 2 GPa than expected from the depth of the 660-km discontinuity. Here we report a new P-V-T EOS of gold based on the inversion of quasi-hydrostatic compression and shock wave data using the Mie-Grüneisen relation and the Birch-Murnaghan-Debye equation. The previously poorly constrained pressure derivative of isothermal bulk modulus and the volume dependence of Grüneisen parameter (q=d lnγ/d ln V) are determined by including both phonon and electron effects implicitly: K′_0T=5.0±0.2 and q=1.0±0.1. This combined with other accurately measured parameters enables us to calculate pressure at a given volume and temperature. At 660-km depth conditions, this new EOS yields 1.0±0.2 GPa higher pressure than Anderson et al.’s EOS which has been used in the multi-anvil experiments. However, after the correction, there still exists a 1.5-GPa discrepancy between the post-spinel boundary measured by multi-anvil studies and the 660-km discontinuity. Other potential error sources, such as thermocouple emf dependence on pressure or systematic errors in spectroradiometry, should be investigated. Theoretical and experimental studies to better understand electronic and anharmonic effects in gold at high P-T are also needed.     

上地幔电导率-温度关系及在大地电磁测深法中应用综述

上地幔热结构的研究对探索地球内部物质状态和岩石圈形成演化过程, 评估自然资源的存储状况, 监控火山区岩浆活动从而降低自然灾害风险等有着重要的意义.本文概述了基于不同岩石物理学高温高压实验得到的上地幔矿物电导率-温度关系.通过分析并对比前人的实验结果, 讨论了不同实验得到的同一矿物的电导率-温度经验关系中参数差异产生的原因.大地电磁测深法(MT)以其探测深度大, 受浅部高阻体影响较小且对低阻体分辨率较高等特点, 在研究壳幔电性结构、热状态和地球动力学机制中得到了广泛的应用.以MT得到的上地幔电阻率模型为基础, 结合岩石物理学实验标定的矿物电导率-温度经验关系, 建立上地幔的温度和熔融百分比模型, 这项工作在研究上地幔热状态中起着不可或缺的作用.笔者总结并讨论了前人以MT方法获得电性结构为基础, 利用电导率-温度经验关系评价上地幔热状态的应用实例, 并对未来的研究工作做了展望, 同时对其可行性做了评估.