Spatial distribution of soil heavy metal pollution estimated by different interpolation methods: Accuracy and uncertainty analysis

Mapping the spatial distribution of contaminants in soils is the basis of pollution evaluation and risk control. Interpolation methods are extensively applied in the mapping processes to estimate the heavy metal concentrations at unsampled sites. The performances of interpolation methods (inverse distance weighting, local polynomial, ordinary kriging and radial basis functions) were assessed and compared using the root mean square error for cross validation. The results indicated that all interpolation methods provided a high prediction accuracy of the mean concentration of soil heavy metals. However, the classic method based on percentages of polluted samples, gave a pollution area 23.54-41.92% larger than that estimated by interpolation methods. The difference in contaminated area estimation among the four methods reached 6.14%. According to the interpolation results, the spatial uncertainty of polluted areas was mainly located in three types of region: (a) the local maxima concentration region surrounded by low concentration (clean) sites, (b) the local minima concentration region surrounded with highly polluted samples; and (c) the boundaries of the contaminated areas. (C) 2010 Elsevier Ltd. All rights reserved.

土壤稳定入渗速率空间插值不确定性研究

以黄土丘陵沟壑第三副区的藉河流域为研究区,根据65个实测点数据,采用普通克里格法、反距离权重法、样条函数法等插值方法,分析了测点数量变化、栅格像元尺寸变化及插值方法的差异对土壤稳定入渗速率空间插值结果的影响,剖析了空间插值中的不确定性。结果表明:(1)参与插值站点越多,所得插值结果不确定性越小;(2)像元尺寸在25~800 m间变化对土壤稳定入渗速率的插值结果影响微弱;(3)不同插值方法对插值结果的精度影响较大,说明插值方法的差异对插值结果的不确定性有较大影响。

陕西省水蚀土壤因子指标插值方法比较研究

基于地统计学原理和半变异理论,采用地理信息系统技术定量研究了陕西省水蚀土壤因子指标的空间变异特征。分析比较了反距离权重法、样条函数法与普通克吕格法对陕西省土壤因子指标空间插值的精度。结果表明,研究区土壤抗冲系数具有中等强度的空间相关性,块金系数为32.29%,而稳渗速率、崩解速率、抗剪强度均表现为强烈的空间相关性,块金系数分别为13.19%,11.61%和12.98%。综合考虑平均相对误差、均方差及插值效果,认为普通克吕格法最好,更能反映土壤参数的空间特征并符合区域水土流失模型对数据的要求。对于普通克吕格法,稳渗速率的Lag步长为30 000 m,半方差理论模型为指数模型;抗冲系数、崩解速率、抗剪强度的Lag步长为55 000,半方差理论模型均为高斯模型。在空间分布上,各指标随土壤类型由北到南呈现明显的地带性规律。

区域水土流失研究中降雨量空间插值方法选取

在区域水土流失模型研究中,空间插值可提供每个计算栅格的气象要素资料。考虑到研究区域降雨与高程相关性很弱,不宜采用梯度距离反比法(GIDS),故采用距离反比法(IDW)和普通克里格法(Kriging),对延安示范区及其周围共50个站点2000—2003年的5—10月逐月降雨量进行插值。交叉验证结果表明:对2种插值方法,二者经过对数变换后平均相对误差(MRE)为8.30%和7.67%,分别比原始数据插值后的MRE下降了23.17%和23.50%,说明插值精度得到了提升,对研究区域某一年逐月降水的插值Kriging方法比IDW方法更加精确。

基于显式分形插值的高保真地震数据重建方法研究

