Next Article in Journal
Spatial and Temporal Distribution of Ecballium elaterium in Almond Orchards
Previous Article in Journal
Avocado Production and Local Trade in the Southern Highlands of Tanzania: A Case of an Emerging Trade Commodity from Horticulture
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soil Physical Properties Spatial Variability under Long-Term No-Tillage Corn

1
College of Agriculture and Human Sciences, Prairie View A&M University, Prairie View, TX 77446, USA
2
Civil and Environmental Engineering, University of California, Merced, CA 95343, USA
3
Division of Agriculture and Natural Resources, University of California, Davis, CA 95618, USA
4
Faculty of Sustainable Design Engineering, University of Prince Edward Island, Charlottetown, PE C1A 4P3, Canada
5
College of Tropical Agriculture and Human Sciences, University of Hawaii at Mānoa, Honolulu, HI 96822, USA
6
Department of Plant and Soil Science, Texas Tech University, Lubbock, TX 79409, USA
*
Author to whom correspondence should be addressed.
Agronomy 2019, 9(11), 750; https://doi.org/10.3390/agronomy9110750
Submission received: 16 September 2019 / Revised: 6 November 2019 / Accepted: 11 November 2019 / Published: 13 November 2019
(This article belongs to the Section Soil and Plant Nutrition)

Abstract

:
Spatial variability of soil physical and hydrological properties within or among agricultural fields could be intrinsically induced due to geologic and pedologic soil forming factors, but some of the variability may be induced by anthropogenic activities such as tillage practices. No-tillage has been gaining ground as a successful conservation practice, and quantifying spatial variability of soil physical properties induced by no-tillage practices is a prerequisite for making appropriate site-specific agricultural management decisions and/or reformulating some management practices. In particular, there remains very limited information on the spatial variability of soil physical properties under long-term no-tillage corn and tropical soil conditions. Therefore, the main objective of this study was to quantify the spatial variability of some selected soil physical properties (soil surface temperature (ST), volumetric water content (θv), soil resistance (TIP), total porosity (θt), bulk density (ρb), organic carbon, and saturated hydraulic conductivity (Ksat)) using classical and geostatistical methods. The study site was a 2 ha field cropped no-tillage sweet corn for nearly 10 years on Oahu, Hawaii. The field was divided into 10 × 10 and 20 × 20 m grids. Soil samples were collected at each grid for measuring ρb, θt, and soil organic carbon (SOC) in the laboratory following standard methods. Saturated hydraulic conductivity, TIP at 10 and 20 cm depths, soil surface temperature, and θv were also measured. Porosity and ρb have low and low to moderate variability, respectively based on the relative ranking of the magnitude of variability drawn from the coefficient of variation. Variability of the SOC, TIP, and Ksat ranges from moderate to high. Based on the best-fitted semivariogram model for finer grid data, 9.8 m and 142.2 m are the cut off beyond which the measured parameter does not show any spatial correlation for SOC, and TIP at 10 cm depth, respectively. Bulk density shows the highest spatial dependence (range = 226.8 m) among all measured properties. Spatial distribution of the soil properties based on kriging shows a high level of variability even though the sampled field is relatively small.

1. Introduction

Knowledge of soil physical, hydrological, chemical, and biological properties is important for many purposes including precision farming, environmental management, and crop growth modeling. Soils are characterized by a high degree of spatial variability because of the combined effects of physical, chemical, and biological processes that operate with different intensities and at different scales [1]. Changes in soil properties and processes, in turn, impact plant growth and the environment [2] and require accurate knowledge of these properties for efficient utilization of agricultural inputs. Reduced tillage practices in agronomic crops such as corn, soybeans, cotton, sorghum, and cereal grains were introduced in the mid-twentieth century to conserve soil and water [3]. In the United States, acreage under conservation tillage has increased from 2% in 1968 to 41% in 2004 [4,5]. A recent study by Claassen et al. [6] showed adoption of different conservation tillage practices varies across crops; it revealed that conservation tillage was used on a majority of wheat (67%, 2017), corn (65%, 2016), and soybeans (70%, 2012) whereas conservation tillage was used on just 40% of cotton acres (2015) [6]. Crops under conservation tillage are grown with minimum soil disturbance where crop residues are left in the field after harvesting of previous crop resulting in at least 30 percent of the soil is covered by crop residues [3]. Conservation tillage practices include no-till/strip-till, ridge-till, and mulch-till [7]. No-tillage, also known as no-till or zero tillage protects the soil surface and reduces soil erosion. However, no-tillage practice can lead to over-compaction of the soil [8]. Several studies assessed the long-term effect of no-tillage on soil compaction; they reported that long-term practice of no-tillage increases soil organic carbon (SOC) and pore connectivity in the surface layer [8,9]. Some studies [10,11] reported the formation of a relatively over compacted layer under no-tillage systems whereas other studies [12,13] did not experience any compacted layer formation under no-tillage conditions. There is evidence that no-tillage practices cause significant improvements in water retention, aggregate stability, temperature, and other properties compared to conventional tillage [14,15]. Dao [16] found that although tillage initially decreases ρb; however, soil under no-till has lower ρb than the soil under traditional tillage practices. In contrast, Hammel [17] found that soils under no-tillage had a higher ρb as compared to conventional tillage in the top 30 cm soil surface. It seems there is site-specific effect of the response of some soil hydraulic properties to tillage treatments [18,19].
Spatial variability of soil properties may be enhanced by tillage and other on-farm management practices [20]. Soil characteristics generally show spatial dependence [21], which means that samples close to each other are more similar than samples farther apart. Classic statistical methods have generally been used to study soil variability in the field; however, geostatistical methods have been used to understand the spatial correlation structure [22] of soil properties. Geostatistics take into account the spatial autocorrelation of a given variable in predicting its value in another unsampled location. Previous studies of spatial variability of soil physical properties within or among agricultural fields have yielded contrasting results across geographical locations, crops, soils, climates, and agricultural practices particularly tillage practices [18]. Although no-tillage has become a successful conservation practice, quantifying spatial variability of soil physical properties induced by no-tillage practices is a prerequisite for making appropriate site-specific agricultural management decisions or reformulating some management practices. In particular, there remains very limited information on the spatial variability of soil physical properties under long-term no-tillage corn under tropical soil conditions. Therefore, the objective of this work was to quantify some of the selected soil physical properties (ST, θv, θt, TIP, ρb, SOC, and Ksat) and describe their spatial dependence using classical and geostatistical methods.

2. Materials and Methods

2.1. Study Area

