Spatial interpolation method based on integrated RBF neural networks for estimating heavy metals in soil of a mountain region

A novel spatial interpolation method based on integrated radial basis function artificial neural networks （IRBFANNs） is proposed to provide accurate and stable predictions of heavy metals concentrations in soil at un- sampled sites in a mountain region. The IRBFANNs hybridize the advantages of the artificial neural networks and the neural networks integration approach. Three experimental projects under different sampling densities are carried out to study the performance of the proposed IRBFANNs-based interpolation method. This novel method is compared with six peer spatial interpolation methods based on the root mean square error and visual evaluation of the distribution maps of Mn elements. The experimental results show that the proposed method performs better in accuracy and stability. Moreover, the proposed method can provide more details in the spatial distribution maps than the compared interpolation methods in the cases of sparse sampling density.

Spatial Interpolation of Daily Precipitation in China:1951-2005

Climate research relies heavily on good quality instrumental data; for modeling efforts gridded data are needed. So far, relatively little effort has been made to create gridded climate data for China. This is especially true for high-resolution daily data. This work, focuses on identifying an accurate method to produce gridded daily precipitation in China based on the observed data at 753 stations for the period 1951-2005. Five interpolation methods, including ordinary nearest neighbor, local polynomial, radial basis function, inverse distance weighting, and ordinary kriging, have been used and compared. Cross-validation shows that the ordinary kriging based on seasonal semi-variograms gives the best performance, closely followed by the inverse distance weighting with a power of 2. Finally the ordinary kriging is chosen to interpolate the station data to a 18 km× 18 km grid system covering the whole country. Precipitation for each 0.5°×0.5° latitude-longitude block is then obtained by averaging the values at the grid nodes within the block. Owing to the higher station density in the eastern part of the country, the interpolation errors are much smaller than those in the west （west of 100°E）. Excluding 145 stations in the western region, the daily, monthly, and annual relative mean absolute errors of the interpolation for the remaining 608 stations are 74%, 29%, and 16%, respectively. The interpolated daily precipitation has been made available on the internet for the scientific community.

A comparative study of spatial interpolation methods fordetermining fishery resources density in the Yellow Sea

Spatial interpolation is a common tool used in the study of fishery ecology, especially for the construction of ecosystem models. To develop an appropriate interpolation method of determining fishery resources density in the Yellow Sea, we tested four frequently used methods, including inverse distance weighted interpolation(IDW), global polynomial interpolation(GPI), local polynomial interpolation(LPI) and ordinary kriging(OK).A cross-validation diagnostic was used to analyze the efficacy of interpolation, and a visual examination was conducted to evaluate the spatial performance of the different methods. The results showed that the original data were not normally distributed. A log transformation was then used to make the data fit a normal distribution. During four survey periods, an exponential model was shown to be the best semivariogram model in August and October 2014, while data from January and May 2015 exhibited the pure nugget effect.Using a paired-samples t test, no significant differences(P>0.05) between predicted and observed data were found in all four of the interpolation methods during the four survey periods. Results of the cross-validation diagnostic demonstrated that OK performed the best in August 2014, while IDW performed better during the other three survey periods. The GPI and LPI methods had relatively poor interpolation results compared to IDW and OK. With respect to the spatial distribution, OK was balanced and was not as disconnected as IDW nor as overly smooth as GPI and LPI, although OK still produced a few 'bull's-eye' patterns in some areas.However, the degree of autocorrelation sometimes limits the application of OK. Thus, OK is highly recommended if data are spatially autocorrelated. With respect to feasibility and accuracy, we recommend IDW to be used as a routine interpolation method. IDW is more accurate than GPI and LPI and has a combination of desirable properties, such as easy accessibility and rapid processing.

Comparing different spatial interpolation methods to predict the distribution of fishes:A case study of Coilia nasus in the Changjiang River Estuary

