Estimating variograms of soil properties by the method-of-moments and maximum likelihood

Variograms of soil properties are usually obtained by estimating the variogram for distinct lag classes by the method‐of‐moments and fitting an appropriate model to the estimates. An alternative is to fit a model by maximum likelihood to data on the assumption that they are a realization of a multivariate Gaussian process. This paper compares the two using both simulation and real data. The method‐of‐moments and maximum likelihood were used to estimate the variograms of data simulated from stationary Gaussian processes. In one example, where the simulated field was sampled at different intensities, maximum likelihood estimation was consistently more efficient than the method‐of‐moments, but this result was not general and the relative performance of the methods depends on the form of the variogram. Where the nugget variance was relatively small and the correlation range of the data was large the method‐of‐moments was at an advantage and likewise in the presence of data from a contaminating distribution. When fields were simulated with positive skew this affected the results of both the method‐of‐moments and maximum likelihood. The two methods were used to estimate variograms from actual metal concentrations in topsoil in the Swiss Jura, and the variograms were used for kriging. Both estimators were susceptible to sampling problems which resulted in over‐ or underestimation of the variance of three of the metals by kriging. For four other metals the results for kriging using the variogram obtained by maximum likelihood were consistently closer to the theoretical expectation than the results for kriging with the variogram obtained by the method‐of‐moments, although the differences between the results using the two approaches were not significantly different from each other or from expectation. Soil scientists should use both procedures in their analysis and compare the results.     

On the Akaike Information Criterion for choosing models for variograms of soil properties

A problem in the application of geostatistics to soil is to find satisfactory models for variograms of soil properties. It is usually solved by fitting plausible models to the sample variogram by weighted least squares approximation. The residual sum of squares can always be diminished, and the fit improved in that sense, by adding parameters to the model. A satisfactory compromise between goodness of fit and parsimony can be achieved by applying the Akaike Information Criterion (AIC). For a given set of data the variable part of the AIC is estimated by where n is the number of experimental points on the variogram, R is the residual sum of squares and p is the number of parameters in the model. The model to choose is the one for which  is least.
The AIC is closely related to Akaike's earlier final prediction error and the Schwarz criterion. It is also equivalent to an F test when adding parameters in nested models.     

Comparing sampling needs for variograms of soil properties computed by the method of moments and residual maximum likelihood

It has been generally accepted that the method of moments (MoM) variogram, which has been widely applied in soil science, requires about 100 sites at an appropriate interval apart to describe the variation adequately. This sample size is often larger than can be afforded for soil surveys of agricultural fields or contaminated sites. Furthermore, it might be a much larger sample size than is needed where the scale of variation is large. A possible alternative in such situations is the residual maximum likelihood (REML) variogram because fewer data appear to be required. The REML method is parametric and is considered reliable where there is trend in the data because it is based on generalized increments that filter trend out and only the covariance parameters are estimated. Previous research has suggested that fewer data are needed to compute a reliable variogram using a maximum likelihood approach such as REML, however, the results can vary according to the nature of the spatial variation. There remain issues to examine: how many fewer data can be used, how should the sampling sites be distributed over the site of interest, and how do different degrees of spatial variation affect the data requirements? The soil of four field sites of different size, physiography, parent material and soil type was sampled intensively, and MoM and REML variograms were calculated for clay content. The data were then sub-sampled to give different sample sizes and distributions of sites and the variograms were computed again. The model parameters for the sets of variograms for each site were used for cross-validation. Predictions based on REML variograms were generally more accurate than those from MoM variograms with fewer than 100 sampling sites. A sample size of around 50 sites at an appropriate distance apart, possibly determined from variograms of ancillary data, appears adequate to compute REML variograms for kriging soil properties for precision agriculture and contaminated sites.     

Vegetal cover to estimate soil erosion hazard in Rhodesia

