Bouguer gravity study of middle section of Tanlu fault

Although Tanlu fault is one of the most important tectonic fault zones and active earthquake belts in eastern China, little is known about its deep structure. In this study, we use the existing Bouguer gravity data to study the middle section of the Tanlu fault zone, which is also known as the Yishu fault zone. Our gravity inversion results indicate that the Moho has an abrupt offset in depth at the Tanlu fault zone and it has a relatively smooth variation away from the fault zone. The crustal structures on both sides are different from each other. Sediment is thin on the west side with an average thickness of less than 5 km, while it is as thick as 6 km on the east side. The thinnest sediment （3-4 km） is at the fault zone. Moho depth increases from 33 to 34 km on east side and from 36 to 38 km on west side. Tanlu fault zone is shown as a wide zone of linear gradient in the Bouguer gravity anomaly.     

应用布格重力异常研究郯庐断裂构造

使用布格重力资料对郯庐断裂带的中段部分（沂沭断裂带）进行了研究. 结果表明, 郯庐断裂带莫霍面及地壳内界面均发生错断,断裂带两侧地壳各界面起伏平缓. 该结果与前人的郯庐断裂带是切穿地壳的深大断裂的认识相一致. 在郯庐断裂带两侧,地壳结构明显不同,西侧沉积层较薄,平均在5 km以下;东侧多数在6 km以上;在断裂带中央沉积物最薄,大约为3～4 km. 断裂带东侧莫霍面埋深浅,大约为33～34 km;西侧莫霍面埋深明显增加,达到36～38 km．反映了莫霍面深度在断裂带附近整体是向西增加的. 郯庐断裂带在重力场分布中则表现为一条宽度较大的线性布格重力异常梯度带.      

THE DEPTH OF MOHO INTERFACE AND CRUSTAL THICKNESS IN SICHUAN-YUNNAN REGION,CHINA

By using moving average method to separate Bouguer gravity anomaly field in Sichuan-Yunnan region, we got the low-frequency Bouguer gravity anomaly field which reflects the undulating of Moho interface. The initial model is obtained after seismic model transformation and elevation correction. Then, we used Parker method to invert the low-frequency Bouguer gravity anomaly field to obtain the depth of Moho interface and crustal thickness in the area. The results show that the Qinghai-Tibet block in the northwest of the study area deepens and thickens from the edge to the interior, with the depth of Moho interface and the crust thickness of about 52~62km and 54~66km, respectively. The depth of Moho interface in Sichuan Basin is about 38~42km. In Sichuan-Yunnan block, the depth of Moho interface is about 42~62km from southeast to northwest. Beneath the West Yunnan block, west of the Red River fault zone, the Moho depth is about 34~52km from south to north. The Longmen Mountains and Red River fault zone are the gradient zone of the Moho depth change. Along the Red River fault zone, the depth difference of Moho interface is increasing gradually from north to south. No obvious uplift is found on the Moho interface of Panzhihua rift valley. The depth of Moho interface distribution in Sichuan and Yunnan is obviously restricted by the collision between the Indian plate and the Eurasian plate and the lateral subduction of the Indo-China peninsula. The mean square error of the depth of Moho interface is less than 1.7km between the result of divisional density interface inversion and artificial seismic exploration. At the same time, we compared the integral with divisional inversion result. It shows that:in areas where there is obvious difference between the crust velocity and density structure in different tectonic blocks, the use of high resolution seismic exploration data as the constraints to the divisional density interface inversion can effectively improve the reliability of inversion results.     

基于变密度模型的位场界面反演

在波数域中用重力反演莫霍面深度时通常假定壳幔密度差为一常数,但这只是一种近似的密度模型,本文采用了密度随深度呈指数变化的变密度模型来反演莫霍面深度,给出了利用指数密度模型在波数域中计算重力异常的正演公式及界面深度的反演公式.利用指数密度模型及重力资料反演了青藏高原莫霍面的深度,分析了莫霍面的特征.结果表明,青藏高原莫霍面呈现出边缘浅、中部深的特点,边缘变化快、梯度大,中间变化梯度趋缓.中心地带的羌塘地体莫霍面深度达74 km,向四周慢慢变浅至67 km左右,边缘地区突然变浅至50km左右.通过常密度模型、变密度模型及地震反演得到的莫霍面的比较,证实变密度模型更适合于莫霍面结构的反演.     

南海海盆三维重力约束反演莫霍面深度及其特征

