地形—均衡补偿重力、大地水准面异常频谱分析

将地形高、地球内部质量异常以及重力、大地水准面展开成球谐级数，依据岩石圈弹性挠曲均衡补偿理论建立地形—均衡补偿重力、大地水准面异常的球谐级数表达式.由此我们可以研究地形—均衡补偿重力、大地水准面异常与球谐级数阶次的关系，以及不同波长地形荷载与岩石圈挠曲补偿的关系，探讨地形—均衡补偿重力、大地水准面的频谱特性.从观测大地水准面异常和自由空气重力异常扣除地形—均衡补偿大地水准面、重力异常，可以得到均衡大地水准面异常和均衡重力异常.均衡大地水准面异常已经消除了浅层物质不均匀的影响，反映的是地球深部物质密度不均匀分布.均衡重力异常显示出中短波长特性，反映的是地壳上地幔物质分布的失衡和物质调整的动力学特征.     

黄海及邻区莫霍面起伏特征

针对地壳构造形成的动力学机制,在广泛收集区域地质、地质构造和地球物理等资料的基础上,着重对重力数据进行了数字化、坐标、重力公式、投影方式和比例尺的统一化处理,进而进行了网格化处理.为将区域异常和局部异常分离开来,以便以对深部构造的系统研究,笔者选取趋势分析法分别对研究区内预处理后的布格重力异常数据(Δg)进行了三维多项式迭代拟合计算,得到了区域布格重力异常二阶、五阶、十阶趋势分析结果.继而,利用二级近似公式迭代法对布格重力异常五阶趋势分析区域异常数据进行了三维运算处理,得到了黄海及周边地区的莫霍面埋深值.分析了黄海及邻区莫霍面起伏特征,并进行了深部构造区划,探讨了深部断裂构造与莫霍面起伏间的成因联系,为深部构造和活动断裂演化的地球动力学研究提供了依据.     

鄂尔多斯地块及邻区重力均衡研究

为更加精细地了解鄂尔多斯及其邻区的地壳均衡状态,本文在考虑地壳密度横向变化的基础上计算了该区域的重力均衡异常。首先,利用研究区内重力异常和真实莫霍面深度进行反演得到壳幔密度差分布,以艾里均衡理论为基础计算了莫霍面的理论深度;其次,通过莫霍深度理论值与真实值对比来获得研究区的重力均衡异常;最后,结合区内地震分布,讨论了重力均衡异常与地震活动的关系。结果表明:研究区内盆地以负均衡异常为主,山地则表现为正均衡异常,鄂尔多斯地块内部基本均衡;研究区内中强地震主要分布在重力均衡异常梯度带上。      

格莱尼改正及其在地球物理数据处理中的意义

格莱尼改正是重力改正的一种,是对球面半径约166.7 km以外地形和均衡造成的重力综合效应的改正,改正后所得格莱尼异常可对应研究区范围内地壳结构特征,具有重要的研究意义.但近年来国内的重力研究中,对格莱尼改正的认识和重要性还不是十分透彻,特别是在大区域研究中,忽略这一远区重力效应可能会对研究结果的解释造成误导.针对这一问题,本文介绍了格莱尼改正的含义、由来与研究历史,指出了在大区域研究中需要考虑格莱尼效应的重要性,以“扇形球壳块”法为例给出了格莱尼改正的计算原理,列举了格莱尼异常在地球物理中的常见应用,并通过重力改正与莫霍面深度反演实例,通过对反演结果的控制点检测与精度评价,指出格莱尼异常反演所得莫霍面埋深更接近地震方法结果,从而表明了格莱尼改正在地球物理研究中的不可忽视性,并指出完全布格改正与格莱尼改正在计算方法上可实现统一.     

南海及邻区莫霍面深度分布特征及其与各构造单元的关系

