SELFE: A semi-implicit Eulerian–Lagrangian finite-element model for cross-scale ocean circulation

Unstructured-grid models grounded on semi-implicit, finite-volume, Eulerian–Lagrangian algorithms, such as UnTRIM and ELCIRC, have enjoyed considerable success recently in simulating 3D estuarine and coastal circulation. However, opportunities for improving the accuracy of this type of models were identified during extensive simulations of a tightly coupled estuary–plume–shelf system in the Columbia River system. Efforts to improve numerical accuracy resulted in SELFE, a new finite-element model for cross-scale ocean modeling. SELFE retains key benefits, including computational efficiency of existing semi-implicit Eulerian–Lagrangian finite-volume models, but relaxes restrictions on grids, uses higher-order shape functions for elevation, and enables superior flexibility in representing the bathymetry. Better representation of the bathymetry is enabled by a novel, “localized” vertical grid that resembles unstructured grids. At a particular horizontal location, SELFE uses either S coordinates or SZ coordinates, but the equations are consistently solved in Z space. SELFE also performs well relative to volume conservation and spurious oscillations, two problems that plague some finite-element models. This paper introduces SELFE as an open-source code available for community use and enhancement. The main focus here is on describing the formulation of the model and on showing results for a range of progressively demanding benchmark tests. While leaving details to separate publications, we also briefly illustrate the superior performance of SELFE over ELCIRC in a field application to the Columbia River estuary and plume.     

Method of Data Reduction and Uncertainty Estimation for Platinum-Group Element Data Using Inductively Coupled Plasma-Mass Spectrometry

A data reduction method is described for determining platinum-group element (PGE) abundances by inductively coupled plasma-mass spectrometry (ICP-MS) using external calibration or the method of standard addition. Gravimetric measurement of volumes, the analysis of reference materials and the use of procedural blanks were all used to minimise systematic errors. Internal standards were used to correct for instrument drift. A linear least squares regression model was used to calculate concentrations from drift-corrected counts per second (cps). Furthermore, mathematical manipulations also contribute to the uncertainty estimates of a procedure. Typical uncertainty estimate calculations for ICP-MS data manipulations involve: (1) Carrying standard deviations from the raw cps through the data reduction or (2) calculating a standard deviation from multiple final concentration calculations. It is demonstrated that method 2 may underestimate the uncertainty estimate of the calculated data. Methods 1 and 2 do not typically include an uncertainty estimate component from a regression model. As such models contribute to the uncertainty estimates affecting the calculated data, an uncertainty estimate component from the regression must be included in any final error calculations. Confidence intervals are used to account for uncertainty estimates from the regression model. These confidence intervals are simpler to calculate than uncertainty estimates from method 1, for example. The data reduction and uncertainty estimation method described here addresses problems of reporting PGE data from an article in the literature and addresses both precision and accuracy. The method can be applied to any analytical technique where drift corrections or regression models are used.     

巢湖形成演变与现代沉积作用

通过地质调查, 钻孔样品分析及历史资料对比, 得出:1.巢湖位于在中新生代形成的巢湖断陷盆地南部, 形成于更新世发育的河谷平原上, 由于全新世受气候及相应海平面变化和长江河道演变等三者共同影响而积水成湖。成湖历史大约2000年左右, 属河成型浅水湖。2.湖底沉积物平均中值粒径5-6ф, 并形成有三角洲沉积, 开阔湖沉积等五种沉积类型, 东、西、中三湖区的沉积物中分別为0.21, 0.37, 0.14cm/a, 全湖年淤积总量为93×l0⁴t/a。3.湖泊污染加重起始于70年代, 磷除来自外源外, 内负荷引起的释放量达220.56t/a, 这对巢湖富营养化亦有着重要的影响。     

贵州半边街锗锌矿床锗的富集特征及其地质意义

锗(Ge)属于战略性关键矿产资源, 是众多高新技术领域的重要原材料。Ge主要以共/伴生形式产出于Zn、Pb和Cu等热液贱金属矿床和煤矿床中, 目前全球Ge资源主要来源于铅锌矿的副产品。近期取得较大找矿突破的半边街锌矿床(7.70Mt矿石, Zn平均品位5.1%)伴生有超大型规模的金属Ge(>900t @110×10-6), 然而Ge超常富集的机制不清。本次研究利用LA-ICPMS对闪锌矿进行了微区原位微量元素组成和Mapping分析, 以揭示Ge等元素的富集特征、赋存状态和替代方式, 探讨Ge的富集机制。分析结果表明, 半边街矿床闪锌矿中高含量(平均值>100×10-6)微量元素有Fe、Pb、Ge和Cd等, 较高含量(平均值100×10-6~10×10-6)的有Mn、Ga、Sn和Tl, 以及较低含量(平均值< 10×10-6)的有Ag、Cu、In和Sb。相比于各类型铅锌矿床, 半边街闪锌矿Ge超常富集(274×10-6~1938×10-6, 平均为1055×10-6, 达到1000倍富集)。该矿床中Ge以类质同象或纳米颗粒形式赋存于闪锌矿中。微量元素摩尔量相关性图揭示Ge替代Zn的方式以Ge4++Pb2+↔3Zn2+为主, 其次为nGe4++Mn2+↔(2n+1)Zn2+和nGe4++2Tl3+↔(2n+3)Zn2+。综合分析认为, 半边街矿床Ge可能来源于相对富Ge的前寒纪基底; 富有机质的盆地流体在Ge的搬运中发挥重要作用; Ge、Cd等稀散元素的共生分异制约了Ge的最终超常富集。因此, 半边街矿床中Ge的超常富集是"源、运、聚"复杂过程多要素耦合作用的结果。     

