Abstract

We assessed the spatial variability of soil physical properties in a clay-loam soil cropped to corn and soybean. The study was conducted at Lincoln University in Jefferson City, Missouri. Soil samples were taken at four depths: 0–10 cm, 10–20, 20–40, and 40–60 cm and were oven dried at 105°C for 72 hours. Bulk density (BDY), volumetric (VWC) and gravimetric (GWC) water contents, volumetric air content (VAC), total pore space (TPS), air-filled (AFPS) and water-filled (WFPS) pore space, the relative gas diffusion coefficient (DIFF), and the pore tortuosity factor (TORT) were calculated. Results showed that, in comparison to depth 1, means for AFPS, Diff, TPS, and VAC decreased in Depth 2. Opposingly, BDY, Tort, VWC, and WFPS increased in depth 2. Semivariogram analysis showed that GWC, VWC, BDY, and TPS in depth 2 fitted to an exponential variogram model. The range of spatial variability () for BDY, TPS, VAC, WFPS, AFPS, DIFF, and TORT was the same (25.77 m) in depths 1 and 4, suggesting that these soil properties can be sampled together at the same distance. The analysis also showed the presence of a strong (≤25%) to weak (>75%) spatial dependence for soil physical properties.

1. Introduction

Characterizing the spatial variability and distribution of soil properties is important in predicting the rates of ecosystem processes with respect to natural and anthropogenic factors [1] and in understanding how ecosystems and their services work [2]. In agriculture, studies of the effects of land management on soil properties have shown that cultivation generally increases the potential for soil degradation due to the breakdown of soil aggregates and the reduction of soil cohesion, water content and nutrient holding capacity [3, 4]. Cultivation, especially when accompanied by tillage, has been reported to have significant effects on topsoil structure and thus the ability of soil to fulfill essential soil functions and services in relation to root growth, gas and water transport and organic matter turnover [57]. Soil properties vary considerably under different crops, tillage type and intensity, fertilizer types and application rates. Consequently, the physical properties of the soil are also affected by many factors that change vertically with depth, laterally across fields and temporally in response to climate and human activity [8]. Since this variability affects plant growth, nutrient dynamics, and other soil processes, knowledge of the spatial variability of soil physical properties is therefore necessary. To study the spatial distribution of soil properties, techniques such as classical statistics and geostatistics have been widely applied [911]. Geostatistics provides the basis for the interpolation and interpretation of the spatial variability of soil properties [9, 1214]. Information on the spatial variability of soil properties leads to better management decisions aimed at correcting problems and at least maintaining productivity and sustainability of the soils and thus increasing the precision of farming practices [1, 15]. A better understanding of the spatial variability of soil properties would enable refining agricultural management practices by identifying sites where remediation and management are needed. This promotes sustainable soil and land use and also provides a valuable base against which subsequent future measurements can be proposed [14]. Despite the importance of this topic in agriculture, the literature is not abundant on the variability of soil physical properties in Central Missouri. Furthermore, existing studies on the spatial variability of soil properties have focused on the top soil (0–20 cm) with less or no studies at deeper soil depths (30–100 cm). The objective of this study was therefore to assess the spatial variability of soil physical properties at various depths (0–10 cm, 10–20, 20–40 and 40–60 cm) in a clay-loam soil cropped to corn and soybean, and determine how knowledge on this variability can affect soil management practices.

2. Materials and Methods

2.1. Experimental Site

The study was conducted at Lincoln University’s Freeman farm in Jefferson City, Missouri. The geographic coordinates of the study site are 38°58′116′′N latitude and 92°10′53′′W longitude. The soil of the experiment site is a Waldron clay-loam (Fine, smectitic, calcareous, mesic Aeric Fluvaquents). The study area is almost flat, with an average slope of 2%. The experimental field was made of 48 plots of 12.19 m width by 21.34 m length each. The 48 plots were arranged in a grid of 4 plots in the width by 12 plots in length as shown in Figure 1. One half of the plots was planted to corn (Zea mays) while the other half was planted to soybean (Glycine max). Soybean and corn plots all received 26.31 kg/ha of nitrogen, 67.25 kg/ha of phosphorus, and 89.67 kg/ha of potassium. Corn plots received 201.75 kg/ha of additional nitrogen in the form of urea.

2.2. Soil Sampling