Spatial-temporal distribution of marine fishes is strongly influenced by environmental factors.To obtain a more continuous distribution of these variables usually measured by stationary sampling designs,spatial interpolation methods(SIMs)is usually used.However,different SIMs may obtain varied estimation values with significant differences,thus affecting the prediction of fish spatial distribution.In this study,different SIMs were used to obtain continuous environmental variables(water depth,water temperature,salinity,dissolved oxygen(DO),p H,chlorophyll a and chemical oxygen demand(COD))in the Changjiang River Estuary(CRE),including inverse distance weighted(IDW)interpolation,ordinary Kriging(OK)(semivariogram model:exponential(OKE),Gaussian(OKG)and spherical(OKS))and radial basis function(RBF)(regularized spline function(RS)and tension spline function(TS)).The accuracy and effect of SIMs were cross-validated,and two-stage generalized additive model(GAM)was used to predict the distribution of Coilia nasus from 2012 to 2014 in CRE.DO and COD were removed before model prediction due to their autocorrelation coefficient based on variance inflation factors analysis.Results showed that the estimated values of environmental variables obtained by the different SIMs differed(i.e.,mean values,range etc.).Cross-validation revealed that the most suitable SIMs of water depth and chlorophyll a was IDW,water temperature and salinity was RS,and p H was OKG.Further,different interpolation results affected the predicted spatial distribution of Coilia nasus in the CRE.The mean values of the predicted abundance were similar,but the differences between and among the maximum value were large.Studies showed that different SIMs can affect estimated values of the environmental variables in the CRE(especially salinity).These variations further suggest that the most applicable SIMs to each variable will also differ.Thus,it is necessary to take these potential impacts into consideration when studying the relationship between the spatial distribution of fishes and environmental changes in the CRE.

The efficacy of Kriging spatial interpolation for filling temporal gaps in daily air temperature data: A case study in northeast China

Quality-controlled and serially complete daily air temperature data are essential to evaluating and modelling the influences of climate change on the permafrost in cold regions. Due to malfunctions and location changes of observing stations, temporal gaps （i.e., missing data） are common in collected datasets. The objective of this study was to assess the efficacy of Kriging spatial interpolation for estimating missing data to fill the temporal gaps in daily air temperature data in northeast China. A cross-validation experiment was conducted. Daily air temperature series from 1960 to 2012 at each station were estimated by using the universal Kriging （UK） and Kriging with an external drift （KED）, as appropriate, as if all the ob-servations at a given station were completely missing. The temporal and spatial variation patterns of estimation uncertainties were also checked. Results showed that Kriging spatial interpolation was generally desirable for estimating missing data in daily air temperature, and in this study KED performed slightly better than UK. At most stations the correlation coefficients （R2） between the observed and estimated daily series were 〉0.98, and root mean square errors （RMSEs） of the estimated daily mean （Tmean）, maximum （Tmax）, and minimum （Tmin） of air temperature were 〈3 ℃. However, the estimation quality was strongly affected by seasonality and had spatial variation. In general, estimation uncertainties were small in summer and large in winter. On average, the RMSE in winter was approximately 1 ℃ higher than that in summer. In addition, estimation uncertainties in mountainous areas with complex terrain were significantly larger than those in plain areas.

Spatial Interpolation Method of the Meteorological Factor in Northeast China Based on GIS

[ Objective ] The research aimed to study the best spatial interpolation method of the meteorological factor in Northeast China. [ Method ] Based on geostatistical analysis tool of the Arclnfo GIS software, several spatial interpolation methods were used to estimate the meteorological fac- tore （annual rainfall and monthly average temperature） in Northeast China, such as inverse distance weighted （IDW）, radial basis function （RBF） and Kriging. Then, the best interpolation method of one meteorological factor was selected. [ Result] For monthly average temperature, Kriging method was better than others. For annual rainfall, precision of the evaluated value with RBF method was higher than that of the IDW and Kriging methods. [Conclusion] There was obvious regional difference of the meteorological factor in Northeast China. Monthly average temperature in south was higher than that in north, and annual rainfall in southeast was more than that in northwest in Northeast China.

