COMPARISON OF INDOOR RADON MAPS OBTAINED WITH DIFFERENT SOFTWARES AND METHODS
|
|
- Elinor Haynes
- 6 years ago
- Views:
Transcription
1 Radon in the Living Environment, 039 COMPARISON OF INDOOR RADON MAPS OBTAINED WITH DIFFERENT SOFTWARES AND METHODS François Tondeur (2), Antoine Kies (1), André Robinet (1), Huichao Zhu (3) (1) Centre Universitaire de Luxembourg (CUNLUX) avenue de la faïencerie 162a, L1511 Luxembourg; (2) Institut Supérieur Industriel de Bruxelles (ISIB) rue Royale 150, B0 Bruxelles, Belgium; (3) Faculté Polytechnique de Mons (FPMs) rue de Houdain 9, B7000 Mons, Belgium; Corresponding author: F. Tondeur, Tel: , Fax: , tondeur@isib.be The use of two mapping methods, kriging and moving average, in mapping the indoor radon risk, is investigated. Both methods are applied to three databases: simulated radon data, data collected in Southern Belgium with low sampling density, and data collected in Luxembourg. We use commercial software (SURFER R 6) for kriging, as well as softwares developed by the authors, especially for the radon case, for kriging and for moving average. Simulated data prove to be very useful in this context. We conclude that kriging as implemented in SURFER R 6 may not be well adapted to radon mapping. It seems easier to obtain reasonable maps with the moving average, though it may be difficult to find the good compromise between showing only significant features, and showing all of them. Key words: indoor radon, risk mapping, software, method INTRODUCTION Indoor radon concentrations may be considered as a stochastic function of spatial or spatio-temporal coordinates. This function is characterised by a combination of deterministic and stochastic factors, including an important local noise. Following Cressie (1993), we may write the concentration measured in a house of coordinates (x,y) as the decomposition: C(x,y) = µ(x,y) + ω(x,y) + η(x,y) + ε(x,y) where µ is the deterministic long-range mean structure, ω is a stochastic function reflecting smooth small-scale variations, η is a stochastic function reflecting microscale variations (i.e. variations whose range is smaller than the typical sampling interval) and ε a measurement noise. Usually, η and ε cannot be distinguished in practice and we shall consider η+ε globally as the noise, unless otherwise specified. Radon mapping always consists in trying to smooth out the noise, in order to extract some kind of average trend of risk. Whether the factors controlling the variations of this average trend are of deterministic or stochastic nature, or a mixture of both, is often not discussed. However, the statistical methods to use in the averaging procedure could be more adapted to one case than to the other. The first aim of the present contribution is to present and discuss the application of two mapping methods to different radon databases: kriging and weighted moving average (Cressie,1993). Kriging 373
2 039 Radon in the Living Environment, assumes that all variations, except the noise, are contained in the stochastic function ω; it gives the best prediction of µ + ω(x,y). The moving average instead is a way to evaluate the variable deterministic mean µ(x,y). The second aim of this paper is to test the possibility to use a commercial software (SURFER R 6, Keckler 1994) in radon mapping. Its results will be compared to those of homemade codes specially designed for the radon problem. Our goal at this stage is not to select a method that can be appropriate to obtain one specified kind of radon risk map. We rather want to make a first evaluation of the behaviour of the two considered methods when applied to the radon risk problems, in order to determine the options and factors that are most important in this context. THE DATABASES Maps are generated for different databases that include a spatial location of the data. We first consider simulated data with medium density (2 data per km 2 ) generated by assuming a specified mean value depending on the co-ordinates and a typical lognormal noise with a GSD of 2.5. Three patterns were generated on km x km squares: (a) a constant average of Bq/m 3, (b) an average with a constant slope from to 200 Bq/m 3 ; (c) a cross-shaped area with an average 200 Bq/m 3 within it and Bq/m 3 out of it. We are grateful to J. Miles (NRPB) who generated these data and kindly enabled us to use them in the present work. Using these data to test the mapping methods and softwares proved to be very enlightening. In the present paper, we only present results obtained for the cross pattern, from which the most important conclusions can be drawn. However the sole discussion on simulated data could be unrealistic if the conclusions were not confirmed when applying the same methods to real indoor data. Therefore we also consider here the indoor radon database of the Ministry of Health of Luxembourg, obtained with alpha-track measurements, which has a medium sampling density (1 to 2 data per km 2 ). This database presently only contains the geometrical mean of data grouped by village. This can be considered as a kind of preliminary smoothing. Luxembourg is quite simple from the geological point of view, being divided into two roughly equivalent areas: Oesling in the North and Gutland in the South. Oesling is a part of the Palaeozoic Ardenne-Eifel Massif, where Lower Devonian is dominant. Gutland is mostly constituted of Mesozoic formations. Radon levels are significantly higher in Oesling (Kies, 1996). For Southern Belgium (the Walloon region), the indoor radon database was collected by ISIB with charcoal canister detectors exposed for 3-4 days. This database was presented earlier (Tondeur and Gerardy, 1994) but has been extended by new measurements. However, it still has a quite low sampling density (0.1 data per km 2 ) except for a few small well-sampled areas. It also has a somewhat larger noise due to the short sampling time associated with charcoal detectors. Southern Belgium is also divided into two main parts. The NNW part is made of Mesozoic and Cenozoic grounds. The Palaeozoic dominates the SSE part. None of the two parts is homogenous, as a few localised Palaeozoic 374
3 Radon in the Living Environment, 039 outcrops are found in the NNW part, whereas Mesozoic formations, connected to Gutland, are found in the extreme South. MAPPING METHODS In most of this paper, we shall consider that radon concentrations are, at least approximately, lognormally distributed, and that it is more adequate to take the natural logarithm of C(x,y) as the data to be mapped. Kriging Kriging is a widely used geostatistical method that assumes the stationarity of the expectation value µ of the stochastic variable ln(c) (or at least its stationarity within the area implied by the kriging prediction for each coordinate) but takes into account the correlations between the data that are contained in the short-range ω(x,y). Any global trend in the data is assumed to be the effect of these correlations. The key of kriging is the variogram model. The variogram is defined as: 2 γ(h) = var(ln C(s+h) ln C(s)) where the vector s=(x,y). The variogram of radon data usually show different features: (a) a non-zero value at zero distance (nugget effect) due to the local stochastic noise ε+η (important in radon data); it can be formally separated into two contributions c ME from measurement errors and c MS from microscale variations related respectively to ε and η. (b) a rapid variation related to short-range correlations; (c) a constant trend (sill) at large distance, that indicates the absence of correlations at that distance. Kriging requires that some mathematical function is used to model the variogram. Without nugget effect, kriging acts as an interpolator that gives the best prediction of the stochastic variable, taking the correlations into account through the variogram model. Including a nugget effect transforms kriging into a smoothing interpolator. Moving average The moving average calculates the mean value of all data within a circle (or any other type of neighbourhood), and takes it as the predicted value at the centre. It implicitly assumes that no correlations are present, but accepts a linear trend within that area (a variant allowing for a quadratic trend is also examined here). This trend is understood as a deterministic variation of the mean. The data can be weighted according to their distance d to the centre. With an inverse-square weight (d - 2 ), the method acts as a pure interpolator. We rather use it here as a smoother, with a quadratic weight decreasing with increasing distance (a - b.d 2 ). 375
4 039 Radon in the Living Environment, The mapping softwares Three mapping softwares are used. The first one is a commercial software: SURFER R 6 [5]. It is a grid-base contour program giving the possibility to generate a grid from the data using several gridding methods. We use here the kriging method, that is proposed as the default method. Additional smoothing is possible after gridding. The factors to be specified in the kriging method are the variogram model and the nugget effect. Different variogram models were tested, finally the linear and Gaussian model are retained. Results presented here use the isotropic linear model γ(h)=g.h limited to h<14 km and G=2/14. A satisfactory treatment of radon data with SURFER proved to be impossible without using high values of the nugget parameters (c ME =1, c MS =0.1), that cause kriging to become a smoothing interpolator. But even the nugget smoothing proved to be insufficient for the high variability of radon data. We had to proceed through a further smoothing to eliminate remaining spikes and better approach the general trends. The smoothing consists in an average of 25 values on the kilometric grid, with double weight for the central value. The same kriging and smoothing parameters are used throughout the present work. A few results obtained with lognormal kriging for the Luxembourg database will be presented, using a software developed for radon mapping at FPMs, with a model fitted to the variogram deduced from the data (Zhu, 1996). Finally, we use a weighted moving geometrical mean software with weight decreasing with the distance w(d)=a-bd 2, developed for radon mapping at ISIB. This software also applies a preliminary weighting of the data that is inversely proportional to the local density of data, in order to avoid the excessive weight that could be attributed to well-sampled areas. The geometrical mean at each selected point is calculated within a circle of radius R containing a fixed number N of data, in order to get the same accuracy at each point. The distance-dependent weight is maximum at the centre and vanishes on the circumference of the circle: w(d)=1-(d/r) 2. This was found to be essential to improve the continuity of the calculated µ(x,y). In a variant of the moving average program, we tried to improve the results by performing a quadratic least-squares fit to the data within the circle, keeping the same weights, the fitted value at the center being retained as the predicted µ(x,y). Only minute variations were observed in the maps, and this option will not be applied in the following paragraphs. RESULTS AND DISCUSSION Cross Figure 1 shows results obtained for the cross-like pattern with the moving average method (with N=50 and with N=200), and with SURFER R 6 (kriging alone, and kriging followed by smoothing). Both the kriging map and the N=50 moving average map show strong remaining spurious fluctuations that are only due to the stochastic clustering of high or low values or to outliers, as by construction the 376
5 Radon in the Living Environment, 039 data do not include short-range variations within the two areas of the map. The moving average with N=50 however gives a slightly less fluctuating map than the kriging map. Increasing N to 200 allows to reproduce reasonably the constant mean expected in each of these two areas. The same result is obtained by adding the smoothing step after kriging. Both maps are now strikingly similar. In this context, the usefulness of kriging is not obvious. Clearly, the nugget smoothing alone is not able to smooth out the noise included in simulated data. On the other hand, it seems that the choice of the method is not really crucial, as two calculations based on quite different approaches finally lead to similar results. Southern Belgium (Walloon region) The maps obtained with the moving average method (N=50 and N=200) and with kriging (SURFER R ) for the Walloon region are given in figure 2. Again pure kriging (not shown) is insufficient for smoothing out the noise. After smoothing on 25 km 2, details that are probably not significant (not related to local geology or housing types) still appear. This is due, of course, to the low density of sampling in this database. As smoothing in SURFER R inactivates a band of territory along the border, further smoothing rapidly reduces the area in which predictions are made. Thus, no satisfactory map is obtained with SURFER R for this database. Better maps are obtained with the moving average method. The map obtained with N=50 seems more regular than with simulated data. However, some of the small structures it displays are probably not significant. Taking instead N=200 largely smoothes out many structures that are probably significant, because they are clearly related to geology (e.g. the hot spot in the middle of the NNW part). This drawback is again clearly due to the low sampling density of the database. It is not possible in this case to display all significant structures and simultaneously to avoid the nonsignificant ones. Luxembourg Here better-looking maps are much easier to obtain, because the sampling density is higher, and because of the preliminary smoothing of the data that are grouped by village. They are given in figure 3. The average number of data per village is 6.6. The moving average maps are now given for 8 and 32 villages respectively, which is roughly equivalent on the average to the 50 and 200 data used above (this explains the labels Lux50 and Lux200 used in the figures). Small but probably significant features are lost when going from 8 to 32 villages with the moving average method. Though the sampling density is here much better, the situation is thus similar to the one observed in the previous section: it is not obvious to find the compromise between only displaying significant structures present in the data and displaying all of them. Pure kriging map without nugget with SURFER R 6 is given in figure 3 for the purpose of comparison with a map (figure 5) obtained by lognormal kriging (Rendu,1979) and a fitted variogram model. The variogram is fitted, without nugget, to the experimental variogram (figure 4). For figures 4 and 5, Luxembourg has been divided into Oesling and Gutland, with two different variogram models. The experimental variograms are isotropic in each region. They are fitted by the model function: 377
6 039 Radon in the Living Environment, where γ(h) = c 1 γ 1 (h) + c 2 γ 2 (h) γ i (h) = 1.5 h/a i h 3 3 /a i, h<a i γ i (h) = 1, h>a i. For zone A (Oesling) c 1 = ,a 1 =5516,c 2 = ,a 2 =19456, and for zone B(Gutland) c 1 = ,a 1 =1624,c 2 = ,a 2 =31064, wit c i in units ln(bq.m-3) and a i in m. Without nugget, kriging is an exact interpolator. This would not be reasonable to use for radon mapping if the data were not grouped into villages. When using kriging in SURFER R with the same nugget as above, a map is obtained that is quite similar to the N=32 moving average map (Lux200). No additional smoothing is necessary, because of the preliminary smoothing made by grouping the data. CONCLUSIONS SURFER R 6 with kriging is not well-adapted to radon mapping. Reasonable results were only obtained with pure kriging for the Luxembourg database, with the data grouped into villages, and with a quite good sampling density. Additional smoothing after kriging strongly improves the results when the sampling density is good, as for simulated data, but it was not possible to obtain a reasonable map with SURFER R for the low-density ISIB database of Southern Belgium. In general, map it is easier and more transparent to obtain a satisfactory with the moving average method. However, when the sampling density is too low, it is hard to find the good compromise between displaying the significant features, and not displaying features that are not significant. Any structure in the data, like local hot spots (or cold spots, but they generally deserve less attention) should be, wherever possible, represented by enough data, the number of 200 data being suggested by this work. When both methods (kriging/smoothing and moving average) give acceptable results, the maps are very similar. Thus the choice of one method rather than the other does not seem to be the crucial point, despite their different statistical basic assumptions. Simulated data may be of great interest in testing mapping methods and softwares for the radon problem, and we strongly recommend their systematic use. Of course the results will be more interesting if all authors use the same set of simulated data. The cross data set provided by Jon Miles is quite well adapted. We believe that it would be even more interesting after including a few smaller local structures (hot and cold spots). REFERENCES [1] Cressie N.A.C., Statistics for spatial data, Wiley,
7 Radon in the Living Environment, 039 [2] Keckler D., SURFER R for Windows, User Guide, Golden Software, 1994 [3] Kies A., Biell A., Rowlinson L., Feider M., Radon survey in the Grand Duchy of Luxembourg- Indoor measurements related to house features, soil, geology and environment, Envir.Int. 1996; 22:S5-8 [4] Rendu J.M., Normal and lognormal estimation, J.INT.Assoc.Math.geol.,1979;11: [5] Tondeur F., Gerardy I., Geographic and Geologic Distribution of the Indoor Radon Risk in Southern Belgium, IRPA regional congress on radiological protection, W. Nimmo-Scott, D.J.Golding editors, Nuclear Technology Publishing, 1994, [6] Zhu H.C., Charlet J.M., Doremus P., Kriging radon concentrations og groundwaters in western Ardennes, Environmetrics, 1996;7:
8 039 Radon in the Living Environment, Cross 50 Cross Krigging + nugget Krigging + nugget + smoothing Figure1: Maps obtained for the simulated data with cross-like pattern. The four maps are obtained respectively with the moving average method with N=50 (cross50) and N=200 (cross200), with kriging in SURFER R 6, and with kriging + smoothing. The variable is the natural logarithm of the radon concentration measured in Bq/m 3. 3
9 Radon in the Living Environment, Wallonie 50 Wallonie Krigging + nugget + smoothing Krigging + nugget + heavy smoothing Figure 2: Maps obtained for the Walloon region, with the moving average method with N=50 and N=200, with kriging + smoothing in SURFER R 6, and with kriging + heavier smoothing. The variable is the natural logarithm of the radon concentration measured in Bq/m
10 039 Radon in the Living Environment, Lux 50 Lux Lux krigging Lux krigging + nugget Figure 3: Maps obtained for Luxembourg, with the moving average method with N=8 (Lux50) and N=32 (Lux200), with kriging without nugget in SURFER R 6, and with kriging with nugget. The variable is the natural logarithm of the radon concentration measured in Bq/m
11 Radon in the Living Environment, 039 Figure 4: Isotropic variograms deduced from experimental data and fitted variogram models for Luxembourg: zone A = Oesling, zone B = Gutland. 383
12 039 Radon in the Living Environment, (m) Figure 5: Map obtained for Luxembourg with lognormal kriging without nugget, the variogram model being fitted to the experimental variogram. The variable is the predicted concentration. 384
INDOOR RADON ON LOESS DEPOSITS IN WALLONIA
INDOOR RADON ON LOESS DEPOSITS IN WALLONIA François TONDEUR, Isabelle GERARDY, Nathalie GERARDY Institut supérieur industriel de Bruxelles, 150 rue Royale, B1000 BRUSSELS, BELGIUM ABSTRACT Pleistocene
More informationQUEST OF PROXIES FOR INDOOR RADON RISK : BELGIAN EXPERIENCE
QUEST OF PROXIES FOR INDOOR RADON RISK : BELGIAN EXPERIENCE FRANÇOIS TONDEUR I S IB - HE2B, 150 R UE ROYALE, 1000 BRUSSELS - BELGIUM 1 Proxy for indoor Rn risk: 1. Context Property well-correlated to indoor
More informationAdvanced analysis and modelling tools for spatial environmental data. Case study: indoor radon data in Switzerland
EnviroInfo 2004 (Geneva) Sh@ring EnviroInfo 2004 Advanced analysis and modelling tools for spatial environmental data. Case study: indoor radon data in Switzerland Mikhail Kanevski 1, Michel Maignan 1
More informationPRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH
PRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH SURESH TRIPATHI Geostatistical Society of India Assumptions and Geostatistical Variogram
More informationAutomatic Determination of Uncertainty versus Data Density
Automatic Determination of Uncertainty versus Data Density Brandon Wilde and Clayton V. Deutsch It is useful to know how various measures of uncertainty respond to changes in data density. Calculating
More informationSURVEY OF THE 137 Cs CONTAMINATION IN BELGIUM AFTER THE CHERNOBYL ACCIDENT BY IN-SITU GAMMA SPECTROMETRY
SURVEY OF THE 137 Cs CONTAMINATION IN BELGIUM AFTER THE CHERNOBYL ACCIDENT BY IN-SITU GAMMA SPECTROMETRY J. Uyttenhove 1, S. Pommé 2, B. Van Waeyenberge 3 1 University Gent, Physics Lab, Krijgslaan 281
More informationRADON RISK MAPPING IN IRELAND
Radon in the Living Environment, 137 RADON RISK MAPPING IN IRELAND S.G. Fennell 1, Y. Pawitan 2, G.M. Mackin 1, J.S. Madden 1 and A.T. McGarry 1 1 Radiological Protection Institute of Ireland, 3 Clonskeagh
More informationTowards a European map of indoor radon
Towards a European map of indoor radon Current issues and forthcoming challenges G. Dubois JRC - European Commission IES - Emissions and Health Unit Radioactivity Environmental Monitoring () group Gregoire.dubois@jrc.it
More informationQuantifying uncertainty of geological 3D layer models, constructed with a-priori
Quantifying uncertainty of geological 3D layer models, constructed with a-priori geological expertise Jan Gunnink, Denise Maljers 2 and Jan Hummelman 2, TNO Built Environment and Geosciences Geological
More informationIndex. Geostatistics for Environmental Scientists, 2nd Edition R. Webster and M. A. Oliver 2007 John Wiley & Sons, Ltd. ISBN:
Index Akaike information criterion (AIC) 105, 290 analysis of variance 35, 44, 127 132 angular transformation 22 anisotropy 59, 99 affine or geometric 59, 100 101 anisotropy ratio 101 exploring and displaying
More informationStatistícal Methods for Spatial Data Analysis
Texts in Statistícal Science Statistícal Methods for Spatial Data Analysis V- Oliver Schabenberger Carol A. Gotway PCT CHAPMAN & K Contents Preface xv 1 Introduction 1 1.1 The Need for Spatial Analysis
More information7 Geostatistics. Figure 7.1 Focus of geostatistics
7 Geostatistics 7.1 Introduction Geostatistics is the part of statistics that is concerned with geo-referenced data, i.e. data that are linked to spatial coordinates. To describe the spatial variation
More informationGeostatistics in Hydrology: Kriging interpolation
Chapter Geostatistics in Hydrology: Kriging interpolation Hydrologic properties, such as rainfall, aquifer characteristics (porosity, hydraulic conductivity, transmissivity, storage coefficient, etc.),
More informationMapping Precipitation in Switzerland with Ordinary and Indicator Kriging
Journal of Geographic Information and Decision Analysis, vol. 2, no. 2, pp. 65-76, 1998 Mapping Precipitation in Switzerland with Ordinary and Indicator Kriging Peter M. Atkinson Department of Geography,
More informationKEYWORDS : Indoor radon mapping, radon risk map, nugget effect, radon potential
The radon noise and its geostatistical implications: risk mapping or mapping at risk? G. Dubois and P. Bossew European Commission DG Joint Research Centre, Institute for Environment and Sustainability,
More information4th HR-HU and 15th HU geomathematical congress Geomathematics as Geoscience Reliability enhancement of groundwater estimations
Reliability enhancement of groundwater estimations Zoltán Zsolt Fehér 1,2, János Rakonczai 1, 1 Institute of Geoscience, University of Szeged, H-6722 Szeged, Hungary, 2 e-mail: zzfeher@geo.u-szeged.hu
More informationSome practical aspects of the use of lognormal models for confidence limits and block distributions in South African gold mines
Some practical aspects of the use of lognormal models for confidence limits and block distributions in South African gold mines by D.G. Krige* Synopsis For the purpose of determining confidence limits
More informationSpatial Data Analysis in Archaeology Anthropology 589b. Kriging Artifact Density Surfaces in ArcGIS
Spatial Data Analysis in Archaeology Anthropology 589b Fraser D. Neiman University of Virginia 2.19.07 Spring 2007 Kriging Artifact Density Surfaces in ArcGIS 1. The ingredients. -A data file -- in.dbf
More informationA Short Note on the Proportional Effect and Direct Sequential Simulation
A Short Note on the Proportional Effect and Direct Sequential Simulation Abstract B. Oz (boz@ualberta.ca) and C. V. Deutsch (cdeutsch@ualberta.ca) University of Alberta, Edmonton, Alberta, CANADA Direct
More informationSPATIAL VARIABILITY MAPPING OF N-VALUE OF SOILS OF MUMBAI CITY USING ARCGIS
SPATIAL VARIABILITY MAPPING OF N-VALUE OF SOILS OF MUMBAI CITY USING ARCGIS RESHMA RASKAR - PHULE 1, KSHITIJA NADGOUDA 2 1 Assistant Professor, Department of Civil Engineering, Sardar Patel College of
More informationTricks to Creating a Resource Block Model. St John s, Newfoundland and Labrador November 4, 2015
Tricks to Creating a Resource Block Model St John s, Newfoundland and Labrador November 4, 2015 Agenda 2 Domain Selection Top Cut (Grade Capping) Compositing Specific Gravity Variograms Block Size Search
More informationI don t have much to say here: data are often sampled this way but we more typically model them in continuous space, or on a graph
Spatial analysis Huge topic! Key references Diggle (point patterns); Cressie (everything); Diggle and Ribeiro (geostatistics); Dormann et al (GLMMs for species presence/abundance); Haining; (Pinheiro and
More informationENGRG Introduction to GIS
ENGRG 59910 Introduction to GIS Michael Piasecki October 13, 2017 Lecture 06: Spatial Analysis Outline Today Concepts What is spatial interpolation Why is necessary Sample of interpolation (size and pattern)
More informationTime-lapse filtering and improved repeatability with automatic factorial co-kriging. Thierry Coléou CGG Reservoir Services Massy
Time-lapse filtering and improved repeatability with automatic factorial co-kriging. Thierry Coléou CGG Reservoir Services Massy 1 Outline Introduction Variogram and Autocorrelation Factorial Kriging Factorial
More informationSpatial Backfitting of Roller Measurement Values from a Florida Test Bed
Spatial Backfitting of Roller Measurement Values from a Florida Test Bed Daniel K. Heersink 1, Reinhard Furrer 1, and Mike A. Mooney 2 1 Institute of Mathematics, University of Zurich, CH-8057 Zurich 2
More informationThe Snap lake diamond deposit - mineable resource.
Title Page: Author: Designation: Affiliation: Address: Tel: Email: The Snap lake diamond deposit - mineable resource. Fanie Nel Senior Mineral Resource Analyst De Beers Consolidated Mines Mineral Resource
More informationSPATIAL ELECTRICAL LOADS MODELING USING THE GEOSTATISTICAL METHODS
19 th International CODATA Conference THE INFORMATION SOCIETY: NEW HORIZONS FOR SCIENCE Berlin, Germany 7-1 November 24 SPATIAL ELECTRICAL LOADS MODELING USING THE GEOSTATISTICAL METHODS Barbara Namysłowska-Wilczyńska
More informationAnomaly Density Estimation from Strip Transect Data: Pueblo of Isleta Example
Anomaly Density Estimation from Strip Transect Data: Pueblo of Isleta Example Sean A. McKenna, Sandia National Laboratories Brent Pulsipher, Pacific Northwest National Laboratory May 5 Distribution Statement
More informationSpatial Analysis II. Spatial data analysis Spatial analysis and inference
Spatial Analysis II Spatial data analysis Spatial analysis and inference Roadmap Spatial Analysis I Outline: What is spatial analysis? Spatial Joins Step 1: Analysis of attributes Step 2: Preparing for
More informationPOPULAR CARTOGRAPHIC AREAL INTERPOLATION METHODS VIEWED FROM A GEOSTATISTICAL PERSPECTIVE
CO-282 POPULAR CARTOGRAPHIC AREAL INTERPOLATION METHODS VIEWED FROM A GEOSTATISTICAL PERSPECTIVE KYRIAKIDIS P. University of California Santa Barbara, MYTILENE, GREECE ABSTRACT Cartographic areal interpolation
More informationWhat is Non-Linear Estimation?
What is Non-Linear Estimation? You may have heard the terms Linear Estimation and Non-Linear Estimation used in relation to spatial estimation of a resource variable and perhaps wondered exactly what they
More informationSpatial Data Mining. Regression and Classification Techniques
Spatial Data Mining Regression and Classification Techniques 1 Spatial Regression and Classisfication Discrete class labels (left) vs. continues quantities (right) measured at locations (2D for geographic
More informationConditional Distribution Fitting of High Dimensional Stationary Data
Conditional Distribution Fitting of High Dimensional Stationary Data Miguel Cuba and Oy Leuangthong The second order stationary assumption implies the spatial variability defined by the variogram is constant
More informationBuilding Blocks for Direct Sequential Simulation on Unstructured Grids
Building Blocks for Direct Sequential Simulation on Unstructured Grids Abstract M. J. Pyrcz (mpyrcz@ualberta.ca) and C. V. Deutsch (cdeutsch@ualberta.ca) University of Alberta, Edmonton, Alberta, CANADA
More information11/8/2018. Spatial Interpolation & Geostatistics. Kriging Step 1
(Z i Z j ) 2 / 2 (Z i Zj) 2 / 2 Semivariance y 11/8/2018 Spatial Interpolation & Geostatistics Kriging Step 1 Describe spatial variation with Semivariogram Lag Distance between pairs of points Lag Mean
More informationA robust statistically based approach to estimating the probability of contamination occurring between sampling locations
A robust statistically based approach to estimating the probability of contamination occurring between sampling locations Peter Beck Principal Environmental Scientist Image placeholder Image placeholder
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry
More informationInfluence of parameter estimation uncertainty in Kriging: Part 2 Test and case study applications
Hydrology and Earth System Influence Sciences, of 5(), parameter 5 3 estimation (1) uncertainty EGS in Kriging: Part Test and case study applications Influence of parameter estimation uncertainty in Kriging:
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Department of Forestry & Department of Geography, Michigan State University, Lansing Michigan, U.S.A. 2 Biostatistics,
More informationUsing linear and non-linear kriging interpolators to produce probability maps
To be presented at 2001 Annual Conference of the International Association for Mathematical Geology, Cancun, Mexico, September, IAMG2001. Using linear and non-linear kriging interpolators to produce probability
More informationChapter 4 - Fundamentals of spatial processes Lecture notes
TK4150 - Intro 1 Chapter 4 - Fundamentals of spatial processes Lecture notes Odd Kolbjørnsen and Geir Storvik January 30, 2017 STK4150 - Intro 2 Spatial processes Typically correlation between nearby sites
More informationCOLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION
COLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION G. MARIETHOZ, PH. RENARD, R. FROIDEVAUX 2. CHYN, University of Neuchâtel, rue Emile Argand, CH - 2009 Neuchâtel, Switzerland 2 FSS Consultants, 9,
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 1 y Kriging Step 1 Describe spatial variation with Semivariogram (Z i Z j ) 2 / 2 Point cloud Map 3
More informationBayesian Transgaussian Kriging
1 Bayesian Transgaussian Kriging Hannes Müller Institut für Statistik University of Klagenfurt 9020 Austria Keywords: Kriging, Bayesian statistics AMS: 62H11,60G90 Abstract In geostatistics a widely used
More informationMapping the Baseline of Terrestrial Gamma Radiation in China
Radiation Environment and Medicine 2017 Vol.6, No.1 29 33 Note Mapping the Baseline of Terrestrial Gamma Radiation in China Zhen Yang, Weihai Zhuo* and Bo Chen Institute of Radiation Medicine, Fudan University,
More informationFinding the Nearest Positive Definite Matrix for Input to Semiautomatic Variogram Fitting (varfit_lmc)
Finding the Nearest Positive Definite Matrix for Input to Semiautomatic Variogram Fitting (varfit_lmc) Arja Jewbali (arja.jewbali@riotinto.com) Resource Estimation Geologist Rio Tinto Iron Ore In resource
More informationCOLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION
COLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION G. MARIETHOZ, PH. RENARD, R. FROIDEVAUX 2. CHYN, University of Neuchâtel, rue Emile Argand, CH - 2009 Neuchâtel, Switzerland 2 FSS Consultants, 9,
More informationSpatial analysis. 0 move the objects and the results change
0 Outline: Roadmap 0 What is spatial analysis? 0 Transformations 0 Introduction to spatial interpolation 0 Classification of spatial interpolation methods 0 Interpolation methods 0 Areal interpolation
More informationPropagation of Errors in Spatial Analysis
Stephen F. Austin State University SFA ScholarWorks Faculty Presentations Spatial Science 2001 Propagation of Errors in Spatial Analysis Peter P. Siska I-Kuai Hung Arthur Temple College of Forestry and
More informationSPATIAL-TEMPORAL TECHNIQUES FOR PREDICTION AND COMPRESSION OF SOIL FERTILITY DATA
SPATIAL-TEMPORAL TECHNIQUES FOR PREDICTION AND COMPRESSION OF SOIL FERTILITY DATA D. Pokrajac Center for Information Science and Technology Temple University Philadelphia, Pennsylvania A. Lazarevic Computer
More informationEstimation of direction of increase of gold mineralisation using pair-copulas
22nd International Congress on Modelling and Simulation, Hobart, Tasmania, Australia, 3 to 8 December 2017 mssanz.org.au/modsim2017 Estimation of direction of increase of gold mineralisation using pair-copulas
More informationCOMPARISON OF DIGITAL ELEVATION MODELLING METHODS FOR URBAN ENVIRONMENT
COMPARISON OF DIGITAL ELEVATION MODELLING METHODS FOR URBAN ENVIRONMENT Cahyono Susetyo Department of Urban and Regional Planning, Institut Teknologi Sepuluh Nopember, Indonesia Gedung PWK, Kampus ITS,
More informationApplying MapCalc Map Analysis Software
Applying MapCalc Map Analysis Software Generating Surface Maps from Point Data: A farmer wants to generate a set of maps from soil samples he has been collecting for several years. Previously, he would
More informationBasics of Point-Referenced Data Models
Basics of Point-Referenced Data Models Basic tool is a spatial process, {Y (s), s D}, where D R r Chapter 2: Basics of Point-Referenced Data Models p. 1/45 Basics of Point-Referenced Data Models Basic
More informationTransactions on Information and Communications Technologies vol 16, 1996 WIT Press, ISSN
Artificial Neural Networks and Geostatistics for Environmental Mapping Kanevsky M. (), Arutyunyan R. (), Bolshov L. (), Demyanov V. (), Savelieva E. (), Haas T. (2), Maignan M.(3) () Institute of Nuclear
More informationEnvironmental Data Mining and Modelling Based on Machine Learning Algorithms and Geostatistics
Environmental Data Mining and Modelling Based on Machine Learning Algorithms and Geostatistics M. Kanevski a, R. Parkin b, A. Pozdnukhov b, V. Timonin b, M. Maignan c, B. Yatsalo d, S. Canu e a IDIAP Dalle
More informationSupplementary material: Methodological annex
1 Supplementary material: Methodological annex Correcting the spatial representation bias: the grid sample approach Our land-use time series used non-ideal data sources, which differed in spatial and thematic
More informationOptimizing Sampling Schemes for Mapping and Dredging Polluted Sediment Layers
This file was created by scanning the printed publication. Errors identified by the software have been corrected; however, some errors may remain. Optimizing Sampling Schemes for Mapping and Dredging Polluted
More informationA full scale, non stationary approach for the kriging of large spatio(-temporal) datasets
A full scale, non stationary approach for the kriging of large spatio(-temporal) datasets Thomas Romary, Nicolas Desassis & Francky Fouedjio Mines ParisTech Centre de Géosciences, Equipe Géostatistique
More informationGeostatistics for radiological characterization: overview and application cases
Geostatistics for radiological characterization: overview and application cases Yvon DESNOYERS a* a GEOVARIANCES, 49bis av. Franklin Roosevelt, 77215 Avon, France *corresponding author: yvon.desnoyers@geovariances.com
More informationExtent of Radiological Contamination in Soil at Four Sites near the Fukushima Daiichi Power Plant, Japan (ArcGIS)
Extent of Radiological Contamination in Soil at Four Sites near the Fukushima Daiichi Power Plant, Japan (ArcGIS) Contact: Ted Parks, AMEC Foster Wheeler, theodore.parks@amecfw.com, Alex Mikszewski, AMEC
More informationThe Proportional Effect of Spatial Variables
The Proportional Effect of Spatial Variables J. G. Manchuk, O. Leuangthong and C. V. Deutsch Centre for Computational Geostatistics, Department of Civil and Environmental Engineering University of Alberta
More information11. Kriging. ACE 492 SA - Spatial Analysis Fall 2003
11. Kriging ACE 492 SA - Spatial Analysis Fall 2003 c 2003 by Luc Anselin, All Rights Reserved 1 Objectives The goal of this lab is to further familiarize yourself with ESRI s Geostatistical Analyst, extending
More informationINDOOR RADON MAPPING FOR NEW YORK STATE
INDOOR RADON MAPPING FOR NEW YORK STATE C. Kunz, J. Green, C. Schwenker, E. Rigilski, and M. Kitto NYS Department of Health, Wadsworth Center, Empire State Plaza Albany, NY ABSTRACT The percent of homes
More informationLecture 5 Geostatistics
Lecture 5 Geostatistics Lecture Outline Spatial Estimation Spatial Interpolation Spatial Prediction Sampling Spatial Interpolation Methods Spatial Prediction Methods Interpolating Raster Surfaces with
More informationNon-gaussian spatiotemporal modeling
Dec, 2008 1/ 37 Non-gaussian spatiotemporal modeling Thais C O da Fonseca Joint work with Prof Mark F J Steel Department of Statistics University of Warwick Dec, 2008 Dec, 2008 2/ 37 1 Introduction Motivation
More information3D Visualization of Geothermal Features
Proceedings World Geothermal Congress 25 Antalya, Turkey, 24-29 April 25 3D Visualization of Geothermal Features * Umran Serpen, * Eylem Kaya and + Niyazi Aksoy * Petroleum and Natural Gas Engineering
More informationAcceptable Ergodic Fluctuations and Simulation of Skewed Distributions
Acceptable Ergodic Fluctuations and Simulation of Skewed Distributions Oy Leuangthong, Jason McLennan and Clayton V. Deutsch Centre for Computational Geostatistics Department of Civil & Environmental Engineering
More informationRegional Variation of Seasonal Behaviour for Indoor Radon Levels
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS006) p.3924 Regional Variation of Seasonal Behaviour for Indoor Radon Levels Burke, Órlaith University of Oxford,
More informationInternational Atomic Energy Agency Learning programme: Radon gas. Module 4: Developing and Implementing a Representative Indoor Radon Survey
International Atomic Energy Agency Learning programme: Radon gas Module 4: Developing and Implementing a Representative Indoor Radon Survey Content Scope of this module Representative radon survey aims
More informationAn Introduction to Pattern Statistics
An Introduction to Pattern Statistics Nearest Neighbors The CSR hypothesis Clark/Evans and modification Cuzick and Edwards and controls All events k function Weighted k function Comparative k functions
More informationMagnetic Case Study: Raglan Mine Laura Davis May 24, 2006
Magnetic Case Study: Raglan Mine Laura Davis May 24, 2006 Research Objectives The objective of this study was to test the tools available in EMIGMA (PetRos Eikon) for their utility in analyzing magnetic
More informationModeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation
Modeling of Atmospheric Effects on InSAR Measurements With the Method of Stochastic Simulation Z. W. LI, X. L. DING Department of Land Surveying and Geo-Informatics, Hong Kong Polytechnic University, Hung
More informationMultiple realizations: Model variance and data uncertainty
Stanford Exploration Project, Report 108, April 29, 2001, pages 1?? Multiple realizations: Model variance and data uncertainty Robert G. Clapp 1 ABSTRACT Geophysicists typically produce a single model,
More informationToward an automatic real-time mapping system for radiation hazards
Toward an automatic real-time mapping system for radiation hazards Paul H. Hiemstra 1, Edzer J. Pebesma 2, Chris J.W. Twenhöfel 3, Gerard B.M. Heuvelink 4 1 Faculty of Geosciences / University of Utrecht
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO1639 Importance of density-compensated temperature change for deep North Atlantic Ocean heat uptake C. Mauritzen 1,2, A. Melsom 1, R. T. Sutton 3 1 Norwegian
More informationDecision-Oriented Environmental Mapping with Radial Basis Function Neural Networks
Decision-Oriented Environmental Mapping with Radial Basis Function Neural Networks V. Demyanov (1), N. Gilardi (2), M. Kanevski (1,2), M. Maignan (3), V. Polishchuk (1) (1) Institute of Nuclear Safety
More informationAssessing the covariance function in geostatistics
Statistics & Probability Letters 52 (2001) 199 206 Assessing the covariance function in geostatistics Ana F. Militino, M. Dolores Ugarte Departamento de Estadstica e Investigacion Operativa, Universidad
More informationSoil uranium, soil gas radon and indoor radon empirical relationships in the UK and other European countries
Soil uranium, soil gas radon and indoor radon empirical relationships in the UK and other European countries Don Appleton (BGS) and Jon Miles (HPA) EUROPEAN GEOGENIC RADON POTENTIAL MAP Wide range of numerical
More informationDrill-Holes and Blast-Holes
Drill-Holes and Blast-Holes Serge Antoine Séguret 1 and Sebastian De La Fuente 2 Abstract The following is a geostatistical study of copper measurements on samples from diamond drill-holes and blast-holes.
More informationSpatial analysis of annual rainfall using ordinary kriging techniques and lognormal kriging in the Cheliff watershed. Algeria.
15 th International Conference on Environmental Science and Technology Rhodes, Greece, 31 August to 2 September 2017 Spatial analysis of annual rainfall using ordinary kriging techniques and lognormal
More informationThis document is a preview generated by EVS
INTERNATIONAL STANDARD ISO 11665-6 First edition 2012-07-15 Measurement of radioactivity in the environment Air: radon-222 Part 6: Spot measurement method of the activity concentration Mesurage de la radioactivité
More informationCombining geological surface data and geostatistical model for Enhanced Subsurface geological model
Combining geological surface data and geostatistical model for Enhanced Subsurface geological model M. Kurniawan Alfadli, Nanda Natasia, Iyan Haryanto Faculty of Geological Engineering Jalan Raya Bandung
More informationExploring the World of Ordinary Kriging. Dennis J. J. Walvoort. Wageningen University & Research Center Wageningen, The Netherlands
Exploring the World of Ordinary Kriging Wageningen University & Research Center Wageningen, The Netherlands July 2004 (version 0.2) What is? What is it about? Potential Users a computer program for exploring
More information5.1 2D example 59 Figure 5.1: Parabolic velocity field in a straight two-dimensional pipe. Figure 5.2: Concentration on the input boundary of the pipe. The vertical axis corresponds to r 2 -coordinate,
More informationGridding of precipitation and air temperature observations in Belgium. Michel Journée Royal Meteorological Institute of Belgium (RMI)
Gridding of precipitation and air temperature observations in Belgium Michel Journée Royal Meteorological Institute of Belgium (RMI) Gridding of meteorological data A variety of hydrologic, ecological,
More informationEco504, Part II Spring 2010 C. Sims PITFALLS OF LINEAR APPROXIMATION OF STOCHASTIC MODELS
Eco504, Part II Spring 2010 C. Sims PITFALLS OF LINEAR APPROXIMATION OF STOCHASTIC MODELS 1. A LIST OF PITFALLS Linearized models are of course valid only locally. In stochastic economic models, this usually
More informationSpatial Inference of Nitrate Concentrations in Groundwater
Spatial Inference of Nitrate Concentrations in Groundwater Dawn Woodard Operations Research & Information Engineering Cornell University joint work with Robert Wolpert, Duke Univ. Dept. of Statistical
More informationAdvances in Locally Varying Anisotropy With MDS
Paper 102, CCG Annual Report 11, 2009 ( 2009) Advances in Locally Varying Anisotropy With MDS J.B. Boisvert and C. V. Deutsch Often, geology displays non-linear features such as veins, channels or folds/faults
More informationA MultiGaussian Approach to Assess Block Grade Uncertainty
A MultiGaussian Approach to Assess Block Grade Uncertainty Julián M. Ortiz 1, Oy Leuangthong 2, and Clayton V. Deutsch 2 1 Department of Mining Engineering, University of Chile 2 Department of Civil &
More informationMineral resource classification: a comparison of new and existing techniques
Mineral resource classification: a comparison of new and existing techniques by D.S.F. Silva* and J. B. Boisvert* Synopsis A survey of 120 recent NI 43-101 technical reports was conducted to evaluate the
More informationThe Kawasaki Identity and the Fluctuation Theorem
Chapter 6 The Kawasaki Identity and the Fluctuation Theorem This chapter describes the Kawasaki function, exp( Ω t ), and shows that the Kawasaki function follows an identity when the Fluctuation Theorem
More informationManuscript of paper for APCOM 2003.
1 Manuscript of paper for APCOM 2003. AN ANALYSIS OF THE PRACTICAL AND ECONOMIC IMPLICATIONS OF SYSTEMATIC UNDERGROUND DRILLING IN DEEP SOUTH AFRICAN GOLD MINES W. ASSIBEY-BONSU Consultant: Geostatistics
More informationAnother Look at Non-Euclidean Variography
Another Look at Non-Euclidean Variography G. Dubois European Commission DG Joint Research Centre Institute for Environment and Sustainability, Ispra, Italy. Email: gregoire.dubois@jrc.it ABSTRACT: Tobler
More informationISO Measurement of radioactivity in the environment Air: radon-222 Part 5: Continuous measurement method of the activity concentration
INTERNATIONAL STANDARD ISO 11665-5 First edition 2012-07-15 Measurement of radioactivity in the environment Air: radon-222 Part 5: Continuous measurement method of the activity concentration Mesurage de
More informationBasic Geostatistics: Pattern Description
Basic Geostatistics: Pattern Description A SpaceStat Software Tutorial Copyright 2013, BioMedware, Inc. (www.biomedware.com). All rights reserved. SpaceStat and BioMedware are trademarks of BioMedware,
More informationA sensitivity and uncertainty analysis. Ministry of the Walloon Region Agricultural Research Centre
Development of an agrometeorological model integrating leaf wetness duration estimation and weather radar data to assess the risk of head blight infection in wheat A sensitivity and uncertainty analysis
More informationRADON MONITORING AT THE JOINT RESEARCH CENTRE ISPRA. S. Verdelocco, D. Walker, P. Turkowsky, and C. Osimani
Radon in the Living Environment, 081 RADON MONITORING AT THE JOINT RESEARCH CENTRE ISPRA S. Verdelocco, D. Walker, P. Turkowsky, and C. Osimani Joint Research Centre Ispra, Radiation Protection Unit, Via
More informationComparison of rainfall distribution method
Team 6 Comparison of rainfall distribution method In this section different methods of rainfall distribution are compared. METEO-France is the French meteorological agency, a public administrative institution
More informationEdge Detection in Computer Vision Systems
1 CS332 Visual Processing in Computer and Biological Vision Systems Edge Detection in Computer Vision Systems This handout summarizes much of the material on the detection and description of intensity
More information