利用南海海盆及周边最新的重力,经过海底地形、沉积层的重力效应改正,并采用岩石圈减薄模型的温度场公式,校正了从张裂边缘到扩张海盆的热扰动重力效应.通过研究区的地震剖面和少量声呐数据得到的莫霍面深度点作为约束,采用基于"起伏界面初始模型"的深度修正量反演迭代公式,反演、计算了研究区的莫霍面深度及地壳厚度.结果表明,海盆区莫霍面深度在8~14 km之间,地壳厚度在3~9 km之间;东部海盆和西南海盆残留扩张中心沿NNE向展布向西南延伸至112°E,莫霍面深度超过12 km,地壳厚度在6 km以上,而西北海盆没有明显的增厚扩张中心;在西南海盆北缘的中沙地块南侧,存在一个近EW向地壳减薄带,地壳厚度在9~10 km;莫霍面深度14 km的等深线和地壳厚度9 km的等值线可指示洋陆边界位置.     

龙门山断裂带地壳密度结构

研究龙门山及邻区地壳密度结构对于认识该地区地震活动性具有重要意义.根据龙门山及邻区( 100°～105°E,28°～33°N)的布格重力异常资料,选取了跨越龙门山断裂带的6条重力测线,在深地震测深资料约束下,使用Geosoft软件分别反演出了龙门山地区地下的沉积层、康拉德界面和莫霍面的深度分布.研究结果表明:龙门山断裂带两侧的地壳结构明显不同,西面高原地区沉积层较薄,大部分为基岩出露;而东边盆地沉积层明显较厚,多在6km以上.莫霍面和康拉德面在两侧均相对平缓,康拉德面从东部的大约24km增加到青藏高原山区的35km左右;莫霍面深度从东部盆地的大约42km增加到西部青藏高原的67km左右.龙门山断裂带整体表现为一条近SN向的陡变重力梯度带,并在其地壳内各界面均发生错断,莫霍面和康拉德面错断距离分别达6 ～ 7km和3～ 5km.该区地壳的这种陡变和不均匀性是导致地震活动性强烈的主要原因之一.     

3D upper crustal density structure of the Deccan Syneclise,Central India

A constrained 3D density model of the upper crust along a part of the Deccan Syneclise is carried out based on the complete Bouguer anomaly data. Spectral analysis of the complete Bouguer gravity anomaly map of the study region suggests two major sources: short wavelength anomalies (<100 km) caused primarily due to the density inhomogeneities at shallow crustal level and long wavelength anomalies (>100 km) produced due to the sources deeper than the upper crust. A residual map of the short wavelength anomalies is prepared from the complete Bouguer anomaly using Butterworth high‐pass filter (100 km cut‐off wavelength). Utilizing the constraints from deep resistivity sounding, magnetotellurics and deep seismic sounding studies, 2.5D density models have been generated along 39 profiles of this region. The mismatch between the calculated response of the a priori 2.5D model with the residual (short wavelength) gravity anomalies is minimized by introducing high‐density intrusive bodies (≥2.81 g/cm³) in the basement. With these 2.5D density models, the initial geometry of our 3D density model, which includes alluvium, Deccan trap, Mesozoic sediment and high‐density intrusive bodies in the basement up to a depth of 7 km (upper crust), is generated. In the final 3D model, Deccan trap extends from 200 m to nearly 1700 m below the 90–150 m thick Quaternary sediment. Further down, the sub‐trappean Mesozoic sediment is present at a depth range of 600–3000 m followed by the basement. The derived 3D density model also indicates six intrusive bodies of density 2.83 g/cm³ in the basement at an average depth of about 4–7 km that best fits the residual gravity anomaly of the study area.     

Crustal Configuration of the Northwest Himalaya Based on Modeling of Gravity Data

The gravity response and crustal shortening in the Himalayan belt are modeled in detail for the first time in the NW Himalaya. The Bouguer gravity anomaly along a ~450-km-long (projected) transect from the Sub-Himalaya in the south to the Karakoram fault in the north across the Indus-Tsangpo Suture Zone is modeled using spectral analysis, wavelet transform and forward modeling. The spectral analysis suggests three-layer interfaces in the lithosphere at 68-, 34- and 11-km depths corresponding to the Moho, the Conrad discontinuity and the Himalayan decollement thrust, respectively. The coherence, admittance and cross spectra suggest crustal shortening because of convergence compensated by lithospheric folding at 536- and 178-km wavelength at the Moho and the upper-crustal level. An average effective elastic thickness of around 31 km is calculated using the coherence method. The gravity data are modeled to demarcate intracrustal to subcrustal regional thrust/fault zones. The geometrical constraints of these faults are obtained in the space scale domain using the wavelet transform, showing good correlation with the major tectonic boundaries. The crustal configuration along the transect shows how the Moho depth increases from 45 to 80 km towards the north with the locus of flexure of the Indian crust beneath the Higher Himalayan zone. The combination of forward modeling and wavelet analysis gives insight into the subsurface extent and geometry of regional structures across the NW Himalaya.     

