A geostatistical approach for areal rainfall statistics assessment

Areal rainfall statistics are more relevant in flood hydrology and water resources management than point rainfall statistics when it comes to help designing dams or hydraulic structures. This paper presents a geostatistically based method to derive the areal statistics from point statistics. Assuming that the distribution models of point rainfall and areal belong to the same class of models and that the rainfall process is stationary, it is shown how the parameters of the areal distribution model can directly be computed from the parameters of the point distribution models in case of a non stationary process, an approximation is derived that yielded good results when applied to a mountainous region in Southern France. The method also allows the computation of the areal reduction factors in a very general form.     

A geostatistical approach for areal rainfall statistics assessment

Areal rainfall statistics are more relevant in flood hydrology and water resources management than point rainfall statistics when it comes to help designing dams or hydraulic structures. This paper presents a geostatistically based method to derive the areal statistics from point statistics. Assuming that the distribution models of point rainfall and areal belong to the same class of models and that the rainfall process is stationary, it is shown how the parameters of the areal distribution model can directly be computed from the parameters of the point distribution models in case of a non stationary process, an approximation is derived that yielded good results when applied to a mountainous region in Southern France. The method also allows the computation of the areal reduction factors in a very general form.     

黑土坡耕地土壤有机质空间变异及其与土壤侵蚀的关系——以黑龙江省鹤山农场为例

选取黑龙江省鹤山农场面积为0.91 km²的典型黑土区的坡耕地作为研究样地。按横纵100 m间隔共采集101个样点,运用地理信息系统和地统计学相结合的方法研究分析0~15 cm土层有机质空间变异及其与土壤侵蚀的关系。结果表明：位于典型黑土区样地的有机质含量集中在3%~5%范围内,均值为4.13%,高于黑龙江省的有机质平均水平。有机质含量空间变异明显,且主要受土壤侵蚀的影响：高侵蚀区对应低有机质区,中度侵蚀区对应中等有机质区,沉积区对应高有机质区。顺坡种植平均坡度2.2°时,每侵蚀1 000 t/km²土壤,有机质含量降低0.8%。土壤有机质空间变异可采用球状模型表达,自相关明显,进一步表明土壤侵蚀导致的再分布。对比分析确定200 m采样间距能够能准确表达该区表层有机质含量的空间特征,为精准施肥提供了采样依据。     

Spatial analysis of land use and shallow groundwater vulnerability in the watershed adjacent to Assateague Island National Seashore,Maryland and Virginia,USA

Spatial relations between land use and groundwater quality in the watershed adjacent to Assateague Island National Seashore, Maryland and Virginia, USA were analyzed by the use of two spatial models. One model used a logit analysis and the other was based on geostatistics. The models were developed and compared on the basis of existing concentrations of nitrate as nitrogen in samples from 529 domestic wells. The models were applied to produce spatial probability maps that show areas in the watershed where concentrations of nitrate in groundwater are likely to exceed a predetermined management threshold value. Maps of the watershed generated by logistic regression and probability kriging analysis showing where the probability of nitrate concentrations would exceed 3 mg/L (>0.50) compared favorably. Logistic regression was less dependent on the spatial distribution of sampled wells, and identified an additional high probability area within the watershed that was missed by probability kriging. The spatial probability maps could be used to determine the natural or anthropogenic factors that best explain the occurrence and distribution of elevated concentrations of nitrate (or other constituents) in shallow groundwater. This information can be used by local land-use planners, ecologists, and managers to protect water supplies and identify land-use planning solutions and monitoring programs in vulnerable areas.     

Optimizing a piezometric network in the estimation of the groundwater budget: a case study from a crystalline-rock watershed in southern India

An estimate of the groundwater budget at the catchment scale is extremely important for the sustainable management of available water resources. Water resources are generally subjected to over-exploitation for agricultural and domestic purposes in agrarian economies like India. The double water-table fluctuation method is a reliable method for calculating the water budget in semi-arid crystalline rock areas. Extensive measurements of water levels from a dense network before and after the monsoon rainfall were made in a 53 km² watershed in southern India and various components of the water balance were then calculated. Later, water level data underwent geostatistical analyses to determine the priority and/or redundancy of each measurement point using a cross-validation method. An optimal network evolved from these analyses. The network was then used in re-calculation of the water-balance components. It was established that such an optimized network provides far fewer measurement points without considerably changing the conclusions regarding groundwater budget. This exercise is helpful in reducing the time and expenditure involved in exhaustive piezometric surveys and also in determining the water budget for large watersheds (watersheds greater than 50 km²).     

