Magma ascent beneath Unzen Volcano, SW Japan, deduced from the electrical resistivity structure

The resistivity structure of Unzen Volcano has been revealed by extensive magnetotelluric surveys since the first eruption on November 17, 1990. This structure comprises a highly resistive surface layer, a low-resistive second layer at several hundred meters depth, interpreted as a water-saturated layer, a resistive third layer, and a low-resistive fourth layer at 10 km depth, possibly related to the deep magmatic activity. The structure has influenced the volcanic activity of Unzen. This activity was characterized by a series of dramatic changes in eruption type: a minor phreatic eruption on November 17, 1990; phreatic eruptions after February 12, 1991, preceded by several weeks of volcanic tremor; phreatomagmatic eruptions after April 9, and dome effusion beginning May 19, 1991. This paper presents a hypothesis in which the top of the magma column rose about 20 m/day, reached the base of the water-saturated layer at the end of January, 1991, and approached the upper boundary of this layer on April 9. Thus, the temporal change of eruption type and associated phenomena are systematically explained by an interaction between magma and groundwater contained in the saturated layer.     

长白山天池火山全新世喷发与岩石地球化学特征

全新世以来长白山天池火山以爆炸式喷发为特点,其喷发物包括空降堆积浮岩、熔结凝灰岩和碎成熔岩,除了极少量粗面质成分（SiO2<65％）,均属于碱流质火山碎屑岩（SiO2＝70%～73％,K2O+Na2O=10%～11%）。火山岩的矿物学（造岩矿物和残余矿物）和主微量元素、SrNd同位素地球化学研究表明,天池火山全新世火山岩与造锥阶段、造盾阶段火山岩是同源（幔源）岩浆演化的产物,全新世火山岩都具有强烈而相似的Eu、Ba、Sr、Ti负异常,指示了岩浆演化主要受斜长石的强烈分离结晶作用控制。     

中国东部大同火山群发育的构造地貌背景

大同火山群是我国东部重要独特地貌类型之一。火山群及其所在的大同盆地位于上地幔隆起带上,受中国东部最突出的构造地貌格局——NNE向雁行右旋剪切拉张断陷带控制。火山群及熔岩溢出带即明显沿着这个右旋剪切拉张破裂带以及与NW压扭性断裂交汇处发育,同时受到断裂多次活动的影响。     

Numerical inversion of deformation caused by pressure sources: application to Mount Etna (Italy)

The interpretation of geodetic data in volcanic areas is usually based on analytical deformation models. Although numerical finite element (FE) modelling allows realistic features such as topography and crustal heterogeneities to be included, the technique is not computationally convenient for solving inverse problems using classical methods. In this paper, we develop a general tool to perform inversions of geodetic data by means of 3-D FE models. The forward model is a library of numerical displacement solutions, where each entry of the library is the surface displacement due to a single stress component applied to an element of the grid. The final solution is a weighted combination of the six stress components applied to a single element-source. The pre-computed forward models are implemented in a global search algorithm, followed by an appraisal of the sampled solutions. After providing extended testing, we apply the method to model the 1993–1997 inflation phase at Mt Etna, documented by GPS and EDM measurements. We consider four different forward libraries, computed in models characterized by homogeneous/heterogeneous medium and flat/topographic free surface. Our results suggest that the elastic heterogeneities of the medium can significantly alter the position of the inferred source, while the topography has minor effect.     

Automatic classification of dome instabilities based on Doppler radar measurements at Merapi volcano, Indonesia: Part I

Monitoring lava dome instabilities is crucial to efficiently monitor active dome building volcanoes. The Doppler radar technique provides a unique opportunity to gather information about the number of instability events occurring at the growing dome and about the dynamic processes that take place during different types of instabilities. So far, three different kinds of processes have been identified: sliding material, gravitational break-offs and explosive outbursts. In addition, Doppler radars provide rain measurements, which can be used to investigate possible correlations between rainfall and dome activity. Two radar systems have been installed at Merapi volcano in October 2001 and January 2005 to continuously monitor dome instabilities. Due to the large number of instability events that occur during times of high activity, manual processing and analysis of instability events is not practical for monitoring purposes. Therefore, an automatic classification system has been developed, which is capable of identifying different kinds of instabilities as well as rainfall. Two different kinds of classifier models have been applied: (1) neural network and (2) K-nearest-neighbour classifier model. Both classify Doppler spectra according to the underlying dynamic process, that is, rain, sliding material, gravitational break-off or explosive outburst. The classifiers are able to identify disturbances, which have no physical source, but are merely artefacts from the radar device itself. Because radar events are sequences of Doppler spectra, a rule set has been defined, which finally determines the event class. All classifiers have been trained and tested on independent data sets to estimate the classification performance. The overall classification rate is about 90 per cent. Discrimination of instabilities and non-volcanic events reaches about 98 per cent accuracy.     

Stable inversions for complete moment tensors

The seismic moment tensors for certain types of sources, such as volcanic earthquakes and nuclear explosions are expected to contain an isotropic component. Some earlier efforts to calculate the isotropic component of these sources are flawed due to an error in the method of Jost & Herrmann. We corrected the method after Herrmann & Hutchensen and found great improvement in the recovery of non-double-couple moment tensors that include an isotropic component. Tests with synthetic data demonstrate the stability of the corrected linear inversion method, and we recalculate the moment tensor solutions reported in Dreger et al. for Long Valley caldera events and Dreger & Woods for Nevada Test Site nuclear explosions. We confirm the findings of Dreger et al. that the Long Valley volcanic sources contain large statistically significant isotropic components. The nuclear explosions have strikingly anomalous source mechanisms, which contain very large isotropic components, making it evident that these events are not tectonic in origin. This indicates that moment tensor inversions could be an important tool for nuclear monitoring.     

