The Gravity Field of the Vogtland and NW Bohemia: Presentation of a Project

The Vogtland and NW Bohemia are characterized geoscientifically by periodically occurrence of swarm earthquakes. The basic geophysical mechanism is not yet sufficiently clarified, just like detail questions to geology in especially the deeper underground. Complex geophysical investigations in the seismoactive region indicate geodynamic phenomena like mass redistribution or stress accumulation and release (Spiák et al., 1998). According to Grünthal (1989) a weakness zone is suggested in the region of the swarm earthquakes. This zone can be caused by fluid-tectonics (Kämpf et al., 1992), a mantle plume (pers. com. J. Svancara, 1999) and/or by the geometry of the geological structures (Neunhöfer & Güth, 1988). A three-dimensional gravimetric model can clear up the underground situation. By means of high-resolution gravimetry a three-dimensional model will be developped for the Vogtland and NW Bohemia region. In the first step a homogeneous Bouguer map of the Vogtland and NW Bohemia was created (fig. 1) containing gravity structures analysed by Ibrmajer & Suk (1989) and Blízkovsky et al. (1985). The used gravimetric data were made available by the Saxonian National Office for Environment and Geology, by the Czech Geological Survey, Prague and by the GGA Hannover. In the context with the interpretation of the deep-seismic profile MVE 90 a two-dimensional gravimetric modeling was carried out (Behr et al., 1994), too. Anomaly-producing source bodies apparently do not offer themselves in a two-dimensional model, because after Jung (1961) the length of a gravimetric source structure must be about four times larger than it's width. The technique of the three-dimensional gravimetric modeling by means of any polyhedrons was developed by Götze (1976, 1984). Gravimetry is a potential method and supplies an infinite number of solutions, so the model has to be developed close to other geoscientific results. The aim is to construct a high-resolution three-dimensional underground model, which includes the upper earth's crust and the deep-seated structures of the middle and lower crust, too. The determination of the mass distribution in the underground supplies contradicting or supporting facts for geodynamic views in the Vogtland and NW Bohemia for example of Bankwitz et al. (1993). The interpretation of the Bouguer map of the Vogtland and a three-dimensional gravimetric model ought to contribute a substantial, also geodynamic part to understand the origin and the emergence of the swarm earthquakes in this region.     

新疆精河6.6级地震周边地区密度构造、均衡异常以及岩石圈挠曲机理

