赣南地区石雷石英闪长岩的成因:岩石化学、副矿物微量元素、锆石U-Pb年代学与Sr-Nd-Hf同位素制约

石雷石英闪长岩是赣南崇-余-犹地区比较特殊的闪长质侵入体。锆石的原位U-Pb定年表明,该岩体侵位于433.5±3.4Ma。全岩主量元素特征上显示出中偏酸性(SiO2=56.92%～64.70%),富Al(Al2O3=14.10%～14.83%),富碱(Alk=6.41%～7.40%)特别是富钾(K2O=3.86%～4.85%),镁、铁含量较高,MgO:3.47%～5.95%,FeOT:5.23%～8.14%以及低磷(P2O5=0.27%～0.4%)的特点;微量元素上主要富集K、Rb、Cs等大离子亲石元素和轻稀土元素,亏损Nb、Ta、Ti、P等高场强元素。磷灰石微量元素特征上显示高度富集稀土元素特别是轻稀土元素的特征;具有Eu的负异常(δEu=0.37～0.45)。ISr位于0.7073～0.7132之间,εNd(t)变化于-8.41～-4.97之间,两阶段钕模式年龄介于1.58～1.86Ga之间,Hf同位素组成相对均一,εHf(t)主要集中变化于-8～-2之间,两阶段Hf模式年龄加权平均为1.77±0.09Ga,这些特征都暗示了该石英闪长质岩体的形成是强烈壳幔相互作用的产物,区内加里东晚期可能发生了局部的岩石圈的减薄。     

Numerical mass balance for soil-moisture transport problems

The Galerkin finite element method coupled with the Crank-Nicolson time advance procedure is often used as a numerical analog for unsaturated soil-moisture transport problems. The Crank-Nicolson procedure leads to numerical mass balance problems which results in instability. A new temporal and spatial integration procedure is proposed that exactly satisfies mass balance for the approximating function used. This is accomplished by fitting polynomials continuously throughout the time and space domain and integrating the governing differential equations. To reduce computational effort, the resulting higher order polynomials are reduced to quadratic and linear piece-wise continuous polynomial approximation functions analogous to the finite element approach. Results indicate a substantial improvement in accuracy over the combined Galerkin and Crank-Nicolson methods when comparing to simplified problems where analytical solutions are available.     

A note on time integration of unsaturated soil-moisture transport

The non-linear soil-moisture diffusivity model can be approximately linearized by using values of diffusivity assumed constant for small intervals of space and time. By a series expansion of the diffusivity function and integrating the resulting series of differential equations with respect to time, an improved numerical model is developed. Results from application of this new approach to a sharp wetting-front soil infiltration problem indicates that a 67% saving in numerical effort is achieved at comparable estimation accuracy levels when using the traditional finite timestep Crank-Nicolson approach.     

Nodal domain integration model of one-dimensional advection—diffusion

The nodal domain integration method is applied to a one-dimensional advection—diffusion mathematical model without a source term. Comparison of the resulting numerical model to the well known Galerkin finite element, subdomain, and finite difference domain models indicates that a single numerical statement can be developed which includes the Galerkin finite element, subdomain, and finite difference models as special cases.     

Mass lumping models of the linear diffusion equation

The nodal domain integration method is used to develop a numerical model of the linear diffusion equation. The nodal domain integration approach is shown to represent an infinity of finite element mass matrix lumping schemes including the Galerkin and subdomain integration versions of the weighted residual method and an integrated finite difference method. Neumann, Dirichlet and mixed boundary conditions are accommodated analogous to the Galerkin finite element method. In order to reduce the overall integrated approximation relative error, a mass matrix lumping formulation is developed which is based on the Crank-Nicolson time advancement approximation. The optimum mass lumping factors are found to be strongly related to the model timestep size.     

西秦岭勉略带陆内构造变形研究

秦岭造山带勉略缝合带是古特提斯洋盆向北俯冲形成的华北与华南最后拼接带。这个主缝合带俯冲-碰撞过程中以由北向南的一系列韧性逆冲推覆构造为特征,形成由前泥盆系、泥盆-石炭系和蛇绿混杂岩等不同构造岩片叠置的复杂构造带,碰撞时代从245Ma一直延续到230Ma左右。最近,作者对勉略缝合带内发育的韧性和脆性左行走滑剪切变形进行了研究,结果表明这些顺造山带的左行韧性走滑剪切变形带的变形时代为223±2Ma,与碰撞后花岗岩所确定的碰撞后构造环境的起始时间(225Ma)一致,显示这些韧性走滑剪切变形带是勉略带陆内变形初期变形产物。亦即华北、扬子大陆碰撞之后很快就转入陆内变形阶段,并且是以顺造山带的侧向走滑位移为主要变形方式。勉略带内顺造山带的脆性左行走滑断层的发育,表明这种顺造山带的侧向位移过程从深部到地壳浅层是一致的。因此,大陆碰撞在直接碰撞之后很快转变为顺造山带的侧向走滑位移为主的陆内变形,这种位移可能表现为两个大陆碰撞后的相对走滑,或是碰撞带中强烈变形部分顺造山带的侧向挤出,从而消减了正向碰撞所造成的地壳缩短和增厚。     

A two-dimensional dam-break flood plain model

A simple two-dimensional dam-break model is developed for flood plain study purposes. Both a finite difference grid and an irregular triangle element integrated finite difference formulation are presented. The governing flow equations are approximately solved as a diffusion model coupled to the equation of continuity. Application of the model to a hypothetical dam-break study indicates that the approach can be used to predict a two-dimensional dam-break flood plain over a broad, flat plain more accurately than a one-dimensional model, especially when the flow can break-out of the main channel and then return to the channel at other downstream reaches.     

西藏丁青弧前蛇绿岩的地球化学特征

丁青蛇绿岩位于班化湖－丁青－怒江蛇绿岩带的东段，其地幔岩出露规模是该带中最大的。本文报道的丁青蛇绿岩主要由地幔橄榄岩、堆晶岩、辉长岩和斜长花岗岩组成。蛇绿岩剖面上覆硅岩中的放射虫化石是早株罗世和晚三叠世诺利克期的，中侏罗统砂岩和砾岩不整合覆盖在蛇绿岩之上，由此确定丁青蛇绿岩是晚三叠－早侏罗世的，在中株罗世之前侵位，丁青蛇绿岩属于玻安岩系，玻安岩的特点是富Si、Mg和大离子亲石元素（LILE），贫高场强元素（Ti、P、Zr、Y、Yb和Nb）。丁青蛇绿岩的堆晶岩、辉长岩和辉绿岩均具“U”型REE分布，暗示丁青玻安岩是由于亏损的地幔源岩和来自消减带的水和流体两组分的混合形成的。丁青玻安岩的地球化学特征类似西太平洋第三纪玻安岩，而明显不同于MORB的地球化学性质， 表明丁青玻安岩应当形成于洋内岛弧的弧前环境，属于弧前蛇绿岩。