西太平洋卫星测高重力场与地球动力学特征

通过多卫星测高数据的综合处理，获得西太平洋卫星测高重力场，进行不同尺度、深度构造动力信息的分离，探讨诸边缘海盆的地球动力学问题。测高大地水准面反映了研究区板块相互作用的特点，其高频成分可以刻画各海盆的构造特征。测高空间重力异常也可刻画陆架构造及盆地分布，由其推算出的海底地形含有大量的海底构造信息。各边缘海盆的莫霍面埋深具有往南变浅的趋势，与菲律宾海各海盆的莫同埋深大致相当，说明岛弧两侧的构造力强度基本相似。大尺度地幔流应力场总体上反映了欧亚板块向东南蠕散和太平洋板块向北西扩张的特点；日本海北侧和南海巽他陆架的中尺度上地幔对流与地幔柱之间有着密切关系，西北菲律宾海的上地幔对流强化了日本－琉球－台湾－菲律宾岛弧的活动强度；上尺度地幔流主要限于软流圈层内部，在各海盆分散，而在冲绳 海槽和马里亚纳海槽则会聚，可与均衡重力异常比。还讨论了大、中、小地幔流体系的特点及相互之间的关系，籍以阐明海盆及槽演化的地球动力学过程。     

南桑威奇俯冲带大地水准面异常及其地幔对流特征分析

以南桑威奇俯冲带为例,根据EGM2008超高阶地球重力场模型、卫星重力数据为基础,利用移去-恢复原理计算了研究区大地水准面,实现了研究区不同场源深度大地水准面异常信息的分离,根据Runcorn模型计算了研究区小尺度地幔对流应力场,并结合天然地震空间展布和前人研究成果,对俯冲带结构特征与地幔对流模式进行了探讨。结果表明:南桑威奇俯冲带具有俯冲倾角较大、地震震级较低、弧前侵蚀明显等典型的马里亚纳型俯冲带特征,俯冲带南北部俯冲深度存在明显差异,中段偏北俯冲深度可达500 km;受到软流圈与上地幔上部物质密度差异的控制,东斯科舍海脊下存在沿海脊轴向南流动强地幔流;俯冲带结构与小尺度地幔对流应力场具有很强的相关性。本研究对于搞清南桑威奇俯冲带深部构造特征,理解俯冲运动、地幔对流方向及其动力控制机制提供了新的研究思路和方法。     

太平洋卫星测高重力场与地球动力学特征

通过多卫星测高数据的综合处理，获得西太平洋卫星测高重力场，进行不同尺度、深度构造动力信息的分离，探讨诸边缘海盆的地球动力学问题。测高大地水准面反映了研究区板块相互作用的特点，其高频成分可以刻画各海盆的构造特征。测高空间重力异常也可刻画陆架构造及盆地分布，由其推算出的海底地形含有大量的海底构造信息。各边缘海盆的莫霍面埋深具有往南变浅的趋势，与菲律宾海各海盆的莫霍面埋深大致相当，说明岛弧两侧的构造动力强度基本相似。大尺度地幔流应力场总体上反映了欧亚板块向东南蠕散和太平洋板块向北西扩张的特点；日本海北侧和南海巽他陆架的中尺度上地幔对流与地幔柱之间有着密切关系，西菲律宾海的上地幔对流强化了日本-琉球-台湾-菲律宾岛弧的活动强度；小尺度地幔流主要限于软流圈层内部，在各海盆分散，而在冲绳海槽和马里亚纳海槽则会聚，可与均衡重力异常类比。还讨论了大、中、小地幔流体系的特点及相互之间的关系，籍以阐明海盆及海槽演化的地球动力学过程。     

测量系统差对大地水准面精度影响分析

从边值理论出发,推导了顾及Stokes核函数改化的观测值系统误差传播公式,基于超高阶重力场模型设计了仿真试验,以定量分析顾及核函数改化的重力异常系统差对大地水准面计算精度影响。仿真试验结果表明,重力异常系统差对大地水准面的影响表现为系统性特征,顾及核函数改化的Stokes积分公式能显著改进大地水准面计算精度,即使重力异常系统误差达3 mGal,大地水准面精度也优于1 cm,这些成果将为大地水准面精化工程的方案优化和质量评估提供有益参考。     

天津港似大地水准面的精化

基于天津港基础控制测量成果与EGM2008模型,采用多面函数拟合法,进行天津港三大港区的似大地水准面精化,并结合已有成果进行精度分析。结果表明,天津港似大地水准面精化的模型检核精度为±7mm,港区内高程异常跨度为0.622m,高程异常偏差由西向东逐渐减小,南北相差不大,高程异常等值线整体沿南北走向。经工程实践验证,基于天津港似大地水准面精化模型,港区内可用D级GNSS高程测量代替三、四等水准测量,可用网络RTK高程控制测量代替四等水准测量。     