用远震接收函数反演福建地区宽频带台站下方莫霍界面深度

收集福建省“九五”数字地震遥测台网中8个宽频带台站的远震波形资料,应用接收函数的研究方法计算各个台站下方的接收函数。采用非线性的反演方法获得这些台站下方的S波速度结构．确定这些台站下方莫霍界面深度的分布情况。分析得到的反演结果,福建地区莫霍面的起伏不大．平均的地壳厚度约为32km。在0～2km之间均存在一层低速层,这与地表覆盖着一层松散的沉积层是相对应的。内陆地区台站附近莫霍界面深度较沿海地区略高,沿海台站的莫霍界面深度北部略高于南部。     

基于接收函数约束的川滇地区莫霍面深度反演研究

本文利用川滇地区宽频地震接收函数结果和WGM2012全球布格重力场模型数据,采用正则化参数和接收函数结果交叉验证得到最优莫霍面参考深度和上下界面密度差,使用基于球坐标系下的快速非线性重力反演方法建立川滇地区莫霍面深度模型.研究结果显示,川滇地区整体莫霍面深度介于30～69km,青藏高原内部地区莫霍面深度大于50km;四川盆地莫霍面深度在36～38km;攀枝花地区莫霍面出现明显的隆起和下凹,变化范围在42～48km;川滇菱形地块莫霍面深度在40～50km;滇西和滇南地块莫霍面深度由南向北逐渐变深,变化范围在38～44km.本文反演莫霍面深度与接收函数结果平均误差为0.18km,与该区域天然地震层析成像、人工地震探测以及重力数据反演结果基本一致,但细节更加丰富,进一步确认了莫霍面在攀西裂谷地区存在隆起,小江断裂带下方存在下凹的特征.该结果可作为精细化川滇地区地壳密度界面模型,为研究该地区岩石圈结构和地质构造演化提供参考.     

南海北部地壳密度结构:基于约束三维重力反演

地壳结构的揭示是研究陆缘伸展机制的基础.尽管在南海北部陆缘已开展了大量地壳尺度的二维地震探测,但目前还存在许多覆盖空白的地区,这些调查所得到的结果无法提供地壳结构的区域视图.为了揭示南海北部的地壳结构,本研究将水深和沉积物厚度信息作为约束条件,对空间重力异常进行区域三维反演,并利用地震研究的结果来衡量反演结果的质量.沿地震测线的密度切片显示,重力反演结果与地震研究结果具有良好的一致性.本研究通过选择两个与地震研究结果最匹配的密度异常等值面分别作为康拉德面和莫霍面,获取了莫霍面深度和上、下地壳的厚度.根据假设的相应初始地壳厚度,本研究进一步计算了全地壳、上地壳和下地壳的拉张因子.通过与已发表的居里面深度比较,发现南海北部大部分地区的居里面深度均位于莫霍面之下,其中西沙海槽的拉张因子βw大于3.5并且缺乏岩浆活动,是地幔橄榄岩蛇纹石化的潜在区域.南海北部陆缘的地壳拉张因子显示其经历了伸展方向为128°和160°的两组张裂运动,分别对应于神弧运动和珠琼运动一幕及二幕,在张裂过程中应力场发生了顺时针旋转.此外,上、下地壳的拉张因子表明北部陆缘普遍存在正向和反向差异伸展,陆架区域表现为反向差异伸展,洋陆过渡带为正向差异伸展,推测这种正向和反向差异伸展可能是由下地壳流动导致的,由地壳厚度差异引起的横向梯度力、软流圈浮力和沉积物负载共同驱动.     

The gravity and isostatic Moho undulations in Qinghai–Tibet plateau

It is our interest to study the Moho depths in the Qinghai–Tibet. An iterative hybrid spectral–classical methodology is applied to invert the gravity data and obtain the 3D variation in Moho depth. The gravity inversion is constrained by results from deep seismic sounding and seismological investigations. The Moho is found between 70 and 75 km depth over most of Tibet. Maximum depths of up to 80 km are found along the margins of the plateau, and shallower depths of 65 km correlate with an important suture running along central Tibet (Bangong Nujiang). At Moho level most of Tibet is isostatically compensated at 90–110%, according to the Airy isostatic model. The Qaidam basin in North-Eastern Tibet and the Tarim basin to the North-West are found to be over-compensated.     

