SEASAT geoid anomalies and the Macquarie Ridge complex

The seismically active Macquarie Ridge complex forms the Pacific-India plate boundary between New Zealand and the Pacific-Antarctic spreading center. The Late Cenozoic deformation of New Zealand and focal mechanisms of recent large earthquakes in the Macquarie Ridge complex appear consistent with the current plate tectonic models. These models predict a combination of strike-slip and convergent motion in the northern Macquarie Ridge, and strike-slip motion in the southern part. The Hjort trench is the southernmost expression of the Macquarie Ridge complex. Regional considerations of the magnetic lineations imply that some oceanic crust may have been consumed at the Hjort trench. Although this arcuate trench seems inconsistent with the predicted strike-slip setting, a deep trough also occurs in the Romanche fracture zone.Geoid anomalies observed over spreading ridges, subduction zones, and fracture zones are different. Therefore, geoid anomalies may be diagnostic of plate boundary type. We use SEASAT data to examine the Macquarie Ridge complex and find that the geoid anomalies for the northern Hjort trench region are different from the geoid anomalies for the Romanche trough. The Hjort trench region is characterized by an oblique subduction zone geoid anomaly, e.g., the Aleutian-Komandorski region. Also, limited first-motion data for the large 1924 earthquake that occurred in the northern Hjort trench suggest a thrust focal mechanism. We conclude that subduction is occurring at the Hjort trench. The existence of active subduction in this area implies that young oceanic lithosphere can subduct beneath older oceanic lithosphere.     

Long-wavelength gravity anomalies and intraplate seismicity

Anomalously high mid-plate and Mid-Oceanic Ridge seismicity apparently characterizes the flanks of the regional free-air gravity and topographic high roughly centered on the postulated Iceland hot spot; anomalously high mid-plate seismicity also occurs around the deep free-air gravity low centered just south of India. Although the correlation between gravity and seismicity is still somewhat inconclusive, it seems to suggest active asthenosphere flow associated with the regional gravity anomalies. Added to other stresses in the plate, the stress transmitted to the plate by this flow would encourage its fracture.     

Axial depth anomalies from 10 to 50° north along the Mid-Atlantic Ridge: correlation with other mantle properties

We present an analysis of intermediate- to long-wavelength (a few hundred to a few thousand kilometers) axial depth anomalies along the Mid-Atlantic Ridge between 10 and 50°N. The maximum depth of the rift valley is chosen as the elevation datum for oceanic crust of zero-age. The large depth anomalies are correlated, for short wavelengths (less than a few hundred kilometers), with some fracture zones irrespective of their offset and, for intermediate to long wavelengths, with mantle properties expressed in (1) excess elevation near triple junctions, (2) trace element and isotope geochemistry of the basaltic liquids emplaced at the ridge axis, and (3) anomalies in the Earth's gravity field. We suggest that the correlations may be explained in two ways: the depth anomalies of intermediate to long wavelength may represent the sites of upwelling and downwelling of the deep mantle; alternatively, the depth anomalies could be due to a regionalization in temperature or mineralogy in the asthenosphere.     

Regional compensation of subducted lithosphere: effects on geoid,gravity and topography from a preliminary model

It is proposed that the anomalous mass of subducted lithospheric slabs is partially compensated by two mechanisms: (a) by stresses transmitted up the slab to the surface lithosphere, causing the trench, and (b) by stresses transmitted through the mantle wedge overlying the slab, which cause a broad shallow depression. Results calculated from a preliminary two-dimensional model show that it can give a first-order account of observed topography, gravity and geoid anomalies. In particular, the model can reconcile the large mass anomaly predicted for the slab by thermal conduction theory with the smaller than expected gravity and geoid anomalies. Preliminary results indicate substantial compensation by mechanism (a), implying rather large stresses at the top of the slab: these could be important for theories of earthquake mechanisms and plate driving mechanisms. Roughly 50% of the possible compensation by mechanism (b) is indicated: the resulting depression, the order of 500 m deep and over 500 km wide, could explain anomalously deep marginal basins and some marine transgressions of continents.     

Isostasy of aseismic tectonic units in the South Atlantic Ocean and geodynamic implications

Mechanical isostasy of the lithosphere in the South Atlantic Ocean was studied using information on gravity anomalies and bathymetry with additional constraints imposed by the altimetric geoid. The isostatic responses (admittances) over the Walvis Ridge, Rio Grande Rise and Trindade Seamount Chain were computed using a three-dimensional algorithm. The eastern Walvis Ridge and the Rio Grande Rise have the same response, which is well explained by an Airy model of isostasy. The other features are regionally supported. A variation in the thickness of the elastic plate was found along the western Walvis Ridge. A high value of the elastic plate thickness (20 km) was found under the Trindade Chain. Geodynamic implications are discussed in the light of these results.     