Comparative Study on Spatial Interpolation Methods for 3D Modeling of Coal and Rock Layers Based on Borehole Data

The patial interpolation of borehole data is an important means of stratigraphic structure to construct a three-dimensional model of coal strata,and the reasonable selection of an effective spatial interpolation method will directly affect the accuracy of three-dimensional modeling of the strata.To select an effective spatial interpolation method and improve the accuracy of 3D modeling of formations,four interpolation methods(the inverse distance weight interpolation algorithm,the local polynomial interpolation algorithm,the radial basis neural network interpolation algorithm and the kriging interpolation algorithm)were compared and analyzed.In particular,the methods of interpolation algorithm,interpolation surface,sample test error,and cross-validation error were used.The experiment of 13-1 seam coal in the Huainan mining area showed the spatial surface interpolation effect of the radial basis neural network interpolation algorithm(RBF)compared with the inverse distance weight interpolation algorithm(IDW),local polynomial interpolation algorithm(LPI)and kriging algorithm.The three interpolation methods have higher accuracy and are more suitable for surface interpolation of coal seams,which is of great significance for improving the accuracy of subsequent 3D modeling of coal seams.

Spatial Interpolation Applied a Crustal Thickness in Brazil

The use of spatial interpolation methods of data is becoming increasingly common in geophysical analysis, for that reason, currently, several software already contain many of these methods, allowing more detailed studies. In the present work four interpolation methods are evaluated, for the crustal thickness data of Brazil tectonic provinces, with the intention of making Moho’s map of the regions. The methods used were IDW, Natural Neighbor, Spline and Kriging. We compiled 257 data that constituted a geographic database implemented in the template Postgree PostGIS and were processed using the tools of interpolation located in the Spatyal Analyst Tools program ArcGIS?9 ESRI. Traditional methods, IDW, Natural Neighbor and Spline, generate artifacts in their results, the effects of aim, not consistent with the behavior of crust. Such anomalies are generated because of mathematical formulation methods added to data compiled gravimetry. The analysis results of geostatistical Kriging are more refined and consistent, showing no specific anormalities, i.e., the crustal thickness variation (thinning and thickening) is introduced gradually. Initial our estimates were separated in four specific blocks. With the approval of new networks (BRASIS, RSISNE and RSIS), the crustal thickness database for Brazil may be amended or supplemented so that new models may be generated more consistently, complementing studies of regional tectonics evolution and seismicity.

A MODIFIED EDGE-ORIENTED SPATIAL INTERPOLATION ALGORITHM FOR CONSECUTIVE BLOCKS ERROR CONCEALMENT

This paper proposes a low-complexity spatial-domain Error Concealment (EC) algorithm for recovering consecutive blocks error in still images or Intra-coded (I) frames of video sequences. The proposed algorithm works with the following steps. Firstly the Sobel operator is performed on the top and bottom adjacent pixels to detect the most likely edge direction of current block area. After that one-Dimensional (1D) matching is used on the available block boundaries. Displacement between edge direction candidate and most likely edge direction is taken into consideration as an important factor to improve stability of 1D boundary matching. Then the corrupted pixels are recovered by linear weighting interpolation along the estimated edge direction. Finally the interpolated values are merged to get last recovered picture. Simulation results demonstrate that the proposed algorithms obtain good subjective quality and higher Peak Signal-to-Noise Ratio (PSNR) than the methods in literatures for most images.

Non-parametric machine learning methods for interpolation of spatially varying non-stationary and non-Gaussian geotechnical properties