银川盆地深地震断层的三维透射成像

为了获得三维地震透射成像技术在活断层探测中的有效性和应用价值的评价,在银川盆地中北部布设了一个三维地震透射台阵,利用该台阵获得的基底初至折射波和莫霍界面反射波资料,采用有限差分反演、时间项反演和连续速度模型反演方法,对台阵区域基底及上地壳结构进行了分析.结果表明：研究区基底呈东西浅、中部深的界面形态,且西陡东缓,最深处大致位于芦花台至西大滩一带,埋深达7 km.芦花台断层、银川—平罗断层、黄河断层在研究区内均表现为北北东走向的速度差异条带,且断层两侧基底及沉积界面埋深存在显著变化.芦花台断层东倾,倾角较陡,延伸至研究区基底之下;银川—平罗断层倾向西,是一条超基底的隐伏断层;黄河断层西倾,延伸深度超过研究区基底.本探测结果证明,初至折射波与莫霍面反射波探测相结合的三维地震透射台阵技术能够给出研究区上地壳三维细结构图像,不仅可以揭示主要断裂的展布位置、浅部空间形态和特征,而且可以揭示断裂向基底之下的延伸状况.     

基于抛物线密度模型的频率域三维界面反演及其在川滇地区的应用

密度界面反演作为了解地球内部结构的一种重要方法,长期以来都是重力学研究的主要内容.本文结合抛物线密度模型及频率域算法的优点,将抛物线密度函数应用于Parker-Oldenburg算法,经过理论推导得到了抛物线密度模型的频率域公式,从而建立了基于抛物线密度模型的三维密度界面重力异常正反演的算法和流程.理论模型数据试验表明本方法快速、有效,适用于大多数浅部比深部增加更快的实际地壳密度.研究中还利用该方法对川滇地区重力异常进行了反演,获得了该区的莫霍面深度分布,并与接收函数研究结果进行对比分析,进一步验证了本文方法的正确性和有效性.     

MAPPING OF CRUSTAL DISCONTINUITIES BY WAVELENGTH FILTERING OF THE GRAVITY FIELD1

A statistical technique, based on the concept of a 1D energy density spectrum of the observed gravity field, has been used to compute ensemble average depths to various horizons containing causative sources of random geometric shape, size, density, etc. The plot of the logarithm of the energy of the observed Bouguer anomaly versus the angular frequency can be approximated, over a certain frequency band, by a linear segment whose slope is related to an average ensemble depth around which a random distribution of numerous anomalous sources exists. Suitable matched filters, based on the computed values of intercepts and slopes of several linear segments approximating the spectrum, have been used to deconvolve the gravity effects associated with the causative sources, occurring around their respective mean depths. The individual deconvolved gravity effects thus separated out have been modelled using the sin x/x method by assuming a fluctuating interface between two formations. The applicability of the present method has been assessed using two observed Bouguer anomaly profiles: one from Ujjain to Mahan, and the other from Jhansi to Mandla where Deep Seismic Sounding (DSS) results are available. The proposed geological crustal models along these two profiles exhibit reasonably good agreement with those obtained from DSS results. A geologically plausible model of the crust in a virgin region has been presented along a Bouguer anomaly profile from Jaipur to Raipur. The following main conclusions have been drawn from the present analysis: (1) The depths to the Moho and Archaean basement interfaces fluctuate between 33.2 and 36.8 km and between 4.6 and 7.0 km respectively. (2) The Narmada-Son Lineament (NSL) does not coincide exactly with the Moho upwarp beneath it. However, this offset is greater in the eastern part of the NSL rather than in the western part. (3) The development of the Satpura horst structure is due to a rise in the Moho interface in a compressional regime. (4) The intrabasement feature (depth from 5 to 12 km) represents a hybrid massif possibly formed due to an admixture of sialic and simatic crust under a tensional regime in the Ujjain-Mahan section.     

3D Crustal and Lithospheric Structures in the Southeastern Mediterranean and Northeastern Egypt