The modern geoid and ancient plate boundaries

The non-hydrostatic geoid is dominated by three large anomalies: an area of high gravity potential in the equatorial Pacific; another stretching from Greenland through Africa to the southwest Indian Ocean; and a semi-continuous low region passing from Hudson's Bay through Siberia to India and on to Antarctica. None of these three high-amplitude (greater than 60 m) and long-wavelength anomalies corresponds to present-day plate boundaries. However, if the modern geoid is plotted over the positions of continents and plate boundaries at 125 Ma B.P. (reconstructed relative to hotspots) a strong correlation emerges. The modern geoidal low corresponds in position to the areas of subduction surrounding the Pacific 125 Ma ago. The geoidal high now centered on Africa is entirely contained within ancient Pangaea, and the equatorial Pacific high overlies the location of the spreading centers preserved in the magnetic anomalies of the central Pacific. The most plausible cause of the large geoidal undulations is lower mantle convection only weakly coupled to plate motions. The correspondence between modern geoid and ancient plate boundaries implies either that the coupling was much more intimate in the past, or that there is a lag of at least 100 Ma in response of the lower mantle to upper mantle conditions.     

青藏地区大地水准面形态及其与构造动力学的关系

利用密合法求解青藏地区大地水准面模型，研究青藏地区大地水准面主要特征及其与该区内部构造和活动的相关性．表明青藏地区大地水准而异常与青藏高原隆起的形状极为相似，大地水准面整体负异常与下地幔相关，形状起伏与上地幔相对应，岩石层反映的大地水准而异常与青藏高原活动构造有较好的对应关系．     

Effects of plate and slab viscosities on the geoid

The effects of plate rheology (strong plate interiors and weak plate margins) and stiff subducted lithosphere (slabs) on the geoid and plate motions, considered jointly, are examined with three-dimensional spherical models of mantle flow. Buoyancy forces are based on the internal distribution of subducted lithosphere estimated from the last 160 Ma of subduction history. While the ratio of the lower mantle/upper mantle viscosity has a strong effect on the long-wavelength geoid, as has been shown before, we find that plate rheology is also significant and that its inclusion yields a better geoid model while simultaneously reproducing basic features of observed plate motion. Slab viscosity can strongly affect the geoid, depending on whether a slab is coupled to the surface. In particular, deep, high-viscosity slabs beneath the northern Pacific that are disconnected from the surface as a result of subduction history produce significant long-wavelength geoid highs that differ from the observation. This suggests that slabs in the lower mantle may be not as stiff as predicted from a simple thermally activated rheology, if the slab model is accurate.     

On the effect of a low viscosity asthenosphere on the temporal change of the geoid—A challenge for future gravity missions

New satellite technology to measure changes in the Earth’s gravity field gives new possibilities to detect layers of low viscosity inside the Earth. We used density models for the Earth mantle based on slab history as well as on tomography and fitted the viscosity by comparison of predicted gravity to the new CHAMP gravity model. We first confirm that the fit to the observed geoid is insensitive to the presence of a low viscosity anomaly in the upper mantle as long as the layer is thin ( 200 km) and the viscosity reduction is less than two orders of magnitude. Then we investigated the temporal change in geoid by comparing two stages of slablet sinking based on subduction history or by advection of tomography derived densities and compared the spectra of the geoid change for cases with and without a low viscosity layer, but about equal fit to the observed geoid. The presence of a low viscosity layer causes relaxation at smaller wavelength and thus leads to a spectrum with relatively stronger power in higher modes and a peak around degrees 5 and 6. Comparing the spectra to the expected degree resolution for GRACE data for a 5 years mission duration shows a weak possibility to detect changes in the Earth’s gravity field due to large scale mantle circulation, provided that other causes of geoid changes can be taken into account with sufficient accuracy. A discrimination between the two viscosity cases, however, demands a new generation of gravity field observing satellites.     

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

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

利用多种地球物理观测资料直接反演地幔对流模型

