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.     

北京连续降雪过程分析

该文对2002年12月份北京出现的1841年以来历史上最长的连续6天的降雪天气进行了合成分析和诊断分析。北半球合成环流形势显示，连续降雪发生在东亚稳定的纬向环流形势下，其上游地区强大经向环流和阻塞高压使中纬度东亚和西太平洋纬向环流得以维持。东亚合成环流的垂直结构表明，连续降雪期间华北地区始终处于锋区中，并存在低空辐合、高空辐散的散度场垂直分布结构。卫星水汽云图显示出有水汽通道向河套倒槽云系输送水汽。华北地区的诊断分析表明，地面和边界层中山东地区分裂高压南侧向北回流的偏南气流是造成降雪的主要水汽通道。变形场的流场结构使水汽在北京附近辐合。垂直剖面展现了华北回流降雪天气近地面层的浅薄的冷空气垫，和暖湿空气回流在冷空气垫上的爬升，并在其上形成一个浅薄的饱和层。     

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.     

Spatially controlled Fe and Si isotope variations: an alternative view on the formation of the Torres del Paine pluton

We present new Fe and Si isotope ratio data for the Torres del Paine igneous complex in southern Chile. The multi-composition pluton consists of an approximately 1 km vertical exposure of homogenous granite overlying a contemporaneous 250-m-thick mafic gabbro suite. This first-of-its-kind spatially dependent Fe and Si isotope investigation of a convergent margin-related pluton aims to understand the nature of granite and silicic igneous rock formation. Results collected by MC-ICP-MS show a trend of increasing δ⁵⁶Fe and δ³⁰Si with increasing silica content as well as a systematic increase in δ⁵⁶Fe away from the mafic base of the pluton. The marginal Torres del Paine granites have heavier Fe isotope signatures (δ⁵⁶Fe = +0.25 ± 0.02 2se) compared to granites found in the interior pluton (δ⁵⁶Fe = +0.17 ± 0.02 2se). Cerro Toro country rock values are isotopically light in both Fe and Si isotopic systems (δ⁵⁶Fe = +0.05 ± 0.02 ‰; δ³⁰Si = ?0.38 ± 0.07 ‰). The variations in the Fe and Si isotopic data cannot be accounted for by local assimilation of the wall rocks, in situ fractional crystallization, late-stage fluid exsolution or some combination of these processes. Instead, we conclude that thermal diffusion or source magma variation is the most likely process producing Fe isotope ratio variations in the Torres del Paine pluton.