岩石－土壤－脐橙系统中稀土元素迁聚特征

在赣南稀土高背景区和稀土低背景区不同地质环境脐橙果同,系统采集岩、土、脐橙植物样,用ICP-MS法测定样品中15种稀土元素的含量,并对稀土元素在岩石、土壤和脐橙植物体内的分布、迁移、累积特征进行了研究.结果表明:a.基岩的稀土元素含量影响(决定)土壤各层位的稀土元素含量;b.土壤环境中稀土元素含量高低直接影响脐橙植物体稀土元素含量,但不完全成正比;c.基岩的稀土元素配分模式基本上贯穿土壤各层位、脐橙植物体各器官;d.岩、土、脐橙三大环节稀土无素的百分含量,总的趋势是轻稀土比重逐渐增高,重稀土比重逐渐降低,但重稀土在脐橙植物体末稍器官比例有所增加;e.岩土脐橙生态系统立体刮面中,稀土元素的垂直方向运动呈向量关系;f.由于土壤-植物壁垒作用和脐橙植物体自身的选择性吸收、控制性积累功能,高、低稀土背景区脐橙果肉的稀土含量均低于国家食品卫生限量标准.     

辽南大莲花泡全新世地层中某些元素分布规律及古环境意义

在辽南大莲花泡全新世地层早、中、晚 3个沉积阶段中 Ca元素含量表现为低 -高 -低变化 ,Al,Mg,Cr,Ni,Cu,V,B,Mn等元素呈高 -低 -高分布。因子分析揭示出控制元素这种分布的主要因素是碎屑成因和生物成因。Ca元素的“稀释剂”作用改变了其他元素的实际含量 ,造成与干凉 -暖湿 -干凉的全新世气候变化不相匹配的假象。     

新疆萨尔托海糜棱岩化石英菱镁岩元素地球化学特征及其对金成矿作用的贡献

新疆萨尔托海石英菱镁岩产在达拉布特蛇绿混杂岩带中,是蛇纹岩在剪切带深部发生热液交代作用的产物。蛇纹岩先转变为滑石片岩,再进一步转变为石英菱镁岩。石英菱镁岩局部发生剪切变形,形成糜棱岩化石英菱镁岩。与蛇纹岩相比,石英菱镁岩的MgO和SiO_2含量降低,Al_2O_3和CaO含量升高;而糜棱岩化石英菱镁岩的MgO含量较石英菱镁岩降低,Al_2O_3和SiO_2含量较石英菱镁岩升高。微量元素地球化学对比研究表明石英菱镁岩继承了蛇纹岩的微量元素特征;而糜棱岩化石英菱镁岩的微量元素含量较蛇纹岩发生了显著变化,稀土元素、高场强元素和自然金的含量明显升高,指示这些元素在剪切变形过程中发生迁移富集。剪切变形伴随强烈的流体/岩石反应,并导致剪切带流体组成和物理化学性质发生改变,从而影响矿物结晶或分解,并控制微量元素的带入迁出。蛇纹岩转变为石英菱镁岩的过程释放Au,可为糜棱岩化石英菱镁岩中的金矿化提供成矿元素。     

基于要素规则分类的城镇地籍数据建库方法研究

城镇地籍数据的有效分类与正确组织是地籍数据库规范建立和正常运行的前提条件。本文提出了基于地籍要素规则的分类方法,按照地籍数据的表达类型,描述了宗地背景图形数据、宗地图形数据、测量控制点数据、属性数据和档案数据的各个要素定义和编码规则。以成都市地籍数据库的建设为实例,根据各种来源的地籍数据分别开展了不同类型地籍数据的快速分类与入库工作。     

九江下蜀黄土的常量元素组成及其物源指示

长江中下游地区下蜀黄土的物质来源一直是个热点问题。本文对江西九江下蜀黄土和甘肃临洮晚更新世以来黄土的常量元素组成进行了研究,并与末次间冰期以来黄土高原其它点位的黄土、长江中下游其它地点的下蜀黄土和长江沉积物的元素组成进行对比,以期探讨九江下蜀黄土的物质来源。结果表明：（1）九江下蜀黄土的地球化学元素以SiO2、Al2O3、Fe2O3为主。与临洮黄土的元素组成相比,九江下蜀黄土具有较高的SiO2、Al2O3、Fe2O3、TiO2含量和较低的CaO、MgO含量。（2）在空间上,下蜀黄土的SiO2/Al2O3比值未表现出黄土高原黄土（北方黄土）的SiO2/Al2O3比值随纬度规律性变化的特点,揭示下蜀黄土的主要源区可能不是北方的内陆荒漠或黄土高原。（3）TiO2/Al2O3-K2O/Al2O3和TiO2/Al2O3-Fe2O3/Al2O3图解显示下蜀黄土与北方黄土存在显著差异,进一步表明下蜀黄土与北方黄土的主要源区可能不同。     

Garnet growth at high- and ultra-high pressure conditions and the effect of element fractionation on mineral modes and composition

In this paper we show that thermodynamic forward modelling, using Gibbs energy minimisation with consideration of element fractionation into refractory phases and/or liberated fluids, is able to extract information about the complex physical and chemical evolution of a deeply subducted rock volume. By comparing complex compositional growth zonations in garnets from high-and ultra-high pressure samples with those derived from thermodynamic forward modelling, we yield an insight into the effects of element fractionation on composition and modes of the co-genetic metamorphic phase assemblage. Our results demonstrate that fractionation effects cause discontinuous growth and re-crystallisation of metamorphic minerals in high pressure rocks. Reduced or hindered mineral growth at UHP conditions can control the inclusion and preservation of minerals indicative for UHP metamorphism, such as coesite, thus masking peak pressure conditions reached in subducted rocks.Further, our results demonstrate that fractional garnet crystallisation leads to strong compositional gradients and step-like zonation patterns in garnet, a feature often observed in high-and ultra-high pressure rocks. Thermodynamic forward modelling allows the interpretation of commonly observed garnet growth zonation patterns in terms of garnet forming reactions and the relative timing of garnet growth with respect to the rock's pressure–temperature path. Such a correlation is essential for the determination of tectonic and metamorphic rates in subduction zones as well as for the understanding of trace element signatures in subduction related rocks. It therefore should be commonplace in the investigation of metamorphic processes in subduction zones.     

