南极南设得兰群岛海域地质新进展:“海洋六号”南极科考初步成果

2016—2017年跨年度执行的中国第33次南极科考“海洋六号”航次, 首次获取了南极半岛南设得兰群岛北部海域、布兰斯菲尔德海峡以及欺骗岛三个区域的17 850 km2范围内多波束全覆盖高精度水深数据, 地形分辨率高于目前公开的国际数据。多道地震结果显示调查区多处存在天然气水合物地球物理标志BSR。地热测量热流值表明南设得兰群岛以北热流值与全球海洋热流值数据吻合, 不存在明显异常, 群岛以南热流值数据偏高。沉积物样品为典型的冰海沉积物, 在南设得兰群岛与南设得兰海沟之间, 沉积物以粉砂-砂为主, 夹杂砾石, 随水深增加, 沉积物粒径变细; 而在布兰斯菲尔德海区表层沉积物分布具有不对称性, 海槽南部陆架-上陆坡以及海槽东北和相邻陆架, 沉积物以砾石、砂质为主; 海槽和下陆坡为粉砂质硅质软泥。沉积物顶空气样品现场测试结果表明, 表层沉积物50 cm左右深度内的甲烷含量没有明显的差异性, 但在岛坡深水区两个测站甲烷含量在270 cm 左右深度时都存在局部高异常, 且两者总甲烷含量的变化趋势有一定的相似性, 可能与水合物分解迁移有关。     

Late Quaternary Slip Rate of the Xiugou Segment,Eastern Kunlun Fault Zone

The Eastern Kunlun fault zone (EKLF) is a large left-lateral strike-slip fault, whose slip rate is meaningful to seismic hazard assessment and geodynamics of the Tibetan Plateau. Previous studies suggested that the late Quaternary average slip rate was stable and uniform (10~13 mm/a) in the central and western segment of the EKLF. But there were a few researches of accurate slip rate in the central segment on the EKLF. Therefore, we focused on an offset and well preserved alluvial fan from Xiugou basin, located in the east of Xidatan-Dongdatan, to make it clear. Moreover, we used high-resolution satellite images and digital elevation model extracted from SPOT7 stereo image pairs to restore the offset alluvial fan, and combined terrestrial cosmogenic nuclides method, including 13 quartz-rich samples from this fan surface, 1 quartz-rich sample from the main active channel bed and 1 ¹⁰Be depth profile from this fan edge to eliminate the ¹⁰Be concentration of inheritance accurately, with 1 optically stimulated luminescence sample to obtain the reliable age of this alluvial fan together. Referring to field observations, this alluvial fan was offset left-laterally by (1 862±103) m, and its age is (76.55±3.20)~(106.37±3.38) ka which can be determined through the actual geologic setting and improving chi-square test. Thus, we used the Monte Carlo method to obtain a left-lateral slip rate of (20.3+3.5/-2.3) mm/a with 68% confidence envelopes since the late Pleistocene in the Xiugou basin. As a result, combining with the results of previous studies, the left-lateral slip rate indicated that the obviously decreasing activity transferred from late Pleistocene to Holocene on the central segment of the EKLF.     

Generation and Development of Structures and Their Controls on Ore Mineralization in the Lehong Zinc-lead Deposit,Northeastern YunnanEI北大核心CSCD

The Lehong zinc-lead deposit is one of the new-found large Zn-Pb deposits in the Sichuan- Yunnan-Guizhou Zn-Pb poly-metallic mineralization area. The ore-bodies of the Lehong deposit are strictly controlled by structures, and the fault tectonites can be divided into four groups which include morbruk rock, cataclastic rock and tectonic breccia. The mechanical properties of the faults suggest that there are five tectonic systems formed in the Caledonian-early Hercynian, Indosinian-early Yanshanian, mid-Yanshanian, late Yanshanian and Himalayan periods, respectively. The Lehong fault and the Qiaojia-Lianfeng fault are two important ore-controlling structures which might have acted as channel ways of the ore fluids. The ores are hosted in the subordinary fracture zone, joint fracture zone and interlayer fault zones in trailing folds of the Lehong fault and the Baobaoshang syncline and Jinjiagou anticline. Hence, the Lehong deposit is a tectonic -controlled ore deposit. These research do not only lay a foundation for studying deposit genesis and ore-finding direction in the Lehong deposit, but also provides important enlightenments for the ore-finding prognosis of lead-zinc deposits in northwestern Guizhou ore concentration district. © 2018, Science Press. All right reserved.     

Geometric-constrained multi-view image matching method based on semi-global optimization

