Spatial analysis of elements in soils from Chengdu City     

Shuyun XIE Zhengyu BAO Bing LI 《中国地球化学学报》2006,25(B08):7-8

Three kinds of spatial analysis methods （geostatistics, concentration-area fractal model and the multifractal analysis called the moment method） were used for almost 50 elements, including heavy metals, disperse elements, rare elements and even others, in 6586 top soil （0-20 cm） samples and 1833 deep soil （150-200 cm） samples from Chengdu metropolitan area of 12400 km^2, southwestern China. The ranges of spatial correlation revealed by variograms are quite different for different kinds of elements in the top and deep soils. The most interest is the fact that the multifractal spectra of environmentally important elements such as Pb, Cr, Cd and Ni in top soils in the metropolitan area show systematic change from those in the deep soils, revealing a strong anthropogenic addition, while Hg, Zn, As, Cu and all common elements show no such kind of addition. In terms of multifractal properties based on the multifractal spectrum curves, those disperse and rare elements show great deviation from other major and trace elements, which is also of great interest.  相似文献   

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

C. R. Dietrich  G. N. Newsam 《Mathematical Geology》1995,27(2):207-228

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.  相似文献   

Hierarchical decomposition of variance with applications in environmental mapping based on satellite images     

Ferenc Csillag  Sándor Kabos 《Mathematical Geology》1996,28(4):385-405

A quadtree-based image segmentation procedure (HQ) is presented to map complex environmental conditions. It applies a hierarchical nested analysis of variance within the framework of multiresolution wavelet approximation. The procedure leads to an optimal solution for determining mapping units based on spatial variability with constraints on the arrangement and shape of the units. Linkages to geostatisiics are pointed out, but the HQ decomposition algorithm does not require any homogeneity criteria. The computer implementation can be parameterized by either the number of required mapping units or the maximum within-unit variance, or it can provide a spectrum of significances of nested ANOVA. The detailed mathematical background and methodology is illustrated by a salt-affected grassland mapping study (Hortobágy, Hungary), where heterogeneous environmental characteristics have been sampled and predicted based on remotely sensed images using these principles.  相似文献   

Geographical Information Systems and Geostatistics for Modelling Radioactively Contaminated Land Areas   总被引：2，自引：0，他引：2  

NAZZARENO?DIODATOEmail author  MICHELE?CECCARELLI 《Natural Hazards》2005,35(2):229-242

The Chernobyl plume contaminated vast lands of Europe with radiocaesium (¹37Cs) in 1986 because of the deposition of radionuclides on the ground by wet and dry deposition processes. Nevertheless, in a nuclear emergency, contamination data may be very sparse and there is need to make rapid and scientifically supported decisions. Here we analyze the rainfall field, an important precursor of the wet deposition, during the passage of the plume. Thus, estimating rainfall spatial variability can help to identify possible contaminated areas and associated risks when rainfall exceeded a given threshold. In this paper, we show that the conditional probabilities of exceeding threshold rainfall values could be spatially assessed using the mutual benefits of linking geostatistical and geographical information system (GIS) to quantify the evaluation of the risk involved in decision making. In particular, the non-parametric geostatistic technique, termed Indicator Kriging (IK), enables one to efficiently estimate the probability that the true value exceeds the threshold values by means of the indicator coding transform. Afterward, GIS has been used to find the areas probably affected by wash-out (probability >0.5 that rainfall is above a certain threshold). The experimental study has been focused on a test site in Beneventan agroecosystem (Southern Italy) to model the spatial uncertainty over a continuous area from sparse rainfall data. This enabled to generate probability maps delineating area potentially affected by to contamination to be monitored after wet deposition of Chernobyl releases.  相似文献   

贵州百花湖表层水中不同形态汞的分布规律   总被引：9，自引：1，他引：9  

侯亚敏  冯新斌  仇广乐  阎海鱼 《湖泊科学》2004,16(2):126-133

均匀布点采集百花湖表层水样44个,对水样不同形态汞含量进行了分析,并现场测定水质参数。采用经典统计学与半变异函数拟合相结合的方法,研究了百花湖表层水中不同形态汞含量的空间变异性。根据半变异函数拟合的结果,在ArcGIS软件的支持下进行克立格插值,以此来揭示百花湖表层水中不同形态汞含量的空间分布规律。研究表明:百花湖表层水中的汞污染主要来自上游,随水体流动汞污染减少;但也有一些区域可能因富营养化而使颗粒态汞和总汞含量增加。  相似文献   

Improving Sequential Simulation with a Structured Path Guided by Information Content   总被引：1，自引：0，他引：1  

Yuhong Liu  Andre Journel 《Mathematical Geology》2004,36(8):945-964

Multiple-point simulation is a newly developed geostatistical method that aims at combining the strengths of two mainstream geostatistical methods: object-based and pixel-based methods. It maintains the flexibility of pixel-based algorithms in data conditioning, while enhancing its capability of reproducing realistic geological shapes, which is traditionally reserved to object-based algorithms. However, the current snesim program for multiple-point simulation has difficulty in reproducing large-scale structures, which have a significant impact on the flow response. To address this problem, we propose to simulate along a structured path based on an information content measure. This structured path accounts for not only the information from the data, but also some prior structural information provided by geological knowledge. Various case studies show a better reproduction of large-scale structures. This concept of simulating along a structured path guided by information content can be applied to any sequential simulation algorithms, including traditional variogram-based two-point geostatistical algorithms.  相似文献   