为了研究南海及邻区莫霍面分布特征及其与边缘海盆、海沟、岛弧、新生代沉积盆地的关系等构造单元的关系,本文通过对研究区的空间重力异常数据进行全布格改正,得到研究区内的布格重力异常,并以近年来的声纳浮标探测与海底地震仪探测剖面所得到的莫霍面深度资料为控制点采用三维带控制点界面反演方法得到了研究区的莫霍面深度图和地壳厚度图.本...     

论华北地区的均衡状态(一)——方法和局部补偿

根据经典均衡的原理,本文重点分析了五种局部均衡补偿模式.在计算方法上利用频域三维重力场的理论公式,快速求得均衡校正值,据此计算了66个模型的均衡重力异常场,注意到对这许多不同参量的模型而言,它们的均衡异常场在形态与分布特征上基本一致,其中计入了地表和上下地壳密度差异分布的 Airy 模型具有最佳的补偿效果,它的均衡面在莫霍界面以下的上地幔中,标准深度50km.从整体上看,华北地台处于亚均衡状态,均衡异常的均值为1810-5m/s2.均衡重力异常的分布表现出明显的块体特征,正均衡异常区主要分布在东部胶辽地块和冀中平原北缘,在汾渭裂谷区存在负异常.模型对比表明,以莫霍界面作为均衡补偿面的模型是不可取的;Airy 模型比 Pratt 模型的补偿效果略好,这同地壳构造以层状为主而侧向变化有限的特征相符.有关复合补偿、均衡重力异常的基本特征和深部构造的关系等结果,将在文章的第二部分发表.      

南海深部构造特征及其地质意义:来自重磁位场反演的认识

莫霍面和居里面是认识深部过程重要的地质与地球物理界面.为了进一步理解南海深部构造活动与洋盆扩张的关系,本文以OBS剖面和深反射地震剖面作为约束,对卫星测高重力异常进行海水、沉积层影响校正,采用最小曲率位场分离方法消除局部密度体的重力影响,获取了反映莫霍面起伏的重力异常,并利用双界面模型重力场快速反演方法计算得到了南海地区莫霍面深度值.通过与居里面起伏的对比研究,发现南海莫霍面和居里面整体均表现为"洋盆浅、周缘深"的菱形特征,两者在洋陆转换区呈现明显的窄梯级带特征,反映了南海扩张期岩石圈的强烈伸展减薄、南北向构造拉张作用等深部构造过程.洋盆莫霍面和居里面的西南向楔形形态是对南海由东向西渐进式扩张的深部构造响应.洋盆南部莫霍面浅于北部,这与扩张中心逐渐向南迁移的特征一致,而洋盆居里面南深北浅的特征则可能与洋盆的简单剪切扩张方式以及洋盆北部的岩浆活动更活跃有关.南海地区莫霍面和居里面呈现交错叠置关系,南、北陆缘表现为明显的深部构造差异,说明南海为非对称式扩张.北部陆缘区居里面深度浅于莫霍面,而洋盆区和南部陆缘区居里面深于莫霍面,这与南、北陆缘性质的差异和南部陆缘复杂的中-新生代俯冲碰撞等构造演化相关,而洋盆区居里面深于莫霍面的现象推测与大洋上地幔橄榄岩蛇纹石化导致的岩石磁性增强有关.     

基于重力水准复测的局部大地水准面变化与强震关系研究

大地水准面是静止海洋面(平均海水面)以及假设这个海洋面在大陆内部延伸形成的一个封闭曲面,它既是地球形状的反映,又是地球内部物质分布和运动的体现.大地水准面被认为是地球重力场的几何表象和重力等位面,也是地球上部圈层的一个物理介面.     

中国及邻区岩石层大地水准面正演模拟

依据地震层析成像、地震测深和大地电磁测深结果的地球内部结构和物性分布,模拟了中国及邻区岩石圈中短波长大地水准面起伏,正演模拟结果显示：地形起伏引起的大地水准成异常巨大,幅值变化可达２１０ｍ,最大值在青藏高原,达１５０余ｍ,沉积层底面起伏产生的大地水准面异常非常小,幅值只有６．５ｍ,表现出短波长特征,岩石层底面志伏生产的大地水准面异常极为平缓,约为２０ｍ。莫霍面产生的大地水准面异常形态与地形起伏引起的大地水准面异常形态相异,异常幅值只有后者的２／３,在青藏高原只有后者一半,岩石层内密度不均匀引起的大地水准面异常的幅值约为７５ｍ。     

