青藏高原东南缘现今地壳形变与多尺度应变率场特征

利用青藏高原东南缘1999~2007年与2011~2017年高精度GNSS监测资料,充分顾及GNSS测站非均匀分布的特性,构建青藏高原东南缘地壳形变多尺度球面小波模型,定量分析汶川强震前后该区域地壳形变与不同空间尺度下地壳应变率场变化特征。研究结果表明,汶川强震前后研究区整体地壳运动具有一定继承性发展特征,均在龙门山、安宁河、则木河与小江断裂处形成明显的速度差异梯度带,且高应变率值也主要聚集在上述主干断裂及附近区域;汶川强震后,研究区域整体地壳运动速率量值,特别是龙门山断裂带西北侧地壳运动速度量值显著增大,羌塘、巴颜喀拉与川滇地块也呈现出加速向南东运移并推动华南块体的趋势。不同尺度下的应变率场反映出不同空间范围下区域应变积累特征,青藏高原东南缘区域在尺度因子q=7时的计算结果是合理的（合理的最大尺度因子）;当尺度因子q=6时,能较好地揭示出区域整体构造活动特性,即清晰地揭示出汶川强震后龙门山断裂处呈现出的显著主压应变、面压缩与最大剪应变率高值特征;当尺度因子q=3~7时（最佳组合尺度因子）,可较好地综合揭示出区域地壳大尺度（整体）形变与局部形变特征。汶川强震后,研究区域主干断裂带,特别是震中及其附近区域地震活动性显著增强。     

豫西南泥湖地区辉长岩地球化学特征及其与成矿关系

豫西南泥湖钼矿集区是以花岗斑岩体为中心的钼多金属矿床系列,区内出露大面积加里东期变辉长岩和燕山期辉长岩。通过分析区内辉长岩的岩石学、年代学、地球化学和成矿元素等方面特征,研究其在成岩成矿过程中的作用。区内加里东期变辉长岩和燕山期辉长岩成矿元素总稀土含量低,稀土元素球粒陨石标准化曲线右倾平缓,具有明显的Eu,Tm元素的异常;同比世界基性岩,自晚元古代上地幔主要成矿元素(除钨元素以外)基本没有明显的富集,主成矿元素钼含量呈明显的亏损,预示燕山期成矿花岗质岩浆没有从上地幔获得主要成矿物质。初步认为大规模的钼多金属成矿带的形成与基性岩浆的底垫作用相关,基性岩浆的底垫作用为下地壳部分熔融形成岩浆房提供了热源。     

招平断裂带中段深部结构构造——来自地震剖面的证据

玲珑花岗岩岩基东侧的招平断裂带是胶西北矿集区的主要控矿构造之一。布设于招平断裂带中段大尹格庄矿区的地震剖面,揭示了招平断裂带中段深部的结构构造特征,地震剖面显示,玲珑岩基内部有多个类似"拱弧构造"的岩浆侵入形成的弧形界面,由此推测存在2~3期的岩浆活动。剖面显示了3期构造组合样式,第1期为以招平断裂带和宋家河断裂为组合倾向SE的铲式断裂,底部似汇入统一的构造滑脱带;第2期为倾向NW的错断了招平断裂带的反倾伸展断裂;第3期为NE向延伸以花状构造为特征的走滑断裂。前两期反映构造发育处于伸展环境,第3期显示为左行走滑。招平断裂带深部延伸存在分支复合现象,形成由2个主裂面夹持的透镜状岩片。地震反射剖面数据的获得对深刻理解该区域中生代岩浆构造演化过程和深部找矿具有重要的借鉴意义。     

山东龙口市大磨曲家金矿地质特征及找矿方向探讨

大磨曲家金矿位于玲珑金矿田东风矿床主脉的NE延伸部位.矿区内已发现矿脉25条.通过对玲珑金矿田的成矿条件分析和东风矿床的综合研究,认为本区有很好的找矿前景.     

Coupling winds to ocean surface currents over the global ocean