Based on the fractal theories, contractive mapping principles as well as the fixed point theory, by means of affine transform, this dissertation develops a novel Explicit Fractal Interpolation Function（EFIF）which can be used to reconstruct the seismic data with high fidelity and precision. Spatial trace interpolation is one of the important issues in seismic data processing. Under the ideal circumstances, seismic data should be sampled with a uniform spatial coverage. However, practical constraints such as the complex surface conditions indicate that the sampling density may be sparse or for other reasons some traces may be lost. The wide spacing between receivers can result in sparse sampling along traverse lines, thus result in a spatial aliasing of short-wavelength features. Hence, the method of interpolation is of very importance. It not only needs to make the amplitude information obvious but the phase information, especially that of the point that the phase changes acutely. Many people put forward several interpolation methods, yet this dissertation focuses attention on a special class of fractal interpolation function, referred to as explicit fractal interpolation function to improve the accuracy of the interpolation reconstruction and to make the local information obvious. The traditional fractal interpolation method mainly based on the randomly Fractional Brown Motion (FBM) model, furthermore, the vertical scaling factor which plays a critical role in the implementation of fractal interpolation is assigned the same value during the whole interpolating process, so it can not make the local information obvious. In addition, the maximal defect of the traditional fractal interpolation method is that it cannot obtain the function values on each interpolating nodes, thereby it cannot analyze the node error quantitatively and cannot evaluate the feasibility of this method. Detailed discussions about the applications of fractal interpolation in seismology have not been given by the pioneers, let alone the interpolating processing of the single trace seismogram. On the basis of the previous work and fractal theory this dissertation discusses the fractal interpolation thoroughly and the stability of this special kind of interpolating function is discussed, at the same time the explicit presentation of the vertical scaling factor which controls the precision of the interpolation has been proposed. This novel method develops the traditional fractal interpolation method and converts the fractal interpolation with random algorithms into the interpolation with determined algorithms. The data structure of binary tree method has been applied during the process of interpolation, and it avoids the process of iteration that is inevitable in traditional fractal interpolation and improves the computation efficiency. To illustrate the validity of the novel method, this dissertation develops several theoretical models and synthesizes the common shot gathers and seismograms and reconstructs the traces that were erased from the initial section using the explicit fractal interpolation method. In order to compare the differences between the theoretical traces that were erased in the initial section and the resulting traces after reconstruction on waveform and amplitudes quantitatively, each missing traces are reconstructed and the residuals are analyzed. The numerical experiments demonstrate that the novel fractal interpolation method is not only applicable to reconstruct the seismograms with small offset but to the seismograms with large offset. The seismograms reconstructed by explicit fractal interpolation method resemble the original ones well. The waveform of the missing traces could be estimated very well and also the amplitudes of the interpolated traces are a good approximation of the original ones. The high precision and computational efficiency of the explicit fractal interpolation make it a useful tool to reconstruct the seismic data; it can not only make the local information obvious but preserve the overall characteristics of the object investigated. To illustrate the influence of the explicit fractal interpolation method to the accuracy of the imaging of the structure in the earth’s interior, this dissertation applies the method mentioned above to the reverse-time migration. The imaging sections obtained by using the fractal interpolated reflected data resemble the original ones very well. The numerical experiments demonstrate that even with the sparse sampling we can still obtain the high accurate imaging of the earth’s interior’s structure by means of the explicit fractal interpolation method. So we can obtain the imaging results of the earth’s interior with fine quality by using relatively small number of seismic stations. With the fractal interpolation method we will improve the efficiency and the accuracy of the reverse-time migration under economic conditions. To verify the application effect to real data of the method presented in this paper, we tested the method by using the real data provided by the Broadband Seismic Array Laboratory, IGGCAS. The results demonstrate that the accuracy of explicit fractal interpolation is still very high even with the real data with large epicenter and large offset. The amplitudes and the phase of the reconstructed station data resemble the original ones that were erased in the initial section very well. Altogether, the novel fractal interpolation function provides a new and useful tool to reconstruct the seismic data with high precision and efficiency, and presents an alternative to image the deep structure of the earth accurately.

Determination of Point of Zero Charge of Tunisian Kaolinites by Potentiometric and Mass Titration Methods

This paper deals with determining points of zero charge of natural and Na+-saturated mineral kaolinites using two methods: (1) acid-base potentiometric titration was employed to obtain the adsorption of H+ and OH- on amphoteric surfaces in solutions of varying ionic strengths in order to determinate graphically the point of zero net proton charge (PZNPC) defined equally as point of zero salt effect (PZSE); (2) mass titration curve at different electrolyte concentrations in order to estimate PZNPCs by interpolation and to compare with those determined by potentiometric titrations. The two methods involved points of zero charge approximately similar for the two kaolinites between 6.5-7.8, comparable to those reported previously and were in the range expected for these clay minerals. The comparison of potentiometric surface titration curves obtained at 25 °C and those published in the literature reveals significant discrepancies both in the shape and in the pH of PZNPCs values.

Evaluation of the Effectiveness of Representative Methods for Determining Young's Modulus and Hardness from Instrumented Indentation Data

The effectiveness of Oliver & Pharr's (O&P's) method, Cheng & Cheng's (C&C's) method, and a new method developed by our group for estimating Young's modulus and hardness based on instrumented indentation was evaluated for the case of yield stress to reduced Young's modulus ratio (sigma(y)/E-r) >= 4.55 x 10(-4) and hardening coefficient (n) <= 0.45. Dimensional theorem and finite element simulations were applied to produce reference results for this purpose. Both O&P's and C&C's methods overestimated the Young's modulus under some conditions, whereas the error can be controlled within +/- 16% if the formulation was modified with appropriate correction functions. Similar modification was not introduced to our method for determining Young's modulus, while the maximum error of results was around +/- 13%. The errors of hardness values obtained from all the three methods could be even larger and were irreducible with any correction scheme. It is therefore suggested that when hardness values of different materials are concerned, relative comparison of the data obtained from a single standard measurement technique would be more practically useful. It is noted that the ranges of error derived from the analysis could be different if different ranges of material parameters sigma(y)/E-r and n are considered.