青藏高原重力场与壳幔结构分析

本文基于重力资料, 分析了青藏高原壳幔结构模型、 高原陆内形变动力学条件、 高原深部物质运动特征及动力学机制。 研究表明, 重力布格异常和自由空间异常除了分别反映大地水准面之下的“剩余”密度信息和大地水准面之上的“附加”密度信息之外, 还可以组合在一起反映壳幔结构的流变学信息。 在整体处于Airy 重力均衡状态下, 如果局部布格异常与空间异常同向减小, 则是弱地壳强地幔的反映; 如果布格异常减小空间异常增大, 则是强地壳弱地幔的反映。 笔者认为, 青藏高原南部多为强地壳弱地幔地段, 东部既有强地壳弱地幔地段, 也有下地壳柔性-上地幔脆性地段, 北部多为弱地壳强地幔地段。 高原南北两侧板块边界的挤压力对高原做功, 重力位能使高原物质向低位势转移, 产生流变变形, 导致南区和北区主要为挤压变形区, 东区主要为构造伸展-侧向挤出区。 由于壳幔结构的差异, 不同地区驱动变形所需位能大小不同。 相同位能条件下, 南部更易于隆升, 东部更易于流变伸展。     

Gravimetric geoid for Egypt implementing Moho depths and optimal geoid fitting approach

The aim of this paper is to investigate the effect of implementing the experimentally determined GEMMA Moho depths (GOCE Exploitation for Moho Modeling and Applications), which are partly seismically estimated, in gravimetric geoid computation in Egypt. The window remove-restore technique has been proposed to avoid the double consideration of the topographic-isostatic masses in the neighbourhood of the computational point. The plate loading theory has been used to model the seismically determined Moho depths. A constant density contrast between the lower crust and the upper mantle has been postulated. The tailored geopotential model EGTGM2014 has been used for the long wavelength contributions of the Earth’s gravity field. A comparison with a geoid computed using the EGM2008 and Airy floating hypothesis has been made. For all cases, a gravimetric geoid for Egypt has been computed using Stokes’ integral in the frequency domain by 1-D FFT technique. The computed geoids are fitted to the GPS-levelling derived geoid using an optimum geoid fitting technique for Egypt introduced by the author. The results show that using the seismically determined Moho depths within the plate loading theory and the EGTGM2014 tailored geopotential model gives a geoid with external accuracy of about 16 cm.     

Transformation between gravimetric and GPS/levelling-derived geoids using additional gravity information

The transformation from the gravimetric to the GPS/levelling-derived geoid using additional gravity information for the covariance function of geoid height differences has been investigated in a test area in south-western Canada. A “corrector surface” model, which accounts for datum inconsistencies, long-wavelength geoid errors, vertical network distortions and GPS errors, has been constructed using least-squares collocation. The local covariance function of geoid height differences is usually obtained from residual values between the GPS/levelling and gravimetric geoid heights after the elimination of all known systematic distortions. If additional gravity data (in the form of gravity anomalies) are available, the covariance function of geoid height differences can be determined by the following steps: (1) transforming the GPS/levelling-derived geoid heights into gravity anomalies; (2) forming differences between the computed in step 1 and given gravity anomalies; (3) determining the parameters of the local covariance function of the gravity anomaly differences; (4) constructing an analytical covariance model for the geoid height differences from the covariance function of the gravity anomaly differences using the parameters derived in step 3. The advantage of the proposed method stems from the great number of gravity data used to derive the empirical covariance function. A comparison with the least-squares adjustment shows that the standard deviation of the residuals of the predicted geoid height differences with respect to the control point values decreases by 2.4 cm.     

