Variation of crustal thickness in the Philippine Sea deduced from three-dimensional gravity modeling

Abstract   Crustal thickness of the northern to central Philippine Sea was gravimetrically determined on the simple assumption of four layers: seawater, sediments, crust and lithospheric mantle, with densities of 1030, 2300, 2800 and 3300 kg/m³, respectively. As for the correction of the regional gravity variation, a 15 km difference of the lithospheric thickness with a density difference of 50 kg/m³ against the asthenosphere below between both sides of the Kyushu-Palau Ridge was taken into consideration. Mantle Bouguer anomalies were calculated on the assumption of constant crustal thickness of 6 km, and then the crustal thickness was obtained by three-dimensional gravity inversion method. The results show occurrence of thin crust areas with a thickness of approximately 5 km in the southern part and at the western margin of the Shikoku Basin and also of thick crust areas in the northwestern and northeastern parts of the Parece Vela Basin. We suggest that these are because of the variation of magma supply at the time of sea floor spreading in the Shikoku and Parece Vela Basins, which is possibly related to the variation of spreading rate and enhanced magmatism near the past arc volcanic fronts. The results further show the occurrence of crust thinner than 5 km in the northeastern part of the West Philippine Basin, of crust thicker than 15 km in the Amami Plateau, the Daito and Oki-Daito Ridges, and also in the northern part of Kyushu-Palau Ridge, whereas the southern part of the Kyushu-Palau Ridge the crust is thicker than 10 km. It was also inferred that small basins in the Daito Ridge province have the thinnest oceanic crust of less than 5 km in the Kita-Daito Basin.     

Regional gravity field modeling based on rectangular harmonic analysis

Regional gravity field modeling with high-precision and high-resolution is one of the most important scientific objectives in geodesy,and can provide fundamental information for geophysics,geodynamics,seismology,and mineral exploration.Rectangular harmonic analysis(RHA)is proposed for regional gravity field modeling in this paper.By solving the Laplace’s equation of gravitational potential in local Cartesian coordinate system,the rectangular harmonic expansions of disturbing potential,gravity anomaly,gravity disturbance,geoid undulation and deflection of the vertical are derived,and so are the formula for signal degree variance and error degree variance of the rectangular harmonic coefficients(RHC).We also present the mathematical model and detailed algorithm for the solution of RHC using RHA from gravity observations.In order to reduce the edge effects caused by periodic continuation in RHA,we propose the strategy of extending the size of computation domain.The RHA-based modeling method is validated by conducting numerical experiments based on simulated ground and airborne gravity data that are generated from geopotential model EGM2008 and contaminated by Gauss white noise with standard deviation of 2 mGal.The accuracy of the 2.5′×2.5′geoid undulations computed from ground and airborne gravity data is 1 and 1.4cm,respectively.The standard error of the gravity disturbances that downward continued from the flight height of 4 km to the geoid is only 3.1 mGal.Numerical results confirm that RHA is able to provide a reliable and accurate regional gravity field model,which may be a new option for the representation of the fine structure of regional gravity field.     

岩浆与构造作用对洋壳厚度的影响——以西北印度洋为例

洋中脊及邻区洋盆的洋壳厚度能很好地反映区域岩浆补给特征,对于研究洋中脊内部及周缘岩浆活动和构造演化过程具有很好的指示意义.西北印度洋中脊作为典型的慢速扩张洋中脊,其扩张过程与周缘构造活动具有很强的时空关系.本文利用剩余地幔布格重力异常反演了西北印度洋洋壳厚度,由此分析区域内洋壳厚度分布和岩浆补给特征.研究发现,西北印度...     

Constraining the crustal structure under the central and western Tian Shan based on teleseismic receiver functions and gravity anomalies