Evidence is presented to show that percent vegetal cover is the major factor determining the erosion hazard from crops and grassland in Rhodesia. Percent vegetal cover is proposed as a satisfactory quantitative expression for inclusion in soil-loss estimation equations. This proposal is considered to have substantial advantages over the cropping management factor of the Universal Soil Loss Equation and to be particularly suited to developing countries. A system of classifying vegetation is suggested to account for the secondary influence on soil loss of variations in plant height, stem and rooting characteristics.     

两种土壤水分特征曲线间接推求方法对黄土的适应性评价

土壤水分特征曲线(SWCC)是模拟土壤水分运动和溶质运移的一个重要参数,利用十壤的基本物理性质来间接推求SWCC的方法已经成为当今土壤物理学领域的研究热点.为了比较两种SWCC间接推求方法——Arya-Paris物理经验方法(简称AP方法)和Tyler-Wheatcraft分形几何方法(简称TW方法)对黄土的适应性,该文分析了黄土高原296组土壤颗粒分布、容重和水分特征曲线等资料,利用简化的Fredlund(Fred3P)模型模拟得到连续的十壤颗粒分布曲线,然后应用AP和TW方法预测出相应吸力下的土壤含水量.研究结果表明,对于黄十性土壤,AP和TW两种方法的预测结果均达到了一定的精度,相比较而言AP方法的预测效果明显优于TW方法,且受质地影响小.     

采用非负矩阵分解方法的秸秆-土壤光谱解混及盐分估测

黄河三角洲地区发展节约型农业,大力推行秸秆还田技术以改善土壤生态环境,当利用遥感技术对部分秸秆覆盖地区进行土壤盐渍化监测时,混合像元制约土壤盐分遥感估测的精度。该研究设置10组不同秸秆覆盖度的盐渍化地表场景,多次光谱测量并取平均值作为其反射光谱。探究非负矩阵分解（Non-negative Matrix Factorization,NMF）方法分离农业秸秆和土壤光谱的有效性,并基于分离出的土壤光谱构建偏最小二乘法土壤盐分估测模型。结果表明：1）秸秆覆盖会使土壤光谱反射率增高,相比于纯土壤光谱在1 730 nm、2 090 nm附近出现吸收谷;2）通过NMF进行混合光谱解混之后,可以有效地将土壤光谱从秸秆覆盖光谱中分离出来,且土壤光谱保留了土壤盐分信息;3）对于不同秸杆覆盖水平下的土壤盐分估测模型而言,与原始混合光谱相比,利用土壤光谱NMF分离数据构建的模型精度普遍提高,建模集决定系数平均提高0.07,均方根误差降低1.21。验证集决定系数、相对分析误差分别平均提高0.07、0.25,均方根误差降低1.22。该研究可为黄河三角洲地区部分秸秆覆盖盐渍土地的近地面遥感估测精度的提高提供方法。     

New method to estimate root biomass in soil through root-derived carbon

Quantification of root biomass through the conventional root excavation and washing method is inefficient. A pot experiment was conducted to estimate root-derived carbon (C) in soil. Spring wheat (Triticum aestivum L. cv. ‘Quantum’) was grown in plastic containers (6 L) filled with sterilized sandy soil in a greenhouse. Plants were enriched with ¹³CO₂ in a glass chamber twice at growth stages GS-37 and GS-59 for 70 min at each time. In one treatment, roots were separated from soil at crop maturity, washed and dried for the determination of biomass. Isotope ratios were then separately analyzed for roots and soil. In a second treatment, roots were thoroughly mixed with the whole soil and representative samples were analyzed for ¹³C abundance at crop maturity. Control plants were untreated with ¹³C, in which roots were separated from soil. The root biomass was calculated based on the root-derived C, which was measured through ¹³C abundance in the soil and root mixed samples. A substantial amount of root-derived C (24%) was unaccounted while separating the roots from soil. Similarly, about 36% of the root biomass was underestimated if conventional root excavation and washing method is used. It has been shown that root biomass can be estimated more accurately from the root-derived C using ¹³C tracer method than the estimates made by the conventional excavation and washing method. We propose this as an alternative method for the estimation of root-derived C in soil, based on which root biomass can be estimated.     