大地水准面高的计算及其在海底地形反演中的应用

为充分挖掘海洋重力数据在反演海底地形中的应用潜力,尝试探索利用大地水准面高反演海底地形的技术途径,并以夏威夷—皇帝海山链拐点所在海区作为反演试验区进行验证。首先采用Belikov列推法计算伴随(缔和)勒让德函数,利用EIGEN-6C4地球重力场模型解算获取了分辨率为1'的大地水准面高格网数值模型;然后通过综合分析反演比例函数和转换函数特点、研究海区大地水准面高与海底地形的相干特性以及大地水准面高本身尺度特征,获得了利用大地水准面高反演海底地形的频段范围;最终以试验海区大地水准面高为数据输入,构建了相应的海底地形模型(BNT模型),并与ETOPO1等海深模型进行比对分析。试验结果表明:BNT模型检核差值在一倍均方差范围检核点数量占比70.60%,相比正态分布更加集中;BNT模型检核精度低于ETOPO1等海深模型;海深模型检核精度随着水深增加不断提升,水深小于1 000 m时,海深模型相对误差出现较大发散现象;计算海域ETOPO1模型精度最高,GEBCO模型和DTU10模型检核精度相当。     

多分辨率数据输入输出法确定大地水准面

为从量化角度分析在采用输入输出法计算局部大地水准面中,输入数据(重力异常)的分辨率对计算结果的影响,以Moritz协方差函数模型估计功率谱密度函数,基于EGM2008模型重力异常计算了3种不同分辨率的输入数据下大地水准面高,仿真实验结果表明,随着分辨率的提升,将有更多的点能够参与拟合理论协方差函数模型系数,从而使得协方差函数模型在表征局部重力场特性上更加细化,大地水准面高的计算精度随分辨率的提升而提高。     

一种区域似大地水准面模型拟合方法分析

GPS技术结合(似)大地水准面模型,可快速高精度的获取测点正常高,因此精化区域(似)大地水准面成为各国家和地区建立现代高程基准的重要任务。利用GPS水准进行区域似大地水准面拟合,通过实验比较分析了二次多项式、三次多项式、多面函数、移动曲面模型,论证了各自模型的特点,实验结果表明:对于地势较平坦的实验区,以上几种方法拟合的区域似大地水准面精度都在4cm以内,都能满足一般工程应用的精度要求,其中移动曲面模型拟合的实验区似大地水准面最佳。     

一种区域似大地水准面模型拟合方法分析

GPS技术结合(似)大地水准面模型,可快速高精度的获取测点正常高,因此精化区域(似)大地水准面成为各国家和地区建立现代高程基准的重要任务。利用GPS水准进行区域似大地水准面拟合,通过实验比较分析了二次多项式、三次多项式、多面函数、移动曲面模型,论证了各自模型的特点,实验结果表明:对于地势较平坦的实验区,以上几种方法拟合的区域似大地水准面精度都在4cm以内,都能满足一般工程应用的精度要求,其中移动曲面模型拟合的实验区似大地水准面最佳。     

基于EGM2008模型的近海高程传递方法研究

研究了基于EGM2008大地水准面模型的高程传递原理,通过青岛地区已知控制点数据计算得到EGM2008大地水准面模型与我国大地水准面差距,利用EGM2008模型计算高程异常的方法,结合GNSS技术实现高程传递。利用青岛市C级GNSS控制网数据,基于EGM2008模型,采用3种不同的拟合方法,建立了区域高程异常残差场,反演高程异常值与已知控制点高程异常值进行比对。利用青岛市D级GNSS点的数据对构建的格网模型的精度进行检核。结果表明:EGM2008大地水准面模型与我国的大地水准面存在22 cm左右的偏差;基于EGM2008模型的高程传递精度可达厘米级,可用于近海高程传递,几种拟合方法精度相当。     

Residual depth and residual geoid anomalies of the South Atlantic ocean

Abstract

We have obtained the residual depth and residual geoid anomalies of the South Atlantic Ocean and interpreted them in terms of upper mantle processes. Starting from the 5’ X 5’ SYNBAPS data, we computed the 1° X 1° mean water depth. We digitized a recent sediment thickness map of the South Atlantic and corrected for the isostatic sediment loading effects. A plot of the corrected basement depth against crustal ages shows a good match to the depth‐age curve of the plate model. We then subtracted the predicted plate model depths from the corrected basement depths and obtained the 1° X 1° residual depth. The residual depth anomalies have positive values over the topographic highs and relative negative values over the ocean basins. A prominent asymmetry is observed between the residual depth over the Argentine Basin and that over the Cape Basin.

