胶东海阳所超镁铁岩中强亲铁元素的地球化学特征及其构造环境探讨

强亲铁元素与亲石元素具有不同的地球化学行为，因此能够从不同的角度为造山带中超镁铁岩的成因及演化提供重要信息。位于苏鲁造山带东北端的胶东海阳所超镁铁岩主要由橄榄岩和辉石岩组成，它们常以团块状赋存于花岗质片麻岩中。虽然前人对这些超镁铁岩已经开展大量岩石学研究，但关于其成因及构造属性仍存在较大争议。本文开展了海阳所超镁铁岩的全岩主微量元素、强亲铁元素及Re-Os同位素的分析工作，结果显示蛇纹石化橄榄岩具有较高的MgO和Fe₂O₃^T含量，较低的Al₂O₃、TiO₂和CaO含量，明显富集流体迁移元素（U、Pb），亏损高场强元素（Zr、Hf），强亲铁元素没有发生明显分异，但Ru显示正异常，表明海阳所蛇纹石化橄榄岩是经历了低-中等程度部分熔融及熔/流体交代作用影响的残余地幔橄榄岩。海阳所辉石岩的主量元素表现出明显的结晶分异特征，稀土元素较原始地幔富集，铂族元素（PGEs）含量较低且发生了明显的分异，表明辉石岩的地幔源区经历过高程度的部分熔融和硫化物的分离。海阳所蛇纹石化橄榄岩的Os同位素地球化学特征表现出大洋亲和性，与辉石岩不具有熔体-残留体的关系。由于该地区发育较深层次的韧性剪切带，蛇纹石化橄榄岩中的橄榄石与辉石表现出韧性变形的特征，同时有辉石岩侵入到橄榄岩的现象，表明该地区的蛇纹石化地幔橄榄岩与辉石岩既不同时，也不同源，因此，暗示了该套岩石组合可能形成于大洋核杂岩（OCC）与洋脊型蛇绿岩（MOR）堆晶岩交互发育环境。     

北大别饶拔寨石榴辉石岩P-T-t轨迹及其构造意义

P-T-t轨迹作为变质岩的重要研究方法之一，对揭示岩石的构造演化过程具有重要意义。北大别饶拔寨镁铁-超镁铁岩形成的构造环境和就位过程长期以来尚存争议。本文通过岩相学观察、矿物化学研究和温压计计算，揭示出饶拔寨石榴辉石岩经历了四个变质演化阶段：1）超高压变质阶段（M1）：主要根据石榴子石中金红石的出溶，单斜辉石中石英的出溶和磷灰石中不透明矿物的出溶，结合前人研究，认为饶拔寨石榴辉石岩经历过超高压变质阶段（≥2.5GPa）；2）高压麻粒岩相阶段（M2）：矿物组合为石榴子石（变斑晶）和单斜辉石（基质）+斜长石（基质），记录的温压条件为T=648~700℃，P=1.47~1.94GPa；3）中压麻粒岩相阶段（M3）：以石榴子石外围发育的主要由斜方辉石+斜长石组成的内圈"白眼圈"为特征，形成的温压条件为T=781~796℃，P=0.92~0.98GPa；4）角闪岩相阶段（M4）：以石榴子石变斑晶周围发育的外圈"白眼圈"为特征，其矿物组合为角闪石+斜长石的后成合晶，形成的温压条件为T=663~685℃，P=0.50~0.58GPa。石榴辉石岩的锆石SIMS U-Pb定年得到了3组不同的交点年龄，分别为208.1~202.1Ma、227.6Ma和817.7Ma。根据锆石包裹体中角闪石+斜长石的组合，推断208.1~202.1Ma代表了角闪岩相变质年龄的下限，227.6Ma则是高压麻粒岩相的下限年龄，而新元古代年龄（817.7Ma）与区域上的变质岩原岩年龄一致，可能代表了其原岩年龄。结合前人研究，饶拔寨石榴辉石岩记录了顺时针的P-T-t轨迹，揭示了板片俯冲（超高压变质）-碰撞-折返（降压升温过程，~227Ma）-抬升（降压降温过程，208~202Ma）的完整过程。     

利用主成分分析改善土地利用变化的遥感监测精度——以珠江三角洲城市用地扩张为例

近年来,珠江三角洲由于经济的快速发展,城市用地急剧增加,利用多时相的遥感图,可以定量地监测这种城市化的现象。但理,由一般的遥感动态监测方法所得的结果往往夸大变化的程度,以及获得一些不合理的结论．该文提出主成分分析的方法来改善遥感动态监测的精度。将该方法应用应用于珠江三角洲发展最快的东莞市,获得了较满意的结果。     

北京云蒙山地区花岗岩穹隆及伸展构造的探讨

北京云蒙山花岗岩为一中生代侵入的花岗岩穹隆,花岗岩穹隆的叶理普遍发育,叶理轨迹基本平行于穹隆的外部边界,并显示出从核部到边部逐渐增强,东南侧明显强于西北侧的特点。变形构造研究显示,花岗岩穹隆的边部及围岩中普遍存在不同层次及不同运动方向的伸展构造。东南侧以具河防口－水峪伸展型韧性剪切带为特征,剪切运动标志显示为从ＮＷ－ＳＥ的正剪切运动,有限应变分析估算其剪切位移量在１０ｋｍ以上,剪切带上部被河防口正     

Research on temporal and spatial distribution, evolutionary character and mechanism of crustal deformation field before and after the Tangshan earthquake

Ｒｅｓｅａｒｃｈｏｎｔｅｍｐｏｒａｌａｎｄｓｐａｔｉａｌｄｉｓｔｒｉｂｕ┐ｔｉｏｎ,ｅｖｏｌｕｔｉｏｎａｒｙｃｈａｒａｃｔｅｒａｎｄｍｅｃｈａｎｉｓｍｏｆｃｒｕｓｔａｌｄｅｆｏｒｍａｔｉｏｎｆｉｅｌｄｂｅｆｏｒｅａｎｄａｆｔｅｒｔｈｅＴａｎｇｓｈａｎ...     

Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.     

Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.     

Regional Seas and Sustainable Development Results of an International Workshop

Reports

Regional Seas and Sustainable Development Results of an International Workshop     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

合肥遥测台网地震速报系统

论述了适用于遥测地震台网大震速报的一套计算机处理系统,系统设计以《全国遥测台网观测技术规范》和《近震分析》理论为依据,实现了数据输入,震相选择,数据处理,地震定位,震级计算,结果输出及贮存等多项功能,达到了快速确定地震三要素,提高速报速度和精度之目的,是一个较完善的,实用的速报系统。