A general approach to estimate soil water content from thermal inertia

Remote sensing is a promising technique for obtaining information of the earth's surface. Remotely sensed thermal inertia has been suggested for mapping soil water content. However, a general relationship between soil thermal inertia and water content is required to estimate soil water content from remotely sensed thermal inertia. In this study, we propose a new model that relates soil thermal inertia as a function of water content. The model requires readily available soil characteristics such as soil texture and bulk density. Heat pulse measurements of thermal inertia as a function of water content on nine soils of different textures were made to generate a universal Kerstan function. Model validation was performed independently in both laboratory and field, and the retrieved soil water contents from the new model were compared with previous models. Laboratory evaluation on an Iowa silt loam showed that the RMSE of the new model was 0.029 m³ m⁻³, significantly less than [Murray, T., Verhoef, A., 2007. Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements. I. A universal approach to calculate thermal inertia. Agric. For. Meteorol. 147, 80–87] model (0.109 m³ m⁻³) and [Ma, A.N., Xue, Y., 1990. A study of remote sensing information model of soil moisture. In: Proceedings of the 11th Asian Conference on Remote Sensing. I. November 15-21. International Academic Publishers, Beijing, pp. P-11-1P-11-5.] model (0.105 m³ m⁻³). Similar results were obtained in a field test on a Chinese silt loam: the RMSE of the new model, [Murray, T., Verhoef, A., 2007. Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements. I. A universal approach to calculate thermal inertia. Agric. For. Meteorol. 147, 80–87] model, and [Ma, A.N., Xue, Y., 1990. A study of remote sensing information model of soil moisture. In: Proceedings of the 11th Asian Conference on Remote Sensing. I. November 15-21. International Academic Publishers, Beijing, pp. P-11-1P-11-5.] model were 0.018, 0.071, and 0.159 m³ m⁻³, respectively. Additionally the model was validated using literature data in which soil thermal properties were estimated from in situ temperature measurements. The mean errors of estimated water content were generally less than 0.02 m³ m⁻³. We concluded that the new model was able to provide accurate water content predictions from soil thermal inertia.     

A method to estimate the water vapour enhancement factor in soil

Enhanced water vapour diffusion under temperature gradients has been proposed as a mechanism to explain the discrepancies between measured and predicted water fluxes in soils. Because of the difficulties in measuring soil vapour diffusion directly, modelling approaches have been used to estimate the vapour enhancement factor (η) by matching theory to measurements. In the method proposed by Hiraiwa & Kasubuchi (2000) , soil thermal conductivity associated with conduction heat transfer (λ_c) is assumed to be equal to the apparent soil thermal conductivity (λ) measured at a low temperature, and η is significantly under‐estimated. In the present study, an improved approach for estimating η is used, in which λ_c is taken as the apparent soil thermal conductivity associated with infinite atmospheric pressure. The λ at infinite atmospheric pressure is estimated by extrapolating λ measurements made at finite air pressures. By subtracting λ_c from measured λ values at a given atmospheric pressure, the contribution of thermal vapour diffusion to heat transfer (λ_v) is obtained and then used to estimate η. In the case of a lysimeter sand, λ_v accounts for 4–25, 8–29 and 13–35% of λ at 3.5, 22.5 and 32.5°C, respectively, at soil water contents greater than 0.02 m³ m⁻³. Thus, the latent heat transfer through vapour diffusion is important even at temperatures as low as 3.5°C. The agreement between predictions from the new method and selected literature values suggests that the improved approach is able to provide accurate estimate of η. The results from this study show that the magnitude of latent heat transfer resulting from thermal vapour diffusion is strongly soil texture‐dependent. Thus, it is important to estimate η on specific soils rather than assuming η from literature values.     

土壤质量指标和评价方法