Soil samples were collected in the middle of each plot after planting and full seeds emergence. Cylindrical cores of 3.15 cm radius and 10 or 20 cm height were used to collect soil samples at four depths: 0–10 cm, 10–20, 20–40 and 40–60 cm, corresponding to depths 1, 2, 3, and 4, respectively. The cylinders of 10 cm height were used for soil samples collection at depths 1 and 2 while the 20 cm height cylinders used for sampling at depths 3 and 4. A total of 576 soil samples were collected as follows: 48 plots × 4 depths × 3 replicates (at the middle of each plot). Collected samples were taken to the laboratory where they were weighed (fresh weight of sample; FWS) then oven dried at 105°C for 72 hrs. The weight was taken after oven drying (dry weight of soil; DWS). Soil physical properties were calculated as follows: Soil bulk density (BDY, g·cm−3) = (DWS/V), where DWS is the dry weight of soil and the volume of cylinder (total volume of soil); Volumetric water content (VWC, cm3·cm−3) = (FWS − DWS)/), with FWS being the fresh weight of soil; gravimetric water content (GWC, g·g−1) = [(FWS − DWS)/DWS] where FWS is the fresh weight of soil; total pore space (TPS, cm3·cm−3) = 1 − (BDY/PDY), where PDY is the soil particle density (taken as 2.65 g cm−3); volumetric air content (VAC, cm3·cm−3) = TPS − VWC; water-filled pore space (WFPS, %) = 100 * (VWC/TPS); air-filled pore space (AFPS, %) = 100 * (VAC/TPS); relative gas diffusion coeffient (Diff., cm2s−1·cm−2·s) = (VAC)2, pore space tortuosity (Tort., m·m−1) = (1/VAC) [16].

2.3. Statistical and Geospatial Analysis

After calculation, data on soil physical properties was first transferred to Statistix 9.0 to compute summaries of simple statistics, then to GS+ (Geostatistics for environmental science) 7.0 for semivariogram analysis. A semivariogram (a measure of the strength of statistical correlation as a function of distance) is defined by the following equation [17]: where is the experimental semivariogram value at a distance interval , is number of sample value pairs within the distance interval , and , and are sample values at two points separated by the distance . Exponential and spherical models were the empirical semivariograms. The stationary models, that is, exponential (2) and spherical model (3) that fitted to experimental semivariograms were defined in the following equations [18]: where is the nugget, is the partial sill, and is the range of spatial dependence to reach the sill . The ratio and the range are the parameters that characterize the spatial structure of a soil property. The relation is the proportion in the dependence zone, and the range defines the distance over which the soil property values are correlated with each other [19]. A low value for the ratio and a high range generally indicate that high precision of the property can be obtained by Parfitt et al. [19]. The classification proposed by Cambardella et al. [14], which considers the degree of spatial dependence as strong when DSD ≤ 25%; moderate when 25 < DSD ≤ 75%; and weak when DSD > 75%, was used in this study to classify the degree of spatial dependence of each soil property.

3. Results and Discussion

3.1. Summaries of Statistics for Soil Physical Properties

Overall, descriptive statistics for soil properties in this study showed moderate to high skewness for some of the properties (Table 1). The highly skewed soil parameters included soil bulk density (BDY), diffusivity (DIFF), and volumetric water content (VWC), whereas total pore space (TPS) was moderately skewed. Air-filled pore space (AFPS) had a low skewness. Highly skewed parameters indicate that these properties have a local distribution; that is, high values were found for these properties at some points, but most values were low [20]. The other soil physical properties were approximately normally distributed on the field. The underlying reason for soil properties being normally or nonnormally distributed may be associated with differences in management practices, land use, vegetation cover, and topographic effects on the variability of soil erosion across the landscape of the field. These factors can be the sources for a large or very small variation of soil properties in some of the samples, which leads to the nonnormal distribution [21]. A wide range of spatial variability was observed for soil physical properties (Table 1). For instance, soil bulk density (BDY) ranged from 1.01 to 1.23 g cm−3 for depth 1, 1.15 to 1.46 g cm−3 for depth 2, 0.96 to 1.19 g cm−3 and 1.04 to 1.18 g cm−3 for depths 3 and 4, respectively (Figure 2). Soil bulk density was also significantly higher in the second depth (1.4 g cm−3) than all the other 3 depths, where it varied between 1.18 g cm−3 and 1.24 g cm−3. The mean value of AFPS was significantly lower in the second depth (26.5 cm3·cm−3) than in all other 3 depths, where it varied from 39.34 to 45.7 cm3·cm−3. Soil pore tortuosity factor (TORT) and water-filled pore space (WFPS) were also significantly higher in the second depth (12.46 cm·cm−1 and 73.46%, resp.). However, the relative gas diffusion coefficient (DIFF), gravimetric water content (GWC), total pore space (TPS), and volumetric air content (VAC) were significantly lower in the second depth (0.02 m2s−1m−2s, 0.21 g·g−1, 0.42 cm3·cm−3, and 0.12 cm3·cm−3, resp.) (Table 1). The variability in soil physical properties is understandable since the soil of this site has a smectite layer (claypan) in the 10–20 cm, which corresponded to our second sampling depth.This layer of smectite is hard and compact, with very low pore space, high mass-volume ratio (bulk density) and high water retention capability (because of their large surface area). As a consequence of the presence of this smectite layer in depth 2, the mean of water-filled pore space (WFPS) was slightly lower in the first depth (54%) than in all four depths. In fact, air predominates the pore space in the first depth and cultivation loosened the soil, thereby allowing the water trapped in the pore space to evaporate. Higher GWC, VWC, and TPS at the lower depths (20–60 cm) mean that crops (especially corn and soybean grown in the field) were able to access water and dissolved nutrients through their roots. In fact, despite the claypan layer (10–20 cm), it has been reported by various researchers that crop roots were able to penetrate into and through this layer of smectitic clay [2224] and that root growth may increase within the claypan layer [23] as a result of plant adaptation to water-limited soil layers. In general, the use of the coefficient of variation (CV) is a common procedure to assess variability in soil properties since it allows comparison among properties with different units of measurement. Overall, the coefficient of variation for all soil physical properties, in the four depths, ranged from 4.83 to 91.61% (Table 1). The pore tortuosity factor (TORT) showed the highest variation while soil bulk density (BDY) showed the least variation. The CV indicated that there was a strong spatial variability of the soil properties investigated. However, to have a better assessment of such spatial variability across the entire field, a geostatistical analysis was used.