This study was conducted at the Waimanalo Research Station (21°20′15′′ N; 157°43′30′′ W) of the College of Tropical Agriculture and Human Resources, the University of Hawaii at Mānoa (Figure 1). The 2 ha experimental field (130 m × 160 m) was continuously cropped with sweet corn crop (Zea mays var. rugosa) under no-tillage for 10 years. The climate of the study area is tropical with a mean annual air temperature of 23 °C, and an average annual precipitation of 938 mm most of which occurs between November and April. The soil at this location has 2% to 6% slope and is classified as Waialua gravelly clay variant (isohyperthermic Pachic Haplustolls). The percentage of sand, silt, and clay at the top 18 cm of soil profile are 8.7, 45.2, and 46.1, respectively [23]. Corn residues were not removed and management traffic patterns were parallel to the crop rows to minimize compaction in planting zone. A non-selective herbicide, Glyphosate (N-phosphonomethyl glycine), was applied when needed to control weeds prior to planting.

2.2. Data Collection

2.2.1. Soil Sample Collection and Laboratory Analysis

The 10 m × 10 m sampling grid was established on the study field (Figure 1). Undisturbed soil core samples (radius, 2.5 cm; height, 6 cm) were collected from the top 10 cm depth of each grid point between 6 December, 2008 and 9 December, 2008 using brass cylinders. These samples were used to determine ρb and θt, following the methods of Grossman and Reinsch [24] and Dane and Hopmans [25], respectively.
In the laboratory, soil core samples were placed in a tray filled with water for saturation after securing the soil from the bottom using nylon mesh. Samples were left in water for approximately 24 h to ensure their complete saturation. The saturated samples were weighed and then oven dried at 105 °C for at least 24 h.
After ρb and θt, analysis, the soil samples of the top 10 cm soil layers were used to determine the SOC using the loss on ignition (LOI) method [26] using a muffle furnace (Model 550 Isotemp Series, Fisher Scientific, Pittsburgh, PA, USA) and following the standard procedure of Schulte and Hopkins [27]. Five grams of air-dry soil samples with a diameter <2 mm were placed in 15 mL ceramic cups, oven-dried at 105 °C for 24 h, cooled in a desiccator for a couple of hours, and then weighed (M105). The samples were then combusted at 550 °C for 5 h [28]. The combusted samples were again transferred to an oven at a temperature of 105 °C for a couple of hours. The samples were then taken out, cooled down in desiccators, and weighed again (M550). The percent LOI was calculated as follows:
LOI ( % ) = M 105 M 550 M 105 × 100
A regression model was established between LOI and SOC using the Hawaii soil SSURGO data from the United States Department of Agriculture Natural Resources Conservation Services (USDA-NRCS). Established regression equation (Figure 2) was used to calculate SOC based on LOI values.

2.2.2. Field Measurement

The soil surface temperature, θv, and soil resistance (TIP) were measured at each intersection point of 10 m × 10 m grid during 17 December, 2008 and 18 December, 2008. Similar measurements were also conducted on a 20 m × 20 m grid on 6 January, 2009 to quantify the temporal variability of ST, θv, and TIP. The soil temperature was measured using the VWR® two-channel thermometer (VWR International, LLC., Radnor, PA, USA). In situ θv was measured using the FIELDSCOUT ™ TDR 300 soil water content sensor, which has an accuracy ± 3.0% volumetric water content with electrical conductivity <2 dS m−1 of (Spectrum Technologies, Inc., Plainfield, IL, USA) in the top 10 cm depth. The soil penetration resistance was measured at 10 cm and 20 cm depths using DICKY-john® soil compaction tester (DICKY-john® Corporation, Auburn, IL, USA). All the above three measurements were conducted in three replicates within few cm and their mean values were used in the data analysis.
Infiltration measurements were conducted on 20 m × 20 m grid over the entire field and a 10 m × 10 m grid on the transect (160,0) using a tension infiltrometer (Soil Measurement Systems, Tucson, AZ, USA) between 6 December, 2008 and 10 December, 2008. The infiltration tests were carried out at two tensions, −15 and −20 cm H2O, with 5 cm tension differences. The falling head of water in the graduated water reservoir was monitored at 2 to 5 min intervals depending on the infiltration rates. The initial readings were taken frequently. At three consecutive intervals of around 5 min, if the rates did not change, steady state conditions were assumed, and the resulting data were used to calculate Ksat using the following equation [29]:
K s a t = Q π   r 2 exp ( α ψ ) [ 1 + 4 π   r α ]
where Q (LT−1) is steady state water flux, r (m) is the radius of perforated disk, α (T−1) is inverse macro-porosity capillary length, ψ (L) is tension at which Q was measured. The inverse macro-porosity capillary length α was estimated using the relationship given by Hussen and Warrick [30] as follows:
α = ln [ Q ( ψ 2 ) / Q ( ψ 1 ) ] [ ψ 2 ψ 1 ]
where ψ1 and ψ2 are the lower and upper bounds of the tension.

2.2.3. Statistical and Geostatistical Analysis

Mean, maximum, minimum, and variance were calculated for each of the properties measured during this study. Frequency distributions were also produced for each soil property and Shapiro-Wilk test [31] was used for normality test using Statistix software package (Analytical Software 2003, Tallahassee, FL, USA).
Spatial variability of ST, θv, TIP, ρb, θt, SOC, and Ksat was described using semivariograms assuming stationarity among similar lag that explains the average dissimilarity between the samples separated by a distance h. Semivariogram values at each lag separation γ(h) were calculated as follows:
γ ( h ) = 1 2 N ( h ) i = 1 N ( h ) [ Z ( x i ) Z ( x i + h ) ] 2
where N(h) is the number of data pairs separated by distance h, Z(xi) is the measured value of the variable at location xi and Z(xi + h) is the measured value of the variable at a distance of h from the location xi.
Experimental semivariogram, cross-variogram, correlogram, and Moran’s I were computed for each property using GS+ [32] software package- version 7 (Gamma Design Software, LLC, Plainwell, Michigan, USA) and plotted with separation distance (h). During the variogram analysis, active lag distance was taken as half of the minimum extent of sampling area to avoid the bordering effect. In the estimation of the semi-variance, as the separation distance gets bigger, the number of samples separated by that distance gets smaller. In this study, all the geostatistical analysis was performed considering the isotropic condition and directional trend was ignored. The computed semivariogram and cross-variogram values were fitted with the theoretical semivariogram models. The best-fit model with the lowest value of residual sum of squares was selected for each soil property. Mathematically, the theoretical semivariogram models used in this study are as follows:
Spherical:
γ ( h ) = C 0 + C [ 3 h 2 A 0 0.5 ( h A 0 ) 3 ] 0     h   A 0 = C 0 + C h   >   A 0
Exponential:
γ ( h ) = C 0 + C [ 1 exp ( h A 0 ) ]
Gaussian:
γ ( h ) = C 0 + C [ 1 exp ( h 2 A 0 2 ) ]
Linear:
γ ( h ) = C 0 + h C A 0
where C0 is the nugget variance and C is the structural variance. C0 + C and A0 are known as the sill and the theoretical range, respectively. In case of the spherical model, effective range is same as the theoretical range, whereas for the exponential and Gaussian models effective range is equal to 3A0 and 3 A0, respectively. Effective range is the distance at which the sill (C0 + C) is within 5% of the asymptote.
Degree of spatial dependence (SD) for each measured property was calculated using the following equation:
SD = C 0 C 0 + C × 100
The spatial dependence of soil properties were classified according to the method detailed by Cambardella et al. [33], which represents the spatial randomness. The spatial dependence is strong if SD < 25%, moderate for SD between 26% and 75% and weak with SD > 75%.
Similar to Pearson’s product moment correlation statistics, Moran’s I [34] statistics for autocorrelation was used to quantify the spatial autocorrelation among the soil properties. In addition to Moran’s I, correlogram was plotted for the analysis of spatial randomness in the measured data. Moran’s I analysis has the index of autocorrelation among group of paired samples separated by increasing lag distance. Similar to correlogram, the values of Moran’s I are in the range of −1 to 1. Positive values of Moran’s I indicate that data values are spatially clustered, whereas negative values of Moran’s I show that neighboring data values are dissimilar [35]. The mathematical form of the index is given by:
I ( h ) = N ( h ) Z i Z i + h Z i + h 2
where I(h) is the Moran’s I for distance interval h, Zi and Zi+h are the measured sample values at the point i and i + h, and N(h) is the number of samples separated by distance h.