We have obtained the 1° X 1° residual geoid of the South Atlantic by subtracting both the long wavelength features and the geoid variations due to the plate cooling from the 1° X 1° Seasat altimeter derived geoid. The long wavelength features are represented by the degree and order 10 geoid of GEM1OC, and the geoid variations due to the plate cooling effects are predicted by the plate model geoid‐age relationship. The residual geoid anomalies also show an asymmetry although weaker than the case of the residual depth, between the Argentine and Cape basins.

By taking the 5° X 5° averages, we removed possible plate bending effects on the depth and geoid anomalies. We then compared those two data sets with respect to the reported hot spot tracks in the South Atlantic. The residual depth generally shows positive values over the hot spot tracks, whereas the residual geoid does not show any sign of them. A prominent asymmetric feature of depth and geoid anomalies is observed over the Argentine and Cape basins. This asymmetry is probably caused by hotter and less dense materials beneath the Cape Basin. Hot spots or other mantle upwellings may be the heat sources.     

Identification of mantle and lithospheric components of the gravity field by isostatic gravity anomalies

All anomalous masses of the Earth are reflected in the free air gravity anomalies and the geoidal undulations. The low viscosity of the asthenosphere significantly reduces the possibility of existence of density inhomogeneities in the layer. This fact provides some physical basis for the separation of the gravity field anomalies. It has been shown by power spectrum analysis of the free air anomalies and gravity field of isostatically compensated model of the lithosphere for the North Atlantic and adjacent areas of America, Europe and Mediterranean, that the attraction of isostatically compensated model is significant for any wave length of the field. It causes significant error in the interpretation if long wavelength constituents of the free air gravity anomalies are considered as a field of deep anomalous masses. The isostatic anomalies und isostatic geoid are free from the influences of isostatically compensated lithosphere. The characteristic feature of the isostatic anomalies power spectrum is a pronounced minimum at the wavelength of about 1000 km. The relative homogeneity of the asthenosphere may explain this minimum. It means that principal density inhomogeneities of the Earth's interior are separated by the asthenospheric layer. Such a minimum has not been observed at the power spectrum of free air anomalies being masked by corresponding wavelength of the field of isostatically compensated lithosphere. Isostatic anomalies that reflect the differences between the real structure of the lithosphere and its isostatically compensated model have wavelengths less than 1000 km. Isostatic anomalies with the wavelength more than 1000 km reflect the attraction of density inhomogeneities situated under the level of isostatic compensation. The basic features of power spectrum of isostatic anomalies are the same for oceanic and continental areas. The method based on Kolmogorov-Wiener filtration which consideres statistical characteristics of the field has been developed to divide the isostatic gravity anomalies into lithosphere and mantle components. For the North Atlantic and adjacent areas the field of mantle inhomogeneities has been determined.     

Relationship between the oceanic geoid and the structure of the oceanic lithosphere

Data from the GEOS 3 and SEASAT Satellites have provided a very accurate geoid map over the oceans. Broad bathymetric features in the oceans such as oceanic swells and plateaus are fully compensated. For these features it can be shown that the geoid anomalies due to the density structure of the lithosphere are proportional to the first moment of the density distribution. Deepening of the ocean basins is attributed to thermal isostasy. The thickness of the oceanic lithosphere increases with age due to the loss of heat to the sea floor. Bathymetry and the geoid provide constraints on the extent of this heat loss. Offsets in the geoid across major fracture zones can also be used to constrain this problem. Geoid-bathymetry correlations show that the Hawaiian and Bermuda swells and the Cape Verde Rise are probably due to lithospheric thinning. A similar correlation for the Walvis Ridge and Agulhas Plateau indicates that these features are probably due to an anomalously light mantle lithosphere.     

Spectral analysis of marine geoid heights and ocean depths: Constraints on models of lithospheric and sublithospheric processes

Cross-spectral analysis has been used to study the relationship between geoid and bathymetry in 16°×16° blocks in the oceans. The admittances resulting from this analysis have been compared with thermomechanical models of the lithosphere and sublithosphere in order to determine modes of topographic compensation in different parts of the oceans. Peak admittances at short wavelengths (<800 km) indicate that loads are supported by the mechanical strength of the lithosphere, while peak admittances at long wavelengths (>800 km) are indicative of lithospheric cooling or dynamic sublithospheric processes. Models of upper mantle convection predict higher admittances at long wavelengths than do models of lithospheric cooling. In most areas the observed admittances can be explained by models of the thermomechanical properties of the lithosphere, but in the eastern Pacific Ocean, the northern Indian Ocean, and over the Cape Verde Rise high long-wavelength admittances are evidence for the existence of upper mantle convection.     

Seismic study of the Ulleung Basin crust and its implications for the opening of the East Sea (Japan Sea)