A Wind stress–Current Coupled System (WCCS) consisting of the HYbrid Coordinate Ocean Model (HYCOM) and an improved wind stress algorithm based on Donelan et al. [Donelan, W.M., Drennan, Katsaros, K.B., 1997. The air–sea momentum flux in mixed wind sea and swell conditions. J. Phys. Oceanogr. 27, 2087–2099] is developed by using the Earth System Modeling Framework (ESMF). The WCCS is applied to the global ocean to study the interactions between the wind stress and the ocean surface currents. In this study, the ocean surface current velocity is taken into consideration in the wind stress calculation and air–sea heat flux calculation. The wind stress that contains the effect of ocean surface current velocity will be used to force the HYCOM. The results indicate that the ocean surface velocity exerts an important influence on the wind stress, which, in turn, significantly affects the global ocean surface currents, air–sea heat fluxes, and the thickness of ocean surface boundary layer. Comparison with the TOGA TAO buoy data, the sea surface temperature from the wind–current coupled simulation showed noticeable improvement over the stand-alone HYCOM simulation.     

用曲面插值方法建立海洋局部地磁场模型

地磁场模型是地磁场的数学表达式,国内外已有多种建立地磁场模型的方法,对于局部地磁场模型来说,应用比较多的是多项式方法和曲面样条函数方法。讨论了基于曲面样条函数建立局部地磁场模型的原理和方法,并在实测数据基础上对模型的精度进行了验证。     

Calibrating the GOCE accelerations with star sensor data and a global gravity field model

A reliable and accurate gradiometer calibration is essential for the scientific return of the gravity field and steady-state ocean circulation explorer (GOCE) mission. This paper describes a new method for external calibration of the GOCE gradiometer accelerations. A global gravity field model in combination with star sensor quaternions is used to compute reference differential accelerations, which may be used to estimate various combinations of gradiometer scale factors, internal gradiometer misalignments and misalignments between star sensor and gradiometer. In many aspects, the new method is complementary to the GOCE in-flight calibration. In contrast to the in-flight calibration, which requires a satellite-shaking phase, the new method uses data from the nominal measurement phases. The results of a simulation study show that gradiometer scale factors can be estimated on a weekly basis with accuracies better than 2 × 10⁻³ for the ultrasensitive and 10⁻² for the less sensitive axes, which is compatible with the requirements of the gravity gradient error. Based on a 58-day data set, scale factors are found that can reduce the errors of the in-flight-calibrated measurements. The elements of the complete inverse calibration matrix, representing both the internal gradiometer misalignments and scale factors, can be estimated with accuracies in general better than 10⁻³.     

Non-stationary covariance function modelling in 2D least-squares collocation

Standard least-squares collocation (LSC) assumes 2D stationarity and 3D isotropy, and relies on a covariance function to account for spatial dependence in the observed data. However, the assumption that the spatial dependence is constant throughout the region of interest may sometimes be violated. Assuming a stationary covariance structure can result in over-smoothing of, e.g., the gravity field in mountains and under-smoothing in great plains. We introduce the kernel convolution method from spatial statistics for non-stationary covariance structures, and demonstrate its advantage for dealing with non-stationarity in geodetic data. We then compared stationary and non- stationary covariance functions in 2D LSC to the empirical example of gravity anomaly interpolation near the Darling Fault, Western Australia, where the field is anisotropic and non-stationary. The results with non-stationary covariance functions are better than standard LSC in terms of formal errors and cross-validation against data not used in the interpolation, demonstrating that the use of non-stationary covariance functions can improve upon standard (stationary) LSC.     

低轨卫星初始状态误差对反演地球重力场的影响

基于动力学法反演地球重力场的基本理论,研究了卫星初始状态向量误差对应用低轨卫星精密轨道数据反演地球重力场的影响。在仅考虑低轨卫星初始状态误差的情况下进行了模拟计算,结果表明:在利用低轨卫星精密轨道数据反演地球重力场时,卫星初始状态向量误差需要重新进行估计;在目前的轨道精度水平下,若不顾及误差方程二次项的影响,反演弧长不宜过长;卫星初始状态速度误差(约1.5mm/s)的影响要大于位置误差(约10 cm)的影响。     

基于能量守恒的星间距离变率与地球重力场频谱关系的建立与分析(英文)

The spectral relationship between range-rate measurements and the gravity potential for low-low satellite-to-satellite tracking mission was established based on the energy conservation theory. Then the performances of satellite separation, the orbital altitude, and the accuracy of range-rate measurements in recovering the earth’s gravity field were simulated and analyzed by this method. Finally, the cumulative geoid errors of the reference mode were obtained by using the configuration parameters of the GRAC...