3. Results and Discussions

3.1. Descriptive Statistics

Daily precipitation, minimum, and maximum temperatures during the sampling period are shown in Figure 3. The average value of Ksat (0.91 cm min−1 or 1.52 × 10−4 m s−1) is lower than the average values reported in similar studies of the long-term effects of no-tillage, for example, Ksat = 3.5 × 10−7 m s−1 for soils with 25% clay reported in Ferreras et al. [36], 4.2 × 10−6 m s−1 for soils with 19.3% clay in Bhattacharyya et al. [37], and 3.54 × 10−6 m s−1 for soils with 26% clay in Azooz and Arshad [38]. Surface temperature has the lowest variability (CV = 1.87) whereas Ksat has the highest (CV = 98.64) among the measured properties (Table 1). In previous works, Tsegaye and Hill [39], and Warrick and Nielsen [40] showed that Ksat variability is normally high and CV values can reach up to 190%; whereas Mulla and McBratney [41] reported the values of CV for Ksat in the range of 48%–352%.
The average θt value of 0.61 cm3 cm−3 (Table 1) is relatively high as compared to the θt values reported by Fares et al. [42] for Waialua gravelly clay variant amended. These high values of θt can be attributed to the lower compaction of the field and the high plant residue and the root activity after harvesting. No-tillage also maintains large pores [43] inside soil aggregates as compared to conventional tillage which destroys soil aggregates as a result of soil tilth. Hill [44] showed greater soil strength for no-till soils compared to conventionally tilled soils at three different locations. Karlen et al. [45] reported an increase in soil aggregate stability, total carbon, microbial activity, and earthworm populations after eliminating tillage in the U.S. Corn Belt region over a 12-year period.
The average soil bulk density of 1.13 g cm3 (Table 1) is slightly lower than that reported in the literature (1.20 g cm−3) for the top 38 cm soil depth of the Waialua gravelly clay variant [23]. This low bulk density value could be partially due to its shallow sampling depth as bulk density tends to increase with depth [17,20]. Hill [44] reported statistically significant difference in bulk density under no-tillage at three locations when compared with conventional tillage. However, Karlen et al., [45] reported no significant difference in bulk density in top 50 cm depth under chisel plow and no-tillage treatments but there was significantly higher total carbon (almost twice as compared to plow and chisel) under no-tillage in the top 10 cm depth. As presented in Table 1, soil penetration resistance increases with depth and the resistance in the deeper depth (20 cm; TIP-20) showed higher spatial variability (CV = 52.54) as compared to the top 10 cm (TIP-10). Compaction in the lower depth tends to be more uniform due to fewer disturbances as compared to the top soil.
Coefficients of variation for bulk density and porosity were lower than the values of these soil physical properties reported by Mulla and McBratney [41]. According to the relative ranking of magnitude of variability, porosity, and bulk density show low and low to moderate variability, respectively [41]. However, variability of soil organic carbon, soil resistance, and Ksat ranges from moderate to high. Frequency distributions of the measured properties show highly skewed and no-normal distribution for majority of soil parameters (Figure 4). Based on Shapiro–Wilk normality test, Ksat data follows a log-normal distribution; this concurs with the findings of other studies, e.g., Vieira et al., [46]; Wierenga et al. [47]; Russo and Bouton [48]. The volumetric water content and soil resistance at 10 cm depth data follow a normal distribution (Table 1). The volumetric water content and soil resistance data distribution varied as a function of the sampling date; this could be due to temporal effects. However, porosity, bulk density, soil surface temperature, and soil organic carbon seem to follow neither a normal nor a log-normal distribution. Soil properties collected based on the 20 m × 20 m grid on 6 January, 2009 were also skewed except for the volumetric water content and soil resistance at 20 cm (Table 1). Parkin and Robinson [49] reported that many soil properties follow skewed distributions are log-normally distributed. Data of soil properties, which did not exhibit normal distributions, according to Shapiro–Wilk test, were log transformed; however, only Ksat data met the stationarity assumption and showed a normal distribution. Cambardella et al. [33] reported similar results for soil organic carbon in central Iowa soils; the log-normal transformation was only able to reduce skewness, data distribution remained kurtotic.

3.2. Variogram Analysis