The Ulleung Basin (Tsushima Basin) in the southwestern East Sea (Japan Sea) is floored by a crust whose affinity is not known whether oceanic or thinned continental. This ambiguity resulted in unconstrained mechanisms of basin evolution. The present work attempts to define the nature of the crust of the Ulleung Basin and its tectonic evolution using seismic wide-angle reflection and refraction data recorded on ocean bottom seismometers (OBSs). Although the thickness of (10 km) of the crust is greater than typical oceanic crust, tau-p analysis of OBS data and forward modeling by 2-D ray tracing suggest that it is oceanic in character: (1) the crust consists of laterally consistent upper and lower layers that are typical of oceanic layers 2 and 3 in seismic velocity and gradient distribution and (2) layer 2C, the transition between layer 2 and layer 3 in oceanic crust, is manifested by a continuous velocity increase from 5.7 to 6.3 km/s over the thickness interval of about 1 km between the upper and lower layers. Therefore it is not likely that the Ulleung Basin was formed by the crustal extension of the southwestern Japan Arc where crustal structure is typically continental. Instead, the thickness of the crust and its velocity structure suggest that the Ulleung Basin was formed by seafloor spreading in a region of hotter than normal mantle surrounding a distant mantle plume, not directly above the core of the plume. It seems that the mantle plume was located in northeast China. This suggestion is consistent with geochemical data that indicate the influence of a mantle plume on the production of volcanic rocks in and around the Ulleung Basin. Thus we propose that the opening models of the southwestern East Sea should incorporate seafloor spreading and the influence of a mantle plume rather than the extension of the crust of the Japan Arc.     

Sio Guyot: A complex volcanic edifice in the western Mid-Pacific Mountains

Sio Guyot, in the westernmost edge of the Mid-Pacific Mountains, is a large, complex volcanic edifice rising to more than 1200 m below sea level. The summit is divided into two flat-topped areas by a WNW-trending sediment-filled trough. Seismic reflection profiles reveal three acoustic units: an upper transparent layer (pelagic cap), a lower opaque layer (reef- and lagoon-derived sediments), and an acoustic (volcanic) basement. Free-air gravity anomalies indicate three eruptive centers or conduits within the main edifice, which apparently was constructed during late Cretaceous time on a broad basement swell or plateau that today is more than 3400 m below sea level (1500 m above regional abyssal depths).     

Equatorial Pacific gravity lineaments: interpretations with basement topography along seismic reflection lines

The central equatorial Pacific is interesting for studying clues to upper mantle processes, as the region lacks complicating effects of continental remnants or major volcanic plateaus. In particular, the most recently produced maps of the free-air gravity field from satellite altimetry show in greater detail the previously reported lineaments west of the East Pacific Rise (EPR) that are aligned with plate motion over the mantle and originally suggested to have formed from mantle convection rolls. In contrast, the gravity field 600 km or farther west of the EPR reveals lineaments with varied orientations. Some are also parallel with plate motion over the mantle but others are sub-parallel with fracture zones or have other orientations. This region is covered by pelagic sediments reaching ~?500–600 m thickness so bathymetry is not so useful for seeking evidence for plate deformation across the lineaments. We instead use depth to basement from three seismic reflection cruises. In some segments of these seismic data crossing the lineaments, we find that the co-variation between gravity and basement depth is roughly compatible with typical densities of basement rocks (basalt, gabbro or mantle), as expected for some explanations for the lineaments (e.g., mantle convection rolls, viscous asthenospheric inter-fingering or extensional deformation). However, some other lineaments are associated with major changes in basement depth with only subtle changes in the gravity field, suggesting topography that is locally supported by varied crustal thickness. Overall, the multiple gravity lineament orientations suggest that they have multiple origins. In particular, we propose that a further asthenospheric inter-fingering instability mechanism could occur from pressure variations in the asthenosphere arising from regional topography and such a mechanism may explain some obliquely oriented gravity lineaments that have no other obvious origin.     

Geophysical interpretation of features in the marine geoid of fennoscandia

An interpretation of the geoid in and around Scandinavia in terms of crustal depth structure has been made. Correlations as high as 0.92 were found between current models of crustal depth and the geoid for marine areas of Scandinavia. The Fennoscandian land uplift together with its corresponding resultant change of geoid were also found to be highly correlated with crustal depth structure. Results of these correlations compare favorably with theoretical models based upon large scale isostatic behavior of the lithosphere. These models indicate that the crustal thickness variations of Scandinavia are compensated generally at depths greater than 100 km. The results indicate that previous attempts to correlate the geoid with the causes of present land uplift have overestimated the remaining isostatic geoid anomaly in Fennoscandia. The application of these results would reduce the isostatic geoid anomaly by as much as 80% for marine areas of Fennoscandia. This may be interpreted as placing the estimated upper mantle viscosities for Fennoscandia closer to 10²¹ Pa s.