A Reanalysis and Reinterpretation of Geodetic and Geological Evidence of Glacial-Isostatic Adjustment in the Churchill Region, Hudson Bay

We review the historical, geological, tide-gauge, GPS and gravimetric evidence advanced in favour of, or against, continuing land uplift around Hudson Bay, Canada. We also reanalyse the tide-gauge and GPS data for Churchill using longer time series than those available to previous investigators. The dependence of the mean rate of relative sea-level change obtained from the tide-gauge record on the length and mid-epoch of the observation interval considered is investigated by means of a newly developed linear-trend analysis diagram. For studying the shorter-period variability of the tide-gauge record, the wavelet transform is used. The mean rate of land uplift obtained from GPS is based on a new analysis using IGS solutions of GFZ. To include the post-glacial land uplift, sea-level indicators from the Churchill region representing the relative sea-level history during the past 8000 years are also used. Finally, the values of the four observables are jointly inverted in terms of mantle viscosity. The optimum values are ~3.2 × 10²⁰ Pa s and ~1.6 × 10²² Pa s for the upper- and lower-mantle viscosities, respectively.     

LBM钻井液在绳索取心钻进中的应用

LBM钻井液是一种低粘度低失水量的高效造浆材料，通过在绳索取心中应用，证明了LBM钻井液具有良好的抑制水敏性地层膨胀与分散的能力．可有效防止粘附卡钻事故及解决钻杆内壁结垢问题．是绳索取心钻进较理想的钻井液材料。     

The Role of History-Dependent Rheology in Plate Boundary Lubrication for Generating One-Sided Subduction

We have developed a two-dimensional dynamical model of asymmetric subduction integrated into the mantle convection without imposed plate velocities. In this model we consider that weak oceanic crust behaves as a lubricator on the thrust fault at the plate boundary. We introduce a rheological layer that depends on the history of the past fracture to simulate the effect of the oceanic crust. The thickness of this layer is set to be as thin as the Earth's oceanic crust. To treat 1-kilometer scale structure at the plate boundary in the 1000-kilometer scale mantle convection calculation, we introduce a new numerical method to solve the hydrodynamic equations using a couple of uniform and nonuniform grids of control volumes. Using our developed models, we have systematically investigated effects of basic rheological parameters that determine the deformation strength of the lithosphere and the oceanic crust on the development of the subducted slab, with a focus on the plate motion controlling mechanism. In our model the plate subduction is produced when the friction coefficient (0.004–0.008) of the modeled oceanic crust and the maximum strength (400 MPa) of the lithosphere are in plausible range inferred from the observations on the plate driving forces and the plate deformation, and the rheology experiments. In this range of the plate strength, yielding induces the plate bending. In this case the speed of plate motion is controlled more by viscosity layering of the underlying mantle than by the plate strength. To examine the setting of the overriding plate, we also consider the two end-member cases in which the overriding plate is fixed or freely-movable. In the case of the freely-movable overriding plate, the trench motion considerably changes the dip angle of the deep slab. Especially in the case with a shallow-angle plate boundary, retrograde slab motion occurs to generate a shallow-angle deep slab.     

青藏高原新生代火车头山碱性及钙碱性两套火山岩的地球化学特征及其物源讨论

藏北羌塘火车头山新生代火山岩可区分为钙碱性及碱性两个不同的系列。钙碱性火山岩主要岩石组合为玄武岩-安山岩-英安岩,其SiO_2介于49%～70%之间,Al_2O_3>10%,Na_2O/K_2O>1;其中玄武岩具平坦型稀土配分型式,LREE/HREE为1.3～1.8,(La/Yb)_N为2.87～4.45,无明显铕异常,δEu为0.96～1.09;该套岩石的Mg~#与SiO_2相关关系以及La/sm-La等亲岩浆元素与超亲岩浆元素协变关系表明,它们应为幔源岩浆经分离结晶演化的产物,其岩石组合类型以及低的Sm/Yb值(Sm/Yb=1.53～5.35)表明它们的原始岩浆应来源于岩石圈地幔尖晶石二辉橄榄岩的局部熔融。本区碱性火山岩为一套典型的钾质岩石系列,主要岩石组合类型为碱玄岩-碱玄质响岩-响岩,其SiO_2介于44%～59%之间,Al_2O_3>14%,Na_2O/K_2O介于0.47～1.51之间;岩石轻稀土强烈富集,LREE/HREE为13.20～15.76,(La/Yb)_N=50.44～91.99;其岩石组合类型以及Mg~#与SiO_2相关关系以及La/Sm-La协变关系同样表明它们为共源岩浆分离结晶演化的产物;然而,其较高的Sm/Yb值(Sm/ Yb=2.63～13.98)表明它们并非地幔橄榄岩直接局部熔融的产物,岩石弱的负Eu异常(δEu=0.77～0.85)以及Th、U的强烈富集和Nb、Ta的相对亏损,又反映了原始岩浆中有显著的地壳物质的贡献;该套钾质碱性系列岩石在La/Co-Th/Co同分母协变图上呈直线型分布,而在La/Co-Sc/Th异分母协变图上呈显著的双曲线分布,从而表明其源区为二源混合型,是青藏高原特殊的壳幔混合层局部熔融的产物,这些特征是新生代青藏高原壳幔层圈物质交换的重要岩石学证据。     