Spatial interpolation has been frequently encountered in earth sciences and engineering.A reasonable appraisal of subsurface heterogeneity plays a significant role in planning,risk assessment and decision making for geotechnical practice.Geostatistics is commonly used to interpolate spatially varying properties at un-sampled locations from scatter measurements.However,successful application of classic geostatistical models requires prior characterization of spatial auto-correlation structures,which poses a great challenge for unexperienced engineers,particularly when only limited measurements are available.Data-driven machine learning methods,such as radial basis function network(RBFN),require minimal human intervention and provide effective alternatives for spatial interpolation of non-stationary and non-Gaussian data,particularly when measurements are sparse.Conventional RBFN,however,is direction independent(i.e.isotropic)and cannot quantify prediction uncertainty in spatial interpolation.In this study,an ensemble RBFN method is proposed that not only allows geotechnical anisotropy to be properly incorporated,but also quantifies uncertainty in spatial interpolation.The proposed method is illustrated using numerical examples of cone penetration test(CPT)data,which involve interpolation of a 2D CPT cross-section from limited continuous 1D CPT soundings in the vertical direction.In addition,a comparative study is performed to benchmark the proposed ensemble RBFN with two other non-parametric data-driven approaches,namely,Multiple Point Statistics(MPS)and Bayesian Compressive Sensing(BCS).The results reveal that the proposed ensemble RBFN provides a better estimation of spatial patterns and associated prediction uncertainty at un-sampled locations when a reasonable amount of data is available as input.Moreover,the prediction accuracy of all the three methods improves as the number of measurements increases,and vice versa.It is also found that BCS prediction is less sensitive to the number of measurement data and outperforms RBFN and MPS when only limited point observations are available.

Reconstructing pre-erosion topography using spatial interpolation techniques： A validation-based approach

Understanding the topographic context preceding the development of erosive landforms is of major relevance in geomorphic research, as topography is an important factor on both water and mass movement-related erosion, and knowledge of the original surface is a condition for quantifying the volume of eroded material. Although any reconstruction implies assuming that the resulting surface reflects the original topography, past works have been dominated by linear interpolation methods, incapable of generating curved surfaces in areas with no data or values out- side the range of variation of inputs. In spite of these limitations, impossibility of validation has led to the assumption of surface representativity never being challenged. In this paper, a validation-based method is applied in order to define the optimal interpolation technique for reconstructing pre-erosion topography in a given study area. In spite of the absence of the original surface, different techniques can be nonetheless evaluated by quantifying their ca- pacity to reproduce known topography in unincised locations within the same geomorphic contexts of existing erosive landforms. A linear method （Triangulated Irregular Network, TIN） and 23 parameterizations of three distinct Spline interpolation techniques were compared using 50 test areas in a context of research on large gully dynamics in the South of Portugal. Results show that almost all Spline methods produced smaller errors than the TIN, and that the latter produced a mean absolute error 61.4% higher than the best Spline method, clearly establishing both the better adjustment of Splines to the geomorphic context considered and the limitations of linear approaches. The proposed method can easily be applied to different interpolation techniques and topographic contexts, enabling better calculations of eroded volumes and denudation rates as well as the investigation of controls by antecedent topographic form over erosive processes.

Comparison of Spatial Interpolation Methods for Gridded Bias Removal in Surface Temperature Forecasts

All numerical weather prediction（NWP） models inherently have substantial biases, especially in the forecast of near-surface weather variables. Statistical methods can be used to remove the systematic error based on historical bias data at observation stations. However, many end users of weather forecasts need bias corrected forecasts at locations that scarcely have any historical bias data. To circumvent this limitation, the bias of surface temperature forecasts on a regular grid covering Iran is removed, by using the information available at observation stations in the vicinity of any given grid point. To this end, the running mean error method is first used to correct the forecasts at observation stations, then four interpolation methods including inverse distance squared weighting with constant lapse rate（IDSW-CLR）, Kriging with constant lapse rate（Kriging-CLR）, gradient inverse distance squared with linear lapse rate（GIDS-LR）, and gradient inverse distance squared with lapse rate determined by classification and regression tree（GIDS-CART）, are employed to interpolate the bias corrected forecasts at neighboring observation stations to any given location. The results show that all four interpolation methods used do reduce the model error significantly,but Kriging-CLR has better performance than the other methods. For Kriging-CLR, root mean square error（RMSE）and mean absolute error（MAE） were decreased by 26% and 29%, respectively, as compared to the raw forecasts. It is found also, that after applying any of the proposed methods, unlike the raw forecasts, the bias corrected forecasts do not show spatial or temporal dependency.