The Tian Shan is a vast range that spans several countries in Asia. Understanding its evolutionary history may provide valuable insights into intracontinental orogenic dynamics. In this study, we explored the crustal characteristics of the Tian Shan and their relationships to the tectonic evolution of the region. A new H-stacking method that combines the P receiver function and gravity anomalies was used to estimate the thickness and ratio of P- to S-wave velocities (v_P/v_S) for 91 broadband seismic stations in the central and western Tian Shan. Our results revealed significant lateral variations in crustal thickness and v_P/v_S. A ~45-km-thick crust and an intermediate-high v_P/v_S (~1.74–1.84) were found in the Kazakh Shield and Tarim Basin, which we interpreted to indicate a mafic crystalline basement and lower crust. The central Tian Shan varied greatly in crustal thickness (40–64 km) and v_P/v_S ratio (1.65–2.00), which may be due to crustal shortening, mafic underplating, and crustal melting. In contrast, we observed a relatively thin crust (42–50 km) with an intermediate v_P/v_S ratio (~1.78) in the western Tian Shan. The differences in the crustal structures between the western and central Tian Shan imply that the Talas-Fergana Fault may be trans-lithospheric.     

龙门山断裂系及邻区地壳重力均衡效应与汶川地震

龙门山断裂系位于松潘-甘孜山地与成都盆地平原之间,地形高差达3500±500 m.根据重力资料给出四川中西部成都盆地、龙门山、邛崃山区、松潘等川西山地(102°E～105°E,30°N～33°N)的地壳厚度(Moho界面深度)分布.为了探究汶川大地震孕育与发生的原因,将其与据艾里(Airy)大陆均衡理论给出的理论均衡地壳厚度进行对比分析,以探讨该地区的地壳均衡状态.研究结果表明:龙门山断裂带系及以西地区深部地壳是处于很不均衡的状态.为了深化认识重力均衡与大地震的相关性,尚需对该地区深部结构与地壳重力均衡效应作精细研究.同时,本文上述研究结果亦有益于震后确定城镇布局和为不同类型建筑物重建选址方面提供深部要素.     

南北地震带南段地壳厚度重震联合最优化反演

重力反演方法是研究地壳结构和物性界面起伏的有效地球物理手段之一.本文收集了南北地震带南段67个已有的固定台站接收函数反演的Moho面深度结果,并使用基于EGM2008重力异常模型计算的布格重力异常,验证了本文提出的重震联合密度界面反演方法的有效性.利用接收函数对台站下方Moho面深度估计作为先验约束,定义了一类评价函数,通过对重力反演算法中尺度因子,平移因子和稳定性因子的最优选择,最小化重力反演结果与接收函数模型之间的差异.结果表明,本文提出的方法,可以有效地同化不同地球物理方法获得的反演模型,且通过重震联合反演可以改进由于对空间分布不均匀的接收函数结果插值可能而引起的误差.本文还通过引入Crust1.0的Moho面深度为初值,同时考虑地壳密度的横向不均匀分布,通过模型之间的联合反演有效改善了地球物理反演模型间的不一致性问题.本文反演得到的最优化Moho面深度模型与已知67个台站位置接收函数模型之间的标准差约1.9 km,小于Crust1.0与接收函数结果模型之间标准差为3.73 km的统计结果.本文研究结果对于同化重震反演结果、精化地壳密度界面模型,都具有十分重要的参考意义.     

喜马拉雅“东构造结”地区的特异重力场与深部地壳结构

跨越中、印、缅三国交界的喜马拉雅“东构造结”地区（92°E～97°E，26°N～30°N）有一半以上的面积尚没有重力测点，是重力数据空白区,故无法直接研究其重力场特征与深部地壳结构（构造）.本文应用卫星重力异常资料作为近似空间重力异常，经计算给出的布格重力异常，其特征与该地区的地形高程呈很好的镜像相关.据此得到该区不同方位的3个地壳深部结构剖面.重力异常反演求得青藏高原地壳厚度>70 km；喜马拉雅造山带为55 km左右；布拉马普特拉河谷盆地为33～35 km；那加山山脉地区为40～45 km，即呈现出3个不同构造单元的展布.同时求得“东构造结"区由高密度的刚性物质构成，在印度洋板块的碰撞、挤压作用下呈向北运移，并插入青藏高原东缘.基于这样的构造格局和深层动力过程，导致了青藏高原东南缘和东北缘的强烈构造运动，大、小地震的频频发生和矿产资源的聚集.     

青藏高原三维地壳结构的重力正演模拟

