Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

The Origin and Evolution of Silica-saturated Alkalic Suites: an Experimental Study

Experimental simulation of incremental crystal fractionationof a hy-normative hawaiite indicates that the spectrum of compositionsfrom mildly alkalic hawaiite to sodic rhyolite found in silica-saturatedalkalic suites of the ocean islands and continental hotspotsand rifts can be produced by fractionation at 9·3 kbarwith bulk water contents >     

基于XML的WebGIS体系结构

针对WebGIS中存在的地理空间数据语义共享问题，在分析XML和GML特点的基础上，论述了基于XML构架WebGIS系统的方案和应用方式，并详析了实现构架的相关技术。     

基于分层的面状地图符号设计

根据面状符号的结构特征和实现机制 ,提出了基于分层的面状地图符号设计方法 ,把面状符号分为若干层 ,每层采用一种填充方式 ,并在基于VC 的Cartosymbol软件中实现。结果表明 ,该方法既便于地理现象的表达 ,又便于计算机实现。     

视电阻率观察中一种抗介质干扰的尝试

埋地金属管线对视电阻率的观测危害甚大。造成观测值大幅改变，跳动，无规律可循。严重时，甚至台站报废，被迫搬迁。通过一些突发事件的启发，结合现场实验，并引证了一些理论计算，笔者摸索出了只须对金属管道作局部微小的改造，就能基本克服其对观测的干扰的方法。该方法在红格台试验半年多来，取得了明显的效果。观测值恢复到安装的水管以前的水平，再无大幅度的跳动。该方法节约经费，简单易行，大大提高了台站在复杂环境中的生存能力。     

Tetrahedral Plots of the Phase Relations for Basalts

The phase relations of quaternary systems are generally represented by projections onto ternary compositional planes. Such projections often obscure relationships that would only be evident in a three-dimensional tetrahedral plot. The tetrahedral plot requires that compositions of the minerals and melts be transformed into Cartesian coordinates. It is shown here how this transformation is carried out. The application is demonstrated by tetrahedral plots of experimental melt compositions of partially molten lherzolite. Furthermore, the plot can be used to evaluate whether or not a particular basaltic composition represents a primary melt. The methods are applicable to any four-component system.