3.2. Spatial Variability of Soil Properties

After computing summaries of simple statistics with Statistix 9.0, data on soil physical properties was transferred to GS+ (geostatistics for environmental science) 7.0 for semivariogram analysis. Semivariogram model fit was determined from the coefficient of determination values, which range from 0 (very poor model fit) to 1 (very good model fit). Table 2 shows soil physical properties which mainly responded to exponential and linear variogram models, with the exponential model providing the best fit. In the 10–20 cm depth, exponential model provided the best fit for BDY , with the spherical model providing very poor model fit. Pore tortuosity also responded to an exponential variogram model in the 20–40 cm depth , although spherical model was noticed. Linear and exponential models were observed in the 40–60 cm depth for TPS , with linear model providing a better fit (Table 2). In general, for all depths, model fit was not very strong with the exception of gravimetric water content and bulk density in the second depth. Overall, the exponential model provided the best fit with about 65% of the physical properties fitting this model. In geostatistical theory, the range of the spatial variability of the semivariogram is the distance between correlated measurements (the minimum lateral distance between two points before the change in property is noticed) and can be an effective criterion for the evaluation of sampling design and mapping of soil properties. The value that the semivariogram model attains at the range (the value on the -axis) is called the sill. The partial sill is the sill minus the nugget [25, 26]. Theoretically, at zero separation distance (lag = 0), the semivariogram value is zero. However, at an infinitesimally small separation distance, the semivariogram often exhibits a nugget effect (the apparent discontinuity at the beginning of many semivariogram graphs), which is some value greater than zero. The nugget effect can be attributed to measurement errors or spatial sources of variation at distances smaller than the sampling interval (or both). Measurement error occurs because of the error inherent in measuring devices. To eliminate this error, multiple samples were taken from each sampling point. Natural phenomena can vary spatially over a range of scales. Variation at microscales smaller than the sampling distances will appear as a part of the nugget effect. Table 2 shows that the spatial correlation (range) of soil properties widely varied from 1 m for volumetric water content (VWC) in depth four to 64 m for gravimetric water content (GWC) in depth 2. However, for the first and second depth (which are agriculturally more important), the range of spatial correlation varied from 3 m for volumetric air content (VAC) in depth 2 to 64 m for GWC in depth 2. Beyond these ranges, there is no spatial dependence (autocorrelation). The spatial dependence can indicate the level of similarity or disturbance of the soil condition. According to López-Granados et al. [27] and Ayoubi et al. [17], a large range indicates that the measured soil property value is influenced by natural and anthropogenic factors over greater distances than parameters which have smaller ranges. Thus, a range of about 64 m for GWC in this study indicates that the measured GWC values can be influenced in the soil over greater distances as compared to the soil parameters having smaller range (Table 2). This means that soil variables with smaller range such as VWC and VAC are good indicators of the more disturbed soils (the more disturbed a soil is, the more variable some soil properties become). The more variable properties have a shorter range of correlation. The different ranges of the spatial dependence among the soil properties may be attributed to differences in response to the erosion—deposition factors, land use-cover, parent material, and human interferences in the study area. The nugget, which is an indication of microvariability was significantly higher for water-filled pore space (WFPS) and air-filled pore space (AFPS) when compared to the others. This can be explained by our sampling distance which could not capture well their spatial dependence. The lowest nugget was for GWC (Table 2). This indicates that GWC had low spatial variability within small distances. Knowledge of the range of influence for various soil properties allows one to construct independent accurate datasets for similar areas in future soil sampling design to perform statistical analysis [17]. This aids in determining where to resample, if necessary, and design future field experiments that avoid spatial dependence. Therefore, for future studies aimed at characterizing the spatial dependency of soil properties in the study area and/or a similar area, it is recommended that the soil properties be sampled at distances shorter than the range found in this study. Cambardella et al. [14] established the classification of degree of spatial dependence (DSD) between adjacent observations of soil property > 75% to correspond to weak spatial structure. In this study, the semivariograms indicated strong spatial dependence (DSD ≤ 25%) for soil physical properties such as bulk density, gravimetric water content, volumetric water content, total pore space, and diffusivity. The rest of the soil physical properties (water-filled pore space, Air-filled pore space, and tortuosity) measured exhibited very weak spatial dependence (DSD > 75%) (Table 2). The strong spatial dependence of the soil properties may be controlled by intrinsic variations in soil characteristics such as texture and mineralogy, whereas extrinsic variations such as tillage and other soil and water management practices may also control the variability of the weak spatially dependent parameters [14].