青藏高原地壳结构的重力学研究主要集中于莫霍面深度反演及单个剖面的重力正演,但是,在亚东-格尔木与下察隅-共和之间没有地震测线的区域,地壳结构的研究受到限制.为了探讨从青藏高原中部到东部的地壳结构特征,本文采用三维重力正演的方法,对青藏高原三维地壳结构进行了重力正演模拟.研究结果显示,地壳各层界面起伏较大,表明高原地壳变形强烈.研究区域中,藏南喜马拉雅地体界面相对较浅,处于界面陡变带,往北进入拉萨地体,界面相对变深.正演剖面中,大的界面起伏与大的断裂构造有关,莫霍面的陡变也暗示着莫霍面存在错断现象.东构造结莫霍面有抬升现象,在玛多附近莫霍面有下凹现象.地壳分层正演显示,32 km以下地壳产生的重力异常占总异常的80%以上,低速层及高速体产生的重力异常最大可占到总异常的10%左右.     

地球重力场恢复中的位旋转效应

分析了地球自转引起的位旋转效应公式中采用近似速度的影响. 对一组GFZ的快速科学轨道、一组TUM的约化动力法轨道以及一组GFZ的事后科学轨道,计算了星历提供的速度与只有地球引力场对卫星产生作用时的卫星速度的差值,其中参考重力场模型分别采用EGM96、EIGEN2和EIGEN_CG01C. 通过比较得出：轨道数据与EIGEN2地球重力场模型的自恰性优于EGM96和EIGEN_CG01C地球重力场模型. 速度差各分量的变化具有很明显的周期性且与卫星轨道的运行周期相吻合. 当要求在卫星轨迹处获得1m2/s2精度的扰动位时,也即要求位旋转效应公式中卫星速度的近似精度小于2mm/s时,GFZ的快速科学轨道、TUM的约化动力法轨道只需要剔除那些速度精度不满足要求的卫星轨迹点;当要求在卫星轨迹处获得0.5m2/s2精度的扰动位时,应当重新估算上述轨道的速度信息,或采用精度更高的GFZ事后科学轨道.     

用Sp震相测定北京地震台附近的地壳厚度

根据CDSN北京地震台宽频带数字地震仪记录到的20个远震记录中的S震相与Sp震相的到时差测定了北京地震台附近地壳厚度值H。得出了以北京地震台为圆心20～45km环形区域内地壳厚度平均值为40．2km左右的结果。并简要的讨论了Sp震相在北京地震台的记录特征。初步得出了Sp震相在北京地震台的可观测范围为11°～40°,S波震相高效率转换为Sp震相的距离约在30°的结果。     

The gravity field and GGOS

The gravity field of the earth is a natural element of the Global Geodetic Observing System (GGOS). Gravity field quantities are like spatial geodetic observations of potential very high accuracy, with measurements, currently at part-per-billion (ppb) accuracy, but gravity field quantities are also unique as they can be globally represented by harmonic functions (long-wavelength geopotential model primarily from satellite gravity field missions), or based on point sampling (airborne and in situ absolute and superconducting gravimetry). From a GGOS global perspective, one of the main challenges is to ensure the consistency of the global and regional geopotential and geoid models, and the temporal changes of the gravity field at large spatial scales. The International Gravity Field Service, an umbrella “level-2” IAG service (incorporating the International Gravity Bureau, International Geoid Service, International Center for Earth Tides, International Center for Global Earth models, and other future new services for, e.g., digital terrain models), would be a natural key element contributing to GGOS. Major parts of the work of the services would, however, remain complementary to the GGOS contributions, which focus on the long-wavelength components of the geopotential and its temporal variations, the consistent procedures for regional data processing in a unified vertical datum and Terrestrial Reference Frame, and the ensuring validations of long-wavelength gravity field data products.     

Dynamic pattern characteristics of fault deformation and gravity field in the development process of Yongdeng M S =5.8 earth quake

ＤｙｎａｍｉｃｐａｔｅｒｎｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｆａｕｌｔｄｅｆｏｒｍａｔｉｏｎａｎｄｇｒａｖｉｔｙｆｉｅｌｄｉｎｔｈｅｄｅｖｅｌｏｐｍｅｎｔｐｒｏｃｅｓｏｆＹｏｎｇｄｅｎｇＭＳ＝５．８ｅａｒｔｈｑｕａｋｅＺＡＩＳＥＮＪＩＡＮＧ...     

由GOCE引力场模型和CNES-CLS2010平均海面高计算的稳态海面地形

利用欧空局发布的三组GOCE引力场模型及CNES-CLS 2010平均海面高数据,计算得到了全球的稳态海面地形,进而得到了全球地转流速度图.在此基础上重点对黑潮进行了对比分析.结果表明:GOCE不同组解的稳定性较好,所计算的稳态海面地形的差异基本在厘米量级内,这间接表明了GOCE引力场模型提供的大地水准面的精度达到了厘米量级.此外,通过将GOCE与GRACE相应结果进行对比发现,GOCE可提供更多的局部信息,特别是对于流速快、水流窄的边界流,如黑潮、墨西哥湾流等,GOCE所得结果更加清晰,速度也更精确.     

应用卫星重力信息对横断山系地区布格重力异常特异分布的纠正

在云南省西部,跨越中、缅两国交界的横断山系地区（97°E~102°E,24°N~30°N）有近一半的面积尚没有重力测点、即重力数据空白区和重力测点稀少的普查级测区.以前的有关文献、图集中所给出对此地区的重力场都是十分模糊的结果与图件.因此应用这些资料无法详细地研究该地区重力场特征与深部地壳结构（构造）.本文应用卫星重力异常资料作为“近似空间重力异常”,经中间层改正后给出的“计算布格重力异常”,其分布特征与该地区的地形高程呈很好的镜像相关.对相应山脉、河谷以及断裂构造都有所反映.特别是在横断山系地区该布格重力异常呈现为近南北的走向.为此,据该“计算布格重力异常”,并选定对该区有代表性的一条重力异常剖面作正反演计算,以得到其地壳深部结构剖面.结果表明,在横断山脉地区的地壳厚度在51~56 km间起伏变化;滇西北云岭山系以及玉龙山区的地壳厚度约在60 km以上. 最后,对所得结果与图件进行了讨论,并提出了几点认识和纠正的建议.     

基于GIS的区域重磁空间信息集成研究

根据区域重磁场的特点,利用其水平一导峰谷值线可确定异常体的线性构造,垂向二阶导数可确定异常体的环形构造的基本原理,在计算机上自动完成重磁数据处理及成图;并基于GIS"空间分析"功能对图形进行模式识别;实现区域重磁构造的自动提取与空间信息的自动集成.通过对我国东昆仑1:100万重力资料的实际应用,表明其结果具有实用性,其信息集成技术可用于改造以往计算机成图并借助透明纸叠置的人工重磁资料解释工作模式.     

基于GIS的区域重磁空间信息集成研究

根据区域重磁场的特点，利用其水平一导峰谷值线可确定异常体的线性构造，垂向二阶导数可确定异常体的环形构造的基本原理，在计算机上自动完成重磁数据处理及成图；并基于GIS"空间分析"功能对图形进行模式识别；实现区域重磁构造的自动提取与空间信息的自动集成.通过对我国东昆仑1:100万重力资料的实际应用，表明其结果具有实用性，其信息集成技术可用于改造以往计算机成图并借助透明纸叠置的人工重磁资料解释工作模式.     

利用GRACE时变重力场反演青藏高原的隆升速率