Semivariogram parameters (range, nugget, and sill) for each soil property with the best-fitted model are shown in Figure 5, and their spatial dependencies are presented in Table 2. Values of semivariogram fitting parameters and models (Figure 5) varied with sampling date and grid size (Table 2). However, irrespective of sampling date and grid size, soil properties show moderate to strong spatial dependence based on the definition of Cambardella et al. [33]. Nugget defines the micro-scale variability and measurement error, whereas partial sill indicates the amount of variation which can be defined by spatial correlation structure. Out of the total variation, the average nugget component is 35%; it ranges between 15% and 50%, which could be due to high micro-scale variation and/or measurement error. Nugget components for soil bulk density (48%), porosity (43%), and Ksat (38%) were higher than the average nugget component (35%). Theoretically, the value of the range can be considered as the scale of distance beyond which the fitted parameters do not show any spatial correlation [46]. For the 10 m × 10 m grid, value of the range varies between 9.8 m for soil organic carbon and 142.2 m for soil resistance at 10 cm depth. The volumetric water content measured at the 20 m × 20 m grid shows the highest spatial dependence (226.8 m) among all measured properties. Higher values of spatial dependence in volumetric water content can be attributed to the rainfall events occurred just before the field measurements (Figure 3). The goodness of fit measured as indicated by R2 and the residual sum of squares (RSS) between measured and modeled semivariogram for the 10 m × 10 m sampling grid were slightly better than those for the 20 m × 20 m grid (Table 2). Among the different semivariogram models tested, the Gaussian model is found as the best fitting model. The spherical model best-fitted volumetric water content and soil resistance data for the 20 m × 20 m grid and volumetric water content and soil resistance combined samples collected on 17 December, 2008 and 18 December, 2008. The spatial correlation of soil organic carbon, porosity, and ln(Ksat) data followed exponential model.
Cokriging is often recommended for spatial interpolation of measured data due to its smaller estimation variance compared to kriging [50] and its ability to accurately predict the property that is difficult or expensive to measure using well-behaved cross-semivariograms with another more easily sampled property [51]. To assess the feasibility of cokriging, cross-semivariograms were established between measured properties. Before calculating the cross-semivariance, the degree of association between the data of every two studied properties was assessed using Pearson’s correlation matrix (Table 3). There is a strong negative correlation (R = −0.72) between bulk density and porosity, and a strong positive correlation (R = 0.84) between soil resistance at 10 cm and 20 cm depth. The cross-semivariance between (i) porosity and bulk density, (ii) bulk density and soil organic carbon, (iii) Ksat and bulk density, and (iv) Ksat and soil resistance were negative. However, the cross-semivariance between soil resistance and bulk density, Ksat, and porosity, Ksat and soil organic carbon were all positive. For the volumetric water content, soil resistance, and soil temperature, the cross-semivariance was sometimes positive and sometimes negative depending on the date and/or the grid size of sampling (Table 4).

3.3. Analysis of Correlogram

Mulla [52] reported a negative cross-semivariance between gravimetric water content and surface temperature, and gravimetric water content and soil resistance with significant variation in range and partial sill. Variation in cross-semivariance within a small field is not surprising since both surface temperature and soil resistance depend on several properties including water content, soil albedo, solar radiation and soil bulk density, organic matter content, and caliche layer. The inconsistent correlation between soil resistance, soil temperature, and volumetric water content could be partly due to the wet condition of the field during the sampling time (Figure 3). In this study, soil temperature was found to be highly correlated with water content as compared to soil resistance and can be used as secondary variable for cokriging. The modeled variogram was poorly correlated with cross-semivariance between soil resistance and volumetric water content at 20 m × 20 m gird (Table 4) indicating its inappropriate scale for explaining the semivariance.
Moran’s I correlogram analysis (Figure 6) for spatial dependence shows no autocorrelation at any lag distance for Ksat. However, soil organic carbon was positively autocorrelated up to lag of 20 m. Porosity, volumetric water content, bulk density, soil resistance, and soil temperature show linear decreasing trend with point of inflection between 40 and 60 m. With lag distance higher than 60 m, the autocorrelation became negative except for porosity. The highest autocorrelation was associated with the 10 m (sampling grid size) lag distance. Most of these properties did not meet the assumption of stationarity, which induces some uncertainty in these results. However, these results can be used as primary information for designing future sampling plan in similar conditions.

3.4. Spatial Map from Kriging

Spatial maps of measured soil properties prepared using the ordinary kriging are presented in Figure 7. The spatial structure analysis based on kriging indicates a spatial variability across the field for volumetric water content, bulk density, porosity, soil organic carbon, soil resistance, and Ksat. Soil temperature did not show a significant spatial variability as it is expected due to small amplitude of diurnal air temperature at the study site and the Hawaiian Islands in general. The volumetric water content is higher in the downslope areas than in the upslope areas. Values of the soil bulk density and porosity were lower and higher, respectively in the lower portion of the study site. This shows the presence of restricting subsurface layer in the northern part of the field that could be due to the extensive continuous corn cultivation in these areas as compared to lower part of the field. We did not observe significant difference in soil resistance at surface (0–10 cm) and subsurface (10–20 cm) layers.

4. Summary and Conclusions

Soil physical and hydrological properties were measured on grid basis from an agricultural field cropped with sweet corn under a long term no-tillage practice. Classic statistics and geostatistics were used to assess the spatial dependence in several soil physical properties (soil surface temperature, volumetric water content, soil resistance, porosity, bulk density, soil organic carbon, and saturated hydraulic conductivity). Surface temperature has the lowest spatial variability (CV = 1.87) whereas Ksat has the highest spatial variability (CV = 98.64). Values of soil porosity and bulk density spatial variability indicators were relatively higher and lower, respectively, as compared to those of similar soils as reported in the literature. Soil resistance increased with depth indicating compaction of the subsurface layer. The degree of spatial dependence varied from strong to moderate for all tested properties.
The variability of volumetric water content, soil temperature, soil resistance, bulk density, porosity, soil organic carbon, and Ksat exhibited spatial dependence of varying range that can be described using different semivariogram models. The range varied from 9.8 m for soil organic carbon to 142.2 m for soil resistance at 10 cm depth. However, volumetric water content measured at 20 m × 20 m grid had the highest range value (226.8 m). Spatial distribution of the soil properties based on the results of the kriging analysis has shown high level of spatial variability across a relatively small area. These observations suggest that soil physical properties are spatially variable and that spatial variability should not be ignored. Thus, assuming uniform soil properties can lead to mismanagement of natural resources. Furthermore, results of this work revealed that since the scale of spatial variability is not constant across the different soil hydrological and physical properties, it is crucial to use the appropriate scale for the specific practical application. For instance, the degree of the spatial dependency for the volumetric water content is about 227 m. Based on this, it is practical to design field irrigation units for this tropical soil at that level where it is safe to apply constant irrigation amounts at that scale. However, for other agricultural practices such as fertilizer’s applications, it might be prudent to use a finer scale, e.g., 10 m, that corresponds to extend of the maximum of spatial dependency of soil organic carbon content which is an important factor that should be considered in any crop fertility program. Further work is needed to fine-tune some of these practical implications of soil hydrological and physical properties spatial variabilities before recommending them for farmers of the region.

Author Contributions

Conceptualization, A.F. and F.A.; methodology, A.F. and S.K.D.; software, M.S., R.A., and S.F.; validation, M.S. and R.A.; formal analysis, M.S., R.A. and A.F.; investigation, M.S., F.A., S.F., S.K.D., A.A., and A.F.; resources, A.F.; data curation, M.S., R.A. and S.F.; writing, R.A., M.S., and A.F.; review and editing, all co-authors contributed to these tasks; visualization, M.S., R.A. and S.F.; supervision, A.F.; project administration, A.F.; funding acquisition, A.F.

