A modified mix-differenced approach for estimating multi-GNSS real-time satellite clock offsets

A modified mixed-differenced approach for estimating multi-GNSS real-time clock offsets is presented. This approach, as compared to the earlier presented mixed-differenced approach which uses epoch-differenced and undifferenced observations, further adds a satellite-differenced process. The proposed approach, based on real-time orbit products and a mix of epoch-differenced and satellite-differenced observations to estimate only satellite clock offsets and tropospheric zenith wet delays, has fewer estimated parameters than other approaches, and thus its implementing procedure is efficient and can be performed and extended easily. To obtain high accuracy, the approach involves three steps. First, the high-accuracy tropospheric zenith wet delay of each station is estimated using mixed-differenced carrier phase observations. Second, satellite clock offset changes between adjacent epochs are estimated using also mixed-differenced carrier phase observations. Third, the satellite clock offsets at the initial epoch are estimated using satellite-differenced pseudorange observations. Finally, the initial epoch clock results and clock offset changes are concatenated to obtain the clock results of the current epoch. To validate the real-time satellite clock results, multi-GNSS post-processing clock products from IGS ACs were selected for comparison. From the comparison, the standard deviations of the GPS, GLONASS, BeiDou and Galileo systems clock results are approximately 0.1–0.4 ns, except for the BeiDou GEO satellites. The root mean squares are about 0.4–2.3 ns, which are similar to those of other international real-time products. When the clock estimates were assessed based on a pseudo-kinematic PPP procedure, the positioning accuracies in the East, North and Up components reach 5.6, 5.5 and 7.6 cm, respectively, which meet the centimeter level and are comparable to the application of other products.     

A MATLAB-based Kriged Kalman Filter software for interpolating missing data in GNSS coordinate time series

GNSS coordinate time series data for permanent reference stations often suffer from random, or even continuous, missing data. Missing data interpolation is necessary due to the fact that some data processing methods require evenly spaced data. Traditional missing data interpolation methods usually use single point time series, without considering spatial correlations between points. We present a MATLAB software for dynamic spatiotemporal interpolation of GNSS missing data based on the Kriged Kalman Filter model. With the graphical user interface, users can load source GNSS data, set parameters, view the interpolated series and save the final results. The SCIGN GPS data indicate that the software is an effective tool for GNSS coordinate time series missing data interpolation.     

Mapping hourly dynamics of urban population using trajectories reconstructed from mobile phone records

Understanding the spatiotemporal dynamics of urban population is crucial for addressing a wide range of urban planning and management issues. Aggregated geospatial big data have been widely used to quantitatively estimate population distribution at fine spatial scales over a given time period. However, it is still a challenge to estimate population density at a fine temporal resolution over a large geographical space, mainly due to the temporal asynchrony of population movement and the challenges to acquiring a complete individual movement record. In this article, we propose a method to estimate hourly population density by examining the time‐series individual trajectories, which were reconstructed from call detail records using BP neural networks. We first used BP neural networks to predict the positions of mobile phone users at an hourly interval and then estimated the hourly population density using log‐linear regression at the cell tower level. The estimated population density is linearly correlated with population census data at the sub‐district level. Trajectory clustering results show five distinct diurnal dynamic patterns of population movement in the study area, revealing spatially explicit characteristics of the diurnal commuting flows, though the driving forces of the flows need further investigation.     

An enhanced approach for surface flow routing over drainage‐constrained triangulated irregular networks

The accuracy and efficiency of the simulations in distributed hydrological models must depend on the proper estimation of flow directions and paths. Numerous studies have been carried out to delineate the drainage patterns based on gridded digital elevation models (DEMs). The triangulated irregular network (TIN) has been increasingly applied in hydrological applications due to the advantages of high storage efficiency and multi‐scale adaptive performance. Much of the previous literature focuses mainly on filling the depressions on gridded DEMs rather than treating the special cases in TIN structures, which has hampered its applications to hydrological models. This study proposes a triangulation‐based solution for the removal of flat areas and pits to enhance the simulation of flow routing on triangulated facet networks. Based on the drainage‐constrained TIN generated from only a gridded DEM by the compound point extraction (CPE) method, the inconsistent situations including flat triangles, V‐shape flat edges and sink nodes are respectively identified and rectified. The optimization algorithm is an iterative process of TIN reconstruction, in which the flat areas are generalized into their center points and the pits are rectified by embedding break lines. To verify the proposed algorithm and investigate the potential for flow routing, flow paths of steepest descent are derived by the vector‐based tracking algorithm based on the optimized TIN. A case study of TIN optimization and flow path tracking was performed on a real‐world DEM. The outcomes indicate that the proposed approach can effectively solve the problem of inconsistencies without a significant loss in accuracy of the terrain model.     