Edge-Oriented Spatial Interpolation for Error Concealment of Consecutive Blocks

This paper proposes a low-complexity spatial-domain error concealment （EC） algorithm for recovering consecutive blocks error in still images or intra-coded （I） frames of video sequences. The proposed algorithm works with the following steps. Firstly the Sobel operator is performed on the top and bottom adjacent pixels to detect the most probable edge direction of current block area. After that one-dimensional （1-D） matching is used on the available block boundaries. Displacement between edge direction candidate and most probable edge direction is taken into consideration as an important factor to improve stability of 1-D boundary matching. Then the corrupted pixels are recovered by linear weighting interpolation along the estimated edge direction. Finally the interpolated values are merged to get last recovered picture. Simulation results demonstrate that the proposed algorithms obtain good subjective quality and higher PSNR than the methods in literatures for most images.

Spatial interpolation of soil nutrients using algebra hyper-curve neural network

Spatial distribution of and interpolation methods for soil nutrients are the basis of soil nutrient management in precision agriculture.For study of application potential and characteristics of algebra hyper-curve neural network(AHCNN)in delineating spatial variability and interpolation of soil properties,956 soil samples were taken from a 50 hectare field with 20 m interval for alkaline hydrolytic nitrogen measurement.The test data set consisted of 100 random samples extracted from the 956 samples,and the training data set extracted from the remaining samples using 20,40,60,80,100 and 120 m grid intervals.Using the AHCNN model,three training plans were designed,including plan AHC1,using spatial coordinates as the only network input,plan AHC2,adding information of four neighboring points as network input,and plan AHC3,adding information of six neighboring points as network input.The interpolation precision of AHCNN method was compared with that of Kriging method.When the number of training samples was big,interpolation precisions of Kriging and AHCNN were similar.When the number of training samples was small,the precisions of both methods deteriorated.Since AHCNN method has no request on data distribution and it is non-linearization of neutron input variables,it is suitable for delineation of spatial distribution of nonlinear soil properties.In addition,AHCNN has an advantage of adaptive self-adjustment of model parameters,which makes it proper for soil nutrient spatial interpolation.After comparison of mean absolute error d,root mean squared error RMSE,and mean relative error%d,and the spatial distribution maps generated from different methods,it can be concluded that using spatial coordinates as the only network input cannot simulate the characteristics of soil nutrient spatial variability well,and the simulation results can be improved greatly after adding neighboring sample points’information and the distance effect as network input.When the number of samples was small,interpolation precision can be improved after properly increasing the number of neighboring sample points.It was also showed that evaluation of interpolation precision using conventional error statistic indexes was defective,and the spatial distribution map should be used as an important evaluation factor.

A new strategy combined HASM and classical interpolation methods for DEM construction in areas without sufficient terrain data

Spatial interpolation is an important method in the process of DEM construction. However, DEMs constructed by interpolation methods may induce serious distortion of surface morphology in areas lack of terrain data. In order to solve this problem, this paper proposes a strategy combining high-accuracy surface modeling(HASM) and classical interpolation methods to construct DEM. Firstly, a triangulated irregular network(TIN) is built based on the original terrain data, and the area of the triangles in the TIN is used to determine whether to add supplementary altimetric points(SA-Points). Then, classical interpolation methods, such as Inverse Distance Weighted(IDW) method, Kriging, and Spline, are applied to assign elevation values to the SA-Points. Finally, the SA-Points are merged with the original terrain data, and HASM is used to construct DEM. In this research, two test areas which are located in Nanjing suburb in Jiangsu Province and Guiyang suburb in Guizhou Province are selected to verify the feasibility of the new strategy. The study results show that:(1) The combination of HASM and classical interpolation methods can significantly improve the elevation accuracy of DEMs compared with DEM constructed by a single method.(2) The process of adding SA-Points proposed in this study can be repeated in many times. For the test areas in this paper, compared with the results with only one execution, the results with more executions are in much more accordance with the actual terrain.(3) Among all the methods discussed in this paper, the one combined HASM and Kriging produce the best result. Compared with the HASM alone, absolute mean error(MAE) and root mean square error(RMSE) of the best result were reduced from 1.29 m and 1.83 m to 0.68 m and 0.45 m(the first test area), and from 0.32 m and 0.38 m to 0.21 m and 0.28 m( The second test area).

A Gaussian process regression-based sea surface temperature interpolation algorithm

The resolution of ocean reanalysis datasets is generally low because of the limited resolution of their associated numerical models.Low-resolution ocean reanalysis datasets are therefore usually interpolated to provide an initial or boundary field for higher-resolution regional ocean models.However,traditional interpolation methods(nearest neighbor interpolation,bilinear interpolation,and bicubic interpolation)lack physical constraints and can generate significant errors at land-sea boundaries and around islands.In this paper,a machine learning method is used to design an interpolation algorithm based on Gaussian process regression.The method uses a multiscale kernel function to process two-dimensional space meteorological ocean processes and introduces multiscale physical feature information(sea surface wind stress,sea surface heat flux,and ocean current velocity).This greatly improves the spatial resolution of ocean features and the interpolation accuracy.The eff ectiveness of the algorithm was validated through interpolation experiments relating to sea surface temperature(SST).The root mean square error(RMSE)of the interpolation algorithm was 38.9%,43.7%,and 62.4%lower than that of bilinear interpolation,bicubic interpolation,and nearest neighbor interpolation,respectively.The interpolation accuracy was also significantly better in off shore area and around islands.The algorithm has an acceptable runtime cost and good temporal and spatial generalizability.

Study on Prediction Model for Spatial Distribution of the Average Temperature Based on GIS——Take Hainan Island as an Example