吉林蛟河地幔橄榄岩捕虏体的Sr-Nd-Hf-Os同位素特征与岩石圈地幔时代

吉林省蛟河市境内大石河新生代玄武岩中含有丰富的地幔橄榄岩包体,详细的岩石学与矿物学研究显示,这些包体的主要岩石类型为尖晶石二辉橄榄岩-方辉橄榄岩,未发现石榴石橄榄岩。岩相学及地球化学资料显示它们都是经历过熔体抽取而形成的岩石圈地幔残留。矿物平衡温度计算发现,本区的这些地幔橄榄岩包体来自地下40～60km 深度,且下部以二辉橄榄岩为主,而上部以贫单斜辉石的二辉橄榄岩和方辉橄榄岩为主,显示明显的岩石圈地幔分层现象。Sr-Nd-Hf 同位素资料反映这些地幔包体均表现为亏损性质,而 Re-Os 同位素资料确定上述岩石圈地幔形成于中元古代,明显老于上覆地壳的新元古宙时代,反映壳幔年龄上的解耦。因此我们推测,该区曾经历过华北克拉通类似的早期岩石圈地幔的整体丢失事件,然后形成于其它地区的中元古宙岩石圈地幔在本区增生。     

高温高压下水的剪切粘滞度的平衡分子动力学计算

本项工作用ＴＩＰ４Ｐ模型对水分子进行了平衡分子动力学模拟。根据Ｇｒｅｅｎ－Ｋｕｂｏ理论,计算了温度从６５２Ｋ到８１４Ｋ,密度从０．５３６０ｇ／ｃｍ３到０．８６３８ｇ／ｃｍ３六个不同温度和密度水的剪切粘滞度,并和实验值进行了对比。在所计算的温度和密度范围,平衡分子动力学计算结果与实验值的平均绝对值差别为１５％,好于使用非平衡分子动力学方法获得的粘滞度值。平衡分子动力学的方法将是获得地学中高温高压下流体输运性质的重要途径。     

Flow of groundwater with variable density and viscosity, Atikokan Research Area, Canada

 Flow of groundwater with variable density and viscosity was simulated at the Atikokan Research Area (ARA) in northwestern Ontario, Canada. An empirical viscosity–concentration equation was modified to include total-dissolved-solids (TDS) data from the ARA. The resulting equation was used successfully to estimate reasonably accurate viscosity values over the expected range of temperature and concentration, in comparison with experimental values derived for sodium chloride solutions. A three-dimensional finite-element code, MOTIF, developed by Atomic Energy of Canada Limited, was used in the simulations. The inclusion of the effects of depth-increasing temperature and TDS-dependent fluid-density distribution, while maintaining only a temperature-dependent viscosity relationship in a simulation, resulted in a more penetrative flow against expected buoyancy effects (i.e., the physics of the system was not honored). Accounting for concentration in the viscosity equation caused water to be less penetrative and more in accordance with the expected physics of the system. A conclusion is that fluid concentration should be considered simultaneously in calculating the density and viscosity of a fluid during modeling of variable-density flow in areas underlain by fluids with high TDS. Results of simulations suggest that both flow directions and magnitudes should be employed simultaneously during the calibration of a model. Large-scale groundwater movement in the ARA may be analyzed with carefully selected vertical no-flow boundaries. By incorporating the geothermal temperature gradient, groundwater recharge increases by 12%; thus, this gradient plays a significant role in groundwater flow at the ARA. Variability in the fluid concentration at the ARA neither decreases nor increases recharge into the groundwater system. The hypothesis that an isolated continuous regional flow system may exist at depth in the ARA is not supported by these simulations. Received, September 1996 Revised, September 1997, February 1998 Accepted, February 1998     

利用空间测地技术结果估计地球平均下地幔粘性(英)

基于IEC－4G冰期后地壳反弹模型,和地球上Laurentia,Fennoscandia,Antarctica,andGreenland四大冰盖最近21000年以来的冰融参数,计算了对地球最大主转动惯量的影响△I33,并进而由现代空间测地技术观测资料分析得到的地球自转非潮汐加速项为约束,估计了地球平均下地幔（670km以下）粘性vLM为（0．9～2．5）·1022Pas,这个结果表明了vLM应具有1022Pas量级．