鄂尔多斯盆地延安地区山西组山2段储层可动流体赋存影响因素分析

山2段是鄂尔多斯盆地延安地区致密砂岩气生产的重要层段,目前对其流体赋存规律认识明显不足。因此在铸体薄片、扫描电镜、X射线衍射、高压压汞、恒速压汞、核磁共振测试的基础上,分析该储层可动流体的赋存特征及其影响因素。结果表明:山2储层主要为岩屑石英砂岩和岩屑砂岩,溶蚀孔和晶间孔为主要储集空间,粒间溶孔显著发育。山2储层可动流体饱和度为34.74%～91.83%,平均为69.94%,T2谱多为双峰态,呈左低右高型。储层孔隙度、渗透率、孔隙类型、孔喉特征及胶结物影响可动流体饱和度。孔隙为可动流体提供主要空间,平均孔隙半径越大,平均喉道半径越大,孔喉比越小,可动流体饱和度越高。硅质含量越高,粒间孔保存越好,可动流体饱和度越高。铁方解石含量越高,孔隙破坏越明显,可动流体饱和越低。高岭石含量较高,长石溶孔及晶间孔发育较好,有利于流体流动,可动流体饱和度也较高。     

北秦岭西段早侏罗世A型花岗岩及其地质意义

本文通过岩相学、岩石地球化学、锆石U?Pb定年和Lu?Hf同位素组成分析等方法,对出露于北秦岭西段宝鸡岩体王家山一带的黑云母花岗岩和其中的包体进行了研究。结果表明,该花岗岩形成时代为187±2 Ma,属于高钾钙碱性—钾玄岩系列岩石,富集Rb、Th、U等大离子亲石元素以及Nb、Zr和Hf等高场强元素,亏损Ba、Sr和Eu,具有高的全岩锆石饱和温度(825℃~838℃),显示A型花岗岩特征,形成于造山后的板内环境,可能为秦岭岩群副变质岩与安山质岩石部分熔融的产物。暗色包体显示塑性流变特征,具有岩浆结构,发育针状磷灰石和具有复杂成分环带的更长环斑结构长石,是幔源岩浆注入酸性岩浆发生混合作用的产物,形成时代为191±2 Ma,其锆石Hf同位素组成变化范围较大,εHf(t)值介于-11.26~-2.51,主要为富集地幔部分熔融产物。综合本文及前人已有研究结果,认为~190 Ma的早侏罗世早期秦岭地区早中生代碰撞造山过程已经结束,区域开始逐渐进入板内伸展构造演化阶段。     

滇西勐兴铅锌矿床热液方解石碳氧同位素、稀土元素特征及其指示意义

勐兴铅锌矿床是滇西保山地块内系列低温热液矿床中的典型代表,但地质科研工作相对薄弱。由于方解石是该矿床中分布广泛且与热液活动息息相关的重要脉石矿物,其形成贯穿整个成矿过程,故对其C-O同位素及稀土元素组成特征进行系统研究,以便探讨矿床成矿流体来源及演化。结果显示,该矿床两阶段方解石的δ~(13)C_(PDB)和δ~(18)O_(SMOW)值与海相碳酸盐岩较为接近,且从阶段Ⅰ→阶段Ⅱ呈微弱下降趋势,表明成矿流体中的C主要来源于围岩碳酸盐岩的溶解,阶段Ⅱ成矿流体在与围岩反应之前可能已经发生CO2热液去气作用;阶段Ⅰ方解石的稀土总量、特征参数及配分模式与结晶灰岩较为一致,而阶段Ⅱ方解石与生物碎屑灰岩较为相似,表明阶段Ⅰ方解石中的REE主要来源于结晶灰岩,而后生物碎屑灰岩对阶段Ⅱ方解石中的稀土贡献较大。以上数据及其矿床地质特征均显示该矿床为MVT型铅锌矿床,结合区域构造演化,认为勐兴矿床的形成与中特提斯洋闭合后腾冲地块与保山地块陆陆碰撞挤压密切相关。     

