Image processing of 2D resistivity data for imaging faults

A methodology to locate automatically limits or boundaries between different geological bodies in 2D electrical tomography is proposed, using a crest line extraction process in gradient images. This method is applied on several synthetic models and on field data set acquired on three experimental sites during the European project PALEOSIS where trenches were dug. The results presented in this work are valid for electrical tomographies data collected with a Wenner-alpha array and computed with an l₁ norm (blocky inversion) as optimization method. For the synthetic cases, three geometric contexts are modelled: a vertical and a dipping fault juxtaposing two different geological formations and a step-like structure. A superficial layer can cover each geological structure. In these three situations, the method locates the synthetic faults and layer boundaries, and determines fault displacement but with several limitations. The estimated fault positions correlate exactly with the synthetic ones if a conductive (or no superficial) layer overlies the studied structure. When a resistive layer with a thickness of 6 m covers the model, faults are positioned with a maximum error of 1 m. Moreover, when a resistive and/or a thick top layer is present, the resolution significantly decreases for the fault displacement estimation (error up to 150%). The tests with the synthetic models for surveys using the Wenner-alpha array indicate that the proposed methodology is best suited to vertical and horizontal contacts. Application of the methodology to real data sets shows that a lateral resistivity contrast of 1:5–1:10 leads to exact faults location. A fault contact with a resistivity contrast of 1:0.75 and overlaid by a resistive layer with a thickness of 1 m gives an error location ranging from 1 to 3 m. Moreover, no result is obtained for a contact with very low contrasts (1:0.85) overlaid by a resistive soil. The method shows poor results when vertical gradients are greater than horizontal ones. This kind of image processing technique should be systematically used for improving the objectiveness of tomography interpretation when looking for limits between geological objects.     

Crust and upper mantle structure of the Ailao Shan-Red River fault zone and adjacent regions

Using arrival data of the body waves recorded by seismic stations, we reconstructed the velocity structure of the crust and upper mantle beneath the southeastern edge of the Tibetan Plateau and the northwestern continental margin of the South China Sea through a travel time tomography technique. The result revealed the apparent tectonic variation along the Ailao Shan-Red River fault zone and its adjacent regions. High velocities are observed in the upper and middle crust beneath the Ailao Shan-Red River fault zone and they reflect the character of the fast uplifting and cooling of the metamorphic belt after the ductile shearing of the fault zone, while low velocities in the lower crust and near the Moho imply a relatively active crust-mantle boundary beneath the fault zone. On the west of the fault zone, the large-scale low velocities in the uppermost mantle beneath western Yunnan prove the influence of the mantle heat flow on volcano, hot spring and magma activities, however, the upper mantle on the eas     

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.     

拉萨地体内低速异常成因的探讨--壳内局部熔融的地震层析证据

青藏高原南部天然地震体波的层析反演图像描绘出拉萨地体内在15～90km深度内低速异常的分布特征。在研究速度结构的基础上，根据熔融程度与纵波波速呈倒数衰减的特性，论证了拉萨地体中低速体的出现与壳内的局部熔融有关。壳内局部熔融则来源于板块俯冲、地壳增厚产生的热而非幔源侵入所造成。同时，也发现了在拉萨以东中下地壳低速体有向东发展的趋势。大面积部分熔融层分布的深度范围约在15～35km处。另外，据藏南地区的SKS波分裂时间延迟很小的特点，推测藏南地幔中可能不存在各向异性，也就是没有足够的热源使各向异性矿物的品格重新定位。这从另一方面佐证了部分熔融仅存在于地壳内而不是地幔中。     

极化干涉SAR森林树高提取算法的对比分析

在对极化干涉SAR森林树高反演的DEM差值算法、相干相位-幅度综合反演算法进行分析的基础上,对基于极化干涉相干优化方法的改进算法进行了探讨。利用黑龙江大兴安岭地区的一对ALOS全极化干涉SAR数据进行实验,并对比分析各算法的反演结果。结果表明,在使用改进的算法进行森林树高反演时可以获取精度较高的反演结果,并且在一定程度上提高了森林树高反演的稳定性,为森林树高反演工作的业务化运行提供一定的依据。