Position and velocity perturbations for the determination of geopotential from space geodetic measurements

Although space geodetic observing systems have been advanced recently to such a revolutionary level that low Earth Orbiting (LEO) satellites can now be tracked almost continuously and at the unprecedented high accuracy, none of the three basic methods for mapping the Earth’s gravity field, namely, Kaula linear perturbation, the numerical integration method and the orbit energy-based method, could meet the demand of these challenging data. Some theoretical effort has been made in order to establish comparable mathematical modellings for these measurements, notably by Mayer-Gürr et al. (J Geod 78:462–480, 2005). Although the numerical integration method has been routinely used to produce models of the Earth’s gravity field, for example, from recent satellite gravity missions CHAMP and GRACE, the modelling error of the method increases with the increase of the length of an arc. In order to best exploit the almost continuity and unprecedented high accuracy provided by modern space observing technology for the determination of the Earth’s gravity field, we propose using measured orbits as approximate values and derive the corresponding coordinate and velocity perturbations. The perturbations derived are quasi-linear, linear and of second-order approximation. Unlike conventional perturbation techniques which are only valid in the vicinity of reference mean values, our coordinate and velocity perturbations are mathematically valid uniformly through a whole orbital arc of any length. In particular, the derived coordinate and velocity perturbations are free of singularity due to the critical inclination and resonance inherent in the solution of artificial satellite motion by using various types of orbital elements. We then transform the coordinate and velocity perturbations into those of the six Keplerian orbital elements. For completeness, we also briefly outline how to use the derived coordinate and velocity perturbations to establish observation equations of space geodetic measurements for the determination of geopotential.     

渤海湾的有机氯农药和多氯联苯

有机氯农药(DDT，BHC)和多氯联苯(PCBs)是人工合成的有机氯化合物。它们的化学性质稳定，在环境中能持久地残留并不易受环境中各种因素的作用而降解，它们在生物体内的累积和生物链中的浓缩已构成对人类和生态系统的潜在危害，早已引起世界各国的关注，并公认为全球性的环境污染物。许多国家的学者对不同海区、各大洋海水、沿岸海水、沉积物以及海洋生物体内的DDT, BHC, PCBs作了广泛深入的调查研究(Elear, 1976; Edwand et al, 1975; Ehrhardt, 1981; Harvey, 1973; Osterroht, 1977;Robert, 1976)。我国对近海海域有机氯农药污染状况作过一些调查，但对多氯联苯未见有过报道。、1980年8月作者在渤海湾采集了26个站位的海水样品，22个站位的沉积物样品和8个毛蚶(Arasubcrenata Lischke)样品，分析了其中DDT, BHC, PCBs的含量水平，并对其在渤海湾中的分布特征作了初步探讨。     

WRF模式三维变分中背景误差协方差估计

利用WRF模式2008年5-10月逐日预报结果,通过NMC方法进行背景误差协方差(B)估计.给出其结构特征,进行单点数值试验,并利用不同B进行1个月的数值模拟试验,检验模拟降水效果.结果表明:通过单点数值试验验证估算的B结构合理.不同的B,资料同化过程差别较大,应用重新统计的B,同化效率更高,目标函数收敛更稳定.模式模...     

塔里木河下游优势草本植物与地下水埋深的关系

对塔里木河下游区域的草本群落进行了调查,获取了其分布状况和群落特征资料,并记录相应的地下水埋深,探讨了区域内距离河流不同处地下水埋深与草本群落特征之间的关系。结果表明:(1)地下水埋深越大,草本群落的生物多样性越小,覆盖度越低,生态结构趋于简单,其中地下水埋深<6 m处草本群落的种类组成和物种多样性相对较高,>7 m处草本植物分布很少,部分地区也受到微地形和地下水埋深的双重影响;(2)在大尺度空间上,草本群落随地下水的分布特点与其他荒漠地区的分布规律基本一致,随着地下水埋深的增加,草本群落的演替次序为:高水位芦苇(Phragmites communis)群落、中水位鹿角草(Glossogyne tenuifolia)群落、中低水位骆驼刺(Alhagi sparsif)群落,不同地下水埋深对应的草本群落之间具有一定的差异性、相似性和连续性;(3)塔里木河下游优势草本物种的生态位宽度均不大,主要是受水环境条件的限制,芦苇与骆驼刺的生态位宽度相对较大,在区域内分布较广。     

白云鄂博赋矿白云岩与微晶丘和碳酸岩墙的碳氧同位素对比研究