Effects of seismic force and pore water pressure on stability of 3D unsaturated hillslopes

Natural Hazards - Three-dimensional (3D) kinematic limit analysis of unsaturated hillslopes is presented in this paper. Different from the traditional two-dimensional (2D) mechanism based on the...     

锡同位素研究进展及其在矿床学中的应用展望

锡同位素是一种新兴的非传统稳定同位素,其在考古学及天体化学中的运用显示出非常大的示踪潜力和价值,然而目前其在地质学(尤其是矿床学)中的研究和应用前景缺乏系统介绍.本文总结分析了世界上目前所发表的主要天然及人工样品的锡同位素数据,发现天然样品中锡同位素组成有较大差异,其中玻璃陨石最富重锡(其δ122/118 Sn值可达2.53‰),而黝锡矿(黄锡矿)最富轻锡(其δ120/116Sn值可达-1.71‰).其中,含锡矿物(如锡石和黝锡矿)中的锡同位素组成变化范围要远远大于全岩样品.地幔及地壳来源的不同岩性或地质体的全岩锡同位素组成有明显差别;锡同位素在一定条件下可以发生分馏,且分馏程度可能远远大于锡同位素的初始值差异.锡石锡同位素对成矿环境非常敏感,其形成时的流体成分、化学反应速率以及物化条件(如温度、盐度、氧逸度、pH值等)等因素均能影响其锡同位素的组成.深部流体(如岩浆来源)结晶的锡石富重锡,而浅部流体(如地层流体)的加入将使锡石富轻锡,因此锡石的锡同位素具有判别不同矿床成因类型的潜力.展望未来,锡同位素的研究有望在以下三方面取得突破:①各地球圈层锡同位素储库数据的精确测定;②矿物原位微区锡同位素的准确快速分析;③热液矿床锡同位素分馏机制的建立.聚焦岩浆-热液演化过程中含锡矿物的锡同位素组成变化,有望揭示含锡流体性质及物化环境,从成矿流体来源、演化、沉淀等角度探讨成矿过程中锡同位素分馏的控制因素及其示踪机制,建立复杂锡成矿系统中的锡同位素演化模型.系统的锡同位素研究可为深入认识多类型锡矿化的"源"、"运""储"过程提供新的思路,为判别有争议锡矿床的成因类型及成矿物质来源提供关键的锡同位素证据,进而为研究大规模锡多金属成矿作用提供全新的视角,具有重要的理论价值及现实意义.     

钙同位素地球化学研究新进展及其在碳酸岩-共生硅酸盐研究中的应用

作为"非传统稳定同位素"家族成员,钙同位素正受到国际地学界日益广泛的关注.钙是主要的造岩元素,也是生物必需的元素.钙在地球各圈层广泛分布,研究钙同位素的地球化学行为将有助于提高我们对各种生物过程和地质过程的认识.钙同位素测定主要采用热电离质谱仪(TIMS)和多接收电感耦合等离子体质谱仪(MC-ICP-MS),分别表示为δ44/40 Ca和δ44/42Ca.目前自然界可观测到的δ44/40Ca变化范围为-2.0‰～6.75‰,约8.7‰.本文系统介绍了近年来钙同位素分析中样品溶解、化学分离、质谱测定以及高温地质过程中的钙同位素分馏及其地质应用等方面的研究成果,尤其对钙同位素在碳酸岩-共生硅酸盐岩研究中的意义、钙同位素组成以及取得的主要认识作了较为全面的介绍.阐述了放射成因40 Ca、部分熔融作用/分离结晶作用、地壳同化作用、古俯冲碳酸盐循环、热液蚀变作用、岩浆起源深度等对碳酸岩、硅酸盐岩的钙同位素组成造成的影响.最后,通过系统对比碳酸岩-共生硅酸盐岩的锂、镁、钙同位素研究成果,认为应该开展多元同位素体系的联合示踪.由于不同同位素体系存在相似性和差异性,而多元同位索体系相结合能有效地加强优势互补,将是同位素地球化学研究发展的一种必然趋势.