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.     

多相糜棱岩内第二相对应变局部化的影响——以秦岭群花岗质糜棱岩为例

糜棱岩韧性变形发生的应变局部化过程,尤其是多相糜棱岩第二相对基质相变形的影响一直是显微构造研究难点.研究表明糜棱岩借助颗粒边界滑移实现多相混合,形成多矿物相集合体.在多相糜棱岩内,第二相在基质相颗粒边界施加齐纳阻力,牵制基质相颗粒边界的迁移速率,破坏基质相颗粒的动态平衡过程,使基质相颗粒位于古应力计对应的颗粒粒度以下,导致基质相整体的表面积增大,促进扩散交换过程,提高了扩散蠕变,降低了基质相位错蠕变和结晶学优选方位(CPO)形成的效率,使变形机制从颗粒粒径不敏感蠕变机制(GSI)过渡为颗粒粒径敏感蠕变机制(GSS).另外,多相糜棱岩内的第二相具有诱导应变局部化的效应,使塑性应变局部化更为强烈,引起物质强度的变化,进而引起岩石变形过程和岩石圈流变行为的改变.选取秦岭群花岗质糜棱岩进行多相矿物糜棱岩定量化研究,结果显示花岗质糜棱岩伴随着云母含量的增多以及各相混合程度的增大,石英的颗粒粒度明显减小,CPO强度显著降低,基质相显微变形受第二相控制逐渐增强.     

渤海湾盆地渤中19-6气田凝析油特征及油气关系

2018年中国东部最大整装凝析气田——渤中19-6气田的发现,一举打开了渤海湾盆地深层天然气勘探的新领域,展现了油型盆地天然气勘探的巨大潜力。渤中凹陷发育多套不同成熟度的烃源岩,位于渤中凹陷西南部的渤中19-6气田凝析油与天然气并存,近年来于该地区的天然气成因及形成条件多有研究,但是对于凝析油的形成则少有研究。本文以渤中19-6气田为例,探讨多套烃源岩发育的复杂油气区凝析油的成因及与天然气的关系。渤中19-6气田6个凝析油样品化学组成及其变化特征精细分析对比表明:该区凝析油以低碳数正构烷烃为主,含有较丰富的环己烷、甲基环己烷等环烷烃和苯、甲苯等低碳数芳烃,金刚烷类化合物较发育,高碳数生物标志化合物含量较低,综合判断凝析油属于混合偏腐殖型高熟原油。结合该区凝析油与天然气对比结果认为,该区凝析油和天然气属于同一套烃源岩在同一阶段所生成的烃类化合物,两者相伴而生,均属于烃源岩原生产物。     

基于Sentinel-1和DEM数据的南岭高植被覆盖区地形线性特征提取方法

在南岭高海拔山脉、植被茂密等地区开展地质调查难度大,人员和物资难以进入,因此需提前根据该区域的地形线性特征合理设计野外调查路线;此外地形线性特征可为区域构造运动信息的认知提供辅助知识.遥感技术以其宏观性、多尺度、多层次的特点成为地质研究和地质勘查的重要手段,然而在高植被覆盖区,常见的光学影像难以穿透植被,而微波雷达遥感则因为较好的植被穿透性,使得其应用成为可能.因此,提出了一种基于雷达卫星Sentinel-1和数字高程模型（digital elevation model,DEM）数据的地形线性特征提取方法.该方法首先利用似然比边缘检测算法提取Sentinel-1影像中的边缘特征,然后对DEM进行线性特征增强以生成山体阴影影像,再利用Canny边缘检测算子提取DEM数据中的山脊线与山谷线等主要线性特征,以DEM提取出的特征线为中心建立缓冲区并与雷达影像的提取结果做相交处理,将得到的线性特征利用道格拉斯-普克算法进行局部直线拟合,最后得到研究区的地形线性特征.结果表明,该方法综合考虑了雷达影像的微观细节信息和DEM数据的宏观趋势信息,在保留主要线性地形特征的同时剔除了伪边缘和噪声点,提取效果较好.      

Variations of methane stable isotopic values from an Alpine peatland on the eastern Qinghai-Tibetan Plateau

Acta Geochimica - Isotopic signature is a powerful tool to discriminate methane (CH4) source types and constrain regional and global scale CH4 budgets. Peatlands on the Qinghai-Tibetan Plateau are...     

Correction: Hydraulic tomography analysis of municipal-well operation data with geology-based groundwater models

Hydrogeology Journal - A Correction to this paper has been published: https://doi.org/10.1007/s10040-021-02351-x     

华南上泥盆统法门阶哲杜斯贝动物群的埋藏学和居群动力学特征*

哲杜斯贝是一种以高丰度和低物种分异度为显著特征、在全球广泛独立分布的小嘴贝类。笔者对华南台地相和台间盆地相5个含哲杜斯贝类上泥盆统剖面进行沉积环境和埋藏特征研究,确定华南哲杜斯贝类的原生沉积环境为台地的正常浅水,这与美国和摩洛哥生长在冷泉环境中的哲杜斯贝不同,表明其既可在浅海环境中生长,也可以在冷泉环境中生长,具有较强的生态可塑性。台地相(原地埋藏)的3个哲杜斯贝种群具有以下特征:(1)均呈现出壳体呈层状紧密堆积、集群生长特点,丰度极高,并且有相当数量的个体呈现出由于拥挤生长而导致的壳体生长畸形现象;(2)缺乏其他底栖后生生物,显示了极低的生物多样性;(3)壳体大小—频率分布和生存曲线揭示了3个哲杜斯贝种群都具有非常低的幼年期死亡率特点,同时同一介壳层内哲杜斯贝个体大小均一,表明绝大多数的个体能够达到或接近自然寿命后死亡。上述证据表明,哲杜斯贝具有密集群居生活的习性,其生活在一个独特的、不适合其他后生底栖生物生存的环境。     