土壤质量是反映土壤保持生物生产力、环境质量以及动植物健康能力的土壤的内在属性。土壤质量评价是根据土壤内部属性对土壤综合状况进行表达。评价必须确定合适的时间和空间尺度,综合量化分析指标和不可量化的描述性指标。评价方法包括定性评价和定量评价,系统化的定性评价仍有其现实意义。定量评价要选择指标的最小数据集(MDS),采用最小数据集(MDS) 土壤变换函数(PTF)的方法将土壤质量指标和土壤功能相联系,对土壤的各个功能因子进行评价,在此基础上利用各种数学评价方法得到综合的量化土壤质量评价结果。     

The use of cation activity ratios to estimate the intensity of soil acidity

Soil solution extracts can be relatively difficult to obtain and to handle. The aim of this study was to investigate whether a simple 0.01 M CaCI, extraction could be used to estimate the intensity factor of soil acidity by employing the Schofield Ratio Law. Activity ratios of H, Ca and A1 were compared in soil solutions and in 0.01 M CaCI, extracts on four soils adjusted to a range of pH values with acid and lime. The ratios pH-1/2pCa (0.01 M CaCl₂) and pH-1/2p(Ca + Mg) (soil solution) were not statistically significantly different. The ratio pH-1/3pAl was less strongly correlated, but again was not significantly different between the soil solution and extractant. The Schofield Ratio Law did not hold over the entire range of data for pCa-+2/3pAl, but the values in the two solutions correlated closely. Cation activity ratios in 0.01 M CaCl₂ can be used to estimate the intensity factors of acidity in a soil, even though the extractant has a higher ionic strength than many soil solutions in variable charge soils. Therefore pH and A1 as measured in 0.01 Cacl₂ extracts may be useful tests for routine diagnostic purposes and soil surveys.     

Can the extent of degradation of soil fungal mycelium during soil incubation be used to estimate ectomycorrhizal biomass in soil?

The extent of degradation of the fungal biomass in forest soil during laboratory incubation was investigated as a measure of ectomycorrhizal (EM) biomass. The method simulates the disappearance of fungal mycelium after root trenching, where the EM fungi, deprived of its energy source (the tree), will start to die off. Incubating a forest humus soil at 25 °C resulted in a decrease in the relative proportion (mol%) of the phospholipid fatty acid 18:2ω6,9 (a fungal marker molecule) within 3-6 months, indicating that fungal biomass was disappearing. Incubation at 5 °C resulted in essentially no change in the amount of 18:2ω6,9. The measurement of ergosterol, another fungal marker molecule, gave similar results. Incubation of different forest soils (pine, spruce and spruce/oak), and assuming that the disappearance of fungal biomass during this period of time was entirely due to EM fungi, resulted in an estimation of EM biomass of between 47 and 84% of the total fungal biomass in these soils. The humus layer had more EM biomass than deeper mineral layers.     

Use of organic carbon and loss-on-ignition to estimate soil organic matter in different soil types and horizons

Summary Loss-on-ignition (LOI) and the organic C content have been used to estimate soil organic matter. Organic matter is often estimated from organic C by applying a factor of 1.724. Several authors have examined the relationship between LOI, used as an estimate of organic matter, and C by simple linear regressions. In the present study, this approach was examined in relation to two sets of data. LOI overestimates organic matter in soils with significant proportions of clay minerals because of bound water, and correcting for bound water gives some LOI: C ratios of less than 1. It is concluded that differences in the nature of the organic matter in different soils and horizons make the simple regression approach unsuitable. More attention needs to be paid to studies of the nature of the organic matter.     

Use of force-deformation curves to estimate Young's modulus and its applications to soil aggregate breakdown

Soil aggregates from a clay and a sandy loam were subjected to uniaxial compression and their force-deformation relationships determined. From these relationships the elastic deformation was identified, and estimates were made of Young's modulus, elastic strain energy and work of fracture. Also estimated were the increases in surface area produced on breakdown and hence surface energies. The results were interpreted in terms of Griffith's crack theory, and apparent crack lengths were calculated. The calculated values were high, but the variability of the data prevented a full assessment of this approach to the study of brittle failure of aggregates.     