Targeting at a reliable image matching of multiple remote sensing images for the generation of digital surface models, this paper presents a geometric-constrained multi-view image matching method, based on an energy minimization framework. By employing a geometrical constraint, the cost value of the energy function was calculated from multiple images, and the cost value was aggregated in an image space using a semi-global optimization approach. A homography transform parameter calculation method is proposed for fast calculation of projection pixel on each image when calculating cost values. It is based on the known interior orientation parameters, exterior orientation parameters, and a given elevation value. For an efficient and reliable processing of multiple remote sensing images, the proposed matching method was performed via a coarse-to-fine strategy through image pyramid. Three sets of airborne remote sensing images were used to evaluate the performance of the proposed method. Results reveal that the multi-view image matching can improve matching reliability. Moreover, the experimental results show that the proposed method performs better than traditional methods.     

Position-domain integrity risk-based ambiguity validation for the integer bootstrap estimator

Integrity monitoring for ambiguity resolution is of significance for utilizing the high-precision carrier phase differential positioning for safety–critical navigational applications. The integer bootstrap estimator can provide an analytical probability density function, which enables the precise evaluation of the integrity risk for ambiguity validation. In order to monitor the effect of unknown ambiguity bias on the integer bootstrap estimator, the position-domain integrity risk of the integer bootstrapped baseline is evaluated under the complete failure modes by using the worst-case protection principle. Furthermore, a partial ambiguity resolution method is developed in order to satisfy the predefined integrity risk requirement. Static and kinematic experiments are carried out to test the proposed method by comparing with the traditional ratio test method and the protection level-based method. The static experimental result has shown that the proposed method can achieve a significant global availability improvement by 51% at most. The kinematic result reveals that the proposed method obtains the best balance between the positioning accuracy and the continuity performance.     

A modified three-dimensional ionospheric tomography algorithm with side rays

The three-dimensional ionospheric tomography (3DCIT) algorithm based on Global Navigation Satellite System (GNSS) observations have been developed into an effective tool for ionospheric monitoring in recent years. However, because the rays that come into or come out from the side of the inversion region cannot be used, the distribution of the rays in the edge and bottom part of the inversion region is scarce and the electron density cannot be effectively improved in the inversion process. We present a three-dimensional tomography algorithm with side rays (3DCIT-SR) applying the side rays to the inversion. The partial slant total electron content (STEC) of side rays in the inversion region is obtained based on the NeQuick2 model and GNSS-STEC. The simulation experiment results show that the algorithm can effectively improve the distribution of GNSS rays in the inversion region. Meanwhile, the iteration accuracy has also been significantly improved. After the same number of iterations, the iterative results of 3DCIT-SR are closer to the truth than 3DCIT, in particular, the inversion of the edge regions is improved noticeably. The GNSS data of the International GNSS Service (IGS) stations in Europe are used to perform real data experiments, and the inversion results show that the electron density profiles of 3DCIT-SR are closer to the ionosonde measurements. The accuracy improvement of 3DCIT-SR is up to 56.3% while the improvement is more obvious during the magnetic storm compared to the case of a calm ionospheric state .     

A new GPS/BDS tropospheric delay resolution approach for monitoring deformation in super high-rise buildings

A new approach for deformation monitoring of super high-rise building using GPS/BDS technology is proposed for the case when prior coordinates are known and the baseline is short but has a large height difference. The approach is based on the ambiguity function method (AFM). Considering that the double-differenced (DD) troposphere delay residual error cannot be ignored, the relative zenith tropospheric delay (RZTD) parameter is introduced into the original AFM equation. Thus, the RZTD and 3D coordinate parameters are together obtained through the modified AFM (MAFM). Due to the low computational efficiency of conventional AFM, an improved particle swarm optimization (IPSO) algorithm is used to search the four optimal parameters X/Y/Z/RZTD and replaces the grid search method. In this study, GPS/BDS deformation monitoring data for buildings with approximately 290 m height difference were used to verify the feasibility of the proposed MAFM. Numerical results show a single-epoch average computation time of approximately 0.3 s, which meets the requirements of near-real-time dynamic monitoring. The average accuracy of the GPS single-epoch RZTD solution is better than 1 cm, the combined GPS/BDS MAFM performance outperforms the GPS-only system, and using multi-epoch observations can further improve the accuracy of the RZTD solution. After RZTD correction, GPS/BDS monitoring precision can be improved, particularly the height dimension, whose precision is improved by approximately 6 cm.     

Analysis of estimated satellite clock biases and their effects on precise point positioning