假定地幔为一个均匀的、粘滞系数为常数、同时均匀分布放射性热源的流体球层,其内部存在的对流则由流体力学3个基本方程:运动方程、能量方程和连续性方程确定.如果假定地幔处于低瑞利数的状态（临界瑞利数1.5倍左右）,那么上述方程中的非线性项可以忽略不计.作为一类可能的模型,本文计算一组用6个边界条件确定6个未知数的线性方程组.这些条件包括板块绝对运动极型场、地球大地水准面异常和地震层析结果提供的地幔密度分布横向不均匀相应的“刚性地球”水准面异常等.模型计算表明:1.地幔中流体运动格局不仅受地幔热动力学参数（瑞利数）控制,而且强烈地受边界条件的影响.2.若不限定下边界为等温边界,则上、下地幔之间并不呈现出活动性明显差异;但是在模型瑞利数加大到一定值时,核-幔边界附近将出现一些局部的小尺度对流环.3.当模型瑞利数从很小增加时,对流格局将发生变化,这些格局可能反应由地幔热动力学参数决定的地幔固有特性.4.当瑞利数为50000和80000时,核-幔边界形变与PcP波得到的结果吻合较好.     

对日本俯冲带与IBM俯冲带俯冲特征的地球物理研究:来自重力与震源分布数据的启示

日本俯冲带与IBM俯冲带位于太平洋板块、菲律宾海板块和欧亚板块三者的交汇地带,是典型的"俯冲工厂"地区,具有重要的研究意义.本文利用震源分布资料与卫星重力数据对日本俯冲带与IBM俯冲带进行了研究.通过空间重力异常反映了俯冲带地区的区域构造形态,在此基础上基于艾利模式计算了均衡异常以反映地壳均衡特征.利用震源分布资料,分别从垂直俯冲带走向与沿俯冲带走向划定了横截剖面(cross-sections)进行了地震提取,讨论了俯冲带地区的Wadati-Benioff带形态特征,并借助于俯冲带地震等深线图直观描述了俯冲带的俯冲形态.在日本俯冲带与伊豆—小笠原俯冲带各选取了一条典型剖面进行了重力2.5D反演,研究了俯冲带地区的壳幔结构特征.研究结果表明,九州—帕劳海脊与IBM岛弧在均衡异常上存在差异,前者已逐渐趋向于地壳均衡.IBM的Wadati-Benioff带存在明显的南北差异,反映出伊豆—小笠原俯冲板片停留在了660km转换带中,而马里亚纳俯冲板片很可能垂直穿过了这一转换带,造成这种南北差异的原因与板块相对运动、岩石圈黏性和年龄差异以及俯冲板片的重力效应等因素有关.在IBM的中部和南部存在板片撕裂现象.日本俯冲带的俯冲洋壳密度随俯冲深度变化较小,洋幔存在一定程度的蛇纹岩化,地幔楔蛇纹岩化作用不典型,海沟处有一范围较小的含水畸变带;伊豆—小笠原俯冲带俯冲洋壳密度随深度增大而明显增大,洋幔蛇纹岩化程度较日本俯冲带低,地幔楔蛇纹岩化作用强烈,板块交汇处存在明显的蛇纹岩底辟.日本俯冲带与IBM俯冲带一线自北向南板片俯冲变陡,两侧板块耦合度降低,与俯冲带两侧的板块运动速率差异有关.     

Comparison of global geopotential models from the champ and grace missions for regional geoid modelling in Turkey

The continuous efforts on establishment and modernization of the geodetic control in Turkey include a number of regional geoid models that have been determined since 1976. The recently released gravimetric Geoid of Turkey, TG03, is used in geodetic applications where GPS-heights need to be converted to the local vertical datum. To reach a regional geoid model with improved accuracy, the selection of the appropriate global geopotential model is of primary importance. This study assesses the performance of a number of recent satellite-only and combined global geopotential models (GGMs) derived from CHAMP and GRACE missions’ data in comparison to the older EGM96 model, which is the underlying reference model for TG03. In this respect, gravity anomalies and geoid heights from the global geopotential models were compared with terrestrial gravity data and low-pass filtered GPS/levelling data, respectively. Also, five new gravimetric geoid models, computed by the Fast Fourier Transform technique using terrestrial gravity data and the geopotential models, were validated at the GPS/levelling benchmarks. The findings were also compared with the validation results of the TG03 model. The tests showed that as it was expected any of the high-degree combined models (EIGEN-CG03C, EIGEN-GL04C, EGM96) can be employed for determining the gravity anomalies over Turkey. In the west of Turkey, EGM96 and EIGEN-CHAMP03S fit the GPS/levelling surface better. However, all the tested GGMs revealed equal performance when they were employed in gravimetric geoid modelling after de-trending the gravimetric geoid model with corrector surface fitting. The new geoid models have improved accuracy (after fit) compared to TG03.     

卫星重力场、地幔对流应力场与板块运动关系的探讨