渤海湾盆地南堡凹陷原油成因类型及其分布规律

南堡凹陷近年油气勘探取得了重大突破,为揭示该区滩海与陆地油气成因及分布规律,对67个原油和油砂样品进行了详细的地球化学研究。南堡凹陷原油具有低硫、中高蜡含量特征,陆相成因特点显著。共分为4种类型原油:(1)陆地沙三段深部层系原油,具有高4-甲基甾烷丰度、低奥利烷丰度与低甲基菲指数值和轻微偏高成熟度特征;(2)陆地Es1—Ed原油与Ng和Nm原油,其特征与第一类原油相反;上第三系原油普遍遭遇生物降解;(3)滩海东营组原油,以较高伽马蜡烷/C30藿烷和高甾烷/藿烷值区别于陆地原油;(4)滩海奥陶系潜山原油,具有低丰度甾萜类生物标志物、相对高丰度孕甾烷、三环萜烷与重排甾烷系列等特征,显示较高成熟度。上述4种成因类型原油指示南堡凹陷陆地、滩海均存在多层导有利烃源岩;同层位中滩海Es1-Ed存在较陆地更为优越的有利烃源岩,南堡滩海地区具有广泛的油气勘探前景。     

闽北近岸海域表层沉积物的风化特征及物源指示

通过对闽北近岸海域的229个表层沉积物样品和周边河流的54个表层沉积物样品进行粒度和常量元素测试,分析了研究区常量元素分布及粒度控制效应,并运用成分变异指数(ICV)和化学蚀变指数(CIA) 2种风化指标,结合(A-CN-K)-(A-CNK-FM)三角图,探讨了其风化特征及物源指示。研究表明,闽北近岸海域从岸向海表层沉积物粒径变粗,底质类型表现为明显的条带状,依次为泥-粉砂-砂质粉砂-砂。分析认为,形成这种条带状的主要原因是其物源及动力机制的不同,近岸主要为长江及周边河流物质输送,分布在50 m水深以浅,而在70 m以深的粗砂主要为残留沉积,50～70 m表现为过渡类型。风化程度研究表明,在空间分布上,从岸向海化学风化程度逐渐增强,70 m以深的异常区是以石英为主的残留沉积,抗风化能力强,黏土矿物含量低,表现为弱化学风化程度。长江、瓯江、闽江沉积物的风化程度表现为随着地理纬度的降低而增强,气候效应明显,研究区与长江沉积物风化趋势更为接近,所经历的形成环境最为相似,受物源效应影响。     

黑龙江省达连河油页岩地球化学特征及成因探讨

论述了达连河含煤段油页岩及油页岩段油页岩的地球化学特征,探讨了两套油页岩成因。这两套油页岩地球化学性质有很多相似之处,主要表现为油页岩都属于中熔或难熔的硅质灰分,SiO₂、Al₂O₃及Fe₂O₃含量很高;无机元素 Mn/Ti小于0.1,Sr/Ba为0.1～0.5,V/Ni为2.6～7.0;有机质类型以腐泥腐殖型为主,正构烷烃OEP值为1.55～3.67,具奇碳数优势,主峰碳为nC₂₃和C₂₉,Pr/Ph比值平均高达2.30,具强烈的姥鲛烷优势,Pr/nC₁₇比值为1.18～3.20,说明这两套油页岩都形成于淡水、近岸及弱氧化-弱还原的湖沼相沉积环境。但含煤段油页岩沉积时期高温、湿润的古气候条件下的喜湿植物面貌,使油页岩母质类型具有较多的富氢组分,这是造成两套油页岩品质差异的主要原因。     

新疆阿合勒－冲乎尔地区泥盆纪相动力学研究

阿合勒－冲乎尔火山沉积盆地是一个以火山岩并伴有陆源和内源沉积物的堆积场所，即该区泥盆纪地层是以事件作用的产物为其主要成分特点。显然，对其古地理的恢复、相的划分采用传统的相序分析方法是不合适的。作者认为要阐明该区泥盆纪的古地理面貌，必须在区域地层研究基础上，综合分析各种相标志和沉积物的特点，寻求它们的物理成因和地质成因解。在这一过程中，相的动力学分析显得尤为重要。本文在上述学术思想指导下，把本区泥盆纪的相分为成因相和环境相两大类，并对其中的微相进行了详细的阐述，从而为把沉积学的原理和方法应用了火山岩的搬运机制和堆积环境研究中提供了先期研究实例，同时为研究和恢复该区以至整个新疆北部泥盆纪的古地理格局打下了坚实的基础。     

围岩蚀变过程中地球化学组份质量迁移计算——以湖南沃溪Au-Sb-W矿床为例

在以往围岩蚀变的质量迁移研究基础上,本文提出了一种判定惰性元素的方法,同时结合等浓度线和加权最小二乘法,给出研究热液围岩蚀变中组分质量迁移的计算方法。在判定惰性元素时,更多地应用了线性回归分析中的方差分析与回归诊断参数与计算方法。这一系列的参数与方法被用于评述沃溪金锑钨矿床褪色化围岩蚀变质量迁移,取得了较好的效果。研究表明,在沃溪矿床围岩蚀变过程中,Ｔｉ表出现了一定的活动性。