Funding

This project was partially supported by USDA National Institute of Food and Agriculture, Hatch Act and Evans Allen Formula Grants at the University of Hawai’i-Manoa and Prairie View A&M University, respectively.

Acknowledgments

James L. Brewbaker and staff of the Waimanalo Research Station are acknowledged for their cooperation during the experiments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Goovaerts, P. Geostatistical tools for characterizing the spatial variability of microbiological and physico-chemical soil properties. Biol. Fertil. Soil. 1998, 27, 315–334. [Google Scholar] [CrossRef]
  2. Ahuja, L.R. Quantifying agricultural management effects on soil properties and processes. Geoderma 2003, 116, 1–2. [Google Scholar] [CrossRef]
  3. Peet, M. Sustainable Practices for Vegetable Production in the South; Focus Publishing, R. Pullins Co.: Newburyport, MA, USA, 1996. [Google Scholar]
  4. Schertz, D.L. Conservation tillage: An analysis of acreage projections in the United States. J. Soil Water Conserv. 1988, 43, 256–258. [Google Scholar]
  5. Conservation Tillage Information Center (CTIC). National Tillage Trends 1990–2004. Available online: http://www.ctic.purdue.edu/media/pdf/Core4/1990-2004data.pdf (accessed on 29 November 2012).
  6. Claassen, R.; Bowman, M.; McFadden, J.; Smith, D.; Wallander, S. Tillage Intensity and Conservation Cropping in the United States; United States Department of Agriculture, Economic Research Service: Washington, DC, USA, 2018; Volume EIB-197.
  7. Baker, N.T. Tillage Practices in the Conterminous United States, 1989–2004—Datasets Aggregated by Watershed: U.S. Geological Survey Data Series 573; Indiana Water Science Center: Indianapolis, IN, USA, 2011; p. 13.
  8. Cavalieri, K.M.V.; Da Silva, A.P.; Tormena, C.A.; Leão, T.P.; Dexter, A.R.; Håkansson, I. Long-term effects of no-tillage on dynamic soil physical properties in a Rhodic Ferrasol in Paraná, Brazil. Soil Till. Res. 2009, 103, 158–164. [Google Scholar] [CrossRef]
  9. Jakab, G.; Madarász, B.; Szabó, J.A.; Tóth, A.; Zacháry, D.; Szalai, Z.; Kertész, A.; Dyson, J. Infiltration and Soil Loss Changes during the Growing Season under Ploughing and Conservation Tillage. Sustainability 2017, 9, 1726. [Google Scholar] [CrossRef]
  10. Oliver, R.; Douzet, J.M.; Scopel, E.; Blanchart, E.; Curmi, P.; Alves Moreira, J.A.; Minette, S.; Guerin, P.; Fortier, M.; Maraux, F. Medium Term Impact of No Tillage on Some Physical Properties of a Brazilian Oxisol of Cerrados (Tropical Humid Savannah of Central Brazil). In Proceedings of the 17th World Congress of Soil Science: Confronting New Realities in the 21st Century, Bangkok, Thailand, 14–21 August 2002; Maraux, F., Ed.; Ksetsart University: Bangkok, Thailand, 2002; pp. 930-1–930-9. [Google Scholar]
  11. Fabrizzi, K.P.; Garcia, F.O.; Costa, J.L.; Picone, L.I. Soil water dynamics, physical properties and corn and wheat responses to minimum and no-tillage systems in the southern Pampas of Argentina. Soil Till. Res. 2005, 81, 57–69. [Google Scholar] [CrossRef]
  12. Logsdon, S.D.; Karlen, D.L. Bulk density as a soil quality indicator during conversion to no-tillage. Soil Till. Res. 2004, 78, 143–149. [Google Scholar] [CrossRef]
  13. Filipovic, D.; Husnjak, S.; Kosutic, S.; Gospodaric, Z. Effects of tillage systems on compaction and crop yield of Albic Luvisol in Croatia. J. Terramech. 2006, 43, 177–189. [Google Scholar] [CrossRef]
  14. Rhoton, F.E. Influence of time on soil response to no-till practices. Soil Sci. Soc. Am. J. 2000, 64, 700–709. [Google Scholar] [CrossRef]
  15. Bilandžija, D.; Zgorelec, Ž.; Kisic, I. Influence of tillage systems on short-term soil CO2 emissions. Hung. Geogr. Bull. 2017, 66, 29–35. [Google Scholar] [CrossRef]
  16. Dao, T.H. Tillage system and crop residue effects on surface compaction of a Paleustoll. Agron. J. 1996, 88, 141–148. [Google Scholar] [CrossRef]
  17. Hammel, J.E. Long term tillage and crop rotation effects on bulk density and soil impedance in northern Idaho. Soil Sci. Soc. Am. J. 1989, 53, 1515–1519. [Google Scholar] [CrossRef]
  18. Strudley, M.W.; Green, T.R.; Ascough, J.C. Tillage effects on soil hydraulic properties in space and time: State of the science. Soil Till. Res. 2008, 99, 4–48. [Google Scholar] [CrossRef]
  19. Jakab, G.; Dobos, E.; Madarász, B.; Szalai, Z.; Szabó, J.A. Spatial and Temporal Changes in Infiltration and Aggregate Stability: A Case Study of a Subhumid Irrigated Cropland. Water 2019, 11, 876. [Google Scholar] [CrossRef]
  20. Iqbal, J.; Thomasson, J.A.; Jenkins, J.N.; Owens, P.R.; Whisler, F.D. Spatial variability analysis of soil physical properties of alluvial soils. Soil Sci. Soc. Am. J. 2005, 69, 1338–1350. [Google Scholar] [CrossRef]
  21. Webster, R.; Oliver, M.A. Geostatistics for Environmental Scientists; John Wiley & Sons: New York, NY, USA, 2001. [Google Scholar]
  22. Isaaks, E.H.; Srivastava, R.M. An Introduction to Applied Geostatistics; Oxford University press: New York, NY, USA, 1989. [Google Scholar]
  23. Ikawa, H.; Sato, H.H.; Chang, A.K.S.; Nakamura, S.; Robello, E.; Periaswamy, S.P. Soils of the Hawaii Agricultural Experiment Station, University of Hawaii: Soil Survey, Laboratory Data and Soil Descriptions (Benchmark Soils Project Technical Report: 4); University of Hawaii: Honolulu, HI, USA, 1985. [Google Scholar]
  24. Grossman, R.B.; Reinsch, T.G. Bulk Density and Linear Extensibility. In Methods of Soil Analysis, Part 4. Physical Methods. Soil Science Society of Ameriac Book Series No.5; Dane, J.H., Topp, G.C., Eds.; SSSA: Madison, WI, USA, 2002; pp. 201–228. [Google Scholar]
  25. Dane, J.H.; Hopmans, J.W. Soil Water Retention and Storage—Introduction. In Methods of Soil Analysis. Part 4. Physical Methods. Soil Science Society of America Book Series No. 5; Dane, J.H., Topp, G.C., Eds.; SSSA: Madison, WI, USA, 2002; pp. 671–674. [Google Scholar]
  26. Konen, M.E.; Jacobs, P.M.; Burras, C.L.; Talaga, B.J.; Mason, J.A. Equations for predicting soil organic carbon using loss-on-ignition for north central U.S. soils. Soil Sci. Soc. Am. J. 2002, 66, 1878–1881. [Google Scholar] [CrossRef]
  27. Schulte, E.E.; Hopkins, B.G. Estimation of Soil Organic Matter by Weight Loss-on-Ignition. In Soil Organic Matter: Analysis and Interpretation, SSSA Special Publication 46, 1996; Magdoff, F.R., Tabatabai, M.A., Hanlon, E.A., Eds.; SSSA: Madison, WI, USA, 1996; pp. 21–31. [Google Scholar]
  28. APHA. Standard Methods for the Examination of Water and Wastewater, 19th ed.; 2540 E. Fixed and Volatile Solids Ignited at 550 °C; American Public Health Association: Washington, DC, USA; pp. 2–57.
  29. Wooding, R.A. Steady infiltration from a shallow circular pond. Water Resour. Res. 1968, 4, 1259–1273. [Google Scholar] [CrossRef]
  30. Hussen, A.A.; Warrick, A.W. Alternative analysis of hydraulic data from disc tension infiltrometers. Water Resour. Res. 1993, 29, 4103–4108. [Google Scholar] [CrossRef]
  31. Shapiro, S.S.; Wilk, M.B. An analysis of variance test for normality (complete samples). Biometrika 1965, 52, 591–611. [Google Scholar] [CrossRef]
  32. Robertson, G.P. GS+: Geostatistics for the Environmental Sciences; Gamma Design Software: Plainwell, MI, USA, 2008. [Google Scholar]
  33. Cambardella, C.A.; Moorman, T.B.; Novak, J.M.; Parkin, T.B.; Karlem, D.L.; Turco, R.F.; Konopka, A.E. Field scale variability of soil properties in central Iowa soil. Soil Sci. Soc. Am. J. 1994, 58, 1501–1511. [Google Scholar] [CrossRef]
  34. Moran, P.A.P. Notes on continuous stochastic phenomena. Biometrika 1950, 37, 17–23. [Google Scholar] [CrossRef]
  35. Bruland, G.L.; Richardson, C.J. Spatial variability of soil properties in created, restored, and paired natural wetlands. Soil Sci. Soc. Am. J. 2005, 69, 273–284. [Google Scholar]
  36. Ferreras, L.A.; Costa, J.L.; Garcia, F.O.; Pecorari, C. Effect of no-tillage on some soil physical properties of a structural degraded Petrocalcic Paleudoll of the southern “Pampas” of Argentina. Soil Till. Res. 2000, 54, 31–39. [Google Scholar] [CrossRef]
  37. Bhattacharyya, R.; Prakash, V.; Kundu, S.; Gupta, H.S. Effect of tillage and crop rotations on pore size distribution and soil hydraulic conductivity in sandy clay loam soil of the Indian Himalayas. Soil Till. Res. 2006, 86, 129–140. [Google Scholar] [CrossRef]
  38. Azooz, R.H.; Arshad, M.A. Soil infiltration and hydraulic conductivity under long-term no-tillage and conventional tillage systems. Can. J. Soil Sci. 1996, 76, 143–152. [Google Scholar] [CrossRef]
  39. Tsegaye, T.; Hill, R.L. Intensive tillage effects on spatial variability of soil physical properties. Soil Sci. 1998, 163, 143–154. [Google Scholar] [CrossRef]
  40. Warrick, A.W.; Nielsen, D.R. Spatial Variability of Soil Physical Properties in the Field. In Applications of Soil Physics; Hillel, D., Ed.; Academic Press: New York, NY, USA, 1980; pp. 319–344. [Google Scholar]
  41. Mulla, D.J.; McBratney, A.B. Soil Spatial Variability. In Soil Physics Companion; Warrick, A.W., Ed.; CRC Press: Boca Raton, FL, USA, 2002; pp. 343–373. [Google Scholar]
  42. Fares, A.; Abbas, F.; Ahmed, A.; Deenik, J.; Safeeq, M. Response of selected Soil Physical and Hydrological Properties to Manure Amendment Rates, Levels, and Types. Soil Sci. 2008, 173, 522–533. [Google Scholar] [CrossRef]
  43. Ankeny, M.D.; Kaspar, T.C.; Horton, R. Characterization of tillage and traffic effects on unconfined infiltration measurements. Soil Sci. Soc. Am. J. 1990, 54, 837–840. [Google Scholar] [CrossRef]
  44. Hill, R.L. Long-Term Conventional and No-Tillage Effects on Selected Soil Physical-Properties. Soil Sci. Soc. Am. J. 1990, 54, 161–166. [Google Scholar] [CrossRef]
  45. Karlen, D.L.; Wollenhaupt, N.C.; Erbach, D.C.; Berry, E.C.; Swan, J.B.; Eash, N.S.; Jordahl, J.L. Long-term tillage effects on soil quality. Soil Till. Res. 1994, 32, 313–327. [Google Scholar] [CrossRef]
  46. Vieira, S.R.; Reynolds, W.D.; Topp, G.C. Spatial Variability of Hydraulic Properties in a Highly Structured Clay Soil. In Proceedings of the International Conference and Workshop on the Validation of Flow and Transport Models for the Unsaturated Zone, Ruidoso, NM, USA, 22–25 May 1988; NM State University: Las Cruces, NM, USA, 1988; pp. 471–483. [Google Scholar]
  47. Wierenga, P.J.; Hills, R.G.; Hudson, D.B. The Las Cruces trench site: Characterization, experimental results, and one-dimensional flow predictions. Water Resour. Res. 1991, 27, 2695–2705. [Google Scholar] [CrossRef]
  48. Russo, D.; Bouton, M. Statistical analysis of spatial variability in unsaturated flow parameters. Water Resour. Res. 1992, 28, 1911–1925. [Google Scholar] [CrossRef]
  49. Parkin, T.B.; Robinson, J.A. Analysis of lognormal data. Adv. Soil Sci. 1992, 20, 193–235. [Google Scholar]
  50. Vauclin, M.; Vieira, S.R.; Vachaud, G.; Nielsen, D.R. The use of cokriging with limited field soil observations. Soil Sci. Soc. Am. J. 1983, 47, 175–184. [Google Scholar] [CrossRef]
  51. Yates, S.R.; Warrick, A.W. Estimating soil water content using cokriging. Soil Sci. Soc. Am. J. 1987, 51, 23–30. [Google Scholar] [CrossRef]
  52. Mulla, D.J. Estimating spatial patterns in water content, matric suction, and hydraulic conductivity. Soil Sci. Soc. Am. J. 1988, 52, 1547–1553. [Google Scholar] [CrossRef]