Contribution of bricks to urban soil properties

Purpose

Bricks are regularly found in urban soils where they can strongly impact soil properties. The purpose of this study is to investigate abundance, especially in the fine earth fraction, and properties of bricks in urban soils, focusing on rooting, plant nutrition and contamination.

Materials and methods

Three different urban soils from the city of Berlin have been studied for their brick contents in the coarse and fine earth fractions by hand sorting. Light (LM) and scanning electron microscopy (SEM) was employed to investigate bricks for proofs of rooting. Third, CEC, pH, EC, C_org, nutrient and contaminant storage and availability have been investigated for bricks and the fine earth fractions of the corresponding soil horizons.

Results and discussion

The fine earth fractions of the investigated soils contain 3 to 5 % of bricks, while the coarse fractions contain up to 50 %. The LM and SEM micrographs made the proof that roots enter brick pores or attach to brick surfaces. Therefore, they can use the water and nutrients stored in bricks and bypass pore system discontinuities between bricks and surrounding soil. The CEC of bricks is grain size dependent and reaches a maximum of 6 cmol_c kg^?1 for particles smaller than 0.063 mm. This dependency is the result of the restricted diffusion into the brick pore system due to the short shaking time in the CEC analysis protocol and of the rising surface with decreasing particle size. From the nutrient storage and availability, we conclude that bricks can better supply plants with K, Mg, Ca and S than the investigated sandy bulk soil.

Conclusions

The nutrient availability from bricks is low compared to control soils, except for Ca and S. Because of the water and nutrient storage, low contamination status and the possible rooting of bricks, they can be used for amelioration of poor sandy soils and for constructed Technosols, preferably employed in small grain sizes.     

Does soil ergosterol concentration provide a reliable estimate of soil fungal biomass?

Our aim was to determine if soil ergosterol concentration provides a quantitative estimate of the soil fungal biomass concentration, as is usually assumed. This was done by comparing soil ergosterol measurements with soil fungal biomass (fungal biomass C) concentrations estimated by microscopic measurements and by the selective inhibition technique linked to substrate-induced respiration (SIR). The measurements were compared in a silty-clay loam soil given a range of previous treatments designed to increase or decrease the soil fungal biomass and so also to change the soil ergosterol concentration. The treatments used were ryegrass amendment, to increase the total and fungal biomass, and CHCl₃-fumigation and the addition of the biocides, captan, bronopol and dinoseb, to decrease both ergosterol and fungal biomass C concentrations. The mineralization of ergosterol following addition to sand innoculated with soil extract, and to a sandy loam soil, was also determined. The added ergosterol was little, if at all, degraded following addition to either sand or the unfumigated or fumigated soil during a 10 d aerobic incubation. Similarly, pesticide addition did not significantly change soil ergosterol concentrations yet the soil fungal biomass C concentration decreased significantly. Thus, the ratio: (soil ergosterol concentration/soil fungal biomass C concentration) was much higher in the pesticide-treated soils than the control soil. Following ryegrass amendment, soil ergosterol concentration increased from about 6-12 μg⁻¹ soil within 5 d and then decreased gradually to about 7 μg g⁻¹ soil by 20 d incubation. Changes in fungal biomass C (measured by direct microscopy) closely mirrored changes in soil ergosterol over this period. However, when the amended soil was fumigated and then incubated for a further 5 d, the initial ergosterol concentration declined from 7 to 5 μg g⁻¹ soil by 20 d incubation (a decline of about 0.4). The comparable decline in fungal biomass C was about eight-fold. Thus the ratio of ergosterol to fungal biomass C increased from 0.005 to about 0.01. There was a significant correlation (r>0.84, P<0.001) between soil ergosterol concentration and fungal biomass measured by either SIR or microscopy. However, three data points played a vital role in the correlation. When these points were excluded the relationship was very poor (r<0.4). Our results therefore suggest that substantial amounts of ergosterol may exist, other than in living cells, for considerable periods, with little, if any mineralization. Thus, these results indicate that ergosterol and fungal biomass C concentrations are not always closely correlated, due to the slow metabolism of ergosterol in recently dead fugal biomass and/or the existence of exocellular ergosterol in soil.     