东昆仑造山带低山头花岗岩体岩石地球化学特征及其地质意义

低山头一带花岗岩体归属东昆仑弧盆系北昆仑岩浆弧带,位于东昆仑成矿带伯喀里克—香日德成矿亚带东昆仑造山带中段,岩石类型主要为石英闪长岩、正长花岗岩及二长花岗岩。为了加强该地区花岗岩体岩石地球化学与成岩成矿背景探讨,对花岗岩体开展了主量元素、微量元素及矿质元素研究。主量元素组成上,石英闪长岩具中硅（57.64%和58.47%）、富钠（Na₂O/K₂O均为2.57）特点,正长花岗岩具高硅（75.45%~75.99%）、富钾（Na₂O/K₂O为0.74~0.94）特点,二长花岗岩具高硅（66.80%~73.45%）、富钠（Na₂O/K₂O为1.50~2.13）特点;花岗岩体铝饱和指数A/CNK<1,为准铝质岩浆岩;碱饱和指数NK/A集中在0.26~0.69之间,属钙碱性岩石;里特曼组合指数σ₄₃在1.18~2.31之间,属钙碱性类型。花岗岩体轻稀土元素相对重稀土元素富集,岩浆分异特征明显,大离子亲石元素（LILE）富集Rb、K、Ba、Th、Sr、Nd,仅正长花岗岩Sr亏损,高场强元素（HFSE）富集Zr、Hf、Ce,而Nb、P、Ti明显亏损,源区物质为壳幔混合物质,属挤压应力环境中同碰撞I型花岗岩。在岩体实测剖面中获得11种元素分析数据,与青海全省、东昆仑成矿带及其亚成矿带平均元素丰度值进行对比,初步划分不同岩性、不同类型、不同时代花岗岩,以及富集的含矿元素为Au、Zn、Y、Pb等。与区域有成矿事实且为I型花岗岩成因进行对比,认为研究区有较好的找矿前景。     

近40年来青藏高原典型高寒草原和湿地蒸散发变化的对比分析

青藏高原地表蒸散发是决定亚洲水塔水储量变化的关键要素.在快速升温背景下,长时间尺度的青藏高原地表蒸散发如何响应气候变化亟需深入探讨.以青藏高原两种典型高寒生态系统(草原和湿地)为研究对象,以野外观测和互补蒸散发模型为研究手段,利用常规气象资料驱动互补蒸散发模型,应用于青藏高原的典型资料稀缺地区,并就模拟结果进行验证评估,揭示了两种典型高寒生态系统近40年的蒸散发变化特征.结果 表明,校正参数后的非线性互补蒸散发模型可较为准确地模拟两种下垫面的蒸散发,亦即该模型在青藏高原资料稀缺区具有较好的应用潜力.1973-2013年,青藏高原典型高寒草原蒸散发呈不显著的增大趋势,而高寒湿地则以2.0 mm/a的速率显著增大.相关分析表明,高寒草原和湿地蒸散发的年际变化主要与水汽压(即空气湿度)有关.阶段性分析发现,1970s至1990s末期,两种生态系统蒸散发皆在波动中逐渐增大;而1997年以后,高寒草原和高寒湿地蒸散发的变化模式表现出明显差异:前者在波动中逐渐减小,后者则持续增大至2000s中期.造成这种差异的原因可归结为高寒湿地受冰川融水的影响,土壤含水量可维持在较高的水平,加之2000s高寒湿地的水汽压和日照时数增大,使得该时段内地表蒸散发仍呈增大之势,亦即上游的冰川融水对下游的湿地蒸散发有重要影响.结果 表明,空间距离较近的两种典型高寒生态系统,由于所受水源补给不同,局地蒸散发对气候变化的响应模式可能有较大差异.     

Integration of zircon and apatite U–Pb geochronology and geochemical mapping of the Wude basalts(Emeishan large igneous province):A tool for a better understanding of the tectonothermal and geodynamic evolution of the Emeishan LIP

Radiogenic isotopic dating and Lu–Hf isotopic composition using laser ablation-inductively coupled plasma-mass spectrometry(LA-ICP-MS)of the Wude basalt in Yunnan province from the Emeishan large igneous province(ELIP)yielded timing of formation and post-eruption tectonothermal event.Holistic lithogeochemistry and elements mapping of basaltic rocks were further reevaluated to provide insights into crustal contamination and formation of the ELIP.A zircon U–Pb age of 251.3±2.0 Ma of the Wude basalt recorded the youngest volcanic eruption event and was consistent with the age span of 251-263 Ma for the emplacement of the ELIP.Such zircons hadε_Hf(t)values ranging from7.3 to+2.2,identical to those of magmatic zircons from the intrusive rocks of the ELIP,suggesting that crust-mantle interaction occurred during magmatic emplacement,or crust-mantle mixing existed in the deep source region prior to deep melting.The apatite U–Pb age at 53.6±3.4 Ma recorded an early Eocene magmatic superimposition of a regional tectonothermal event,corresponding to the Indian–Eurasian plate collision.Negative Nb,Ta,Ti and P anomalies of the Emeishan basalt may reflect crustal contamination.The uneven Nb/La and Th/Ta values distribution throughout the ELIP supported a mantle plume model origin.Therefore,the ELIP was formed as a result of a mantle plume which was later superimposed by a regional tectonothermal event attributed to the Indian–Eurasian plate collision during early Eocene.