Figure 1. A diagram and map of the study area at the University of Hawaii at Mānoa Waimanalo Research Station; it shows sampling grids (Empty circles—10 m × 10 m sampling grids, circle with red dots—20 m × 20 m sampling grids).
Figure 1. A diagram and map of the study area at the University of Hawaii at Mānoa Waimanalo Research Station; it shows sampling grids (Empty circles—10 m × 10 m sampling grids, circle with red dots—20 m × 20 m sampling grids).
Agronomy 09 00750 g001
Figure 2. Relationship between loss on ignition (LOI) and soil organic carbon for the Waialua gravelly clay variant soil in Waimanalo, HI.
Figure 2. Relationship between loss on ignition (LOI) and soil organic carbon for the Waialua gravelly clay variant soil in Waimanalo, HI.
Agronomy 09 00750 g002
Figure 3. Measured values of daily precipitation, and minimum (Tmin) and maximum (Tmax) air temperatures at the study site between 1 December 2008 and 31 January 2009.
Figure 3. Measured values of daily precipitation, and minimum (Tmin) and maximum (Tmax) air temperatures at the study site between 1 December 2008 and 31 January 2009.
Agronomy 09 00750 g003
Figure 4. Frequency distribution of selected soil properties sampled at the 10 m × 10 m grid on 17–18 December 2008.
Figure 4. Frequency distribution of selected soil properties sampled at the 10 m × 10 m grid on 17–18 December 2008.
Agronomy 09 00750 g004
Figure 5. Semivariance values (circle) and the best fitting semivariogram models (solid line) for volumetric water content (θv), soil resistance at 10 cm depth (TIP-10), bulk density, soil resistance at 20 cm depth (TIP-20), porosity, soil temperature, soil organic carbon, and log (Ksat).
Figure 5. Semivariance values (circle) and the best fitting semivariogram models (solid line) for volumetric water content (θv), soil resistance at 10 cm depth (TIP-10), bulk density, soil resistance at 20 cm depth (TIP-20), porosity, soil temperature, soil organic carbon, and log (Ksat).
Agronomy 09 00750 g005
Figure 6. Moran’s I correlograms for volumetric water content (θv), soil resistance at 10 cm depth (TIP-10), bulk density, soil resistance at 20 cm depth (TIP-20), porosity, soil temperature, soil organic carbon, and log (Ksat) studied at 10 m × 10 m grid on 17–18 December 2008.
Figure 6. Moran’s I correlograms for volumetric water content (θv), soil resistance at 10 cm depth (TIP-10), bulk density, soil resistance at 20 cm depth (TIP-20), porosity, soil temperature, soil organic carbon, and log (Ksat) studied at 10 m × 10 m grid on 17–18 December 2008.
Agronomy 09 00750 g006
Figure 7. Maps for soil temperature (ST), volumetric water content (θv), bulk density (ρb), porosity (θt), saturated hydraulic conductivity (Ksat), soil organic carbon (SOC) content (%), soil resistance at 10 cm depth (TIP-10), and soil resistance at 20 cm depth (TIP-20).
Figure 7. Maps for soil temperature (ST), volumetric water content (θv), bulk density (ρb), porosity (θt), saturated hydraulic conductivity (Ksat), soil organic carbon (SOC) content (%), soil resistance at 10 cm depth (TIP-10), and soil resistance at 20 cm depth (TIP-20).
Agronomy 09 00750 g007
Table 1. Descriptive statistics and Shapiro–Wilk normality test results (grouped by sampling date and grid size) for measured soil physical properties: soil temperature (ST), volumetric water content (θv), soil resistance at 10 cm depth (TIP-10), soil resistance at 20 cm depth (TIP-20), bulk density (ρb), porosity (θt), soil organic carbon (SOC), and saturated hydraulic conductivity (Ksat).
Table 1. Descriptive statistics and Shapiro–Wilk normality test results (grouped by sampling date and grid size) for measured soil physical properties: soil temperature (ST), volumetric water content (θv), soil resistance at 10 cm depth (TIP-10), soil resistance at 20 cm depth (TIP-20), bulk density (ρb), porosity (θt), soil organic carbon (SOC), and saturated hydraulic conductivity (Ksat).
Sampling DateVariableMeanVarCV
(%)
MaxMinNShapiro–Wilk Normality Test
10 m × 10 m Grid
12/17/2008ST (°C)23.290.191.8724.0022.03960.000
θv (%)59.4789.7515.9376.2735.00960.085
TIP-10 (kg cm−2)5.508.3652.5412.890.70960.111
TIP-20 (kg cm−2)8.0412.1443.3616.870.94960.235
12/18/2008ST (°C)23.301.054.3926.8021.101300.000
θv (%)67.7242.949.6880.4044.131300.003
TIP-10 (kg cm−2)6.395.3936.3213.360.941300.255
TIP-20 (kg cm−2)8.968.4532.4518.281.641300.000
12/6/2008θt (cm3 cm−3)0.610.009.700. 720. 322230.035
ρb (g cm−3)1.130.0212.751.450.772150.001
SOC (%)1.610.8256.234.520.322280.000
20 m × 20 m Grid
1/6/2009ST (°C)23.562.186.2725.0020.701280.000
θv (%)58.0983.6515.7587.8036.971280.425
TIP-10 (kg cm−2)7.6010.4442.5114.061.521280.010
TIP-20 (kg cm−2)11.8313.0030.4920.863.661280.664
12-6/12-10 2008Ksat (cm min−1)0.910.8098.643.820.02670.000/0.348 *
* Log transformed data, Var = variance, CV = coefficient of variation, Max = maximum, Min = minimum, N = number of samples.
Table 2. Summary of the semivariogram models for: soil temperature (ST), volumetric water content (θv), soil resistance at 10 cm depth (TIP-10), soil resistance at 20 cm depth (TIP-20), bulk density (ρb), porosity (θt), soil organic carbon (SOC), and saturated hydraulic conductivity (Ksat).
Table 2. Summary of the semivariogram models for: soil temperature (ST), volumetric water content (θv), soil resistance at 10 cm depth (TIP-10), soil resistance at 20 cm depth (TIP-20), bulk density (ρb), porosity (θt), soil organic carbon (SOC), and saturated hydraulic conductivity (Ksat).
Sampling DateVariablesModel ‡C0C+C0A0SD †R2RSS
10 m ×10 m Grid
12/17 and 12/18/2008θvSph.18.8883.54123.80S0.9850.50
12/17 and 12/18/2008TIP-10Sph.3.587.4871.50M0.970.280
12/17 and 12/18/2008TIP-20Sph.4.0211.3467.30M0.961.398
12/6/2008θtExp.0.00170.004050.90M0.960.00
12/6/2008ρbGau.0.0100.02152.40M0.980.00
12/6/2008SOCExp.0.1320.8939.80S0.940.003
12/17 and 12/18/2008STGau.0.200.8949.40S1.000.00
12/17/2008θvGau.23.10140.5058.30S0.9931.80
12/17/2008TIP-10Gau.5.0513.2854.30M0.990.15
12/17/2008TIP-20Gau.6.3417.1043.50M0.981.22
12/17/2008STGau.0.130.2650.50M0.950.00
12/18/2008θvExp.27.8755.7582.00M0.9111.60
12/18/2008TIP-10Exp.3.459.43142.20M0.930.21
12/18/2008TIP-20Gau.5.4411.9572.20M0.990.09
12/18/2008STGau.0.271.8067.90S0.980.02
20 m × 20 m Grid
1/6/2009θvSph.28.40127.80226.80S0.85417.00
1/6/2009TIP-10Sph.0.3611.7252.00S0.869.44
1/6/2009TIP-20Sph.1.9214.6554.00S0.8315.18
1/6/2009STGau.0.072.2011.70S0.002.63
12/6 and 12/10/2008Ksat*Exp.0.5171.36632.5M0.670.131
† S = Strong spatial dependence; M = Moderate spatial dependence. ‡ Sph. = Spherical; Exp. = Exponential; Gau. = Gaussian. * Log transformed data; C0 = nuget; C + C0 = sill; A0 = range in m; R2 = coefficient of determination; RSS = residual sum of squares.
Table 3. Pearson correlation between selected parameters: saturated hydraulic conductivity (Ksat), soil temperature (ST), porosity (θt), bulk density (ρb), volumetric water content (θv), soil resistance at 10 cm depth (TIP-10), soil resistance at 20 cm depth (TIP-20), and soil organic carbon (SOC).
Table 3. Pearson correlation between selected parameters: saturated hydraulic conductivity (Ksat), soil temperature (ST), porosity (θt), bulk density (ρb), volumetric water content (θv), soil resistance at 10 cm depth (TIP-10), soil resistance at 20 cm depth (TIP-20), and soil organic carbon (SOC).
Ksat
(cm min−1)
θt
(cm3 cm−3)
ρb
(g cm−3)
θv
(%)
TIP-10
(kg cm−2)
TIP-20
(kg cm−2)
SOC
(%)
Ksat (cm min−1)1
θt (cm3 cm−3)0.051
ρb (g cm−3)−0.04−0.711
θv (%)0.10−0.150.241
TIP-10 (kg cm−2)−0.29 *0.21−0.13−0.171
TIP-20 (kg cm−2)−0.25 *−0.030.12−0.030.81 **1
SOC (%)−0.040.13−0.150.09−0.08−0.151
* Significant at the 0.05 level, ** Significant at the 0.01 level.
Table 4. Summary of cross-variogram model.
Table 4. Summary of cross-variogram model.
Sampling DateVariables †Model ‡C0C + C0A0R2RSS
10 m × 10 m Grid
12/17 and 12/18/2008TIP-10/θvGau.−0.01−3.3049.100.920.74
12/17 and 12/18/2008TIP-20/θvGau.−0.94−3.1654.500.820.80
12/6/2008θt/ρbGau.0.00−0.0143.800.980.00
SOC/ρbGau.−0.013−0.0568.800.830.00
SOC/θtGau.0.0070.0273.600.860.00
TIP-10/ρbGau.0.000.7177.400.830.08
12/17/2008ST/θvGau.−0.001−1.4748.400.670.90
12/17/2008TIP-10/θvGau.−0.69−16.2161.700.968.99
12/17/2008TIP-20/θvGau.−1.48−16.0778.400.8934.22
12/18/2008ST/θvGau.0.0013.0151.700.980.22
12/18/2008TIP-10/θvGau.−0.01−2.9078.800.984.12
12/18/2008TIP-20/θvGau.−0.01−4.9467.100.9417.55
20 m × 20 m Grid
1/6/2009TIP-10/θvGau.4.636.6193.540.0468.60
1/6/2009TIP-20/θvGau.3.566.3893.540.1053.82
Ksat/ρbExp.0.00−0.036.200.040.00
KsattGau.0.000.0117.000.300.00
Ksat/SOCGau.0.0010.32139.500.300.03
Ksat/TIP-10Lin.−0.047−1.16893.430.6860.409
Ksat/TIP-20Gau.−0.633−3.106188.10.4030.425
† Saturated hydraulic conductivity (Ksat), soil temperature (ST), porosity (θt), bulk density (ρb), volumetric water content (θv), soil resistance at 10 cm depth (TIP-10), soil resistance at 20 cm depth (TIP-20), and soil organic carbon (SOC). ‡ Sph. = Spherical; Exp. = Exponential; Gau. = Gaussian; Lin. = Linear; C0 = nuget; C + C0 = sill; A0 = range in m; R2 = coefficient of determination; RSS = residual sum of squares.