Topography of the Moho and Conrad discontinuities in the Kyushu district,Southwest Japan

The first P-arrival time data from local earthquakes are inverted for two-dimensional variation of the depths to the Conrad and Moho discontinuities in the Kyushu district, southwest Japan. At the same time, earthquake hypocenters and station corrections are determined from the data. The depths to the discontinuities are estimated by minimizing the travel time residuals of first P-arrival phases for 608 earthquakes observed at 57 seismic stations. In the land area of Kyushu, the Conrad and Moho discontinuities are located within the depth ranges of 16–18 and 34–40 km, respectively. The Conrad discontinuity is not as largely undulated as the Moho discontinuity. The depth to the Moho is deep along the east coast of Kyushu, and the deepest Moho is closely related to markedly low velocity of P wave. We regard the deepest Moho as reflecting the Kyushu–Palau ridge subducting beneath the Kyushu district, together with the Philippine Sea slab. In western Kyushu, the shallow Moho is spreading in the north–northeast–south–southwest direction in the Okinawa trough region. Based on the presence of low-velocity anomaly in three-dimensional velocity structure and seismogenic stress field of shallow crustal earthquakes, the shallow Moho is interpreted as being due to lower crustal erosion associated with a small-scale mantle upwelling in the Okinawa trough region. The velocity discontinuity undulation basically has insignificant effect on hypocenter determination of the local earthquakes, but the Moho topography makes changes in focal depths of some upper mantle earthquakes. The depth variation of the Moho discontinuity has a good correlation with the Bouguer gravity anomaly map; i.e., the shallow Moho of western Kyushu and the deep Moho of eastern Kyushu closely correlate with the positive and negative Bouguer gravity anomalies, respectively.     

What compensates the topography of southern Norway? Insights from thermo-isostatic modeling

The origin of the high topography of the Norwegian Mountains is currently much debated. Several geophysical studies show that the uppermost mantle below southern Norway has anomalously low velocities as compared to other parts of the Baltic Shield. This study aims to shed lights on the structure of the lithospheric mantle below southern Norway by adapting and further refining a method based on isostatic and thermal equilibrium to compute temperature, temperature-related density and synthetic S-wave velocity in stable continental domains. The one-dimensional steady-state heat equation is used with topographic, Moho depth, crustal density and surface heat flow data. A condition of local isostasy is assumed and geoid undulations are used to constrain the range of possible lithosphere models.Results derived from this method suggest a thickening of the thermal lithosphere below southern Norway from west to east. The western part is found to have higher temperatures, lower densities and lower synthetic S-wave velocities than the eastern part, compatible with results from a recent P-wave travel time residual study. Comparison of the synthetic shear-velocity profiles beneath southwestern Norway with velocity profiles inverted from Rayleigh wave dispersion data suggests that the higher temperatures associated with a thinner lithosphere can explain parts of the seismic low-velocity anomaly.The inferred lithospheric structure is sensitive to uncertainties in the crustal input model, but the main features remain undisturbed by changes in the input data. The results show that the lithosphere of southwestern Norway can be in local isostatic equilibrium, if it is thinner and warmer than the lithosphere of eastern Norway. The present-day high topography may therefore be partially sustained by lower densities in the mantle lithosphere.     

Spatial and Spectral Analysis of Refined Gravity Data for Modelling the Crust–Mantle Interface and Mantle-Lithosphere Structure