Synchronizing seismic networks with ambient noise

The technique of retrieving Green's functions by pairwise correlation of a random wavefield is used to estimate the stability of the timing system of a seismic network. The method is demonstrated with data from a network consisting of stations with four independent clocks and GPS-receivers. For a two years period, daily measurements of the timing error of each of the four systems are obtained. These measurements can be used to correct the timing of the data in periods without GPS reception. The accuracy of this method can exceed the precision of the internal station clocks already after a single day of uncorrected drift.     

2002~2005年长白山火山气体的释放特征与地球化学异常研究——多源高光谱遥感证据

长白山火山区是我国具有潜在喷发危险的活火山,在2002~2005年火山活动性增强出现了岩浆房扰动。利用卫星遥感技术具有观测范围大、观测时间长且连续的优势。因此,本文利用对流层污染探测仪（MOPITT）和大气红外探测仪（AIRS）高光谱遥感数据提取了2002~2005年长白山天池火山区CO总量、O₃总量、水汽总量和地表温度异常信息,讨论了高光谱遥感气体地球化学异常信息与火山活动性之间的关系,并对2002~2005年长白山天池火山区火山活动性进行了研究。结果表明,高光谱遥感数据观测到的气体地球化学（CO、O₃、水汽）异常与地震、形变监测结果以及地面流体（CO₂、He、H₂）观测结果相一致,表明MOPITT和AIRS高光谱遥感卫星观测到的气体异常变化较好地反映了2002~2005年大规模的深部岩浆局部扰动。在2002~2005年火山活动期间,CO总量、O₃总量、水汽总量和地表温度均出现了显著异常且异常出现时段相应准偏差显著增大,反映了气体逸出量在时间上具有不均一性,可能与火山、地震活动过程中地应力的作用和变化有关。从气体异常持续的时间以及地面观测结果综合分析,2002~2005年岩浆房扰动并没有产生长时间的地幔物质流的上升和迁移。在火山活动性增强的时间段内,地表温度出现异常低值,这可能与太平洋板块俯冲过程中引起的断裂拉张增强有关。此外,火山活动过程中逸出的气体进入大气圈产生大气化学反应也会导致高光谱遥感所观测到的气体地球化学异常。研究结果为2002~2005年长白山火山活动性的研究提供了来自高光谱遥感数据的气体地球化学新证据,也对高光谱分辨率遥感数据在火山活动规律的研究以及火山监测中的应用具有一定意义。     

Evolution of Mount Fuji,Japan: Inference from drilling into the subaerial oldest volcano,pre‐Komitake

A buried, old volcanic body (pre‐Komitake Volcano) was discovered during drilling into the northeastern flank of Mount Fuji. The pre‐Komitake Volcano is characterized by hornblende‐bearing andesite and dacite, in contrast to the porphyritic basaltic rocks of Komitake Volcano and to the olivine‐bearing basaltic rocks of Fuji Volcano. K‐Ar age determinations and geological analysis of drilling cores suggest that the pre‐Komitake Volcano began with effusion of basaltic lava flows around 260 ka and ended with explosive eruptions of basaltic andesite and dacite magma around 160 ka. After deposition of a thin soil layer on the pre‐Komitake volcanic rocks, successive effusions of lava flows occurred at Komitake Volcano until 100 ka. Explosive eruptions of Fuji Volcano followed shortly after the activity of Komitake. The long‐term eruption rate of about 3 km³/ka or more for Fuji Volcano is much higher than that estimated for pre‐Komitake and Komitake. The chemical variation within Fuji Volcano, represented by an increase in incompatible elements at nearly constant SiO₂, differs from that within pre‐Komitake and other volcanoes in the northern Izu‐Bonin arc, where incompatible elements increase with increasing SiO₂. These changes in the volcanism in Mount Fuji may have occurred due to a change in regional tectonics around 150 ka, although this remains unproven.     

中国东北地区远震P波走时层析成像研究

利用中国东北流动和固定台网的234个宽频带地震仪记录的远震波形数据,采用波形相关方法拾取了57251个有效相对走时残差数据,进一步采用FMTT(Fast Marching Teleseismic Tomography)层析成像的方法,反演获取了研究区下方深达800 km的P波速度结构.结果显示:在长白山下方发现有一个高速异常结构,这可能就是俯冲到欧亚大陆板块下方的太平洋板块,由于板块的部分下沉,使得板块的形状并没有呈现出明显的板片状.长白山、阿尔山、五大连池火山下方都有低速异常体,长白山和阿尔山下的低速异常向下延伸至地幔转换带,可能与其上部的火山形成有关.五大连池火山下方的低速异常向下延伸至200 km左右,不同埋深的低速异常结构可能意味着五大连池与长白山和阿尔山有着不同的成因.松辽盆地呈现以高速异常为主导高低速异常混合分布的特性,暗示松辽盆地可能有岩石圈拆沉的过程,盆地南部下方的低速异常与长白山和阿尔山下的低速异常有连通性,可能是下地幔热物质上涌的一个通道.