Interpolation by regularized spline with tension: I. Theory and implementation

Bivariate and trivariate functions for interpolation from scattered data are derived. They are constructed by explicit minimization of a general smoothness functional, and they include a tension parameter that controls the character of the interpolation function (e.g., for bivariate case the surface can be tuned from a membrane to a thin steel plate), Tension can be applied also in a chosen direction, for modeling of phenomena with a simple type of anisotropy. The functions have regular derivatives of all orders everywhere. This makes them suitable for analysis of surface geometry and for direct application in models where derivatives are necessary. For processing of large datasets (thousands of data points), which are now common in geosciences, a segmentation algorithm with a flexible size of overlapping neighborhood is presented. Simple examples demonstrating flexibility and accuracy of the functions are presented.On leave from the Department of Physical Geography and Cartography, Comenius University, Mlynská dolina, Bratislava, Czechoslovakia.On leave from the Institute of Physics, Dúbravská cesta 9, Bratislava, Czechoslovakia.     

南极海冰、冰穴和冰川冰及其对水团形成和变性的作用

本文利用前人的成果及笔者1992/1993年的南极海冰观测和收集的资料以及水文观测资料数据阐述了南极海冰的特性,特别是南极海冰过程、冰穴以及冰川冰对南极水团(南极表层水、南极底层水、南极陆架水、南极中层水以及南极冰架水)的形成和变性所起的特殊作用。 南极海冰覆盖面积的年际变化,夏季最大年份是最小年份的2倍多,冬季年间变化较小,最大仅为20％;但其季节变化非常大,冬季平均覆盖面积通常是夏季的5倍。南极海冰对大气-海洋间相互作用有重大影响,特别是深海洋区中冬季的结冰和发育造成的垂向对流、夏季的融化是形成南极表层水(含南极冬季水和南极夏季表层水),进而形成南极中层水的主要原因;南极陆架区的的海冰兴衰过程是形成南极陆架水的直接原因,它与变性南极绕极深层水混合并受到冰川冰的进一步冷却作用,成为形成南极底层水的主要水团;南极冰架底部的冷却、融化和冰架以下水体的结冰作用形成的高盐对流过程产生的南极冰架水,亦是形成南极底层水的贡献者。 冰穴是70年代以来卫星观测的重大发现。对其形成和对大气、海洋的影响作用尚不完全清楚,初步的研究成果表明,冰穴中产生的热盐对流对南极水团的形成、变性、大洋深层的翻转以及海洋-大气间的热量传输和气体交换起有非常重要的作用。     

相律的应用:蒸发结晶过程自由度

在连续原理和浓度限制条件的基础上,分析了卤水蒸发过程的自由度,提出了“蒸发结晶过程自由度”这一概念。发现蒸发结晶过程自由度在整个过程中不是1便是0,在1和0之间交替变化,与选择的组分数和析出固相数目无关。蒸发结晶过程自由度为1,物理意义是在不引起新相产生旧相不消失的前提下过程连续地变化(水分的连续蒸失),几何意义是指在适当的坐标系中的一条线,一条直线或曲线;0的物理意义是指新相的产生和旧相的即将消失或者前一段的连续变化过渡到后一段的连续变化,几何意义是指一点,两蒸发阶段的交点或者过程的终点。简单地讨论了过程的单向性和连续件。     

张广才岭南部杨木岗组的厘定及物源分析

本文通过岩石组合特征和区域对比，将张广才岭南部西蛤拉河子—大锅盔一带分布的浅变质地层重新厘定为杨木岗组。为了确定杨木岗组的形成时代和沉积物源，进行了碎屑锆石U- Pb年代学和微体古生物地层学研究。锆石大多数呈自形—半自形晶，显示典型振荡岩浆生长环带，暗示其岩浆成因。该地层中测得的两组碎屑锆石U- Pb产生多组谐和年龄，其中PM010- TW样品56个测点最小峰值(谐和)年龄为307 Ma，DB02- TW样品51个测点最小峰值(谐和)年龄为275 Ma；覆盖在杨木岗组之上的中生代二浪河组安山岩的定年结果为181. 1±0. 9 Ma，表明杨木岗组形成于早二叠世晚期。杨木岗组中获取疑源类化石组合出现了新元古代晚期—早寒武世和奥陶纪地层常见分子，结合碎屑锆石年龄结果，反映杨木岗组沉积时周围存在早古生代和中—新元古代地质体。碎屑沉积岩Al2O3/TiO2平均值为24. 44，稀土元素球粒陨石标准化曲线具有轻稀土富集、重稀土稳定和负Eu异常特征，结合碎屑锆石的年龄频数可以看出，确定杨木岗组的沉积物主要来源于沉积盆地周围的晚古生代早期中酸性火成岩，次要物源由沉积盆地周边的早古生代地质体和近地表的中—新元古代地质体提供。佳木斯地块为松嫩- 张广才岭地块上晚古生代地层的形成提供了部分物源，暗示佳木斯地块与松嫩- 张广才岭地块于早二叠世之前已完成拼合。     

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

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

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

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

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