Crustal and lithospheric thicknesses of the southeastern Mediterranean Basin region were determined using 3D Bouguer and elevation data analysis. The model is based on the assumption of local isostatic equilibrium. The calculated regional and residual Bouguer anomaly maps were employed for highlighting both deep and shallow structures. Generally, the regional field in the area under study is considered to be mainly influenced by the density contrast between the crust and upper mantle. Use of the gravity and topographic data with earthquake focal depths has improved both the geometry and the density distribution in the 3-D calculated profiles. The oceanic-continental boundary, the basement relief, Moho depth and lithosphere-asthenosphere boundary maps were estimated. The results point to the occurrence of thick continental crust areas with a thickness of approximately 32 km in northern Egypt. Below the coastal regions, the thickness of crust decreases abruptly (transition zone). An inverse correlation between sediment and crustal thicknesses shows up from the study. Furthermore, our density model reveals the existence of a continental crustal zone below the Eratosthenes Seamount block. Nevertheless, the crustal type beneath the Levantine basin is typically oceanic; this is covered by sedimentary sequences more than 14 km thick. The modeled Moho map shows a depth of 28–30 km below Cyprus and a depth of 26–28 km beneath the south Florence Rise in the northern west. However, the Moho lies at a constant shallow depth of 22–24 km below the Levantine Basin, which indicates thinning of the crust beneath this region. The Moho map reveals also a maximum depth of about 33–35 km beneath both the northern Egypt and northern Sinai, both of which are of the continental crust. The resulting mantle density anomalies suggest important variations of the lithosphere-asthenosphere boundary (LAB) topography, indicating prominent lithospheric mantle thinning beneath south Cyprus (LAB ~90 km depth), followed by thickening beneath the Eratosthenes seamount, Florence Rise, Levantine Basin and reaching to maximum thickness below Cyprian Arc (LAB ~115–120 km depth), and further followed by thinning in the north African margin plate and north Sinai subplate (LAB ~90–95 km depth). According to our density model profiles, we find that almost all earthquakes in the study area occurred along the western and central segments of the Cyprian arc while they almost disappear along the eastern segment. The active subduction zone in the Cyprian Arc is associated with large negative anomalies due to its low velocity upper mantle zone, which might be an indication of a serpentinized mantle. This means that collision between Cyprus and the Eratosthenes Seamount block is marked by seismic activity. Additionally, this block is in the process of dynamically subsiding, breaking-up and being underthrusted beneath Cyprus to the north and thrusted onto the Levantine Basin to the south.     

川西高原重磁异常特征与构造背景分析

川西高原位于青藏高原东缘,是我国大陆地壳构造变形及地震活动最强烈的区域.利用最新重力、航磁资料,通过异常分析和反演计算,研究了该区鲜水河断裂、理塘断裂、金沙江断裂的重磁异常特征、莫霍面特征、居里面特征,分析得出了这些断裂的深部地质结构与构造背景.计算表明:川西高原莫霍面东南浅、西北深,地壳厚度在43~63km之间.居里面特征表现为条带状,深度在17~23km之间.其中,鲜水河断裂带对应莫霍面深度梯度带,居里面为高低起伏圈闭.理塘断裂带北段莫霍面局部隆坳相间,南段莫霍面逐渐抬升,居里面呈现由西向东加深的梯度带.金沙江断裂带,居里面形成局部抬升,深部可能存在高温地热异常源.综合分析认为,川西高原地壳结构主要特点为:增厚的下地壳,热-塑性变形的中地壳,脆性变形的上地壳.     

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.     

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

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

大兴安岭造山带及两侧邻区莫霍面深度与泊松比分布特征

利用2009～2016年内蒙古自治区数字地震台网宽频带固定地震台站的远震波形数据,采用接收函数H-k算法获得23个基岩台站下方的莫霍面深度和泊松比,同时,收集并筛选出277个已有探测台阵和流动台站的接收函数研究结果,综合分析给出大兴安岭造山带及两侧邻区莫霍面深度、泊松比的分布特征。研究表明,研究区域的莫霍面在整体上呈现自东向西逐渐加深的特征,莫霍面深度为25.0～42.3km,平均约为33.5km。莫霍面最浅的区域为松辽盆地(深度为27.0～35.0km),最深的区域为大兴安岭重力梯级带以西地区(深度为41.0～42.3km)。研究区域泊松比为0.19～0.33,平均值为0.26,大于全球大陆地壳的平均值。泊松比高值异常区集中在火山岩区及具有较厚沉积层的盆地。台站所处位置的海拔与莫霍面深度之间具有较强的正相关性,艾里补偿模式在研究区成立,莫霍面起伏与区域地形地貌特征间具有显著的镜像关系。大兴安岭地区的莫霍面深度与泊松比间存在显著的反相关关系,而在松辽盆地及周缘地区未发现明显的规律性,这也意味着松辽盆地在构造演化过程中经历了更为复杂的地壳改造过程。