We analyse spatial and spectral characteristics of various refined gravity data used for modelling and gravimetric interpretation of the crust–mantle interface and the mantle-lithosphere structure. Depending on the purpose of the study, refined gravity data have either a strong or weak correlation with the Moho depths (Moho geometry). The compilation of the refined gravity data is purely based on available information on the crustal density structure obtained from seismic surveys without adopting any isostatic hypothesis. We demonstrate that the crust-stripped relative-to-mantle gravity data have a weak correlation with the CRUST2.0 Moho depths of about 0.02. Since gravitational signals due to the crustal density structure and the Moho geometry are subtracted from gravity field, these refined gravity data comprise mainly the information on the mantle lithosphere and sub-lithospheric mantle. On the other hand, the consolidated crust-stripped gravity data, obtained from the gravity field after applying the crust density contrast stripping corrections, comprise mainly the gravitational signal of the Moho geometry, although they also contain the gravitational signal due to anomalous mass density structures within the mantle. In the absence of global models of the mantle structure, the best possible option of computing refined gravity data, suitable for the recovery/refinement of the Moho interface, is to subtract the complete crust-corrected gravity data from the consolidated crust-stripped gravity data. These refined gravity data, that is, the homogenous crust gravity data, have a strong absolute correlation of about 0.99 with the CRUST2.0 Moho depths due to removing a gravitational signal of inhomogeneous density structures within the crust and mantle. Results of the spectral signal decomposition and the subsequent correlation analysis reveal that the correlation of the homogenous crust gravity data with the Moho depths is larger than 0.9 over the investigated harmonic spectrum up to harmonic degree 90. The crust-stripped relative-to-mantle gravity data correlate substantially with the Moho depths above harmonic degree 50 where the correlation exceeds 0.5.     

蒙古及周边地区重力异常和地壳不均匀体分布

基于全球EGM2008自由空气重力异常模型,本文计算了蒙古及周边地区的布格重力异常和AiryHeiskanen均衡重力异常.在此基础上,本文采用Crust 1.0地壳模型为参考,通过重力正演方法,对蒙古及周边地区不同深度地壳密度结构模型的重力异常进行了计算,并对得到的正演布格重力异常与实际重力异常进行了对比和分析.研究结果表明:蒙古西部杭爱山地区与阿尔泰山地区的构造变形差异性明显,现今均衡重力异常中杭爱山周边没有明显的均衡异常高值区,而阿尔泰山地区西南方向存在均衡重力异常高值分布,分析与新构造运动密切相关;Crust 1.0模型给出的壳幔横向密度不均匀体分布对于计算Moho面起伏引起的重力异常作用明显;Crust 1.0给出的地壳内界面变形可以反映深大活动断裂的深部构造变形.研究结果对于认识蒙古东西部构造特征差异,以及现今西部活动断裂的地球物理场特征具有参考意义,也可以为进一步应用Crust 1.0模型为参考开展三维密度结构反演提供一定帮助.     

Expressions for the Global Gravimetric Moho Modeling in Spectral Domain

We apply a newly developed numerical method to improve the Moho geometry by the implementation of gravity data. This method utilizes expressions for the gravimetric forward and inverse modeling derived in a frequency domain. Methods for a spectral analysis and synthesis of the gravity field and crust density structures are applied in the gravimetric forward modeling of the consolidated crust-stripped gravity disturbances, which have a maximum correlation with the (a priori) Moho model. These gravity disturbances are obtained from the Earth’s gravity disturbances after applying the topographic and stripping gravity corrections of major known anomalous crust density structures; in the absence of a global mantle model, mantle density heterogeneities are disregarded. The isostatic scheme applied is based on a complete compensation of the crust relative to the upper mantle density. The functional relation is established between the (unknown) Moho depths and the complete crust-stripped isostatic gravity disturbances, which according to the adopted isostatic scheme have (theoretically) a minimum correlation with the Moho geometry. The system of observation equations, which describes the relation between spherical functions of the isostatic gravity field and the Moho geometry, is defined by means of a linearized Fredholm integral equation of the first kind. The Moho depths are determined based on solving the gravimetric inverse problem. The regularization is applied to stabilize the ill-posed solution. This numerical procedure is utilized to determine the Moho depths globally. The gravimetric result is presented and compared with the seismic Moho model. Our gravimetric result has a relatively good agreement with the CRUST2.0 Moho model by means of the RMS of differences (of 3.5 km). However, the gravimetric solution has a systematic bias. We explain this bias between the gravimetric and seismic Moho models by the unmodelled mantle heterogeneities and uncertainties in the CRUST2.0 global crustal model.     