This study provides a systemic analysis to identify the biases in estimated satellite clocks and illustrates their effects in precise point positioning (PPP). First, the precise satellite clock estimation method considering pseudorange and carrier phase hardware delays is derived. Two methods for satellite clock estimation are compared, and their equivalency is discussed. The results show that apart from the well-known constant code hardware biases, the time-variant phase hardware biases are also absorbed by the estimated clocks. Also, the satellite clocks contain biases caused by modeling errors. To analyze the effects of these biases, they are grouped into initial clock biases (ICBs) and time-dependent biases (TDBs). Then, a detailed analysis of the impact of the biases on PPP-based troposphere and coordinate estimates is conducted. The experimental analysis demonstrates that TDBs affect positioning and tropospheric estimates, and their impacts are more significant in the static mode. The ICBs affect coordinate accuracy, zenith total delay mean bias, and its standard deviations only at the millimeter level for kinematic and static PPP, which is negligible. However, the ICBs affect the convergence period for both static and real kinematic PPP, and the magnitude of their impact largely depends on data quality. Note that satellites clocks are generally estimated with the P1/P2 and L1/L2 ionospheric-free combinations, and that hardware-specific parts of ICBs and TDBs cancel if users employ the same type of observables as the clock providers. Otherwise, the effects of biases cannot be ignored, especially for triple-frequency applications. Also, modeling-specific parts of ICBs and TDBs are significant in real-time clocks, which also affect user applications. Our conclusion is applicable for understanding the effects of these biases.     

Integrated compression of vehicle spatio‐temporal trajectories under the road stroke network constraint

With fast growth of all kinds of trajectory datasets, how to effectively manage the trajectory data of moving objects has received a lot of attention. This study proposes a spatio‐temporal data integrated compression method of vehicle trajectories based on stroke paths coding compression under the road stroke network constraint. The road stroke network is first constructed according to the principle of continuous coherence in Gestalt psychology, and then two types of Huffman tree—a road strokes Huffman tree and a stroke paths Huffman tree—are built, based respectively on the importance function of road strokes and vehicle visiting frequency of stroke paths. After the vehicle trajectories are map matched to the spatial paths in the road network, the Huffman codes of the road strokes and stroke paths are used to compress the trajectory spatial paths. An opening window algorithm is used to simplify the trajectory temporal data depicted on a time–distance polyline by setting the maximum allowable speed difference as the threshold. Through analysis of the relative spatio‐temporal relationship between the preceding and latter feature tracking points, the spatio‐temporal data of the feature tracking points are all converted to binary codes together, accordingly achieving integrated compression of trajectory spatio‐temporal data. A series of comparative experiments between the proposed method and representative state‐of‐the‐art methods are carried out on a real massive taxi trajectory dataset from five aspects, and the experimental results indicate that our method has the highest compression ratio. Meanwhile, this method also has favorable performance in other aspects: compression and decompression time overhead, storage space overhead, and historical dataset training time overhead.     

华北克拉通东缘金岭杂岩体岩浆源区及构造背景:来自岩相学、岩石地球化学及年代学的证据

金岭杂岩体由细粒黑云角闪闪长岩、中细粒黑云角闪闪长岩、中细粒黑云角闪二长闪长岩和细粒角闪二长闪长岩组成,是鲁西地区典型矽卡岩型富铁矿(金岭铁矿)控矿岩体。本次研究对细粒黑云角闪闪长岩和细粒角闪二长闪长岩进行了锆石LA‐ICP‐MS U‐Pb定年,其结果分别为129.2±3.2 Ma和132.8±1.2 Ma,表明该杂岩体的侵位时代为早白垩世。样品SiO2、K2O和Na2O含量分别介于54.17%~63.73%、1.92%~4.76%和3.10%~5.41%之间,K2O/Na2O为0.58~0.94,A/CNK为0.60~0.93,具有轻稀土富集,重稀土亏损的右倾型稀土配分模式,轻、重稀土分馏程度中等((La/Yb)N=9.94~23.49),Eu异常不明显(δEu=0.84~1.10),具中—弱负Ce异常(δCe=0.56~0.92)。样品以富集大离子亲石元素(Ba、K、U、Pb等)、亏损高场强元素(Nb、Ta、Ti)以及高Sr/Y为特征。金岭杂岩体为燕山晚期岩浆活动产物,属准铝质高钾钙碱性系列,岩浆主要来源于富集岩石圈地幔的部分熔融,并在岩浆上升侵位的过程中有地壳物质的同化混染。燕山晚期华北克拉通在古太平洋板块俯冲后后撤引起的板内伸展环境下,增厚陆壳减薄阶段,岩浆上侵就位形成金岭杂岩体。