Fast numerical solution of the linearized Molodensky problem

 When standard boundary element methods (BEM) are used in order to solve the linearized vector Molodensky problem we are confronted with two problems: (1) the absence of O(|x|⁻²) terms in the decay condition is not taken into account, since the single-layer ansatz, which is commonly used as representation of the disturbing potential, is of the order O(|x|⁻¹) as x→∞. This implies that the standard theory of Galerkin BEM is not applicable since the injectivity of the integral operator fails; (2) the N×N stiffness matrix is dense, with N typically of the order 10⁵. Without fast algorithms, which provide suitable approximations to the stiffness matrix by a sparse one with O(N(logN) ^s ), s≥0, non-zero elements, high-resolution global gravity field recovery is not feasible. Solutions to both problems are proposed. (1) A proper variational formulation taking the decay condition into account is based on some closed subspace of co-dimension 3 of the space of square integrable functions on the boundary surface. Instead of imposing the constraints directly on the boundary element trial space, they are incorporated into a variational formulation by penalization with a Lagrange multiplier. The conforming discretization yields an augmented linear system of equations of dimension N+3×N+3. The penalty term guarantees the well-posedness of the problem, and gives precise information about the incompatibility of the data. (2) Since the upper left submatrix of dimension N×N of the augmented system is the stiffness matrix of the standard BEM, the approach allows all techniques to be used to generate sparse approximations to the stiffness matrix, such as wavelets, fast multipole methods, panel clustering etc., without any modification. A combination of panel clustering and fast multipole method is used in order to solve the augmented linear system of equations in O(N) operations. The method is based on an approximation of the kernel function of the integral operator by a degenerate kernel in the far field, which is provided by a multipole expansion of the kernel function. Numerical experiments show that the fast algorithm is superior to the standard BEM algorithm in terms of CPU time by about three orders of magnitude for N=65 538 unknowns. Similar holds for the storage requirements. About 30 iterations are necessary in order to solve the linear system of equations using the generalized minimum residual method (GMRES). The number of iterations is almost independent of the number of unknowns, which indicates good conditioning of the system matrix. Received: 16 October 1999 / Accepted: 28 February 2001     

Two-dimensional numerical approach for the vibration isolation analysis of thin walled wave barriers in poroelastic soils

This paper is concerned with the vibration isolation efficiency analysis of total or partially buried thin walled wave barriers in poroelastic soils. A two-dimensional time harmonic model that treats soils and structures in a direct way by combining appropriately the conventional Boundary Element Method (BEM), the Dual BEM (DBEM) and the Finite Element Method (FEM) is developed to this aim. The wave barriers are impinged by Rayleigh waves obtained from Biot’s poroelasticity equations assuming a permeable free-surface. The suitability of the proposed model is justified by comparison with available previous results. The vibration isolation efficiency of three kinds of wave barriers (open trench, simple wall, open trench-wall) in poroelastic soils is studied by varying their geometry, the soil properties and the frequency. It is found that the efficiency of these wave barriers behaves similarly to these in elastic soils, except for high porosities and small dissipation coefficients. The efficiency of open trench-wall barriers can be evaluated neglecting their walls if they are typical sheet piles. This does not happen with walls of bigger cross-sections, leading in general to efficiency losses. Likewise, increasing the burial depth to trench depth ratio has a negative impact on the efficiency.     

多金属结核生长过程中元素通量的变化规律

本文主要研究太平洋ＣＣ区多金属结核ＣＣＡ３３，ＣＣＢ２９，ＣＣＣ５０，ＣＣＣ１４生长过程中锰、铁、镍、铜的通量（ＦＭｎ），ＦＦｅ，ＦＮｉ，ＦＣｕ）的变化规律．研究表明：１．结核生长过程中Ｍｎ，Ｆｅ，Ｎｉ和Ｃｕ的通量都是波动性的；２．ＦＭｎ与ＦＦｅ，ＦＮｉ，ＦＣｕ均呈正相关性；３．结核的元素通量与结核的成因类型和生长阶段密切相关，并记录了海底缺氧事件；４．结核的元素通量存在区域性差异，是沉积环境不同的反映。     

Sea state characterisation of the test site of an offshore wave energy plant

The Archimedes Wave Swing (AWS) pilot plant, rated at 2 MW, was tested offshore at Póvoa de Varzim in northern Portugal in 2004, gathering a large team of engineers and scientists. The Mechanical Engineering Department of Instituto Superior Técnico was involved, following previous work done for the project. The present paper describes one of the studies conducted. It aims to characterise the sea state at the actual location of the pilot plant using the available pressure sensors. Two approaches were performed: a first one purely based in linear wave theory, neglecting the presence of the device, and a second one, based on the results from a Boundary Element Method (BEM) code named AQUADYN, originally developed at École Centrale de Nantes (France), which allowed a detailed quantification of the effects of the presence of the plant on the wave profile directly above it. Comparisons with a Datawell Waverider buoy located at a certain distance from the plant also played an important role in the current study.