本文通过方解石和白云石的碳和氧同位素分析,对比研究了白云鄂博赋矿白云岩、黑脑包微晶丘、北京西山微晶丘、宽沟北正常沉积灰岩和白云鄂博碳酸岩墙,从而探讨了赋矿白云岩的成因及其与超大型Fe-Nb-REE矿床的成因关系。结果为:①黑脑包腮林忽洞群顶部微晶丘和北京西山寒武系顶部微晶丘碳酸盐的δ~(13)C值都在0±2‰左右,δ~(18)O值为18.3‰～25.1‰,均具有典型海相沉积碳酸盐岩的特点;②白云鄂博东矿采场δ~(13)C值为-7.9‰～-1.1‰,δ~(18)O值为9.1‰～20.9‰;矿区东西两端δ~(13)C值-7.9‰～-0.6‰,δ~(18)O值8.6‰～25.7‰;均介于地幔流体与典型沉积碳酸盐岩之间。部分赋矿白云岩样品中白云石与方解石之间的碳氧同位素分馏△~(13)C和△~(18)O值均小于0‰,表明其受到过次生蚀变作用,低δ~(18)O值白云石样品所对应的负△~(18)O值反映了地幔镁质流体对沉积碳酸盐岩的强烈交代作用;③矿区—富稀土碳酸岩墙的δ~(13)C值为-7.2‰～-4.7‰,δ~(18)O价值为11.9‰～16.4‰,表明其碳酸岩岩浆并非原始地幔来源,而可能与俯冲板块携带的沉积碳酸盐岩与地幔流体在深部的高温混合熔融有关。碳酸岩墙中白云石与方解石之间的碳和氧同位素分馏均小于0‰,说明该碳酸岩墙中的白云石与方解石并非同成因矿物,至少其中之一为     

野外型断层活动连续监测系统

介绍野外型断层活动连续监测系统的组成、布设及主要技术指标。     

大兴安岭北部新林石英二长岩-花岗岩的地球化学特征及其成因

对额尔古纳地块东南部新林镇附近的侵入岩进行了锆石U-Pb年龄测定、Lu-Hf同位素和岩石地球化学研究。结果表明,岩体形成于141.1~123.8 Ma,是早白垩世岩浆活动的产物。新林石英二长岩-花岗岩高硅较富铝贫镁,属于准铝质、高钾钙碱性系列;具有中等-弱负Eu异常(δEu=0.53~0.72),轻重稀土分异比较明显[(La/Yb)_N=12.42~25.38];相对富集大离子亲石元素Rb、Th、U和LREE,高场强元素Ti、Nb、Ta、P和HREE表现出亏损。新林石英二长岩-花岗岩具有高Sr、低Yb及低Y的特征,且相对富K,属于C型埃达克岩,形成于加厚下地壳的部分熔融。Lu-Hf同位素研究显示其ε_(Hf)(t)为2.1~6.2,二阶段模式年龄(T_(DM2))为1025~786 Ma,认为其岩浆源区为新元古代增生的地壳物质。     

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

苏北平原MIS 3阶段海侵事件的沉积记录

MIS 3阶段的气候与环境是古气候与古环境研究的热点,该阶段中国沿海地区的海侵过程是沿海平原环境演化的重要事件。根据苏北平原宝应By1钻孔的 14 C年代结果和沉积特征、微体生物组合特点,确认苏北平原MIS 3阶段存在明显的海侵事件,潮滩与泻湖相沉积层的厚度约7m,深度在By1孔19~26m,海相层位于长江三角洲地区的典型末次冰盛期的硬粘土层之下,14 C年代在距今26~39ka之间,属于MIS 3阶段的沉积,岩芯沉积物的沉积构造以水平的薄层砂质粉砂沉积为主。沉积物中大量的海相有孔虫壳体在研究的岩芯段中是最丰富的; 此外,还发现咸水环境下的Sinocytheridea impressa (Brady),Spinileberis pulchra与Spinileberis furuyaensis,Neomonoceratina chenae Zhao,Keijella bisanensis (Okubo)等形类壳体; 该段岩芯中的植物残体、碳屑也很丰富,有机质含量丰富,显示比较封闭的弱动力环境,指示了潮滩、泻湖沉积环境。     

红层地区砂泥岩互层状斜坡岩体流变特性研究

本文在对侏罗系砂溪庙组砂，泥岩进行的系统流变试验的基础上，详细研究了砂，泥岩的流变特性及其相应的长期强度，以及砂，泥岩互层状组成的斜坡碉体流变特征，特别是模拟研究了风化深度较厚的泥岩产生的粘塑性流动对斜坡应力场及斜坡变形破坏所起的不可忽视的作用。