Indicator Simulation Accounting for Multiple-Point Statistics   总被引：7，自引：0，他引：7  

Julián M. Ortiz  Clayton V. Deutsch 《Mathematical Geology》2004,36(5):545-565

Geostatistical simulation aims at reproducing the variability of the real underlying phenomena. When nonlinear features or large-range connectivity is present, the traditional variogram-based simulation approaches do not provide good reproduction of those features. Connectivity of high and low values is often critical for grades in a mineral deposit. Multiple-point statistics can help to characterize these features. The use of multiple-point statistics in geostatistical simulation was proposed more than 10 years ago, on the basis of the use of training images to extract the statistics. This paper proposes the use of multiple-point statistics extracted from actual data. A method is developed to simulate continuous variables. The indicator kriging probabilities used in sequential indicator simulation are modified by probabilities extracted from multiple-point configurations. The correction is done under the assumption of conditional independence. The practical implementation of the method is illustrated with data from a porphyry copper mine.  相似文献   

Application of kriging with external drift to estimate hydraulic conductivity from electrical-resistivity data in unconsolidated deposits near Montalto Uffugo, Italy   总被引：2，自引：0，他引：2  

S. Troisi  C. Fallico  S. Straface  E. Migliari 《Hydrogeology Journal》2000,8(4):356-367

In the framework of various studies to characterize the aquifer at the groundwater experimental field near Montalto Uffugo, Italy, the present work estimates the spatial distribution of hydraulic conductivity of the unconsolidated deposits that underlie the area, by applying the geostatistical technique of kriging with external drift to electrical-resistivity and hydraulic-conductivity data. The reliability of the estimation method was tested by implementing a model, based on the method of cells, that simulates groundwater flow, with the estimated values of hydraulic conductivity. The results obtained indicate that the estimation method used has an acceptable degree of reliability.  相似文献   

Optimization of mobile radioactivity monitoring networks     

G. B.M. Heuvelink  Z. Jiang  S. De Bruin  C. J.W. Twenhöfel 《International journal of geographical information science》2013,27(3):365-382

In case of a nuclear accident, decision makers rely on high-resolution and accurate information about the spatial distribution of radioactive contamination surrounding the accident site. However, the static nuclear monitoring networks of many European countries are generally too coarse to provide the desired level of spatial accuracy. In the Netherlands, authorities are considering a strategy in which measurement density is increased during an emergency using complementary mobile measuring devices. This raises the question, where should these mobile devices be placed? This article proposes a geostatistical methodology to optimize the allocation of mobile measurement devices, such that the expected weighted sum of false-positive and false-negative areas (i.e. false classification into safe and unsafe zones) is minimized. Radioactivity concentration is modelled as the sum of a deterministic trend and a zero-mean spatially correlated stochastic residual. The trend is defined as the outcome of a physical atmospheric dispersion model, NPK-PUFF. The residual is characterized by a semivariogram of differences between the outputs of various NPK-PUFF model runs, designed to reflect the effect of uncertainty in NPK-PUFF meteorological inputs (e.g. wind speed, wind direction). Spatial simulated annealing is used to obtain the optimal monitoring design, in which accessibility of sampling sites (e.g. distance to roads) is also considered. Although the methodology is computationally demanding, results are promising and the computational load may be considerably reduced to compute optimal mobile monitoring designs in nearly real time.  相似文献   

Spatio-temporal analysis of tree height in a young cork oak plantation     

L. Sedda  P.M. Atkinson  M.R. Filigheddu  G. Cotzia  S. Dettori 《International journal of geographical information science》2013,27(7):1083-1096

Cork oak is one of the most valuable natural forest genera in the Mediterranean basin. Modelling cork oak growth has been a challenge for foresters in recent years because of strong site and genetic influences, below-ground competition, management regimes and age effects. Because cork productivity is related to forest height, which is, in turn, related directly with site characteristics, an increase in the accuracy of height prediction implies improved productivity estimation. A Bayesian maximum entropy (BME) geostatistical model was applied to characterize the space–time pattern of height of young cork oak in a forest stand from central Sardinia in the years 2000, 2002, 2003, 2006 and 2008. Cork oak height maps were produced for each of the 5 years. The main goals were to analyse and interpret through time (i) the changes in spatial correlation and (ii) the changes in spatial distribution of cork oak height. The plantation was characterized by an increasing spatial dependence through time, whereas the temporal range was 2 years. Cork oak height was significantly correlated with wind speed (reduced by a neighbouring forest) in all the years implying a single trend. The correlations were larger for 2006 and 2008 than for previous years. Three other environmental variables (shade, elevation and slope) were less significant and their influence restricted to 2 years only. This research has several implications for the management of cork oak in the young phase.  相似文献