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     

A new view of Italian seismicity using 20 years of instrumental recordings

In this paper, we show the seismicity of the past 20 years that occurred in Italy and surrounding regions. Hypocentral locations have been obtained by using P- and S-wave arrival times from the INGV national and several regional permanent seismic networks. More than 48,000 events, selected from an original data set of about 99,780, are used to reconstruct the most complete seismic picture of the Italian region so far. The seismicity distribution allows inference on seismotectonics of this complex region of subduction versus continental collision. Our results clearly reveal the geometry of the Adria and the Ionian subduction and a continuous normal fault belt in the upper crust, following the Apennines mountain range. The depth of the seismogenic layer is computed from the cut-off of seismicity at depth and shows large variations along and across the seismic active regions. Earthquakes are generated by the different velocity of slab retreat and the subsequent asthenospheric upwelling.     

Deep structure of the northeastern Japan arc and its implications for crustal deformation and shallow seismic activity

Seismic tomography studies in the northeastern Japan arc have revealed the existence of an inclined sheet-like seismic low-velocity and high-attenuation zone in the mantle wedge at depths shallower than about 150 km. This sheet-like low-velocity, high-attenuation zone is oriented sub-parallel to the subducted slab, and is considered to correspond to the upwelling flow portion of the subduction-induced convection. The low-velocity, high-attenuation zone reaches the Moho immediately beneath the volcanic front (or the Ou Backbone Range) running through the middle of the arc nearly parallel to the trench axis, which suggests that the volcanic front is formed by this hot upwelling flow. Aqueous fluids supplied by the subducted slab are probably transported upward through this upwelling flow to reach shallow levels beneath the Backbone Range where they are expelled from solidified magma and migrate further upward. The existence of aqueous fluids may weaken the surrounding crustal rocks, resulting in local contractive deformation and uplift along the Backbone Range under the compressional stress field of the volcanic arc. A strain-rate distribution map generated from GPS data reveals a notable concentration of east–west contraction along the Backbone Range, consistent with this interpretation. Shallow inland earthquakes are also concentrated in the upper crust of this locally large contraction deformation zone. Based on these observations, a simple model is proposed to explain the deformation pattern of the crust and the characteristic shallow seismic activity beneath the northeastern Japan arc.     

临震前地壳变形不动点现象与机制研究

在综合分析国内外地形变前兆研究现状的基础上,结合汶川8.0级地震前GPS、应变、重力等研究结果,提出临震前震源区域可能出现变形不动点现象。地壳变形过程中的不动点现象是临震前区域应力场进入临界状态的重要标志,是地震潜在危险区域划定的重要依据。借助不动点理论,给出汶川地震前变形不动点集合的演化,结合大区域地震活动空区与远场显著形变异常分布,探讨变形不动点现象发生的动力学背景,验证了汶川地震前地壳变形的不动点现象是大区域应力场有序运动的结果。     

伟晶岩型锂矿成矿构造背景研究

锂是银白色的、自然界最轻的金属,广泛应用于纺织、陶瓷、冶金、核能、新能源等领域,有“高能金属”、“21世纪的能源金属”之称。锂矿床按成因类型可分为内生和外生两种类型,内生型可具体分为花岗伟晶岩型、花岗岩型、云英岩型和岩浆热液型,外生型则包括盐湖型及地下卤水型（赵一鸣等,2004）。伟晶岩型锂矿主要开采锂辉石、透锂云母、锂云母、磷铝锂石等,是锂矿资源的重要来源。     

变分同化框架通过背景误差协方差构建动力平衡约束的研究进展

通过背景误差协方差构建动力平衡约束是变分同化框架设计的重要环节。它不仅帮助实现变量间的协同分析,提高观测使用效率,还能改善变分极小化问题的性状。本文在系统梳理通过背景误差协方差引入动力平衡约束的方式、流程的基础上,对求解目前全球和有限区域变分同化系统普遍采用的准地转平衡和静力平衡约束的共性问题和存在的不足作了归纳总结。分析了求解准地转平衡约束的三类方案:动力平衡方程方案、统计方案和动力-统计相结合方案的优缺点。对照比较了不同垂直离散方案下求解静力平衡约束时遇到的欠定问题的表现以及解决途径。最后,展望了基于背景误差协方差构建动力平衡约束在赤道等特殊地区、高分辨率同化系统、以及集合-变分混合同化系统发展中面临的挑战和机遇。