青藏高原隆升对中国、亚洲乃至世界的气候都有着重要影响,研究青藏高原地壳隆升速率具有重大意义.本文利用2004—2015年期间高覆盖度的卫星重力数据,通过去除陆地储水的重力效应获得地壳隆升引起的重力变化速率,基于直立长方体垂直运动与重力变化的关系模型反演了该区域的地壳隆升速率分布.研究结果表明在300 km的空间尺度下青藏高原隆升速率分布具有不均匀的特点,表现为以冈底斯山—唐古拉山—鲜水河断裂带为界线,其两侧的速率差异较大.位于界线以南,沿喜马拉雅推覆构造带的区域平均隆升速率为2.01±0.87 mm·a^-1,其中西侧的印度板块与东侧的缅甸板块隆升速率分别为~2.43 mm·a^-1、~2.89 mm·a^-1,位于两板块之间的区域隆升速率为~0.69 mm·a^-1;位于界线以北,除了天山区域和华北板块的隆升速率为~1 mm·a^-1,其他区域隆升现象不明显,其速率为~0 mm·a^-1.我们发现存在两条均穿过正断裂带区域的隆升速率梯度带,其中一条为从加德满都到塔里木盆地,其恰好穿过青藏高原内部的正断裂带,另一条为从那加山到四川盆地,其恰好穿过大理正断裂带.本文反演的青藏高原隆升速率与以往观测到的GPS结果有很好的一致性,为青藏高原隆升、地壳增厚等科学问题提供理论支持.

     

Dynamic pattern characteristics of fault deformation and gravity field in the development process of Yongdeng M S=5.8 earthquake

In this paper, the fault deformation abnormality, dynamic evolution features of gravity and vertical deformation field in the seismogenic process of the Yongdeng, Gansu Province earthquake on July 22, 1995 are studied primarily. There appeared α-β-γ tri-stage anomaly at three sites near the epicenter, and there appeared anomalies of step and sudden jump at more than 10 sites in outer region since 1993. The high value area before shock, coseismic effect and process of recovery aftershock were monitored by portable gravity data. Data reflects the changing process of fault movement from the quasi-linear to the nonlinear in the near source region during seismogenic development of the Yongdeng earthquake and evolution of gravity field from heterogeneity of seismogenic term to quasi-homogeneity of postseismic term. There exists close relationship between strong earthquake and dynamic evolution of regional stress-strain field. Considering all above, the experience and lessons in this medium-short term prediction test are summarized.     

冲绳海槽及其邻域地球物理场与地壳结构特征

分析了冲绳海槽及其邻域的重、磁场特征，同时对研究区域的布格重力数据进行了上延拓及二次垂直导数计算．根据分析与计算可知，在冲绳海槽及其周边区域，重磁异常的主要走向为NNE向，重、磁场的分布特征表明冲绳海槽及邻域地球物理场变化复杂，变化幅度大．参考地震解释结果，利用重力资料，应用调和级数法对研究区的莫霍界面进行了反演计算，得到了该地区莫霍界面深度分布，计算表明冲绳海槽及邻域地区莫霍界面深度在28-4km之间变化，地壳厚度自西向东逐渐减薄，显示出陆壳-过渡壳-洋壳的转变过程．     

利用GRACE时变重力场反演青藏高原的隆升速率

青藏高原隆升对中国、亚洲乃至世界的气候都有着重要影响,研究青藏高原地壳隆升速率具有重大意义.本文利用2004—2015年期间高覆盖度的卫星重力数据,通过去除陆地储水的重力效应获得地壳隆升引起的重力变化速率,基于直立长方体垂直运动与重力变化的关系模型反演了该区域的地壳隆升速率分布.研究结果表明在300 km的空间尺度下青藏高原隆升速率分布具有不均匀的特点,表现为以冈底斯山—唐古拉山—鲜水河断裂带为界线,其两侧的速率差异较大.位于界线以南,沿喜马拉雅推覆构造带的区域平均隆升速率为2.01±0.87 mm·a^-1,其中西侧的印度板块与东侧的缅甸板块隆升速率分别为~2.43 mm·a^-1、~2.89 mm·a^-1,位于两板块之间的区域隆升速率为~0.69 mm·a^-1;位于界线以北,除了天山区域和华北板块的隆升速率为~1 mm·a^-1,其他区域隆升现象不明显,其速率为~0 mm·a^-1.我们发现存在两条均穿过正断裂带区域的隆升速率梯度带,其中一条为从加德满都到塔里木盆地,其恰好穿过青藏高原内部的正断裂带,另一条为从那加山到四川盆地,其恰好穿过大理正断裂带.本文反演的青藏高原隆升速率与以往观测到的GPS结果有很好的一致性,为青藏高原隆升、地壳增厚等科学问题提供理论支持.