依据EIGEN-6C4重力模型和ETOPO1高程模型数据,围绕新疆精河6.6级地震展开岩石圈均衡与挠曲机理研究,得到如下结论：（1）震中附近的布格与自由空气重力异常分别为-221和-92mGal(10~(-5 )m·s~(-2)）,震中位于重力异常高梯度带上;（2）震中周边地区地壳厚度约为50km,密度结构总体变化平缓,东西方向地壳厚度变化较小,但自南向北地壳厚度逐渐变薄,精河6.6级地震初始破裂发生在上中地壳分界面附近;（3）震中附近岩石圈承载的垂向构造应力为20MPa左右,震中位于岩石圈垂向构造应力极大值附近的高梯度带上;（4）地震周边地区岩石圈有效弹性厚度最优解为26km,加载比最优解为F_1=1,F_2=F_3=0,表明该区域岩石圈相对坚硬,且导致岩石圈变形的初始加载全部来自地表．     

青藏高原东南缘多尺度重力场及构造含义

本文基于EGM2008重力场模型研究了青藏高原东南缘均衡重力异常和多尺度的布格重力异常特征,以鲁甸和景谷地震为例,认识其深部构造环境和动力学过程,为该区域的构造运动和地震孕育环境研究提供依据。结果表明,研究区布格重力异常和均衡重力异常与地质构造格局相关性较好,川滇地块剧烈的区域布格重力异常和非均衡状态与其强烈的地壳变形、断裂及地震活动密切相关。强震多分布在断裂带两侧重力异常的过渡地带和高梯度带,断裂带两侧横向和垂向的显著介质密度差异是强震孕育的深部构造背景。布格重力异常和均衡重力异常揭示的鲁甸、景谷震源区深浅差异性的重力异常特征,暗示鲁甸和景谷地震孕震环境的不同。     

巴颜喀拉地块东部及邻区重力均衡与岩石圈有效弹性厚度

基于SIO(Scripps Institute of Oceanography)最新全球重力和高程模型,计算了巴颜喀拉地块东部及邻区的布格重力异常、均衡重力异常、岩石圈有效弹性厚度及荷载比.结合大地热流、地震速度结构、地震活动和断裂构造分布等,分析了地壳均衡状态和岩石圈有效弹性厚度、地质构造单元间的差异及与地震活动的相关性特征.研究结果表明,该区域布格重力变化范围约为-500～0mGal(1mGal=10~(-5)m·s~(-2),下同),在巴颜喀拉块体东部区域形成弧形重力梯度带,近年来的中强地震活动频发于该梯度带不同部位,应与其应力依次释放有关;均衡重力异常结果表明,其变化范围约为-80～+100mGal,且大部分区域处于±20mGal以内的被认为处于重力均衡的状态,重力非均衡(正或负)多出现于块体边界带附近,地震多发生在靠近块体边界的均衡重力异常(正或负,主要为正)区域内;巴颜喀拉地块东部及邻区岩石圈有效弹性厚度(T_e)为10～65km,不同构造单元之间T_e空间分布差异明显,较低的T_e值出现在龙门山构造带附近,T_e值为20km左右,岩石圈荷载加载比为0.5～0.8,表明现今的岩石圈挠曲状态主要由莫霍面加载形成.进一步分析表明,巴颜喀拉地块东部挤压增生与横向流动同时发生,是造成该区域地震发生与重力均衡异常高值重合、岩石圈有效弹性厚度和大地热流值较低的主要原因.本文获得的地壳均衡特征及岩石圈有效弹性强度结果,加深了对巴颜喀拉东部及邻区岩石圈构造演化过程的认识.     

Geodynamic modelling of the recent stress and strain field in the Vogtland swarm earthquake area using the finite-element-method

The study of several swarm earthquake regions led us to the assumption that swarm quakes mostly occur in combination with fluid migration and often in the neighbourhood of volcanoes. Observations of the strong carbon dioxide flux and investigations of Weinlich et al. (1999) about recent volcanism in the Vogtland/ Western Bohemia are the linking part to global observations of swarm earthquakes. This idea and the theoretical modelling of Yamashita (1999) were the starting point to develop a model using the finite-element-method with the software package ABAQUS concerning stress field, pore pressure, temperature and deformation in the Vogtland swarm quake area. The model is intended to reveal the interaction between fluid migration, temperature changes and deformation in this area. Its horizontal dimensions are 55 km×60 km and 30 km in depth. The main geological units include the Marianske Lazne Fault Zone and the Eger Rift. We focus on that part of the Vogtland/Western Bohemia region where these two fault zones cross. The model takes into consideration the stress field, Mohr–Coulomb failure, thermal stresses as well as creep and poro-elastoplastic rheologies. The results show that the rate of deformation and the stress accumulation caused only by the regional stress field is not high enough to generate the swarm earthquakes. But periodic change of pore pressure (which includes fluid migration) in combination with temperature changes has a severe influence on the rate of deformation and is in the dimension to cause swarm quakes.     

渭河盆地及邻区重力异常小波多尺度分解与解释

基于EGM2008重力场模型计算获得了渭河盆地及邻区布格重力异常。采用小波多尺度分解方法对布格重力异常进行了4阶小波逼近和小波细节分解,同时基于平均径向对数功率谱方法定量化地计算出1~4阶小波细节和小波逼近所对应的场源平均埋深。结合区域地质和地震资料,对获得的重力场结果进行分析,得到如下结论:①鄂尔多斯地块、渭河盆地、秦岭造山带3个一级构造单元的布格重力异常之间存在明显差异;构造区内部重力异常也存在横向的显著差异。布格重力异常的走向、规模、分布特征与二级构造区及主要的断裂具有一定的对应关系。②渭河盆地及邻区布格重力异常1~4阶细节对应4~23 km不同深度的场源信息,鄂尔多斯地块南缘东、西部的地壳结构存在明显的差异;渭河盆地凹陷、凸起构造区边界清晰,断裂边界与重力异常边界具有较好的一致性;秦岭造山带重力异常连贯性不好,东、西部重力异常变化特征表现出明显的差异。③渭河盆地及邻区布格重力异常分布与莫霍面埋深具有非常明显的镜像关系。渭河盆地及邻区地震主要分布在六盘山—陇县—宝鸡断裂带、渭河断裂与渭南塬前断裂交汇处、韩城断裂与双泉—临猗断裂交汇处。渭河盆地及邻区重力异常主要由中上地壳剩余密度体所影响,这可能是该区地震以浅源地震为主的主要原因。     

Optimal Design of a Gravity Survey Network and its Application to Delineate the Jabera-Damoh Structure in the Vindhyan Basin,Central India

Fractal dimension analysis was carried out for optimal designing of 2-D gravity survey network and to determine an optimum range of gridding interval to generate least aliased Bouguer anomaly maps. As a test case, this method has been successfully applied to the Jabera-Damoh region of the Vindhyan Basin, which is considered as a potential hydrocarbon bearing area. In particular, we aim to delineate accurately the lateral extent of a possible hydrocarbon bearing structure. To achieve this aim, fractal dimension of survey network was computed using 2-D distributions of observation points in the planning phase of the survey so that the optimum station spacing for gravity survey can be obtained. A range of optimum gridding interval for the gravity data set was suggested using the box-counting method of fractal dimension determination. Bouguer anomaly maps of the region are prepared utilizing the optimum gridding interval. For the first time, these anomaly maps clearly outline the gravity evidence of an anomalous rifted structure, which is bounded by parallel faults on either side. This structure is interpreted as a favorable basin for the occurrence of hydrocarbons. Another finding of this study has been the delineation of an apparently small ridge-like structure running east-west, dividing this basin in two parts. A subsurface geological model along a profile across the Jabera structure has also been presented.     

Gravimetric study of a potential mineral deposit in the Itapororoca region, Brazil

A semi-detailed gravity survey was carried out over an area of 650 km² localized in the Eo-Neoproterozoic coastal zone of Paraiba State where 548 new gravity stations were added to the existing database. Gravity measurements were made with a LaCoste and Romberg model G meter with a precision of 0.04 mGal. The altitude was determined by barometric levelling with a fixed base achieving a 1.2 m measure of uncertainty, corresponding to an overall accuracy of 0.24 mGal for the Bouguer anomaly. The residual Bouguer map for a 7th degree regional polynomial showed a circumscribed negative anomaly coincident with a localized aero-magnetic anomaly and with hydro-thermally altered outcrops, near the city of Itapororoca. The 3D gravity modelling, constrained by geologic mapping was interpreted as a low density, fractured and/or altered material with a most probable volume of approximately 23 km³, extending to about 8,500 m depth. This result is in accordance with a volcanic body associated with hydrothermal processes accompanied by surface mineralization evidence, which may be of interest to the mining industry.     

Interpretation of gravimetric and magnetic anomalies in the Cameros Basin (North Spain): combination of deep and shallow sources

Outstanding potential field anomalies (gravimetric and magnetic) in the Cameros Basin (N Spain) follow a WNW-ESE trend, parallel to the geological structures resulting from Mesozoic extension and Tertiary basin inversion. The positive Bouguer gravity anomaly (15 mGal) is interpreted as the result of a strong contrast between the density of Tertiary rocks of the foreland basin and the Paleozoic and Mesozoic rocks, combined with crustal thickening in the Iberian Chain with respect to the Ebro Basin. The dipolar magnetic anomaly, slightly shifted to the south with respect to the relative maximum of the Bouguer anomaly, can be interpreted as related to volcanic rocks within the basement, which are linked to Triassic rifting as witnessed by outcrops of basalts along the basin margins. An exhaustive analysis of rock properties (density, magnetic susceptibility and remanence) and basin geometry from other sources (seismic reflection profiles) allow to constrain variations in crustal thickness and the location of large-scale basement faults.     

维西—贵阳剖面重力异常与地壳密度结构特征

维西—贵阳剖面位于青藏高原东南缘,为青藏高原物质往东南逃逸、东构造结侧向挤压及华南地块北西西向推挤作用的重要地段.利用剖面观测的重力与GPS定位数据,结合区域背景重力场、地质构造及深部地球物理成果,反演研究剖面较为细化的地壳密度结构特征.观测研究表明:剖面布格重力异常总幅差变化达190×10-5 m·s-2,具"斜N"分段变化特征,从西往东呈上升(维西至攀枝花,水平梯变大)—下降(攀枝花至会泽,水平梯变较大)—上升(会泽至贵阳,水平梯变较小)态势;高程与布格重力异常比值的趋势性转折部位为康滇地轴核心和小江断裂带东侧,可能与先存构造或新生构造发育有关;剖面地壳密度结构可分上、中和下三层结构,各层底界面平均埋深分别约20km、35km和51km,金沙江—红河断裂带和鲜水河—小江断裂带为地壳结构相对简单与复杂的过渡带;地壳厚度西深东浅,可能是东构造结的侧向挤压所致;下地壳厚度变化相对较大,可能对地壳增厚起主要作用;华坪—攀枝花附近的Moho面隆起和上地壳高密度体的存在暗示上地幔往上底侵作用,对青藏高原物质向南东逃逸和东构造结的侧向挤压均起到一定阻挡作用;中地壳下伏有限低密度薄层有利于其上物质的南东逃逸和顺时针旋转,有利于其下物质受喜马拉雅东构造结作用下往东向运移.     

Location of secondary faults from cross-correlation of the second vertical derivative of gravity anomalies1

A method is described to locate secondary faults, which can be difficult to identify on the Bouguer gravity map. The method is based on cross-correlation between the theoretical anomaly due to a vertical step and the second vertical derivative of the Bouguer anomaly. Faults are located from the closed maxima and minima on the cross-correlation contour map calculated for two perpendicular directions. One-dimensional model computations show that the magnitude of the extremum of the cross-correlation is related to the depth to the top of the hanging wall and the throw of the fault. Application of the method to the Bouguer gravity map of the former mouth of the Yellow River in the Shengli Oilfield area near the Bo Hai Sea shows the effectiveness of the method.     

Seismohydrological effects related to the NW Bohemia earthquake swarms of 2000 and 2008: Similarities and distinctions

For more than 20 years, seismohydrological investigations have been undertaken at the mineral aquifer system of Bad Brambach (Vogtland, Germany). Two strong swarm earthquake series in 2000–2001 and 2008–2009 at the Nový Kostel epicentre (Czech Republic, 10 km E of BB) have enabled for the first time a comparison of seismological and groundwater hydraulic features in a semi-quantitative way. In spite of their similar spatial distribution in 2001 and 2008, the earthquake foci of each swarm migrated differently through time, horizontally as well as in depth. The seismic energy of the 2008–2009 events was released predominantly within 1 month, in contrast to 2000–2001 when it occurred over 3 months. The main distinctive features of each are seen in the hydraulic pressure anomalies which accompanied the earthquake swarms: number, shape, and progression (duration) of the anomalies. The comprehensive hydraulic data, with high temporal resolution, suggest that fluid triggering dominated not only the earthquake initiating phases. In particular, the long-lasting seismicity of the 2008–2009 swarm can be attributed to a continued triggering of weak earthquakes by over-pressured deep fluids. Here, the remaining static strain was obviously not sufficient to generate strong earthquakes as at the beginning of the earthquake swarm periods. Furthermore, the enduring high fluid pressure in 2009 could also indicate a continuation of the long-term gas flow increase observed at several gas outlets in the Vogtland/NW Bohemia region between 1998 and 2008. However, it is not possible at present to derive a systematic relationship between anomaly occurrence and seismic activity, as generally proposed in the context of earthquake prediction discussion.     

Gravity anomaly, lithospheric structure and seismicity of Western Himalayan Syntaxis

A compiled gravity anomaly map of the Western Himalayan Syntaxis is analysed to understand the tectonics of the region around the epicentre of Kashmir earthquake of October 8, 2005 (Mw = 7.6). Isostatic gravity anomalies and effective elastic thickness (EET) of lithosphere are assessed from coherence analysis between Bouguer anomaly and topography. The isostatic residual gravity high and gravity low correspond to the two main seismic zones in this region, viz. Indus–Kohistan Seismic Zone (IKSZ) and Hindu Kush Seismic Zones (HKSZ), respectively, suggesting a connection between siesmicity and gravity anomalies. The gravity high originates from the high-density thrusted rocks along the syntaxial bend of the Main Boundary Thrust and coincides with the region of the crustal thrust earthquakes, including the Kashmir earthquake of 2005. The gravity low of HKSZ coincides with the region of intermediate–deep-focus earthquakes, where crustal rocks are underthrusting with a higher speed to create low density cold mantle. Comparable EET (∼55 km) to the focal depth of crustal earthquakes suggests that whole crust is seismogenic and brittle. An integrated lithospheric model along a profile provides the crustal structure of the boundary zones with crustal thickness of about 60 km under the Karakoram–Pamir regions and suggests continental subduction from either sides (Indian and Eurasian) leading to a complex compressional environment for large earthquakes.     

青藏高原东北缘重力场深部结构及其动力学特征

本文采用欧拉反褶积、场源参数成像(SPI)、场源边界提取(SED)、莫霍面反演、地壳三维可视化等多源方法,对青藏高原东北缘地区的布格重力场进行反演与分析,深入研究该地区的深部结构与变形特征,探讨区域深部孕震环境及动力学机制.研究表明,青藏高原东北缘的布格重力场整体呈负异常值,具有明显的分区性,表现出鄂尔多斯盆地异常值相对偏高、阿拉善块体次之、青藏高原块体极低的特点,其中海源断裂系形成了一条宽缓的弧形重力梯度条带,梯度值达1.2 mGal·km^-1.欧拉结果显示,鄂尔多斯盆地相比于青藏高原块体而言,场源点具有较强的均一性,场源强度值高(密度值高)且深度稳定在25~32 km范围内,而高原块体的中下地壳尺度广泛分布着低密度异常体.SPI图可知,海源弧形断裂系位于“浅源异常”弧形区,反映其地壳较为活跃,易发生中强地震.SED图揭示青藏高原地壳向东北扩展,经过几大断裂系的调节后运动矢量向东或东南转化,SED与GPS、SKS运动特征大致相同,说明地表-地壳-地幔的运动特征有着较强的一致性.青藏高原东北缘地区壳幔变形是连贯的,加之莫霍面由北向南、由东向西是逐渐加深的,因此属于垂向连贯变形机制,不符合下地壳管道流动力学模式.区域形成了似三联点构造格局,其中海源弧形断裂系的深部地壳结构复杂,高低密度异常体复杂交汇,是青藏高原、阿拉善、鄂尔多斯三大块体相互作用的重要枢纽,其运动学特征总体为中段走滑尾端逆冲,而断裂系正处于大型的弧形莫霍面斜坡带之上,具备强震的深部孕震环境,因此大尺度的运动调节与深部孕震条件共同促使了该地区中强震的多发.     

Gravity variations along the Southeast Bohol Ophiolite Complex (SEBOC), Central Philippines: Implications on Ophiolite Emplacement

Abstract The basement complex of Bohol Island consists of the Southeast Bohol Ophiolite Complex (SEBOC), Cansiwang Melange and Alicia Schist. The SEBOC is a complete, but dismembered ophiolite with outcrops generally trending northeast– southwest and dipping north-west. The harzburgite units of the SEBOC are almost always observed to be thrusted onto the Cansiwang Melange, which in turn is thrusted onto the Alicia Schist. Bouguer gravity values on Bohol range from about +60 mGal in the west to +120 mGal in the east, in the region to the north-east of the SEBOC outcrops. Based on the present distribution of the SEBOC units and their thrust fault relationship with the Cansiwang Melange and Alicia Schist, it is proposed that the SEBOC was emplaced by onramping towards the south-eastward direction. However, the orientation of the Bouguer highs suggests that the thrusting direction of the ophiolite units is towards the south-west and not towards the south-east.     

重震反演中国东北地壳上地幔三维密度结构

本文利用重力和地震P波到时数据反演得到了中国东北地区地壳上地幔三维密度结构.与单一的重力或地震反演相比,重震反演一方面有效地克服了重力反演结果垂向分辨率低的问题,另一方面也提高了地震反演结果的可靠性.结果显示:中国东北地区的地壳及上地幔剩余密度异常分布与构造单元具有明显的相关性,造山带对应低密度异常,盆地对应高密度异常;区域内火山下方有明显的低密度体存在,可能是由于太平洋板块俯冲进入上地幔并部分滞留,在滞留板块深部脱水和软流圈热物质共同作用下产生了上涌岩浆,喷发后形成了火山.     

背景噪声面波与布格重力异常联合反演:山西断陷带三维地壳结构

讨论了利用面波与布格重力异常联合反演三维地壳速度结构的新方法,并利用该方法联合反演获得山西断陷带地壳S波速度结构.通过建立速度与密度之间的经验关系,利用非线性迭代反演方法获得最终速度模型.结果显示,联合反演获得的速度模型可以同时提高对面波及重力数据的观测拟合程度,而面波单独反演得到的速度模型则无法很好的拟合重力观测数据.相比较,联合反演速度模型中的大同火山区中下地壳的低速异常幅值小于面波单独反演模型中低速异常体的幅值.联合反演速度模型结果揭示,吕梁山地区在中下地壳存在低速异常,并且和北部的大同火山区低速异常相连接,说明可能导致新生代以来大同火山区岩浆活动的上地幔构造活动(上地幔局部上涌,地幔柱)可能对山西断陷带的形成和构造活动起到了一定的控制作用,并且导致了吕梁山地区中下地壳的低速异常.     

喜玛拉雅“东构造结”地区特异重力场的探讨

跨越中、印、缅三国交界的青藏高原东南的喜玛拉雅“东构造结”地区（92°E~97°E，26°N~30°N）一半以上的面积尚没有重力测点，是重力数据空白区,故无法直接研究其重力场特征与深部地壳结构（构造）.本文分析了卫星重力异常的特性，提出应用卫星重力异常作为近似空间重力异常，并作布格改正后，得到的布格重力异常具有与该地区地形高程呈镜像相关的特征，可用以研究深部地壳结构.据三条重力剖面计算得到该地区三个地壳深部结构剖面的结果，给出青藏高原地壳厚度>70 km；喜马拉雅造山带为55 km左右；布拉马普特拉河谷盆地为33~35 km；那加山山脉地区为40~45 km，显示出三者为三个不同的构造单元.同时给出布拉马普特拉构造单元为相对高密度的刚性物质构成，随着印度洋板块向北运移，在碰撞、挤压下，插入青藏高原东南缘一带.导致该地带的强烈构造运动，和频发大、小地震.最后提出了几点认识和建议.     

临汾强震区地壳结构及发震背景研究

以诸城—宜川深地震测深速度剖面为约束, 对沿该剖面得到的高精度重力数据进行拟合, 并对临汾强震区平面布格重力异常进行处理, 得到了该地区地壳密度结构及平面重力异常分布． 利用上述结果分析了临汾强震区的地壳结构及构造环境, 结合前人相关研究成果, 认为临汾强震区地壳中存在塑性相对较强的介质, 洪洞和临汾两次历史地震皆发生在其与周边弹性介质的转换边界上． 另外, 临汾凹陷南北两侧局部构造环境存在差异, 在区域应力场作用下, 导致了洪洞地震和临汾地震的发生. 两次地震在发震时间、 地点和震级等地震要素上有所不同．      

鲜水河断裂带南北构造差异性的地球物理分析

对鲜水河断裂带重磁异常进行向上延拓,通过计算观测面高度异常与延拓后异常之间的相关系数得出最佳向上延拓高度,该延拓高度所对应的延拓结果即为研究区构造背景产生的异常值。向上延拓结果显示布格重力异常值沿鲜水河断裂带自北西向南东逐渐增大,反映出下地壳底边界沿该方向呈升高趋势,可能由青藏高原地壳软弱物质 “东向逃逸” 所致。化极后ΔT磁异常延拓结果表明鲜水河断裂带南东段的康定—石棉以东为强磁性刚性基底。以鲜水河断裂带为分界的不同地块之间基底岩石及地层物性的不同是断裂带南北磁性差异的主导原因。对鲜水河断裂带两侧各约50 km范围内地震的震源深度进行统计,经投影至剖面及线性拟合求出鲜水河断裂带的三维几何形状。结果表明鲜水河断裂带总体倾向南西,倾角近乎直立,范围约为57°—88°。