本文采用GEM10C的前50阶系数,计算了全球自由空气重力异常和地幔对流产生的岩石层底部应力场.通过对全球重力异常、应力场的分析,着重讨论了板块运动与地幔对流间的关系.对各板块在地幔对流作用下受力的定性分析和定量计算表明,板块的运动主要由低阶应力场决定,板块底部多个地幔对流体（而不是一个以洋中脊为对流上涌中心的单个地幔对流）的共同作用,推动了现代板块的运动.     

从克拉通破坏到板块动力模型的研究

华北克拉通破坏动力机制研究导致了全球动力系统及“板块动力模型”研究.板块运动最有可能的动力是地幔物质流动,但由于地幔物质流动的成因至今尚未查明,所以板块构造学说研究仍处于“运动模型”阶段,而没有进入“动力模型”阶段.如果地幔密度异常是驱动地幔物质流动的成因,那么就有可能基于重力学方法以“板块动力模型”的形式建立地幔密度异常驱动模式;软流圈中可能存在着动力特性不同的区块,地幔密度正异常代表物质盈余、区块内的物质要向区块外移动,地幔密度负异常代表物质亏损、区块外的物质要向区块内移动.本项目采用重力和地震资料相结合研究地球的整体分层,根据重力大地水准面联合地震波速度结构反演求解地幔密度异常,再根据地球正常密度假说和板块运动重力学机制的观点并与现有“板块运动模型”相结合,分析地幔密度异常动力区块,初步建立“基于重力学机制的板块动力模型”;为最终建立多学科机制的“全球板块动力模型”,迈出重要一步.     

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

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

An integrated geodynamic model of the Nankai subduction zone and neighboring regions from geophysical inversion and modeling

I investigate large-scale deep crustal structures of the Nankai subduction zone and neighboring region using regional magnetic and gravity anomalies, heat flow measurements, and earthquake hypocenters. It is found that ages, dip angles, and geothermal states of the subducting slab have direct influences on mantle wedge serpentinization. The weakest serpentinization observed in the Nankai forearc region is associated with the youngest downgoing plate of the Shikoku Basin. Conspicuous gravity anomalies identified in the forearc region are coincidental spatially with magnetic anomalies after the reduction to the pole, a mathematical procedure that helps relocate magnetic sources and boundaries, and allows us to more easily interpret magnetic data. It is argued that these patches of magnetic and gravity anomalies are caused by the same sources of anomalous density and magnetization, and are linked directly to preexisting structures such as magnetic anomalies and their boundaries in the subducting oceanic crust. Since the gravity and magnetic anomaly patches are found to be closely related to interplate seismogenic behaviors in the Nankai subduction zone, I suggest that major magnetic boundaries in the Shikoku Basin are likely weak places for slab tears that trigger seismic segmentations along the subduction zone.     

大地水准面异常在湖南地区的地球物理解释

以湖南地区为例,利用超高阶地球重力位模型EGM2008计算了研究区的重力大地水准面,并采用棱柱体公式和球体公式相结合的方法分别进行了完全地形改正和Airy-Heiskanen局部均衡改正,得到布格大地水准面和均衡大地水准面.对三种大地水准面进行不同波长分量的分离处理,得到包含不同深度异常信息的剩余大地水准面,并结合其他地球物理资料对研究区进行了详细的地球物理解释.结果表明,剩余重力大地水准面可以有效地反映出研究区内的深部构造特征,如深大断裂带分布、构造块体位置、上地幔密度横向分布等,但对地壳内异常结构反映不明显;研究区岩石圈密度变化相对平缓,厚度由东向西增加;根据剩余均衡大地水准面及研究区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.     

罗斯海海域大地水准面异常的成因研究

本文利用数值法求解瞬时地幔对流问题以模拟大地水准面异常.利用两个较新的S波速度异常层析模型SEMUCB_WM1和TX2019slab,将其转换为密度异常作为控制方程的浮力驱动项;采取的黏度结构模型中,上下地幔的黏度比为1∶50.为了研究地幔不同结构对罗斯海海域大地水准面异常的影响,分别提取上、下地幔的密度异常正/负值,作为对流控制方程的输入项,计算相应的模拟大地水准面异常.将模拟大地水准面异常与观测值进行对比,发现罗斯海海域的大地水准面异常主要来自下地幔及上地幔的负密度(波速)异常,下地幔正密度异常对该区域大地水准面负异常也有一定的贡献.本文认为,地幔密度负异常在罗斯海海域大地水准面异常的形成中占据主导作用,地幔对流的动力学效应对该区域大地水准面异常的形成影响较弱.