东北地区重力均衡异常特征的初步研究

本文利用地面实测重力资料和地形高程资料,采用普拉特-海福特(Pratt-ttayford)重力均衡理论模型,取1°×1°方格网,通过使用现成改正表格查取改正值与个别计算点用理论公式计算作校核的方法,计算了我国东北地区75个计算点的均衡重力异常值;并对局部第四系覆盖较厚地区作了第四系密度改正;在此基础上,构制了我国东北N39°—49°,E121°—131°大部分地区的均衡重力异常图;结合区域布格重力异常和区域空间重力异常特征以及莫霍界面的起伏特点作了对比分析和讨论     

欧亚地区均衡残差大地水准面和上地幔强度

首先计算了欧亚地区均衡残差大地水准面．基于地幔热对流的内负荷理论和最新全球层析成像结果，探讨了欧亚地区中波长均衡残差大地水准面的地球动力学意义．研究结果表明，中波长均衡残差大地水准面主要受上地幔粘滞度和岩石层强度的影响，进而得出欧亚地区一些古老地盾和构造稳定地区的上地幔与年轻山脉及构造活动地区的上地幔结构存在着差异．这个差异主要是占老地盾和构造稳定地区，如波罗的海地盾、中西伯利亚地台、东欧等区域，冷却的上地幔已穿透地幔较深，上地幔与岩石层之间耦合较好;而年轻山脉和构造活动区，如帕米尔、天山、贝加尔活动带、青藏高原、日本海周围地区，在上地幔可能存在着热物质即粘滞度很低的软流层，上地幔与岩石层耦合程度较差，甚至有可能解耦．从欧亚地区上地幔属性的差异，可以解释该地区的一些地球动力学问题．     

A regional strategy for geothermal exploration with emphasis on gravity and magnetotellurics

As part of the resource evaluation and exploration program conducted by Los Alamos Scientific Laboratory for the national Hot Dry Rock (HDR) Geothermal Program, a regional magnetotelluric (MT) survey of New Mexico and Arizona is being performed. The MT lines are being located in areas where the results of analysis of residual gravity anomaly maps of Arizona and New Mexico, integrated with other geologic and geophysical studies indicate the greatest potential for HDR resources.The residual gravity anomalies are derived by applying the concept of predicting gravity anomalies from topography. This can be accomplished by employing reductions similar to those used in some isostatic investigations, in which a regional topographic surface is used as the Bouguer reduction datum. The datum is derived by comparison of various harmonics of Bouguer anomalies and elevations of stations. Topography can be used to predict Bouguer anomalies because of isostatic compensation; the resultant anomalies can be considered high frequency residual anomalies or isostatic anomalies corrected for regional compensation. Such maps have been produced for Arizona, New Mexico, west Texas, and Chihuahua, Mexico.The main objective of the MT project is to produce a regional geoelectric contour map of the pervasive deep electrical conductor within the crust and/or upper mantle beneath the Colorado Plateau and the adjacent Basin and Range Province and Rio Grande Rift. The MT survey consists of 200 sites along several long profiles with site spacing of 15–20 km. Pre-existing available MT data are being integrated with the new data. After the data are processed, a one-dimensional inversion is applied to the sounding curve and used as a starting point for 2-D modeling. Such a project and ultimate map will be of major value in studying the regional geophysics and tectonics of the southwest United States as they now apply to HDR resources in particular and geothermal resources in general.Electrical conductivity anomalies of large areal extent are of particular interest in geothermal exploration. Correlation analysis of large conductive anomalies with other geophysical, geological, and geotectonic data is being performed. Preliminary analysis of the data has suggested several major regions of anomalously shallow high electrical conductivity. Among these is the Aquarius area of northwest Arizona which is the site of a longwavelength residual anomaly low, which when modeled and correlated with other geophysical data can be shown to be possibly related to low density and high temperature in the crust at depths of 20 km or less. Preliminary analysis of MT data indicates the possible existence of a mid-crustal high electrical conductivity anomaly in this same region.