Influence of Electrostatic Interaction on Fibrinogen Adsorption on Gold Studied by Imaging Ellipsometry Combined with Electrochemical Methods

Imaging ellipsometry was combined with electrochemical methods for studying electrostatic interactions of protein and solid surfaces. The potential of zero charge for gold-coated silicon wafer/solution interfaces wad determined by AC impedance method. The potential of the gold-coated silicon wafer was controlled at the potential of zero charge, and the adsorption of fibrinogen on the potential-controlled and non-controlled surfaces was measured in real time at the same time by imaging ellipsometry The effect of electrostatic interaction was studied by comparing the difference between the potential of controlled adsorption and the Potential of noncontrolled adsorption. It was shown that the rate of fibrinogen adsorption on the potentiostatic surface was faster than that on the nonpotentiostatic surface. The electrostatic influence on fibrinogen adsorption on the gold-coated silicon wafer was weak, so the hydrophobic interaction should be the major affinity.

The Optical Methods on Measuring Surface Deformation and Surface Wave in The Thermal Capillary Convection

An optical diagnostic system consisting of the Michelson interferometer with the image processor has been developed for the study of the kinetics of the thermal capillary convection. The capillary convection, surface deformation, surface wave and the velocity field in a rectangular cavity with different temperature's sidewalls have been investigated by optical interference method and PIV technique. In order to calculate the surface deformation from the interference fringe, Fourier transformation is used to grating analysis. The quantitative results of the surface deformation and surface wave have been calculated from the interference fringe pattern.

Assessment of Microscale Test Methods of Peeling and Splitting along Surface of Thin-Film/Substrates

Peel test methods are assessed through being applied to a peeling analysis of the ductile film/ceramic substrate system. Through computing the fracture work of the system using the either beam bend model (BB model) or the general plane analysis model (GPA model), surprisingly, a big difference between both model results is found. Although the BB model can capture the plastic dissipation phenomenon for the ductile film case as the GPA model can, it is much sensitive to the choice of the peeling criterion parameters, and it overestimates the plastic bending effect unable to capture crack tip constraint plasticity. In view of the difficulty of measuring interfacial toughness using peel test method when film is the ductile material, a new test method, split test, is recommended and analyzed using the GPA model. The prediction is applied to a wedge-loaded experiment for Al-alloy double-cantilever beam in literature.

On the Gas-Liquid Two-Phase Flow Regimes and their Recognition Methods

气液两相流体系是一个复杂的多变量随机过程体系,流型的定义、流型过渡准则和判别方法等方面的研究是多相流学科目前研究的重点内容。本文就与气液两相流流型及其判别有关的研究状况进行了回顾和评述,力图反映近年来气液两相流流型及其判别问题研究的状态和趋势。

Practical Criterion for Universal Laws and Practical Methods for Construction of Co-variance Sets

给出相对论力学中普遍定律的实用判别法和协变集的实用构造法,还给出实现非普遍定律的“可导出性”的一种实用方法．

Methods of Obtaining Instantaneous Modulus of Viscoelastic Solids Using Displacement-Controlled Instrumented Indentation with Axisymmetric Indenters of Arbitrary Smooth Profiles

We derive a relationship between the initial unloading slope, contact depth, and the instantaneous relaxation modulus for displacement-controlled indentation in linear viscoelastic solids by a rigid indenter with an arbitrary axisymmetric smooth profile. While the same expression is well known for indentation in elastic and in elastic–plastic solids, we show that it is also true for indentation in linear viscoelastic solids, provided that the unloading rate is sufficiently fast. When the unloading rate is slow, a “hold” period between loading and unloading can be used to provide a correction term for the initial unloading slope equation. Finite element calculations are used to illustrate the methods of fast unloading and “hold-at-the-maximum-indenter-displacement” for determining the instantaneous modulus using spherical indenters.

Molecular/Cluster Statistical Thermodynamics Methods To Simulate Quasi-Static Deformations At Finite Temperature

The rapid evolution of nanotechnology appeals for the understanding of global response of nanoscale systems based on atomic interactions, hence necessitates novel, sophisticated, and physically based approaches to bridge the gaps between various length and time scales. In this paper, we propose a group of statistical thermodynamics methods for the simulations of nanoscale systems under quasi-static loading at finite temperature, that is, molecular statistical thermodynamics (MST) method, cluster statistical thermodynamics (CST) method, and the hybrid molecular/cluster statistical thermodynamics (HMCST) method. These methods, by treating atoms as oscillators and particles simultaneously, as well as clusters, comprise different spatial and temporal scales in a unified framework. One appealing feature of these methods is their "seamlessness" or consistency in the same underlying atomistic model in all regions consisting of atoms and clusters, and hence can avoid the ghost force in the simulation. On the other hand, compared with conventional MD simulations, their high computational efficiency appears very attractive, as manifested by the simulations of uniaxial compression and nanoindenation. (C) 2008 Elsevier Ltd. All rights reserved.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.