Conditional Simulation with Patterns

An entirely new approach to stochastic simulation is proposed through the direct simulation of patterns. Unlike pixel-based (single grid cells) or object-based stochastic simulation, pattern-based simulation simulates by pasting patterns directly onto the simulation grid. A pattern is a multi-pixel configuration identifying a meaningful entity (a puzzle piece) of the underlying spatial continuity. The methodology relies on the use of a training image from which the pattern set (database) is extracted. The use of training images is not new. The concept of a training image is extensively used in simulating Markov random fields or for sequentially simulating structures using multiple-point statistics. Both these approaches rely on extracting statistics from the training image, then reproducing these statistics in multiple stochastic realizations, at the same time conditioning to any available data. The proposed approach does not rely, explicitly, on either a statistical or probabilistic methodology. Instead, a sequential simulation method is proposed that borrows heavily from the pattern recognition literature and simulates by pasting at each visited location along a random path a pattern that is compatible with the available local data and any previously simulated patterns. This paper discusses the various implementation details to accomplish this idea. Several 2D illustrative as well as realistic and complex 3D examples are presented to showcase the versatility of the proposed algorithm.     

Hydrogeological modeling of radionuclide transport in low permeability media: a comparison between Boom Clay and Ypresian Clay

Deep low-permeability clay layers are considered as suitable environments for disposal of high-level radioactive waste. In Belgium, the Boom Clay is the reference host formation and the Ypresian Clay an alternative host formation for research and safety and feasibility assessment of deep disposal of nuclear waste. In this study, two hydrogeological models are built to calculate the radionuclide fluxes that would migrate from a potential repository through these two clay formations. Transport parameter heterogeneity is incorporated in the models using geostatistical co-simulations of hydraulic conductivity, diffusion coefficient and diffusion accessible porosity. The calculated radionuclide fluxes in the two clay formations are compared. The results show that in the Ypresian Clay larger differences between the fluxes through the lower and the upper clay boundary occur, larger total output radionuclide amounts are calculated and a larger effect of parameter heterogeneity on the calculated fluxes is observed, compared to the Boom Clay.     

乌鲁木齐市米东区农田土壤重金属含量的空间分布特征

利用地理信息系统（GIS）及地统计学方法,对乌鲁木齐米东区农田土壤中重金属（Hg、Cu、Zn、Pb、Ni、Cd及Cr）的含量进行空间变异性分析。结果表明：除Cu、Zn、Pb和Hg超过土壤背景值外,7种重金属的平均含量均未超过国家环境质量二级标准（GB15618-1995）;7种重金属均具有较好的空间变异结构,可以用指数模型、球状模型和高斯模型拟合,且具有不同程度的块金效应;Cd、Cr具有强烈的空间自相关性,Pb、Ni、Cu、Zn和Hg属中等强度空间相关, 说明其含量受外源污染的影响较大;采用普通克里格插值法得出7种重金属的空间分布图,除Cd空间分布规律不明显外,其它6种重金属均存在显著的空间分布规律。     

基于协同克里金方法的重力梯度全张量三维约束反演

随着重力梯度全张量测量技术的日趋成熟,重力梯度全张量数据的三维反演技术日益受到重视与关注.全张量数据反演与重力数据反演一样仍然面临着严重的多解性问题.本文将基于地质统计学的协同克里金方法应用于重力梯度全张量数据三维反演,建立了密度约束下的多变量协同克里金联合反演方程,以降低反演的多解性.模型试验表明密度信息的加入能够有效降低反演的多解性,提高反演结果的分辨率,尤其是纵向分辨率能够得到显著提高.最后对美国德克萨斯州一个岩丘区所获得实际资料的应用表明了本文方法的实用性.