Share and Cite

MDPI and ACS Style

Awal, R.; Safeeq, M.; Abbas, F.; Fares, S.; Deb, S.K.; Ahmad, A.; Fares, A. Soil Physical Properties Spatial Variability under Long-Term No-Tillage Corn. Agronomy 2019, 9, 750. https://doi.org/10.3390/agronomy9110750

AMA Style

Awal R, Safeeq M, Abbas F, Fares S, Deb SK, Ahmad A, Fares A. Soil Physical Properties Spatial Variability under Long-Term No-Tillage Corn. Agronomy. 2019; 9(11):750. https://doi.org/10.3390/agronomy9110750

Chicago/Turabian Style

Awal, Ripendra, Mohammad Safeeq, Farhat Abbas, Samira Fares, Sanjit K. Deb, Amjad Ahmad, and Ali Fares. 2019. "Soil Physical Properties Spatial Variability under Long-Term No-Tillage Corn" Agronomy 9, no. 11: 750. https://doi.org/10.3390/agronomy9110750

APA Style

Awal, R., Safeeq, M., Abbas, F., Fares, S., Deb, S. K., Ahmad, A., & Fares, A. (2019). Soil Physical Properties Spatial Variability under Long-Term No-Tillage Corn. Agronomy, 9(11), 750. https://doi.org/10.3390/agronomy9110750

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop