基于地形属性使用回归, 克里格和人工神经网络方法估算饱和导水率空间分布

Several methods,including stepwise regression,ordinary kriging,cokriging,kriging with external drift,kriging with varying local means,regression-kriging,ordinary artificial neural networks,and kriging combined with artificial neural networks,were compared to predict spatial variation of saturated hydraulic conductivity from environmental covariates.All methods except ordinary kriging allow for inclusion of secondary variables.The secondary spatial information used was terrain attributes including elevation,slope gradient,slope aspect,profile curvature and contour curvature.A multiple jackknifing procedure was used as a validation method.Root mean square error (RMSE) and mean absolute error (MAE) were used as the validation indices,with the mean RMSE and mean MAE used to judge the prediction quality.Prediction performance by ordinary kriging was poor,indicating that prediction of saturated hydraulic conductivity can be improved by incorporating ancillary data such as terrain variables.Kriging combined with artificial neural networks performed best.These prediction models made better use of ancillary information in predicting saturated hydraulic conductivity compared with the competing models.The combination of geostatistical predictors with neural computing techniques offers more capability for incorporating ancillary information in predictive soil mapping.There is great potential for further research and development of hybrid methods for digital soil mapping.     

基于不同地表曲面模型预测土壤有机碳含量

Local terrain attributes,which are derived directly from the digital elevation model,have been widely applied in digital soil mapping.This study aimed to evaluate the mapping accuracy of soil organic carbon (SOC) concentration in 2 zones of the Heihe River in China,by combining prediction methods with local terrain attributes derived from different polynomial models.The prediction accuracy was used as a benchmark for those who may be more concerned with how accurately the variability of soil properties is modeled in practice,rather than how morphometric variables and their geomorphologic interpretations are understood and calculated.In this study,2 neighborhood types (square and circular) and 6 representative algorithms (Evans-Young,Horn,Zevenbergen-Thorne,Shary,Shi,and Florinsky algorithms) were applied.In general,35 combinations of first-and second-order derivatives were produced as candidate predictors for soil mapping using two mapping methods (i.e.,kriging with an external drift and geographically weighted regression).The results showed that appropriate local terrain attribute algorithms could better capture the spatial variation of SOC concentration in a region where soil properties are strongly influenced by the topography.Among the different combinations of first-and second-order derivatives used,there was a best combination with a more accurate estimate.For different prediction methods,the relative improvement in the two zones varied between 0.30％ and 9.68％.The SOC maps resulting from the higher-order algorithms (Zevenbergen-Thorne and Florinsky) yielded less interpolation errors.Therefore,it was concluded that the performance of predictive methods,which incorporated auxiliary variables,could be improved by attempting different terrain analysis algorithms.     

Integrating properties of soil map delineations into ordinary kriging

Stratification of a region based on soil map delineations followed by within-stratum interpolation is sometimes used to combine soil maps and spatial interpolation. However, not all delineations are equally suitable to subdivide an area into precisely located mutually exclusive strata. This paper proposes a flow-path to characterize the nature of soil map delineations and a methodology to integrate the properties of map delineations into ordinary kriging.Four types of delineations were distinguished based on three criteria: the nature of transition (discontinuous or gradual), the mapping accuracy, and the structure of the within-unit spatial variation. For each type of delineation the ordinary kriging algorithm was modified to integrate its properties in the interpolation.As a test case, the sand content of the topsoil in the province of West-Flanders (Belgium) was mapped, using independent test data for validation. Inside the mapping units and at delineations identified as gradual transitions, our procedure, termed ordinary kriging integrating properties of map delineations (OKPD) , performed similarly to stratified ordinary kriging (SOK). However, close to the delineations identified as inaccurately mapped discontinuities the mean square prediction error of OKPD was 0.64 times that of SOK. Moreover, near these delineations, the prediction variance was largely underestimated by SOK (relative variance = 5.1), whereas OKPD produced a more realistic value (relative variance = 1.5).     

A world dataset of derived soil properties by FAO–UNESCO soil unit for global modelling

A standardized dataset of derived soil properties for the 106 soil units considered on FAO-UNESCO's 1:5 million scale Soil Map of the World is presented. It was derived from a statistical analysis of the 4353 soil profiles held in the WISE (World Inventory of Soil Emission) database, which was developed at the International Soil Reference and Information Centre (ISRIC) for the geographic quantification of soil factors that control processes of global change. Median values are presented by soil unit for selected soil properties including: pH(H₂O); organic carbon content; cation exchange capacity; sum of exchangeable Ca²⁺, Mg²⁺ and K⁺; exchangeable sodium percentage; bulk density; total porosity; available water capacity; soil drainage class; and gravel content class. Medians for these variables are presented both for the topsoil (0–30 cm) and subsoil (30–100 cm), where applicable. The data set can help to refine ratings for soil quality in global environmental models pending the availability of comprehensive georeferenced databases on soil and terrain resources such as SOTER, the World Soil and Terrain Database. In a Geographical Information System (GIS) it can be linked to the units shown on the digital Soil Map of the World through the legend code.     

Kriging and mapping of copper,lead, and mercury contents in surface soil in Shenzhen area

Two-dimensional ordinary block kriging as an optimal interpolation technique was applied to produce regular grids of predicted estimates of copper, lead, and mercury contents in surface soil in the Shenzhen area for mapping purposes. The kriging analysis was based on theoretical variograms calculated from measured data of eighty three top soil samples. The appropriateness of the whole procedure of kriging analysis was evaluated by both cross-validation and kriging standard deviation mapping. The similarity in area variation between the distribution of soil parent material and the mapped results of copper and lead is demonstrated. The possible influence of agricultural contamination on the spatial distribution pattern of mercurcy is also discussed.     

Small scale digital soil mapping in Southeastern Kenya

Digital soil mapping techniques appear to be an interesting alternative for traditional soil survey techniques. However, most applications deal with (semi-)detailed soil surveys where soil variability is determined by a limited number of soil forming factors. The question that remains is whether digital soil mapping techniques are equally suitable for exploratory or reconnaissance soil surveys in more extensive areas with limited data availability. We applied digital soil mapping in a 13,500 km² study area in Kenya with the main aim to create a reconnaissance soil map to assess clay and soil organic carbon contents in terraced maize fields. Soil spatial variability prediction was based on environmental correlation using the concepts of the soil forming factors equation. During field work, 95 composite soil samples were collected. Auxiliary spatially exhaustive data provided insight on the spatial variation of climate, land cover, topography and parent material. The final digital soil maps were elaborated using regression kriging. The variance explained by the regression kriging models was estimated as 13% and 37% for soil organic carbon and clay respectively. These results were confirmed by cross-validation and provide a significant improvement compared to the existing soil survey.     

Water quality implications of dairy slurry applied to cut pastures in the northeast USA

Abstract. Nitrate nitrogen (NO₃-N) leaching from animal production systems in the northeast USA is a major non-point source of pollution in the Chesapeake Bay. We conducted a study to measure NO₃-N leaching from dairy slurry applied to orchardgrass ( Dactylis glomerata L., cv. Pennlate) using large drainage lysimeters to measure the direct impact of four rates of slurry (urine and faeces) N application (0, 168, 336, 672 kg N ha⁻¹ yr⁻¹) on NO₃-N leaching on three soil types. We then used experimentally-based relationships developed earlier between stocking density and NO₃-N leaching loss and leachate NO₃-N concentration to estimate the added impact of animal grazing. Nitrate N leaching losses from only dairy slurry applied at the 0, 158, 336, and 672 kg N ha⁻¹ yr⁻¹ rates were 5.85, 8.26, 8.83, and 12.1 kg N ha⁻¹ yr⁻¹, respectively with corresponding NO₃-N concentrations of 1.60, 2.30, 2.46, and 3.48 mg l⁻¹. These NO₃-N concentrations met the 10 mg l⁻¹ US EPA drinking water standard. However, when a scenario was constructed to include the effect of NO₃-N leaching caused by animal grazing, the NO₃-N drinking water standard was calculated to be exceeded.     

利用数字高程模型改进高山灰岩坑土壤pH值预测

Among spatial interpolation techniques,geostatistics is generally preferred because it takes into account the spatial correlation between neighbouring observations in order to predict attribute values at unsampled locations.A doline of approximately 15 000 m 2 at 1 900 m above sea level (North Italy) was selected as the study area to estimate a digital elevation model (DEM) using geostatistics,to provide a realistic distribution of the errors and to demonstrate whether using widely available secondary data provided more accurate estimates of soil pH than those obtained by univariate kriging.Elevation was measured at 467 randomly distributed points that were converted into a regular DEM using ordinary kriging.Further,110 pits were located using spatial simulated annealing (SSA) method.The interpolation techniques were multi-linear regression analysis (MLR),ordinary kriging (OK),regression kriging (RK),kriging with external drift (KED) and multi-collocated ordinary cokriging (CKmc).A cross-validation test was used to assess the prediction performances of the different algorithms and then evaluate which methods performed best.RK and KED yielded better results than the more complex CKmc and OK.The choice of the most appropriate interpolation method accounting for redundant auxiliary information was strongly conditioned by site specific situations.     

县域土壤质量数字制图方法比较

土壤质量研究几乎涵盖土壤研究的所有领域,土壤质量制图理论与方法是土壤质量研究的一项重要研究内容。该研究以北京市密云县为研究区,基于土壤质量评价最小数据集和指数和法计算的土壤质量指数,探究了在地学模型支持下区域土壤质量数字制图方法。研究设计了5种区域土壤质量数字制图方法,并比较了不同方法的空间数字制图精度。结果显示,目前广泛使用的基于参评指标空间插值结果的土壤质量数字制图方法精度最低、工序较繁琐,且无法反映研究区景观高度异质的特点;而基于计算后的土壤质量指数（soil quality index,SQI）,借助于地统计学方法的土壤质量数字制图方法相对比较科学合理,其中又以基于计算后的SQI和回归克里格法预测效果最好,均方根误差最小,仅为0.01897,相对于基于参评指标空间插值结果的土壤质量数字制图方法,精度相对提高率最大,达到50%以上。综合考虑空间制图精度、工序的繁简程度,在该研究设计的5种方法中基于计算的SQI和回归克里格法最佳,该法避免了地统计插值在景观高度异质区的应用局限性,预测结果与实际最为相符。     

Soil carbon stocks and projected changes according to land use and management: a case study for Kenya

Abstract. Soil organic (SOC) and inorganic carbon (SIC) stocks of Kenya were determined using four different methods to provide baseline data. The assessments used an updated version of the 1:1 M soil and terrain database for the country. Estimates for national SOC stocks to 1 m depth ranged from 3452 to 3797 Tg C. The findings highlight the need for comprehensive databases of soil and terrain data of good quality that consider more than one representative profile per soil component. The 95% confidence limits for the median, area-weighted SOC content were largest in the humid highlands (15.4–15.7 kg C m⁻²) and smallest in the hot arid zone (4.4–4.5 kg C m⁻²). Conversely, for SIC these values were largest in the arid zone (4.3–4.5 kg C m⁻²) and smallest in high rainfall areas (<0.1 kg C m⁻²). Many croplands in Kenya have been over-exploited, resulting in nutrient depletion and loss of organic matter. The SOC gains considered ecologically and technically feasible upon improved management of croplands were estimated at 5.8–9.7 Tg C over the next 25 years. This corresponds to an estimated annual mitigation potential of 5–9% of Kenya's CO₂-C emissions from fossil fuels, cement manufacturing and land use change for 1990.     

Effects of the nitrification inhibitor dicyandiamide on potassium, magnesium and calcium leaching in grazed grassland

Abstract. Leaching of calcium (Ca), potassium (K) and magnesium (Mg) from urine patches in grazed grassland represents a significant loss of valuable nutrients. We studied the effect on cation loss of treating the soil with a nitrification inhibitor, dicyandiamide (DCD), which was used to reduce nitrate loss by leaching. The soil was a free-draining Lismore stony silt loam (Udic Haplustept loamy skeletal) and the pasture was a mixture of perennial ryegrass ( Lolium perenne ) and white clover ( Trifolium repens ). The treatment of the soil with DCD reduced Ca²⁺ leaching by the equivalent of 50%, from 213 to 107 kg Ca ha⁻¹ yr⁻¹ on a field scale. Potassium leaching was reduced by 65%, from 48 to 17 kg K ha⁻¹ yr⁻¹. Magnesium leaching was reduced by 52%, from 17 to 8 kg Mg ha⁻¹ yr⁻¹. We postulate that the reduced leaching loss of these cations was due to the decreased leaching loss of nitrate under the urine patches, and follows from their reduced requirement as counter ions in the drainage water. The treatment of grazed grassland with DCD thus not only decreases nitrate leaching and nitrous oxide emissions as reported previously, but also decreases the leaching loss of cation nutrients such as Ca²⁺, K⁺ and Mg²⁺.     

A review of farm-scale nutrient budgets for organic farms as a tool for management of soil fertility

Abstract. On organic farms, where the importation of materials to build/maintain soil fertility is restricted, it is important that a balance between inputs and outputs of nutrients is achieved to ensure both short-term productivity and long-term sustainability. This paper considers different approaches to nutrient budgeting on organic farms and evaluates the sources of bias in the measurements and/or estimates of the nutrient inputs and outputs. The paper collates 88 nutrient budgets compiled at the farm scale in nine temperate countries. All the nitrogen (N) budgets showed an N surplus (average 83.2 kg N ha^–1 yr^–1). The efficiency of N use, defined as outputs/inputs, was highest (0.9) and lowest (0.2) in arable and beef systems respectively. The phosphorus (P) and potassium (K) budgets showed both surpluses and deficits (average 3.6 kg P ha^–1 yr^–1, 14.2 kg K ha^–1 yr^–1) with horticultural systems showing large surpluses resulting from purchased manure. The estimation of N fixation and quantities of nutrients in purchased manures may introduce significant errors in nutrient budgets. Overall, the data illustrate the diversity of management systems in place on organic farms, and suggest that used together with soil analysis, nutrient budgets are a useful tool for improving the long-term sustainability of organic systems.     

Analysis of Potassium Uptake by Rice Roots Treated with Aluminum Using a Positron Emitting Nuclide, 38K

We investigate the effect of Al on K⁺ uptake by rice roots. Potassium-38 (³⁸K), a positron emitting nuclide (the half-life: 7.61 min), was used to trace K⁺ behavior. When a rice root was treated with 10μM Al for 24 h, the uptake of ³⁸K in the root was increased in the range of 1 to 2 cm from the root tip compared with that of the control sample. Because the root continued to grow without showing any damage of plasma membrane during the Al treatment, it was suggested that the ³⁸K uptake was not occurred through diffusion into the cells. The uptake of ³⁸K in both treatments, with/without Al, was decreased by VO₄^3- (inhibitor of H⁺-ATPase on plasma membrane) and DNP (H⁺ ionophore) treatment, which suggested that the K⁺ uptake was performed through an active transport, such as H⁺:K⁺ transport or H⁺ gradient promoted by an Al treatment.     

The estimation of mineralization, immobilization and nitrification in nitrogen-15 field experiments using computer simulation

A combination of mathematical analysis and computer simulation, using parameters readily measured in a nitrogen-15 field experiment, is employed to determine rates of mineralization, immobilization and nitrification under a growing crop. The procedure also yields the proportion of crop nitrogen uptake occurring as ammonium and nitrate.
When applied to -results from grass lysimeters receiving 250 or 900 kg N ha^–1 a^–1 as ammonium nitrate, the analysis suggested that at 250 kgN ha^–1 a^–1 64–66% of crop nitrogen uptake was as ammonium; at 900 kg N ha^–1 a^–1 the figure was 43–49%. Nitrification at 250kgNha^–1 was only 13–19kgN ha^–1 over 160d while at 900 kg N ha^–1 between 191 and 232 kg N ha^–1 were nitrified.
The results suggested that the apparent inhibition of nitrification in grassland soils may simply reflect poor substrate competition by nitrifying bacteria. Finally, there was a suggestion that mineralization/immobilization was lower at the high fertilizer rate.     

A coherent geostatistical approach for combining choropleth map and field data in the spatial interpolation of soil properties

Information available for mapping continuous soil attributes often includes point field data and choropleth maps (e.g. soil or geological maps) that model the spatial distribution of soil attributes as the juxtaposition of polygons (areas) with constant values. This paper presents two approaches to incorporate both point and areal data in the spatial interpolation of continuous soil attributes. In the first instance, area-to-point kriging is used to map the variability within soil units while ensuring the coherence of the prediction so that the average of disaggregated estimates is equal to the original areal datum. The resulting estimates are then used as local means in residual kriging. The second approach proceeds in one step and capitalizes on: 1) a general formulation of kriging that allows the combination of both point and areal data through the use of area-to-area, area-to-point, and point-to-point covariances in the kriging system, 2) the availability of GIS to discretize polygons of irregular shape and size, and 3) knowledge of the point-support variogram model that can be inferred directly from point measurements, thereby eliminating the need for deconvolution procedures. The two approaches are illustrated using the geological map and heavy metal concentrations recorded in the topsoil of the Swiss Jura. Sensitivity analysis indicates that the new procedures improve prediction over ordinary kriging and traditional residual kriging based on the assumption that the local mean is constant within each mapping unit.     

Using elevation to aid the geostatistical mapping of rainfall erosivity

This paper addresses the issue of incorporating a digital elevation model into the mapping of annual and monthly erosivity values in the Algarve region (Portugal). Besides linear regression of erosivity against elevation, three geostatistical algorithms are introduced: simple kriging with varying local means (SKlm), kriging with an external drift (KED) and colocated cokriging. Cross validation indicates that the straightforward linear regression, which ignores the information provided by neighboring climatic stations, yields the largest prediction errors in most situations. Smaller prediction errors are produced by SKlm and KED that both use elevation to inform on the local mean of erosivity; kriging with an external drift allows one to assess the relation between the two variables within each kriging search neighborhood instead of globally as for simple kriging with varying local means. The best results are generally obtained using cokriging that incorporates the secondary information directly into the computation of the erosivity estimate. The trade-off cost is the inference and modeling of three direct and cross semivariograms.     

Sequential extraction methods for soil organic nitrogen

(pp. 825–831)
This study was carried out to clarify the effects of soil nitrate before cultivation and amounts of basal-dressed nitrogen on additional N application rate and yields of semi-forced tomato for three years from 1998 to 2000. The amounts and timing of additional N dressing were determined based on diagnosis of petiole sap nitrate. The top-dressing was carried out with a liquid fertilizer when the nitrate concentration of a leaflet's petiole sap of leaf beneath fruit which is 2–4 cm declined below 2000 mg L⁻¹.
For standard yield by the method of fertilizer application based on this condition, no basal-dressed nitrogen was required when soil nitrate before cultivation was 150 mg kg⁻¹ dry soil or higher in the 0–30 cm layer; 38 kg ha⁻¹ of basal-dressed nitrogen, which corresponds to 25% of the standard rate of fertilizer application of Chiba Prefecture, was optimum when soil nitrate before cultivation was 100150 mg kg⁻¹ dry soil; 75 kg ha⁻¹ of basal-dressed nitrogen, which corresponds to 50% of the standard, was optimum when soil nitrate before cultivation was under 100 mg kg⁻¹ dry soil. A standard yield was secured and the rate of nitrogen fertilizer application decreased by 49–76% of the standard by keeping the nitrate concentration of tomato petiole sap between 1000–2000 mg L⁻¹ from early harvest time to topping time under these conditions.     

对比景观土壤性质普通克里格法和回归克里格法的比较研究

The accuracy between ordinary kriging and regression kriging was compared based on the combined consideration of sample size, spatial structure, and auxiliary variables (terrain indices and electromagnetic induction surveys) for a variety of soil properties in two contrasting landscapes (agricultural vs. forested). When spatial structure could not be well captured by point-based observations (e. g., when the ratio of sample spacing over correlation range was > 0.5), or when a strong relationship existed between target soil properties and auxiliary variables (e. g., their R² was > 0.6), regression kriging (RK) was more accurate for interpolating soil properties in both landscapes studied. Otherwise, ordinary kriging (OK) was better. Soil depth and wetness condition did not appear to affect the selection of kriging for soil moisture interpolation, because they did not significantly change the ratio of sample spacing over correlation range and the relationship with the auxiliary variables. Because of a smaller ratio of elevation change over total study area (E/A = 1.2) and multiple parent materials in the agricultural land, OK was generally more accurate in that landscape. In contrast, a larger E/A ratio of 6.8 and a single parent material led to RK being preferable in the steep-sloped forested catchment. The results from this study can be useful for selecting kriging for various soil properties and landscapes.     

Potassium balances for arable soils in southern England 1986–1999

Abstract. The behaviour of potassium (K) in a range of arable soils was examined by plotting the change in exchangeable K of the topsoil (Δ K_ex) at the end of a 3–5 year period against the K balance over the same period (fertilizer K applied minus offtake in crops, estimated from farmers' records of yield and straw removal). Based on the assumption that values for offtake per tonne of crop yield used for UK arable crops MAFF 2000) are valid averages, 10–50% of Δ K_ex was explained by the balance, relationships being stronger on shallow/stony soils. Excess fertilizer tended to increase K_ex and reduced fertilization decreased it, requiring between 1.2 and 5.4 kg K ha⁻¹ for each mg L⁻¹Δ K_ex. However, merely to prevent K_ex falling required an extra 20 kg K ha⁻¹ yr⁻¹ fertilizer on Chalk soils and soils formed in the overlying Tertiary and Quaternary deposits, despite clay contents >18%. Whereas, on older geological materials, medium soils needed no extra K and clays gained 17 kg K ha⁻¹ yr⁻¹. It is unlikely that the apparent losses on some soil types are anomalies due to greater crop K contents. Theory and the literature suggest leaching from the topsoil as a major factor; accumulation in the subsoil was not measured. Recommendations for K fertilization of UK soils might be improved by including loss or gain corrections for certain soil types.     

Adsorption of sulphate and fluoride by variable charge soils

The adsorption of sulphate and fluoride by three variable charge soils was studied. Adsorption increased with increase in the amount of sulphate added at constant _pH, and decreased with increase of _pH.
The ratios of the amount of released OH⁻to that of the adsorbed SO^2-₄ at pH 5–0 were 0.12, 0.14 and 0.20 for the three soils, respectively, much lower than the corresponding OH⁻/F⁻ ratios which ranged from 0.3 to 1.0. For a ferric acrisol the OH⁻ released accounted, on average, for only 15% of the SO^2-₄ adsorbed, leaving more than 60% to be explained by the decrease in positive charge and the increase in negative charge carried by the soil.