Effect of Land Use Types on the Spatial Prediction of Soil Nitrogen
|
|
- Jonathan Hudson
- 6 years ago
- Views:
Transcription
1 Effect of Land Use Types on the Spatial Prediction of Soil Nitrogen Mingkai Qu Department of Resource and Environmental Information, College of Resources and Environment, Huazhong Agricultural University, Wuhan , China Weidong Li 1 and Chuanrong Zhang Department of Geography and Center for Environmental Sciences and Engineering, University of Connecticut, Storrs, Connecticut Shanqin Wang 1 Department of Resource and Environmental Information, College of Resources and Environment, Huazhong Agricultural University, Wuhan , China Abstract: Mapping the spatial distribution of soil nutrient contents from sample data has received much attention in the recent decade. Accurately mapping soil nutrients purely based on sample data, however, is difficult due to the sparsity and high cost of samples. Land use types usually influence the contents of soil nutrients at the local level and it is desirable to integrate such information into predictive mapping. The area-and-point kriging (AAPK) method, which was proposed recently, may provide an interpolation technique for such purposes. This study mapped the soil total nitrogen (TN) distribution of Hanchuan County, China, using AAPK with sample data (consisting of 402 points) and land use information. Ordinary kriging (OK) and residual kriging (RK) were compared to evaluate the performance of AAPK. Results showed that: (1) land use types had important impacts on the spatial distribution of soil TN; (2) measured data at 135 validation locations had stronger correlation with the data predicted by AAPK than by RK and OK, and the mean error and root mean square error with AAPK were lower than with RK and OK; and (3) AAPK generated smaller error variances than RK and OK did. This suggests that AAPK represents an effective method for increasing the interpolation accuracy of soil TN. It should be pointed out that some of the land use polygons used in this study are very large and complex, which might impact the effectiveness of AAPK in improving the prediction accuracy. Segmenting them into simple smaller areas might be helpful. INTRODUCTION Nitrogen (N) is an important nutrient in soil, a basic resource for maintaining the Earth s ecosystems, and a primary restrictive factor for crop production (Wu et al., 2009). To improve crop production, nitrogenous fertilizers have been widely applied 1 Corresponding authors; sqwang@mail.hzau.edu.cn; weidongwoody@gmail.com 397 GIScience & Remote Sensing, 2012, 49, No. 3, p Copyright 2012 by Bellwether Publishing, Ltd. All rights reserved.
2 398 q u e t a l. to farmlands. Often, the amount of soil total nitrogen (TN) exceeds the requirements of plant growth. High N fertilization rates generally result in low N use efficiency and high N loss (Li and Zhang, 1999). While effective use of N can improve crop yields, excessive use of nitrogenous fertilizers leads to negative impacts on surrounding environments, especially the aquatic environment (Carpenter et al., 1998; Smith et al., 2001; Lu et al., 2007). Therefore, the sustainable use of soils and protection of the environment require a better understanding of soil TN content in fields and its spatial variability. Geostatistics is a set of powerful spatial statistical techniques that has been widely used to characterize the spatial variability of soil properties and map their spatial distributions (Odeh et al., 1995; Ferguson et al., 1998; Burgess and Webster, 1980). Because the quality of mapping of soil properties affects the performance of site-specific soil management, how to more accurately map the target soil variable with limited sample data becomes an important study issue. In order to obtain more accurate results, an estimation model such as kriging is needed to combine sparse but accurate sample data with related abundant but inaccurate or qualitative information as the input data. One can thus use inexpensively measured data of auxiliary variables to improve the kriging estimation of a target variable that may be expensive to measure. The information available for mapping continuous soil attributes often includes field point sample data and related categorical maps (e.g., land use, soil type, or geological maps). Thus, the issue of combining data measured on different spatial supports (e.g., point and map unit) has been an important research topic in soil science. Residual kriging (RK) has been often used to integrate categorical information that can be used to inform the local mean of the random function (Goovaerts and Journel, 1995; Hengl et al., 2004; Liu et al., 2006; Zhang et al., 2010). In this approach, variography and kriging are conducted only on the stationary residuals (i.e., differences between point measurements and their means in different mapping units), and the final estimates are composed of kriged residual results and local means. By taking into account the variation components of soil (or land use) type effect and residual, some improvement in kriging interpolation accuracy can be expected. However, because RK proceeds in two steps it cannot guarantee that the final map of kriging estimates will honor the areal data: the average of the interpolated values within each area typically does not equal the areal datum. Recently, a new kriging method called area-and-point kriging (AAPK) emerged for combining areal and point data in geostatistical interpolation (Liu and Journel, 2009; Goovaerts, 2010). This method allows combination of both point and areal data through use of area-toarea, area-to-point, and point-to-point covariances in the kriging system. It ensures the coherence of the prediction so that the average of interpolated values within each mapping unit is equal to the original areal datum. Moreover, this approach capitalizes on the availability of GIS to discretize polygons of irregular shape and size and the knowledge of the point-support variogram model that can be inferred directly from point measurements, thereby eliminating the need for a deconvolution procedure (Goovaerts, 2008). Testing case studies demonstrated that AAPK generally has higher prediction accuracy than ordinary and traditional residual kriging based on the assumption that the local mean is constant within each mapping unit and the AAPK variance is a more accurate indicator of the magnitude of prediction errors (Goovaerts, 2010).
3 l a n d u s e t y p e s a n d s o i l n i t r o g e n 399 In China s Jianghan Plain, the major crops are rice and wheat, which are cultivated under different land use regimes. While wheat grows only in dry farmlands, rice requires ample water to flood the fields. So rice fields are also called water farmlands. These two different land use types usually cause differences in many soil properties such as the contents of soil nitrogen and organic matter. In this paper, we used AAPK to map the spatial distribution of soil TN contents in Hanchuan County, China through the integration of both soil sample data and information on land use type. The objectives are to: (1) investigate the effect of land use types on soil TN contents; (2) explore the capability of AAPK in improving the mapping accuracy of the spatial distribution of soil TN with land use type data as auxiliary information; and (3) compare the performance of AAPK, RK, and ordinary kriging (OK) in soil TN predictive mapping. The ultimate objective is to suggest a more appropriate soil nutrient mapping method for precision farming and environmental management. Study Area and Land Use METHODS AND MATERIALS The study was conducted in Hanchuan County, an agricultural region in central China. The study area is located in the Jianghan Plain of central Hubei Province, bounded by the longitudes and East, and the latitudes and North, with an area of 1659 km 2. It belongs to the northern subtropical monsoonal climate zone, with a temperate-humid climate throughout the year and four distinct seasons. The average annual temperature is 16.1 C and the mean annual precipitation is approximately 1198 mm. The land use data were collected from the county s local agricultural department. In this region, the major land use types are paddy fields and dry farmlands. We further edited the land use map and simply classified the land use into four types (or categories) paddy fields, dry farmlands, water bodies, and other land use types. The spatial distribution of the land use types is shown in Figure 1. Paddy fields (i.e., rice farmlands) and dry farmlands use very different forms of cultivation, which usually have strong impacts on soil properties. Other land use types include all locations without tillage, such as woodlands, wild lands, and village and town areas. Water bodies include rivers, lakes, and ponds. This land use type was first separated as a GIS layer and then overlaid on the kriged maps of soil TN, so it was not considered in the interpolation processes of soil TN by different kriging methods. Paddy fields are the dominant land use type of the county, which is one of the most developed cropproducing bases in Hubei Province. Soil Sampling and Lab Analysis In present study, the topsoil sampling points consist of prediction points (n = 402) and validation points (n = 135). Among the 402 prediction samples, 215 samples were taken from paddy fields, 130 samples from dry farmlands, and the remaining samples were taken from other land use types (Table 1). The validation samples contain all of the three land use groups presented in the prediction samples, although they were collected with consideration of randomness and homogeneity in the area. All samples
4 400 q u e t a l. Fig. 1. Soil sample locations and distribution of land use types. were taken in fall after harvest and before the next cropping season in order to avoid the effect of fertilization during crop cultivation. When sampling, soils were collected in the top layers (0 15 cm), 6 8 points at each site over an area of approximately 0.01 ha, and then mixed and divided into parts of 1 2 kg each, which were delivered to a laboratory for analysis. Sample locations were recorded using a hand-held global position system (GPS). All samples were air-dried at room temperature (20 22 C). After stones or other debris were removed, samples were then sieved to ensure the soil particles were less than 2 mm in diameter. Portions of each sample (about 100 g) were ground in an agate grinder and sieved through mm mesh. TN was determined using the Kjeldahl method with H 2 SO 4 + H 2 O 2 digestion (Kim, 2005). Variogram and Ordinary Kriging Geostatistics provides the methods to predict values at unsampled locations from values at sampled locations by taking into account the spatial correlation of sampled points. It can minimize the variance of estimation errors and investigation costs (Ferguson et al., 1998; Saito et al., 2005). The variogram the spatial measure for kriging is an effective tool for evaluating spatial variability (Boyer et al., 1991; Cahn
5 l a n d u s e t y p e s a n d s o i l n i t r o g e n 401 Table 1. Soil TN Content (g kg 1 ) Statistics a for Different Land Use Types Land use type N Range Minimum Maximum Mean SD Skew Kurt CV Total Paddy field Dry farmland Other land use type a SD = standard deviation; skew = skewness; kurt = kurtosis; CV = coefficient of variation (%). et al., 1994). A variogram describes the spatial autocorrelation structure of a continuous variable and provides some insight into possible factors that affect data distribution (Webster and Oliver, 1990). Spatial patterns of soil attributes following intrinsic stationarity can be described using the following experimental variogram: N( h) 1 γˆ ( h) = [ z( x (1) 2N( h) i ) ( x i + h) ] 2 i = 1 where N(h) is the number of data pairs separated by distance h, and z(x i ) and z(x i + h) are the measured values for regionalized variables Z(x) at the locations of x i and x i + h, respectively. In this study, no apparent anisotropy was found for any studied variable through experimental variograms. So all experimental variograms were in the isotropic form, and were fitted using basic mathematical models, such as the spherical and exponential models, for kriging interpolation. Ordinary kriging (OK) was chosen to create the spatial distribution map of soil TN contents, with the maximum search radius being set to the autocorrelation range of the corresponding variable. OK is expressed as a linearly weighted average of observations in the neighborhood of the unsampled location x 0 : n ỹ( x 0 ) = λ i y( x i ) i = 1 where ỹ( x 0 ) is the predicted value at location x 0, y(x i ) is the measured value for soil property at position x i, and λ i is the weight of the corresponding datum, obtained from the ordinary kriging system with the constraint n λ i = 1. Here n is the number of i = 1 sample data in the neighborhood. The interpolated grid had a resolution of 400 m 400 m. Readers may refer to Goovaerts (1997) for a more technical description of kriging and the variogram. Residual Kriging The spatial variability of the TN data is partially caused by the complex distribution of land use types in the study region, which increases the uncertainty of TN prediction. To reduce this uncertainty, in residual kriging (RK) we divided the TN content z(x kj ) of every sample into two parts the mean value within the land use type to which the sample belongs and the corresponding residual r(x kj ): (2)
6 402 q u e t a l. z( x kj ) = mt ( k ) + r( x kj ) (3) where x kj is the location of the sample z(x kj ) and t k is the land use type that x kj belongs to. Therefore, the variance σ z 2 of original of z(xkj ) can be given as the sum of the two components σ s 2 between different land use types and σr 2 within a land use type: σ z = σ s + σ r which shows the influence of land use types on spatial variability of TN and its variation within a land use type, respectively (Liu et al., 2006). Then, the residual r(x kj ) can be treated as a new stationary regionalized variable to be used in spatial interpolation using RK. In this study, the residual data were interpolated using OK and the variogram model of residual data. Area-And-Point Kriging In order to incorporate data with different spatial supports, 2 both point and area data must be considered simultaneously in the kriging system. In this study, point data refer to soil TN sampling data and area data refer to the mean TN contents of sample points within a specific land use type area (map unit). The arithmetic of area-and-point kriging is described in the following (Liu, 2007; Liu and Journel, 2009; Goovaerts, 2010). Consider the problem of estimating the value of a continuous attribute z at any location u within a study area A. The information available consists of a set of point data collected at n discrete locations u α {z(u α ); α = 1,..., n}, supplemented by a set of K areal data {z(v k ); k = 1,..., K} recorded for mapping units v k of various sizes and shapes. Both point and areal data can be simultaneously incorporated into the prediction using the area-and-point kriging (AAPK) estimator, defined as (4) n( u) z AAPK ( u) = λ α ( u)z( u α ) + α = 1 n( u) + K λ k ( u)z( v k ) k = n( u) + 1 (5) where n(u) and K are the number of the surrounding point and areal data, respectively. Point observations are typically selected based on their distance to the interpolation node u while areal data are chosen according to adjacency rules; for example, all polygons adjacent to the polygon including u are used in the estimation. The kriging weights are the solution of the following ordinary kriging system n( u) + k λ j ( u)cx ( i, x j ) + μ( u) = Cx ( i, u), i = 1,, n( u) + K j = 1 n( u) + k λ j ( u) = 1 j = 1 (6) 2 Spatial support here refers to the area represented by a datum, which may be a point, block, pixel, or polygon.
7 l a n d u s e t y p e s a n d s o i l n i t r o g e n 403 where μ(u) is the Lagrange multiplier, x i = u i if i n(u) and x i = v i otherwise. The quantity C( x i, x j ) represents a point-to-point, point-to-block, or block-to-block covariance depending on the indices i and j. As in traditional block kriging, the block-to-point covariance is approximated by the average of the point support covariance C(h) computed between the location u and a set of P k points discretizing the block v k. A similar procedure is used for the block-to-block covariance C( v k, v k ' ) = Cov{ z( v k ), z( v k' )} and involves averaging C(h) computed between any two points that discretize the blocks v k and v k'. The prediction variance associated with the AAPK estimator is computed as 2 σ AAPK n( u) n( u) + K ( u) = C( 0) λ α ( u)c( u α, u) λ k ( u)c ( v k, u) μ( u) α = 1 Interpolation Method Evaluation k = n( u) + 1 Kriging prediction variance can provide information concerning the confidence associated with the kriging estimates, and interpreting such maps is an important step toward quantifying reliability in spatial estimation (Olea, 1999). The larger the prediction variance, the lower is the reliability of the estimate. In this study, we compared the reliability in TN estimation by the methods of AAPK, RK, and OK based on maps of prediction variance. We also compared the Pearson s correlation coefficient (r) between the predicted and measured values of the 135 validation samples as well as the mean error (ME) and root mean square error (RMSE) for the three methods. The ME and RMSE are calculated by the following equations: (7) ME = N { zx ( N i ) z ( x i )} i = 1 (8) RMSE = N zx ( N { i ) z ( x i )} 2 i = 1 where N is the number of validation points, and z*(x i ) and z(x i ) are the measured and predicted values of the validation points, respectively. The ME provides a measure of bias, and the RMSE provides a measure of accuracy. Greater r and lower ME and RMSE values indicate higher prediction accuracy. In this study, Software version 2.1 of S-GeMS (Remy et al., 2009) was used to perform all of the geostatistical computations. Statistical analyses were conducted in SPSS 13.0 for Windows (SPSS, Chicago). Statistics of Soil TN Data RESULTS AND DISCUSSION The descriptive statistics for sampled TN contents are given in Table 1. The mean TN contents of all samples is 1.46 g kg 1. Among the 402 samples, the lowest and (9)
8 404 q u e t a l. Fig. 2. Histograms of (A) original soil TN data and (B) residuals after removing the mean values for each land use type. highest TN contents are 0.31 and 3.04 g kg 1, respectively. The coefficient of variation (CV) of all samples is 37.14%. When all samples are classified into three land use groups, their average TN contents in descending order are 1.62 g kg 1 for paddy fields, 1.38 g kg 1 for dry farmlands and 1.25 g kg 1 for the other land use type. This means that land use types indeed affect the contents of soil TN. Among the CV values of the three land use types, those associated with other land use type and paddy fields are the highest (41.80%) and lowest (33.82%), respectively. Both the original TN content data and the residual data (i.e., after the mean values for different land use types are removed from sample data) fit approximately to normal distributions (Fig. 2). The TN contents are found to be correlated with land use types. Such a characteristic is consistent with the previous studies (Wang et al., 2009). In general, paddy fields have higher TN contents as a result of more mineral and organic matter input and slow decomposition of organic matter, whereas dry farmlands have lower TN contents because of less fertilizer input, fast organic matter decomposition, and not being regarded as the primary land use type for agricultural production in this county. This result implies that categorical land use information should be a valuable input in kriging interpolation to improve spatial prediction accuracy. Spatial Prediction of Soil TN The variogram provides a description of the spatial autocorrelation structure of the variable under study and some insights into possible processes affecting the spatial distribution of the variable (Paz González et al., 2001). Experimental variograms and fitted models for original TN content data and the residual data are presented in Figure 3. The parameters indicate that the experimental variogram of the original TN data (for use in OK and AAPK) is well fitted by an exponential model, and that of the residuals (for use in RK) is a good fit to a spherical model. Because of the removal
9 l a n d u s e t y p e s a n d s o i l n i t r o g e n 405 Fig. 3. Experimental variograms of (A) original data and (B) residuals for soil TN with fitted models.
10 406 q u e t a l. of mean values, the experimental variogram of the residuals is quite different from that of the original TN data (Fig. 3). The sill and range parameters for the residuals are both smaller than those of the original TN data, because the structural variances derived from land use patterns are eliminated with the removed mean values of mapping units. The C 0 /(C 0 +C) ratio is usually used as a criterion to define the spatial autodependence of a variable. Ratio values lower than 25% and higher than 75% correspond to strong and weak spatial dependencies, respectively, while ratio values between 25% and 75% indicate moderate spatial dependence (Cambardella et al., 1994). Generally, strong spatial dependence of soil properties may be attributed to intrinsic factors and weak spatial dependence may be attributed to extrinsic factors (Cambardella et al., 1994). The C 0 /(C 0 +C) ratio of the variogram model for original TN data is 48.7%, exhibiting moderate spatial auto-dependency which may be attributed to both intrinsic factors such as other soil properties and extrinsic factors such as human activities (e.g., land use). The TN spatial distribution maps interpolated using AAPK, RK, and OK are presented in Figure 4. All maps show similar general trends, with higher TN values appearing in the northwest region and lower values mainly occurring in the mid-south region of the county. Among the three predicted distribution maps of soil TN, the AAPK- and OK-interpolated maps are obviously smoother than the RK-interpolated map. This smoothing effect is a widely known characteristic of kriging interpolation that may cause low values to be overestimated and high values to be underestimated (Lark and Webster, 2006). Although the AAPK interpolated map is smoother than others, it may more accurately display the lower or higher contents as expected on different land use types, in particular when compared with the OK-interpolated map. For example, at the top-left corner of the study area where no sample data exist, while the land use type implies a lower TN content, the AAPK-interpolated map indeed shows a lower value but the OK interpolated map indicates a higher value (see Figs. 1 and Fig. 4). The-RK interpolated map is the least smooth one among the three maps. However, it looks rather messy, and there is no guarantee that the coherence constraint is honored once the kriged residuals are added to the trend estimates. Such a constraint can only be imposed through the joint incorporation of field point data and areal map data as done in AAPK (Goovaerts, 2010). From Figure 4, it also can be seen that the differences between the general spatial patterns predicted by the three methods are relatively minor. This means that the advantages of AAPK over RK and OK are not much evident in this case study. One reason may be that the effect of land use types on soil TN contents is not sufficiently large in this case (but it can be large for the contents of some available nutrients). The other reason, we speculate, should be that some of the land use polygons are too large and complex in shape due to the spatial connection of many farmland pieces under the same land use types. 3 For this situation, segmenting a complex-shaped large land use polygon into simple-shaped small areas might be helpful to improving the effectiveness of AAPK. 3 For example, numerous pieces of rice fields divided by narrow ridges are not spatially separated by other land use types, thus forming complex-shaped large polygons. See Figure 1 for this situation.
11 l a n d u s e t y p e s a n d s o i l n i t r o g e n 407 Fig. 4. Soil TN distribution maps generated by (A) AAPK, (B) RK, and (C) OK. Variance and Accuracy Analysis The prediction variance of RK is computed as the sum of the OK variance for residuals and the variance associated with the trend model. From the kriging prediction variance maps (Fig. 5), one can find that the prediction variance of soil TN by AAPK is apparently smaller than that by RK and OK throughout the study area. The mean kriging prediction variances of TN by AAPK, RK, and OK are , , and , respectively. While the kriging variance is a measure of confidence in estimates and the sampling configuration only, and it is independent of the data values (Goovaerts, 1997), for different kriging methods a lower kriging variance may imply the kriging method can globally estimate more accurately (or with more confidence) than the kriging methods with higher kriging variance values. From this point, RK seems the worst estimator here. The correlation coefficients r, the ME values, and the RMSE values between the observed and predicted data at the 135 validation locations for the three kriging options are shown in Figure 6. The r values are 0.81, 0.67, and 0.55 for AAPK, RK, and OK, respectively. The ME and RMSE values are and 0.29 g kg 1 for AAPK, and 0.36 g kg 1 for RK, and and 0.41 g kg 1 for OK, respectively. The larger
12 408 q u e t a l. Fig. 5. Maps of prediction variances generated by (A) AAPK, (B) RK, and (C) OK. r value and the smaller ME and RMSE values for AAPK means it generates a more accurate prediction map than RK does. OK is apparently the least accurate estimator here. CONCLUSION Land use types were found to have some impacts on soil TN contents in the study area, because the mean TN content in paddy fields is apparently higher than that in dry farmlands. The recently emerged AAPK was used to interpolate soil TN so that the influence of land use types can be incorporated. Interpolated results showed that the soil TN contents were generally lower in the mid-south region and higher in the northwest region in the county. Comparison with OK and RK was also conducted to evaluate the merits of the AAPK method. The predictive map of soil TN contents by AAPK with the incorporation of land use data indicated an improvement in prediction accuracy over those by RK and OK, as represented by prediction variance, correlation coefficient, ME, and RMSE values, although the differences were not visually very significant in prediction maps of soil TN contents. It seems that the occurrence of complex-shaped large polygons of land use types impacted the effectiveness of
13 l a n d u s e t y p e s a n d s o i l n i t r o g e n 409 Fig. 6. Scatterplots between observed soil TN data and predicted data for validation using different interpolation methods: (A) AAPK; (B) RK; (C) OK.
14 410 q u e t a l. AAPK in improving prediction accuracy. Segmenting those large complex polygons into smaller simple-shaped areas might help effectively demonstrate the capability of AAPK and might be considered in further study. In general, land use types indeed have important impacts on the spatial distribution of soil TN, and incorporation of such categorical information is necessary for more accurately mapping soil properties. It is also concluded that AAPK is an effective prediction method for increasing the interpolation accuracy of soil TN through integration of both point and area data. Such a conclusion should also be applicable to other soil properties. The interpolated maps of soil TN contents may provide useful information for effective management of soil nutrients and the environment. ACKNOWLEDGMENTS Support from the National Natural Science Foundation of China under Grants and is greatly appreciated. REFERENCES Boyer, D. G., Wright, R. J., Feldhake, C. M., and D. P. Bligh, 1991, Soil Spatial Variability in Steeply Sloping Acid Soil Environment, Soil Science, 161: Burgess, T. M. and R. Webster, 1980, Optimal Interpolation and Isarithmic Mapping of Soil Properties: I. The Semivariogram and Punctual Kriging, Journal of Soil Science, 31: Cahn, M. D., Hummel, J. W., and B. H. Brouer, 1994, Spatial Analysis of Soil Fertility for Site-Specific Crop Management, Soil Science Society of American Journal, 58: Cambardella, C. A., Moorman, T. B., Nocak, J. M., Parkin, T. B., Karlen, D. L., Turco, R. F., and A. E. Konopka, 1994, Field-Scale Variability of Soil Properties in Central Iowa Soils, Soil Science Society of America Journal, 58: Carpenter, S. R., Caraco, N. F., Correll, D. L., Howarth, R. W., Sharpley, A. N., and V. H. Smith, 1998, Non-point Pollution of Surface Waters with Phosphorus and Nitrogen, Ecological Applications, 8: Ferguson, C. C., Darmendrail, D., Freier, K., Jensen, B. K., Jensen, J., Kasamas, H., Urzelai, A., and J. Vegter, 1998, Better Methods for Risk Assessment, in Risk Assessment for Contaminated Sites in Europe, Vol. 1. Scientific Basis, Nottingham, UK: LQM Press, Goovaerts, P., 1997, Geostatistics for Natural Resources Evaluation, New York, NY: Oxford University Press. Goovaerts, P., 2008, Kriging and Semivariogram Deconvolution in Presence of Irregular Geographical Units, Mathematical Geosciences, 40: Goovaerts, P., 2010, Combining Areal and Point Data in Geostatistical Interpolation: Applications to Soil Science and Medical Geography, Mathematical Geosciences, 42: Goovaerts, P. and A. G. Journel, 1995, Integrating Soil Map Information in Modeling the Spatial Variation of Continuous Soil Properties, European Journal of Soil Science, 46:
15 l a n d u s e t y p e s a n d s o i l n i t r o g e n 411 Hengl, T., Heuvelink, G. B. M., and A. Stein, 2004, A Generic Framework for Spatial Prediction of Soil Variables Based on Regression-Kriging, Geoderma, 120: Kim, H. T., 2005, Soil Sampling, Preparation, and Analysis, Boca Raton, FL: CRC, Lark, R. M. and R. Webster, 2006, Geostatistical Mapping of Geomorphic Variables in the Presence of Trend, Earth Surface Processes and Landforms, 31: Li, Y. and J. B. Zhang, 1999, Agricultural Diffuse Pollution from Fertilizers and Pesticides in China, Water Science and Technology, 39: Liu, T. L., Juang, K. W., and D. Y. Lee, 2006, Interpolating Soil Properties using Kriging Combined with Categorical Information of Soil Maps, Soil Science Society of American Journal, 70: Liu, Y., 2007, Geostatistical Integration of Linear Coarse Scale Data and Fine Scale Data, Ph.D. dissertation, Stanford University. Liu, Y. and A. G. Journel, 2009, A Package for Geostatistical Integration of Coarse and Fine Scale Data, Computers and Geosciences, 35: Lu, P., Su, Y., Niu, Z., and J. Wu, 2007, Geostatistical Analysis and Risk Assessment on Soil Total Nitrogen and Total Soil Phosphorus in the Dongting Lake Plain Area, China, Journal of Environmental Quality, 36: Odeh, I. O. A., McBratney, A. B., and D. J. Chittleborough, 1995, Further Results on Prediction from Terrain Attributes: Heterotopic Cokriging and Regression- Kriging, Geoderma, 67: Olea, R. A., 1999, Geostatistics for Engineers and Earth Scientists, New York, NY: Kluwer Academic Publishers. Paz González, A., Taboada Castro, M. T., and S. R. Vieira, 2001, Geostatistical Analy sis of Heavy Metals in a One-Hectare Plot under Natural Vegetation in a Serpentine Area, Canadian Journal of Soil Science, 81: Remy, N., Boucher, A., and J. Wu, 2009, Applied Geostatistics with SGeMS: A User s Guide, New York, NY: Cambridge University Press. Saito, H., McKenna, A., Zimmerman, D. A., and T. C. Coburn, 2005, Geostatistical Interpolation of Object Counts Collected from Multiple Strip Transects: Ordinary Kriging versus Finite Domain Kriging, Stochastic Environmental Research and Risk Assessment, 19: Smith, K. A., Jackson, D. R. and T. J. Pepper, 2001, Nutrient Losses by Surface Runoff Following the Application of Organic Manures to Arable Land: 1. Nitrogen, Environment Pollution, 112: Wang, Y. Q., Zhang, X. C., and C. Q. Huang, 2009, Spatial Variability of Soil Total Nitrogen and Soil Total Phosphorus under Different Land Uses in a Small Watershed on the Loess Plateau, China, Geoderma, 150: Webster, R. and M. A. Oliver, 1990, Statistical Methods in Soil and Land Resource Survey, London, UK: Oxford University Press. Wu, C., Wu, J., Luo, Y., Zhang, L., and S. D. DeGloria, 2009, Spatial Estimation of Soil Total Nitrogen Using Cokriging with Predicted Soil Organic Matter Content, Soil Science Society of America Journal, 73: Zhang, Z., Yu, D., Shi, X., Warner, E., Ren, H., Sun, W., Tan, M, and H. Wang, 2010, Application of Categorical Information in the Spatial Prediction of Soil Organic Carbon in the Red Soil Area of China, Soil Science and Plant Nutrition, 56:
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationEcological Informatics
Ecological Informatics 13 (2013) 99 105 Contents lists available at SciVerse ScienceDirect Ecological Informatics journal homepage: www.elsevier.com/locate/ecolinf Assessing the risk costs in delineating
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 informationASPECTS REGARDING THE USEFULNESS OF GEOGRAPHICALLY WEIGHTED REGRESSION (GWR) FOR DIGITAL MAPPING OF SOIL PARAMETERS
Lucrări Ştiinţifice vol. 52, seria Agronomie ASPECTS REGARDING THE USEFULNESS OF GEOGRAPHICALLY WEIGHTED REGRESSION (GWR) FOR DIGITAL MAPPING OF SOIL PARAMETERS C. PATRICHE 1, I. VASILINIUC 2 1 Romanian
More informationThe Study of Soil Fertility Spatial Variation Feature Based on GIS and Data Mining *
The Study of Soil Fertility Spatial Variation Feature Based on GIS and Data Mining * Chunan Li, Guifen Chen **, Guangwei Zeng, and Jiao Ye College of Information and Technology, Jilin Agricultural University,
More informationTransiogram: A spatial relationship measure for categorical data
International Journal of Geographical Information Science Vol. 20, No. 6, July 2006, 693 699 Technical Note Transiogram: A spatial relationship measure for categorical data WEIDONG LI* Department of Geography,
More informationSpatial variability of soil organic matter and nutrients in paddy fields at various scales in southeast China
Environ Geol (28) 53:1139 1147 DOI 1.17/s254-7-91-8 ORIGINAL ARTICLE Spatial variability of soil organic matter and nutrients in paddy fields at various scales in southeast China Xingmei Liu Æ Keli Zhao
More informationInvestigation of Monthly Pan Evaporation in Turkey with Geostatistical Technique
Investigation of Monthly Pan Evaporation in Turkey with Geostatistical Technique Hatice Çitakoğlu 1, Murat Çobaner 1, Tefaruk Haktanir 1, 1 Department of Civil Engineering, Erciyes University, Kayseri,
More informationABSTRACT INTRODUCTION. J. Soil Sci. Soc. Sri Lanka, Vol. 23, 2011
J. Soil Sci. Soc. Sri Lanka, Vol. 23, 2011 ORIGINAL PAPER ACCEPTED FOR PUBLICATION IN JOURNAL OF SOIL SCIENCE SOCIETY OF SRI LANKA SPATIAL VARIABILITY OF SOIL TEXTURE, ORGANIC CARBON AND CATION EXCHANGE
More informationGLOBAL SYMPOSIUM ON SOIL ORGANIC CARBON, Rome, Italy, March 2017
GLOBAL SYMPOSIUM ON SOIL ORGANIC CARBON, Rome, Italy, 2-23 March 207 Spatial Distribution of Organic Matter Content in Plough Layer of Balikh Flood Plain Soils Northeastern of Syria Hussam H. M. Husein
More informationSoil Moisture Modeling using Geostatistical Techniques at the O Neal Ecological Reserve, Idaho
Final Report: Forecasting Rangeland Condition with GIS in Southeastern Idaho Soil Moisture Modeling using Geostatistical Techniques at the O Neal Ecological Reserve, Idaho Jacob T. Tibbitts, Idaho State
More informationKriging method for estimation of groundwater resources in a basin with scarce monitoring data
36 New Approaches to Hydrological Prediction in Data-sparse Regions (Proc. of Symposium HS. at the Joint IAHS & IAH Convention, Hyderabad, India, September 9). IAHS Publ. 333, 9. Kriging method for estimation
More informationSpatiotemporal Analysis of Environmental Radiation in Korea
WM 0 Conference, February 25 - March, 200, Tucson, AZ Spatiotemporal Analysis of Environmental Radiation in Korea J.Y. Kim, B.C. Lee FNC Technology Co., Ltd. Main Bldg. 56, Seoul National University Research
More informationSpatial Distribution and Uncertainty Assessment of Potential Ecological Risks of Heavy Metals in Soil Using Sequential Gaussian Simulation
This article was downloaded by: [University of Connecticut] On: 27 January 2014, At: 06:59 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
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 informationSimulating the spatial distribution of clay layer occurrence depth in alluvial soils with a Markov chain geostatistical approach
ENVIRONMETRICS Environmetrics 2010; 21: 21 32 Published online 27 March 2009 in Wiley InterScience (www.interscience.wiley.com).981 Simulating the spatial distribution of clay layer occurrence depth in
More informationGeostatistics: Kriging
Geostatistics: Kriging 8.10.2015 Konetekniikka 1, Otakaari 4, 150 10-12 Rangsima Sunila, D.Sc. Background What is Geostatitics Concepts Variogram: experimental, theoretical Anisotropy, Isotropy Lag, Sill,
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 informationComparison of Geographically Weighted Regression and Regression Kriging for Estimating the Spatial Distribution of Soil Organic Matter
Comparison of Geographically Weighted Regression and Regression Kriging for Estimating the Spatial Distribution of Soil Organic Matter Ku Wang Department of Geographical Science, Minjiang University, 350108,
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 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 informationPedometric Techniques in Spatialisation of Soil Properties for Agricultural Land Evaluation
Bulletin UASVM Agriculture, 67(1)/2010 Print ISSN 1843-5246; Electronic ISSN 1843-5386 Pedometric Techniques in Spatialisation of Soil Properties for Agricultural Land Evaluation Iuliana Cornelia TANASĂ
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 informationModeling Spatial Pattern of Precipitation with GIS and Multivariate Geostatistical Methods in Chongqing Tobacco Planting Region, China
Modeling Spatial Pattern of Precipitation with GIS and Multivariate Geostatistical Methods in Chongqing Tobacco Planting Region, China Xuan Wang,2, Jiake Lv,2,*, Chaofu Wei 2, and Deti Xie 2 College of
More informationSimulating Spatial Distributions, Variability and Uncertainty of Soil Arsenic by Geostatistical Simulations in Geographic Information Systems
Open Environmental Sciences, 2008, 2, 26-33 26 Open Access Simulating Spatial Distributions, Variability and Uncertainty of Soil Arsenic by Geostatistical Simulations in Geographic Information Systems
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 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 informationRegional-scale modelling of the spatial distribution of surface and subsurface textural classes in alluvial soils using Markov chain geostatistics
Soil Use and Management, September 2008, 24, 263 272 doi: /j.475-2743.2008.0065.x Regional-scale modelling of the spatial distribution of surface and subsurface textural classes in alluvial soils using
More informationAssessment of Three Spatial Interpolation Models to Obtain the Best One for Cumulative Rainfall Estimation (Case study: Ramsar District)
Assessment of Three Spatial Interpolation Models to Obtain the Best One for Cumulative Rainfall Estimation (Case study: Ramsar District) Hasan Zabihi, Anuar Ahmad, Mohamad Nor Said Department of Geoinformation,
More informationSpatial Variability of Selected Chemical Characteristics of Paddy Soils in Sawah Sempadan, Selangor, Malaysia
Malaysian Journal of Soil Science Vol. 14: 27-39 (2010) ISSN: 1394-7990 Malaysian Society of Soil Science Spatial Variability of Selected Chemical Characteristics of Paddy Soils in Sawah Sempadan, Selangor,
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 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 informationTypes of spatial data. The Nature of Geographic Data. Types of spatial data. Spatial Autocorrelation. Continuous spatial data: geostatistics
The Nature of Geographic Data Types of spatial data Continuous spatial data: geostatistics Samples may be taken at intervals, but the spatial process is continuous e.g. soil quality Discrete data Irregular:
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 informationTypes of Spatial Data
Spatial Data Types of Spatial Data Point pattern Point referenced geostatistical Block referenced Raster / lattice / grid Vector / polygon Point Pattern Data Interested in the location of points, not their
More informationSPATIAL VARIABILITY OF AVAILABLE NUTRIENTS AND SOIL CARBON UNDER ARABLE CROPPING IN CANTERBURY
SPATIAL VARIABILITY OF AVAILABLE NUTRIENTS AND SOIL CARBON UNDER ARABLE CROPPING IN CANTERBURY Weiwen Qiu, Denis Curtin and Mike Beare The New Zealand Institute for Plant & Food Research Limited, Private
More informationAES BIOFLUX Advances in Environmental Sciences - International Journal of the Bioflux Society
AES BIOFLUX Advances in Environmental Sciences - International Journal of the Bioflux Society Providing a soil fertility map using geographic information system, geostatistical techniques and fuzzy logic
More informationImproving Spatial Data Interoperability
Improving Spatial Data Interoperability A Framework for Geostatistical Support-To To-Support Interpolation Michael F. Goodchild, Phaedon C. Kyriakidis, Philipp Schneider, Matt Rice, Qingfeng Guan, Jordan
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 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 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 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 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 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 informationGeostatistical analysis of surface soil texture from Zala county in western Hungary
Geostatistical analysis of surface soil texture from Zala county in western Hungary K. Adhikari *,**, A. Guadagnini **, G. Toth * and T. Hermann *** * Land Management and Natural Hazards Unit, Institute
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 informationOn dealing with spatially correlated residuals in remote sensing and GIS
On dealing with spatially correlated residuals in remote sensing and GIS Nicholas A. S. Hamm 1, Peter M. Atkinson and Edward J. Milton 3 School of Geography University of Southampton Southampton SO17 3AT
More informationRelationships between Soil salinity and geopedological units in Saveh plain, Iran
Available online at www.scholarsresearchlibrary.com Annals of Biological Research, 2012, 3 (5):2292-2296 (http://scholarsresearchlibrary.com/archive.html) ISSN 0976-1233 CODEN (USA): ABRNBW Relationships
More informationGEOINFORMATICS Vol. II - Stochastic Modelling of Spatio-Temporal Phenomena in Earth Sciences - Soares, A.
STOCHASTIC MODELLING OF SPATIOTEMPORAL PHENOMENA IN EARTH SCIENCES Soares, A. CMRP Instituto Superior Técnico, University of Lisbon. Portugal Keywords Spacetime models, geostatistics, stochastic simulation
More informationSPATIAL VARIABILITY OF NITROGEN SUPPLY ASSESSED USING SOIL AND PLANT BIOASSAYS
SPATIAL VARIABILITY OF NITROGEN SUPPLY ASSESSED USING SOIL AND PLANT BIOASSAYS Weiwen Qiu 1, Paul Johnstone 2, Dirk Wallace 1, Nathan Arnold 2, Bruce Searle 2, Joanna Sharp 1, Mike Beare 1 and Denis Curtin
More informationBasics in Geostatistics 2 Geostatistical interpolation/estimation: Kriging methods. Hans Wackernagel. MINES ParisTech.
Basics in Geostatistics 2 Geostatistical interpolation/estimation: Kriging methods Hans Wackernagel MINES ParisTech NERSC April 2013 http://hans.wackernagel.free.fr Basic concepts Geostatistics Hans Wackernagel
More informationResearch Article County-Scale Spatial Variability of Macronutrient Availability Ratios in Paddy Soils
Applied and nvironmental oil cience, Article ID 689482, 1 pages http://dx.doi.org/1.1155/214/689482 Research Article County-cale patial Variability of Macronutrient Availability Ratios in Paddy oils Mingkai
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 informationGEOSTATISTICAL ANALYSIS OF SPATIAL DATA. Goovaerts, P. Biomedware, Inc. and PGeostat, LLC, Ann Arbor, Michigan, USA
GEOSTATISTICAL ANALYSIS OF SPATIAL DATA Goovaerts, P. Biomedware, Inc. and PGeostat, LLC, Ann Arbor, Michigan, USA Keywords: Semivariogram, kriging, spatial patterns, simulation, risk assessment Contents
More informationBayesian Markov Chain Random Field Cosimulation for Improving Land Cover Classification Accuracy
Math Geosci (2015) 47: 148 DOI 10.1007/s11004-014-9553-y Bayesian Markov Chain Random Field Cosimulation for Improving Land Cover Classification Accuracy Weidong Li Chuanrong Zhang Michael R. Willig Dipak
More informationA Geostatistical Approach to Predict the Average Annual Rainfall of Bangladesh
Journal of Data Science 14(2016), 149-166 A Geostatistical Approach to Predict the Average Annual Rainfall of Bangladesh Mohammad Samsul Alam 1 and Syed Shahadat Hossain 1 1 Institute of Statistical Research
More informationGIS
sanaein@gmail.com Astaraei @ferdowsi.um.ac.ir mghaemi @gmail.com negarsiabi @gmail.com. GIS ok ( D-K) RK CO-k RS GIS lesch Corwin Isaaks Srivastava Hengl Laslettet Webster Burgess Goovaerts kriging cokriging.
More informationSurvey on Application of Geostatistical Methods for Estimation of Rainfall in Arid and Semiarid Regions in South West Of Iran
Survey on Application of Geostatistical Methods for Estimation of Rainfall in Arid and Semiarid Regions in South West Of Iran Ebrahim Hosseini Chegini Mohammad Hossein Mahdian Sima Rahimy Bandarabadi Mohammad
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 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 informationVisualizing Spatial Uncertainty of Multinomial Classes in Area-class Mapping
Visualizing Spatial Uncertainty of Multinomial Classes in Area-class Mapping Weidong Li Chuanrong Zhang* Department of Geography, Kent State University, Kent, OH 44242 *Chuanrong Zhang Email: zhangc@uww.edu
More informationANALYSIS OF DEPTH-AREA-DURATION CURVES OF RAINFALL IN SEMIARID AND ARID REGIONS USING GEOSTATISTICAL METHODS: SIRJAN KAFEH NAMAK WATERSHED, IRAN
JOURNAL OF ENVIRONMENTAL HYDROLOGY The Electronic Journal of the International Association for Environmental Hydrology On the World Wide Web at http://www.hydroweb.com VOLUME 14 2006 ANALYSIS OF DEPTH-AREA-DURATION
More informationThis article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and
This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and education use, including for instruction at the authors institution
More informationGeostatistical Applications for Precision Agriculture
M.A. Oliver Editor Geostatistical Applications for Precision Agriculture Springer Contents 1 An Overview of Geostatistics and Precision Agriculture 1 M.A. Oliver 1.1 Introduction 1 1.1.1 A Brief History
More informationGeog 210C Spring 2011 Lab 6. Geostatistics in ArcMap
Geog 210C Spring 2011 Lab 6. Geostatistics in ArcMap Overview In this lab you will think critically about the functionality of spatial interpolation, improve your kriging skills, and learn how to use several
More informationComparison of Methods for Deriving a Digital Elevation Model from Contours and Modelling of the Associated Uncertainty
Comparison of Methods for Deriving a Digital Elevation Model from Contours and Modelling of the Associated Uncertainty Anicet Sindayihebura 1, Marc Van Meirvenne 1 and Stanislas Nsabimana 2 1 Department
More informationSAMPLING DESIGNS OVER TIME BASED ON SPATIAL VARIABILITY OF IMAGES FOR MAPPING AND MONITORING SOIL EROSION COVER FACTOR
SAMPLING DESIGNS OVER TIME BASED ON SPATIAL VARIABILITY OF IMAGES FOR MAPPING AND MONITORING SOIL EROSION COVER FACTOR Guangxing Wang W53 Turner Hall 112 S. Goodwin Ave. University of IL Urbana, IL 6181,
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 informationTraps for the Unwary Subsurface Geoscientist
Traps for the Unwary Subsurface Geoscientist ashley.francis@sorviodvnvm.co.uk http://www.sorviodvnvm.co.uk Presented at SEG Development & Production Forum, 24-29 th June 2001, Taos, New Mexico, USA 24-29
More informationMarkov Chain Modeling of Multinomial Land-Cover Classes
Markov Chain Modeling of Multinomial Land-Cover Classes Chuanrong Zhang Department of Geography and Geology, University of Wisconsin, Whitewater, Wisconsin 53190 Weidong Li Department of Geography, University
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 informationProgress Report Year 2, NAG5-6003: The Dynamics of a Semi-Arid Region in Response to Climate and Water-Use Policy
Progress Report Year 2, NAG5-6003: The Dynamics of a Semi-Arid Region in Response to Climate and Water-Use Policy Principal Investigator: Dr. John F. Mustard Department of Geological Sciences Brown University
More informationSpatio-temporal dynamics of the urban fringe landscapes
Spatio-temporal dynamics of the urban fringe landscapes Yulia Grinblat 1, 2 1 The Porter School of Environmental Studies, Tel Aviv University 2 Department of Geography and Human Environment, Tel Aviv University
More informationInterpolation and 3D Visualization of Geodata
Marek KULCZYCKI and Marcin LIGAS, Poland Key words: interpolation, kriging, real estate market analysis, property price index ABSRAC Data derived from property markets have spatial character, no doubt
More informationFluvial Variography: Characterizing Spatial Dependence on Stream Networks. Dale Zimmerman University of Iowa (joint work with Jay Ver Hoef, NOAA)
Fluvial Variography: Characterizing Spatial Dependence on Stream Networks Dale Zimmerman University of Iowa (joint work with Jay Ver Hoef, NOAA) March 5, 2015 Stream network data Flow Legend o 4.40-5.80
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 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 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 informationGeostatistical Analysis of Soil Properties to Support Spatial Sampling in a Paddy Growing Alfisol
Tropical Agricultural Research Vol. 22 (1): 34-44 (2010) Geostatistical Analysis of Soil Properties to Support Spatial Sampling in a Paddy Growing Alfisol V.G.D. Nayanaka, W.A.U. Vitharana 1* and R.B.
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 informationEstimating threshold-exceeding probability maps of environmental variables with Markov chain random fields
Stoch Environ Res Risk Assess (2010) 24:1113 1126 DOI 10.1007/s00477-010-0389-9 ORIGINAL PAPER Estimating threshold-exceeding probability maps of environmental variables with Markov chain random fields
More informationPrediction of Soil Properties Using Fuzzy Membership
Prediction of Soil Properties Using Fuzzy Membership Zhu, A.X. 1,2 ; Moore, A. 3 ; Burt, J. E. 2 1 State Key Lab of Resources and Environmental Information System, Institute of Geographical Sciences and
More informationBAYESIAN MODEL FOR SPATIAL DEPENDANCE AND PREDICTION OF TUBERCULOSIS
BAYESIAN MODEL FOR SPATIAL DEPENDANCE AND PREDICTION OF TUBERCULOSIS Srinivasan R and Venkatesan P Dept. of Statistics, National Institute for Research Tuberculosis, (Indian Council of Medical Research),
More informationSpatial analysis of soil properties using GIS based geostatistics models
Model. Earth Syst. Environ. (2016) 2:107 DOI 10.1007/s40808-016-0160-4 SHORT COMMUNICATION Spatial analysis of soil properties using GIS based geostatistics models Pravat Kumar Shit 1 Gouri Sankar Bhunia
More informationDROUGHT RISK EVALUATION USING REMOTE SENSING AND GIS : A CASE STUDY IN LOP BURI PROVINCE
DROUGHT RISK EVALUATION USING REMOTE SENSING AND GIS : A CASE STUDY IN LOP BURI PROVINCE K. Prathumchai, Kiyoshi Honda, Kaew Nualchawee Asian Centre for Research on Remote Sensing STAR Program, Asian Institute
More informationGeostatistical Analyst for Deciding Optimal Interpolation Strategies for Delineating Compact Zones
International Journal of Geosciences, 2011, 2, 585-596 doi:10.4236/ijg.2011.24061 Published Online November 2011 (http://www.scirp.org/journal/ijg) 585 Geostatistical Analyst for Deciding Optimal Interpolation
More informationOverview of Statistical Analysis of Spatial Data
Overview of Statistical Analysis of Spatial Data Geog 2C Introduction to Spatial Data Analysis Phaedon C. Kyriakidis www.geog.ucsb.edu/ phaedon Department of Geography University of California Santa Barbara
More informationThe Comparative Analysis of Spatial Structure of Ji Wheat 22 Yield Based on Different Stochastic Samplings
The Comparative Analysis of Spatial Structure of Ji Wheat 22 Yield Based on Different Stochastic Samplings Yujian Yang * and Xueqin Tong S & T Information Engineering Technology Center of Shandong Academy
More informationINVESTIGATING CLIMATE CHANGE IMPACTS ON SURFACE SOIL PROFILE TEMPERATURE (CASE STUDY: AHWAZ SW OF IRAN)
INVESTIGATING CLIMATE CHANGE IMPACTS ON SURFACE SOIL PROFILE TEMPERATURE (CASE STUDY: AHWAZ SW OF IRAN) Kazem Hemmadi 1, Fatemeh Zakerihosseini 2 ABSTRACT In arid and semi-arid regions, warming of soil
More informationThe Relationship between Vegetation Changes and Cut-offs in the Lower Yellow River Based on Satellite and Ground Data
Journal of Natural Disaster Science, Volume 27, Number 1, 2005, pp1-7 The Relationship between Vegetation Changes and Cut-offs in the Lower Yellow River Based on Satellite and Ground Data Xiufeng WANG
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 informationDEM-based Ecological Rainfall-Runoff Modelling in. Mountainous Area of Hong Kong
DEM-based Ecological Rainfall-Runoff Modelling in Mountainous Area of Hong Kong Qiming Zhou 1,2, Junyi Huang 1* 1 Department of Geography and Centre for Geo-computation Studies, Hong Kong Baptist University,
More informationAn internet based tool for land productivity evaluation in plot-level scale: the D-e-Meter system
An internet based tool for land productivity evaluation in plot-level scale: the D-e-Meter system Tamas Hermann 1, Ferenc Speiser 1 and Gergely Toth 2 1 University of Pannonia, Hungary, tamas.hermann@gmail.com
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 informationSTUDY AND COMPARISON OF SOME GEOSTATISTICAL METHODS FOR MAPPING CATION EXCHANGE CAPACITY (CEC) IN SOILS OF NORTHERN IRAN
1. Farrokh ASADZADEH, 1,2. Ali AKBARZADEH, 1. Ali Asghar ZOLFAGHARI, 1. Ruhollah TAGHIZADEH MEHRJARDI, 1. Mojtaba MEHRABANIAN, 3. Hadi RAHIMI LAKE, 4. Mohammad Ali SABETI AMIRHANDEH STUDY AND COMPARISON
More informationMETHODOLOGY WHICH APPLIES GEOSTATISTICS TECHNIQUES TO THE TOPOGRAPHICAL SURVEY
International Journal of Computer Science and Applications, 2008, Vol. 5, No. 3a, pp 67-79 Technomathematics Research Foundation METHODOLOGY WHICH APPLIES GEOSTATISTICS TECHNIQUES TO THE TOPOGRAPHICAL
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 informationSpatial Variability of Some Soil Properties in Paddy Fields (Case Study: Siyahkal, Guilan Province)
AGRICULTURAL COMMUNICATIONS, 2017, 5(1): 7 15. Spatial Variability of Some Soil Properties in Paddy Fields (Case Study: Siyahkal, Guilan Province) NAFISEH YAGHMAEIAN MAHABADI *1 AND ZAHRA AMIRI 2 1 Department
More informationAN OPERATIONAL DROUGHT MONITORING SYSTEM USING SPATIAL INTERPOLATION METHODS FOR PINIOS RIVER BASIN, GREECE
Proceedings of the 13 th International Conference on Environmental Science and Technology Athens, Greece, 5-7 September 2013 AN OPERATIONAL DROUGHT MONITORING SYSTEM USING SPATIAL INTERPOLATION METHODS
More informationUsing 7Be to Trace Temporal Variation of Interrill and Rill Erosion on Slopes
Available online at www.sciencedirect.com Procedia Environmental Sciences 11 (2011) 1220 1226 Using 7Be to Trace Temporal Variation of Interrill and Rill Erosion on Slopes Gang Liu 1, Qiong Zhang 3, Mingyi
More information