3.3. Spatial Distribution of Soil Properties across the Field

Interpolated maps portraying the distribution of soil physical properties in various depths are shown in Figure 3 for soil gravimetric (GWC) and volumetric (VWC) contents and water-filled pore space (WFPS). Gravimetric water content showed a good spatial distribution across the field with the highest values located around the southwestern portion of the field. Volumetric water content also showed good spatial distribution across the field with high values located in the northern, central, and southwestern portions of the field. Water-filled pore has a distribution similar to that of volumetric water content. The other soil properties, however, showed very poor spatial distribution in the field. This is most probably due to their poor sill , model fit and coefficient of determination . Even though the spatial variability was not very pronounced, there were areas on the field that had slightly higher values of these physical properties than the rest of the field. In general, bulk density, total pore space, volumetric air content, Air-filled pore space, diffusivity, and tortuosity were very high in the field even though they did not exhibit very distinguishable variability. This lack of visible spatial variability is supported because the sampling distance (range) is 26 m for these properties.

3.4. Implications of Spatial Variability of Soil Physical Properties on Soil Management

Results of this study indicated that the spatial variability of soil water content (GWC and VWC) was high. This can be explained, among many other reasons, by soil type (clay-loam) which was able to hold more water. But with intensive tillage, this soil water content could be adversely affected. Studies have shown that tillage practices can alter soil physical properties and consequently the hydrological behavior of agricultural fields, especially when a similar tillage system has been practiced for a long period [15, 2831]. Tillage intensity has also considerable effects on spatial structure and spatial variability of soil properties [15, 30]. Therefore, this study can help determine site-specific soil management and decision making. To do so, the spatial variability of soil properties developed through kriging will be an important tool. Different ranges of spatial dependence were noticed in the field. The different ranges of the spatial dependence among the soil properties may be attributed to differences in response to the erosion—deposition factors, land use-cover, parent material, and human interferences in the study area. The different ranges can also be used in future studies to determine the sampling distance of different soil physical properties on the field. Also, the sill can help determine where the variability or change in soil property stops. This will be useful especially for the irrigation purposes. Generally, with farmers facing the decision of whether or not to till and the intensity of tillage, a spatial variability study can help in this decision making. Maps produced in this study can also be used for irrigation purposes as they can clearly indicate which portion of the field needs irrigation (soil water content). To do this, soil water content information can be collected and analyzed geospatially to produce field maps. The process can be repeated frequently to obtain up-to-date soil water content information. To avoid frequent destructive sampling for water content analysis, equipment that allows insitu measurements such as TDR methods and water mark sensors can be used. Since a different range of spatial dependence among soil properties shows differences in response to human interferences and land use-cover, this will help reduce human activities that increase soil bulk density and cause soil compaction like the use of heavy equipment. It can also serve as a reference for the type of crop to be grown (cover crops for erosion susceptible areas).

4. Conclusion

We assessed the spatial variability of soil physical properties in a clay-loam soil cropped to corn and soybean. Results showed that soil physical properties either decreased or increased sharply in the second depth (due to the presence of a smectite layer) before leveling up or dropping off, but without reaching the first depth value in either case. In addition, depending on soil physical property, maps produced by kriging showed either good or poor spatial distribution. The semivariogram analysis showed the presence of a strong (≤25%) to weak (>75%) spatial dependence of soil properties. Our understanding of the behavior of soil properties in this study provides new insights for soil site-specific management in addressing issues such as “where to place the proper interventions” (tillage, irrigation, and crop type to be grown).

Acknowledgment

This research is part of a regional collaborative project supported by the USDA-NIFA, Award no. 2011-68002-30190, Cropping Systems Coordinated Agricultural Project: Climate Change, Mitigation, and Adaptation in Corn-based Cropping Systems. (Project website: http://sustainablecorn.org/).