[ Objectivel The research aimed to study prediction model for spatial distribution of the average temperature based on GIS. [ Method l Average temperature over the years as research object, based on Ordinary Kriging （OK）, Inverse Distance Weight （ IDW）, SPLINE and Mixed In- terpolation （MLR）, monthly temperature data from 1979 to 2008 at 18 long-term meteorological observation stations in Hainan Island were conduc- ted spatial grid treatment. Via contrasts and analyses on different interpolation methods, the optimum interpolation method for average temperature over the years in Hainan Island was selected. [ Resuitl By error analyses of the four interpolation methods for average temperature in recent 30 years in Hainan Island, it was found that accuracy was MLR 〉 IDW 〉 OK 〉 SPLINE. Spatial interpolation effect of MLR was the best for average temperature in Hainan Island. Spatial distribution of the average temperature in Halnan Island had obvious south-high-north-low latitudinal zonality and vertical zonality of gradually declining as altitude rise. In addition, temperature along coast was slightly higher than that in inland. Lapse rate of the temperature in each month in Hainan Island was 0.38 -0.85℃/100 m, and lapse rate of the annual average temperature was about 0.74 ℃/ 100 m. In different areas, lapse rate of the temperature as altitude was different at different time. [ Condusion] The research provided basis for ob- taining continuous distribution situation of the agricultural meteorological factor and establishing accurate prediction model of the spatial distribution in Hainan Island.

Mapping the Genetic Diversity of Castanea sativa: Exploiting Spatial Analysis for Biogeography and Conservation Studies

The current distribution of forest tree species is a result of natural or human mediated historical and contemporary processes. Knowledge of the spatial distribution of the diversity and divergence of populations is crucial for managing and conserving genetic resources in forest tree species. By combining tools from population genetics, landscape ecology and spatial statistics, landscape genetics thus represents a powerful method for evaluating the geographic patterns of genetic resources at the population level. In this study, we explore the possibility of combining genetic diversity data, spatial statistic tools and GIS technologies to map the genetic divergence and diversity of 31 Castanea sativa populations collected in Spain, Italy, Greece, and Turkey. The IDW technique was used to interpolate the diversity values and divergence indices as expected hetereozygosity (He), allelic richness (Rs), private allelic richness (PRs), and membership values (Q) of each population to different clusters. Genetic diversity maps and a synthetic map of the spatial genetic structure of European chestnut populations were produced. Spatial coincidences between landscape elements and statistically significant genetic discontinuities between populations were investigated. Evidence is provided of the significance of cartographic outputs produced in the study and on their usefulness in managing genetic resources.

Spatial Variability of Soil Carbon to Nitrogen Ratio and Its Driving Factors in Ili River Valley,Xinjiang,Northwest China

Soil carbon to nitrogen(C/N) ratio is one of the most important variables reflecting soil quality and ecological function,and an indicator for assessing carbon and nitrogen nutrition balance of soils.Its variation reflects the carbon and nitrogen cycling of soils.In order to explore the spatial variability of soil C/N ratio and its controlling factors of the Ili River valley in Xinjiang Uygur Autonomous Region,Northwest China,the traditional statistical methods,including correlation analysis,geostatistic alanalys and multiple regression analysis were used.The statistical results showed that the soil C/N ratio varied from 7.00 to 23.11,with a mean value of 10.92,and the coefficient of variation was 31.3%.Correlation analysis showed that longitude,altitude,precipitation,soil water,organic carbon,and total nitrogen were positively correlated with the soil C/N ratio(P < 0.01),whereas negative correlations were found between the soil C/N ratio and latitude,temperature,soil bulk density and soil p H.Ordinary Cokriging interpolation showed that r and ME were 0.73 and 0.57,respectively,indicating that the prediction accuracy was high.The spatial autocorrelation of the soil C/N ratio was 6.4 km,and the nugget effect of the soil C/N ratio was 10% with a patchy distribution,in which the area with high value(12.00–20.41) accounted for 22.6% of the total area.Land uses changed the soil C/N ratio with the order of cultivated land > grass land > forest land > garden.Multiple regression analysis showed that geographical and climatic factors,and soil physical and chemical properties could independently explain 26.8%and 55.4% of the spatial features of soil C/N ratio,while human activities could independently explain 5.4% of the spatial features only.The spatial distribution of soil C/N ratio in the study has important reference value for managing soil carbon and nitrogen,and for improving ecological function to similar regions.

Cokriging of Soil Cation Exchange Capacity Using the First Principal Component Derived from Soil Physico-Chemical Properties

As soil cation exchange capacity （CEC） is a vital indicator of soil quality and pollutant sequestration capacity,a study was conducted to evaluate cokriging of CEC with the principal components derived from soil physico-chemical properties.In Qingdao,China,107 soil samples were collected.Soil CEC was estimated by using 86 soil samples for prediction and 21 soil samples for test.The first two principal components （PC1 and PC2） together explained 60.2% of the total variance of soil physico-chemical properties.The PC1 was highly correlated with CEC （r=0.76,P0.01）,whereas there was no significant correlation between CEC and PC2 （r=0.03）.The PC1 was then used as an auxiliary variable for the prediction of soil CEC.Mean error （ME） and root mean square error （RMSE） of kriging for the test dataset were-1.76 and 3.67 cmolc kg-1,and ME and RMSE of cokriging for the test dataset were-1.47 and 2.95 cmolc kg-1,respectively.The cross-validation R2 for the prediction dataset was 0.24 for kriging and 0.39 for cokriging.The results show that cokriging with PC1 is more reliable than kriging for spatial interpolation.In addition,principal components have the highest potential for cokriging predictions when the principal components have good correlations with the primary variables.