Validating the use of Cs measurements to estimate rates of soil redistribution by wind

Wind erosion has degraded over one-half billion hectares of land worldwide. ¹³⁷Cesium (¹³⁷Cs) has been used as a tracer to study long-term rates of soil redistribution by water and, to a lesser extent, by wind. Early studies assumed that the decline in ¹³⁷Cs activity for a potentially eroded soil relative to that for an uneroded soil was linearly proportional to soil loss. More recently, models have emerged that consider the effects of soil cultivation and the particle surface area-dependent partitioning of ¹³⁷Cs on soils. We investigated the partitioning of ¹³⁷Cs in wind-eroded sediments and with soil surface samples sieved into contiguous ranges of particle sizes. We also compared the ¹³⁷Cs activities and stratification of several adjacent soils with known wind erosion and deposition histories. Finally, we tested ¹³⁷Cs-based soil loss models with measured data from sites with documented histories. ¹³⁷Cs activities and mean particle diameters of aeolian samples agreed well with the ¹³⁷Cs activities and respective mean diameters of the sieved surface soil samples. Good agreement between model estimations and measured data indicated that ¹³⁷Cs models developed to estimate soil redistribution by water were also applicable to soil redistribution by wind provided that the models contained an appropriate particle size correction parameter.     

Sample volume effect on determination of soil salinity

Abstract

Five pedons from a 200 hectare area were sampled to study the effect of the soil sample volume on determination of salinity. Four sample sizes, 5, 50, 500 and 5000 g, were used, with five to ten replicates each, for the top 25 cm of the profiles. These samples were prepared and analyzed for soil salinity, soluble cations and anions in 1: 2 (soil: water) extracts. Data revealed that the variations due to the different sample sizes were highly significant. Highest values of EC, cations and anions were obtained when using 5‐g samples. The confidence intervals tended to be lowest for the 5000‐g samples in most of the EC, cation and anion determinations. The 500‐g sample size was awkward and apparently the 5‐g sample size was too small.     

Role of soil properties in sewage sludge toxicity to soil collembolans

Soil properties are one of the most important factors explaining the different toxicity results found in different soils. Although there is knowledge about the role of soil properties on the toxicity of individual chemicals, not much is known about its relevance for sewage sludge amendments. In particular little is known about the effect of soil properties on the toxicity modulation of these complex wastes. In addition, in most studies on sewage sludges the identity of the main substances linked to the toxicity and the influence of soil properties on their bioavailability remains unknown.In this study, the toxicity of a sewage sludge to the soil collembolan Folsomia candida was assessed in nine natural soils from agricultural, grassland and woodland sites, together with the OECD soil. Correlations between the relative toxicity of sludge for collembolans in the different soils and their physical and chemical soil properties were assessed in order to identify the main compounds responsible for the effects observed. Furthermore, the relationships between the toxic effects to collembolans and water-soluble ions released by sludge, pH and electric conductivity were also assessed, together with the modulating effects of soil properties.Sludge toxicity was directly linked to the water extractable ammonium, which explained most of the mortality of the collembolans, and part of the inhibition of reproduction. For the last endpoint, nitrite also contributed significantly to the inhibition observed. The varied levels in water extractable ammonium in the different soils at equal dosages seem to be, in turn, modulated by some soil properties. Higher organic carbon contents were associated with lower toxicity of sludge, both for survival and reproduction, probably related to its higher ammonium sorption capacity. In addition, for reproduction, increasing the C/N ratio and pH appeared to increase the toxicity, probably due to both the greater difficultly in nitrification and the known unsuitability of alkaline soils for this species.