- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Applied and Environmental Soil Science
Volume 2011 (2011), Article ID 198737, 7 pages
Field Scale Studies on the Spatial Variability of Soil Quality Indicators in Washington State, USA
1USDA-Agricultural Research Service, 215 Johnson Hall, Washington State University, Pullman, WA 99164-6421, USA
2USDA-Agricultural Research Service, 1224 Airport Road, Beaver, WV 25813, USA
Received 21 October 2010; Accepted 26 January 2011
Academic Editor: Marco Trevisan
Copyright © 2011 Jeffrey L. Smith and Jonathan J. Halvorson. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Arable lands are needed for sustainable agricultural systems to support an ever-growing human population. Soil quality needs to be defined to assure that new land brought into crop production is sustainable. To evaluate soil quality, a number of soil attributes will need to be measured, evaluated, and integrated into a soil-quality index using the multivariable indicator kriging (MVIK) procedure. This study was conducted to determine the spatial variability and correlation of indicator parameters on a field scale with respect to soil quality and suitability for use with MVIK. The variability of the biological parameters decreased in the order of respiration > enzyme assays and qCO2 > microbial biomass C. The distribution frequency of all parameters except respiration were normal although the spatial distribution across the landscape was highly variable. The biological parameters showed little correlation with each other when all data points were considered; however, when grouped in smaller sections, the correlations were more consistent with observed patterns across the field. To accurately assess soil quality, and arable land use, consideration of spatial and temporal variability, soil conditions, and other controlling factors must be taken into account.
The challenge of feeding 9 billion people by the year 2050 is intimidating. Multiple strategies are needed to meet this challenge. Strategies include reducing human population growth, decreasing protein consumption, increasing crop and animal production, and increasing the agricultural land base for production. While some food staples (crops and livestock) are increasing, others are static or decreasing . Currently food production for 6 billion people occurs on 13% of the global land surface .
The current strategical focus has been on increasing yields and increasing the agricultural land base [3, 4]. The challenge to moving these strategies forward is to evaluate (1) new land for the ability to produce crops and (2) the soil’s resilience and resistance to degrade over time from agriculture management. Thus, the concept of soil quality needs to be developed to evaluate and manage land developed for increased crop production [5, 6]. In a global context, soil quality affects not only soil productivity but is also a significant factor governing environmental quality, human and animal health, and food safety and quality . Soil quality of current and future agricultural land is of similar importance to humankind as air and water quality; thus, it is apparent that simply protecting soil quality by slowing soil degradation or maintaining the current level of soil health will not provide the soil quality that will be needed for future generations. Soil quality must be improved as well.
To properly assess soil quality, appropriate soil indicators or properties must be identified. These indicators must be quantified on a local and landscape basis as a means for making small-scale and regional management decisions. Indicators proposed to asses soil quality are diverse and include chemical, physical, and biological variables as well as descriptive terms . Once the indicators are identified, methods need to be developed to integrate the indicators into a soil quality index. This integration procedure will require the use of information from all indicators as well as their interactions . Several years ago, a procedure was developed to spatially evaluate soil quality. The procedure, multiple variable indicator kriging (MVIK), categorizes numerous variables based on specific criteria producing a probability that an “area” or “soil” meets qualified standards . For example, criteria can be developed and parameters measured and integrated to evaluate if a land area is suitable for development of irrigated agriculture. The MVIK method has been used to map the distribution of soil nutrients , define zones of soil pollutants , evaluate chemical health risks of groundwater , evaluate soil quality indices [6, 13], classify soil degradation in agricultural lands , and assess groundwater quality for irrigated agriculture .
The underlying principle of MVIK is the spatial variability of indicator parameters across the landscape. In this study, we have chosen several chemical and microbiological parameters that have been shown to be important in agriculture and are considered good indicators of soil quality [8, 9]. These indicators, not meant to be inclusive, include microbial biomass, respiration, metabolic quotient (qCO2), dehydrogenase, and phosphatase enzyme activity as well as pH, organic carbon (%C), and electrical conductivity (EC). We explore the spatial variability and relationships of indicators across a small landscape that was intensively sampled. Knowing the range and spatial variability criteria for MVIK, we can develop a system for the evaluation of suitable land for crop production.
2. Materials and Methods
2.1. Soil Collection
Soil samples were taken from a slightly north-sloping agricultural field in Southeastern Washington State, USA. The field had previously been cropped to winter wheat (Triticum aestivum L.). The soil is a Palouse silt loam classified as a fine-silty mixed, superactive, Pachic Ultic Haploxeroll having an average organic C content of 1.5% and a total N content of 0.15%.
The climate of this region is characterized by cold wet winters and hot dry summers with 80% of the annual 500 mm of precipitation being received between October and March. Soil samples (220) were collected in October (Fall) from an area of approximately 0.5 ha (50 × 110 m). The sampling design was a regular 10 × 10 m grid with smaller sampling distances randomly placed throughout the larger grid. At each sampling location, approximately 200 g of soil was collected from the top 10 cm, placed in a plastic bag and stored at 4°C until analyzed.
The 220 soil samples were analyzed for total organic C, pH, electrical conductivity (EC), microbial biomass C (SIRC), basal respiration, phosphatase (PNP), dehydrogenase (TPF), and calculated metabolic quotient (qCO2). Analyses were done in triplicate.
2.2. Biological Assays
Microbial biomass was assayed using the substrate induced respiration (SIR) method developed by Anderson and Domsch . Triplicate 10 g (dry weight) samples of soil were weighed into 40 ml glass vials, brought to 20% moisture (w/w), covered and kept in the dark for 7 days at °C. After the preincubation, each sample was amended with a saturating level of 600 ug glucose (240 mg C) g−1 soil, bringing the soil moisture content to 25% (w/w) (−30 kPa). Water only was added to each control sample also bringing the soil to 25% (w/w) (−30 kPa). All tubes were flushed with hydrated air and capped with a septum. Headspace CO2 was measured from each vial at 3 and 24 h by injecting 0.2 ml of headspace in to a gas chromatograph. Microbial biomass carbon (SIRC) was calculated using the equation published by Anderson and Domsch . Soil respiration (RESP) rate was calculated as the average rate of CO2 production per hour from the control samples (only water added).
Soil dehydrogenase and phosphatase activity was measured using the modified technique of Bolton et al. . The substrate for dehydrogenase activity was 2,3,5-triphenyltetrazolium (3% w/v) and p-nitrophenol phosphate for phosphatase activity. Enzyme assays were made in triplicate on moist soil (10 g dry weight) and reported for dehydrogenase as 10−5 μmol triphenylformazan g−1 min−1 (TPF) and for phosphatase as 10−2 μmol p-nitrophenol g−1 min−1 (PNP).
Metabolic quotient (qCO2) is defined as the basal respiration per unit of microbial biomass in units of μg CO2-C μg−1 biomass-C h−1 and was calculated from the control sample respiration rate at the end of the incubation period divided by the SIRC biomass .
2.3. Chemical Analysis
Total organic C was measured in triplicate on each of the 220 soil samples by a wet oxidation method . The pH and electrical conductivity (EC) measurements were made on saturated pastes of soil (1 : 1 w/w, soil : water).
2.4. Statistical Methods
Summary statistics were calculated for the 220 samples including mean, median, and standard deviation (SD), coefficient of variation (%CV), and skewness. In addition, Pearson correlation coefficients were calculated for untransformed data. We also constructed maps of the spatial distributions for each individual variable over the 0.5 ha area from untransformed data using a weighted least squares algorithm. Every point at an unsampled location is calculated by a weighted quadratic multiple regression on all the points .
3.1. Chemical and Biological Indicators
Univariate statistics for measured soil properties are shown in Table 1. For the chemical properties the mean and median values were similar and the skewness values were low to moderate (EC). The relatively high variability of total C (%CV = 39.6) on this landscape was related to spatial patterns. We observed lower values on the eastern 25% of the field (see Figure 2). However, in general, the %CVs for chemical properties on a landscape area are typical .
All biological variables, except microbial respiration, also showed similar mean and median values with low to moderate skewness (Table 1). The sample distribution for respiration exhibited significant skewness (4.35) due to a few high outlier values. This outlier effect is confirmed by the 75th percentile value (1.76 mg kg−1) being close to the mean value of 1.65 mg kg−1 (data not shown) and also reflected in the CV of 71%. Microbial biomass C ranged from 227 to 1694 mg C kg−1 soil with a CV of 31%, typical for this biological parameter . Values for metabolic quotient (qCO2) ranged over an order of magnitude, from 0.12 to 2.61 × 10-3 μg CO2-C μg−1 biomass-C h−1, with a CV of 42%. The enzyme assays, phosphatase and dehydrogenase, showed similar variation among samples with CVs of 33 and 36%, respectively.
3.2. Spatial Mapping and Correlation
The distribution of electrical conductivity (EC), dehyrogenase (TPF), and respiration (RESP) over the 0.5 ha sampling area are depicted in Figure 1. Most sample respiration rates were near the overall mean of 1.65 mg-C kg−1 h−1 (Figure 1, RESP). However, we also observed a zone of higher values in the 60–70 m Easting area and a few other elevated values at about 85–100 m Easting. In contrast, dehyrogenase activity was distributed throughout the field in zones of high and low activity (Figure 1, TPF), much different than what might have been predicted from the univariate statistics (Table 1) [23, 24]. Highest TPF activity was observed from the 0 m to 20 m Easting, the portion of the field that also had the lowest respiration values. Dehyrogenase activity was lower than the mean at random locations mostly in the north half of the field between 30 and 100 m Easting. Values of EC, an important soil variable since soluble salts may affect microbial and enzyme activity [25, 26] were generally similar to the overall field mean except in the 110 m easting area and in a few locations randomly distributed across the field (Figure 1, EC). The trend in EC was generally opposite of TPF and, in some areas, seemed to be similar in pattern to RESP.
Figure 2 shows the spatial distribution of total C (%C), pH, and microbial biomass (SIRC). Microbial biomass was characterized by a fairly uniform spatial distribution with many values close to the field mean (Figure 2 SIRC). However, we also observed a significant area of the field with values 50% greater than the mean (green). Soil pH was less evenly spatially distributed across the field than SIRC (Figure 2 pH). We observed a zone with relatively high values of pH between 60 and 100 m Easting. Lower values of pH in the eastern edge of the field also coincided with a decrease in SIRC and %C (Figure 2 %C). Total C was distributed with highest concentrations in the west half of the field and lowest concentrations in the east. However, the grading from high to low concentration was not smooth, and the field appears to be divided into two patches. Values near the overall field average (green) occurred mostly in a relatively narrow band with a NW-SW orientation, near the center of the field.
Table 2 provides the Pearson correlation coefficients () for the measured soil properties. We calculated for two spatial scales, full field or subfield. Full field calculations represent an overall value based on all 220 data points. Subfield calculations were based on the data from each of four 27.5 m W to E subsections of the field from 0 to 110 m. Analyses of the entire field generally showed poor correlation’s between variables except between RESP and qCO2 and SIRC and TPF (Table 2). Transformation (log10) the data did not improve the overall correlation’s with the only increase due to normalizing RESP.
Separate analyses of the 27.5 m subfields revealed areas of relatively high and low correlation obscured at the scale of the entire field. For example, the full field correlation coefficient of SIRC and RESP was 0.15 (, ); however, these two variables were significantly and positively correlated in two subfields lying between 0 and 55 m Easting, section of the field, (0.48, , and .54, , ), respectively. In a more extreme example, the overall field correlation of SIRC and PNP was 0.02 (, ), while correlation coefficients in the subfields varied: near 0 in the 0 to 27.5 m section of the field, slightly positive in the 27.5 to 55 m subfield, slightly negative in the 55.0 to 82.5 m subfield, and significantly positive in the 82.5 to 110 m subfield.
By accounting for this intrafield heterogeneity, the correlations between variables were more similar to the qualitative, visual analysis of Figures 1 and 2. For example, RESP and EC exhibited some positive correlation (0.29, , ), in the 0 to 27.5 m Easting subfield (Table 2, Figure 1). However, in the subfield located from 27.5 to 55.0 m, where EC was increasing, the correlation coefficient decreased to 0.10 (, ) (Figure 1). In the east half of the field (55+ m Easting), correlation between RESP and EC became increasingly negative (−0.14, , and −0.18, , ) due to opposing spatial gradients shown in Figure 1. Soil microbial biomass was uncorrelated with %C over the entire field (Table 2, Figure 2); however, SIRC was significantly correlated with pH in all except in the 55.0 to 82.5 m subfield (Table 2) also suggested by Figure 2 which shows SIRC as constant over the landscape and pH as increasing.
The objectives of developing a method to evaluate soil quality include (1) the ability to monitor changes in the soil over time and be able to quantify the direction and rate of change, (2) evaluation of the time rate of change in soil quality due to specific management in both the short and long term, and (3) the ability to detect patterns of soil quality at different spatial and temporal scales. To meet these objectives the evaluation methodology, such as MVIK, must be based on sensitive indicators of change in soil quality that are stable and predictable with respect to other soil properties such that comparisons between different locations and between different time periods can be made.
Several biological indicators of soil quality have been proposed because they are thought to be integrative variables, sensitive to changes in soil degradation or improvement [27–30]. However, we are interested in mapping soil quality on the landscape level to monitor change or implement remediation, and there are several problems associated with the use of only biological indicators over field to landscape scales.
First, biological indicators are likely to be more variable in time and space than chemical or physical parameters and may also be influenced by these parameters . Second, univariate statistics of biological indicators may be of limited use in the context of evaluating soil quality, because they only characterize a system at the overall scale and do not account for spatial patterns . Finally, even if we account for the spatial and temporal variability of biological indicators, the precise relationship between the amount of a particular biological indicator and it’s functional characteristics has yet to be determined .
These problems can become critical when using the MVIK procedure to evaluate specific criteria. The MVIK procedure uses cutoff values to determine if a criteria has been met, thus if indicators vary greatly in time, the cutoff value may not be stable and land evaluation may change.
The coefficient of variation for the biological variables measured in this study ranged from 31% for microbial biomass-C and 30% to 40% for the enzyme assays to 71% for soil respiration. Rochette et al.  found respiration measurements in a 1 ha field (average 50 sampling points) to have a CV ranging from 25% to 69% over the growing season. Bonmati et al.  found the CV for enzymes ranged from 28 to 60% in a small field (<0.1 ha). This is in contrast to many chemical and physical parameters that typically have a CV in the range of 5% to 25% on a spatial basis. While mean values or the relative ranking of variability may not change, actual estimates of statistical dispersion (e.g., standard deviation) will be affected by the size of the sampling unit especially for soil properties with strongly skewed distributions .
Understanding the spatial distribution of chemical and physical indicators of soil quality is likely to be important for corroborating biological data and necessary to help to explain the spatial variability of biological parameters. For example, Figure 2 shows areas of high and low concentrations of three variables occurring simultaneously in a single field. These distinct areas maybe caused by secondary or tertiary factors affecting the primary variable of interest. Greater correlation between soil quality indicators in some sections of the field (Table 2) suggests the factors controlling these variables change over the landscape. For example, unlike the rest of the field, SIRC and pH were uncorrelated in the 55 to 82 m Easting section of the field (Table 2). Since the values of SIRC were distributed, on the average, similarly across the entire field, the lower correlation is likely due to increasing pH values in this section of the field seen in Figure 2. The gradient of pH across the field may have been related to patterns of NO3 and soil moisture content, which decrease from W to E across the field (data not presented). These two soil variables, individually or in combination, may have influenced the increase in pH, and thus its correlation with microbial biomass. Thus, successful use of biological indicators for soil quality analysis will require an understanding their spatial variability as well as that of any underlying factors that might affect this variability.
Our data further illustrate that while the overall univariate distribution of a biological parameter may be normal for a field-size sample unit (Table 1), the variable may be distributed in zones of high and low concentrations within the field (i.e., TPF, Figure 1). Thus the full field mean value may be of little value for guiding management decisions. It would be impossible to manage this field to increase its biological potential based on univariate statistics alone. In order to be useful, we are required to sample in a manner that will allow us to identify and delineate the significant spatial patterns within the field. Contour maps for each soil parameter of interest are desirable but require more sampling and analysis meaning greater expense.
Two other aspects of soil quality indicators can affect land evaluation even using the MVIK procedure. Firstly, evaluating soil quality based on the quantity of a biological indicator assumes that laboratory analyses accurately represent conditions in the field. Secondly, a critical assumption is that the relationship between the quantity of a biological indicator and its activity is well understood and predictable. Yet, little attempt has been made to distinguish between these two facets of soil quality.
Thus based on this study we suggest that soil sampling for soil quality analysis using biological indicators should be conducted when the soils are between 50 to 70% of field capacity (−30 kPa), during mild temperature regimes (Fall or late Spring) and that enough samples be taken that contour maps can be developed. Other chemical data should be collected at the same time, such as pH and electrical conductivity, to provide collaborative information for the soil quality assessment. Biological indicators of soil quality are spatially variable; however, they are also sensitive indicators, thus with proper analysis and interpretation they will enhance our development of a soil quality index.
That biological indicators of soil quality should vary across the landscape or during the course of a season is neither novel nor surprising. This is why we developed the MVIK procedure to smooth the variability into a probably that can be used to define good soil quality or good arable land. However, the more stable the data used in the MVIK procedure the better the prediction of arable land. With the spatial variability information on an area basis we can use procedures such as MVIK to evaluate the probably that a new land base will be sustainable for agricultural production.
- H. C. J. Godfray, J. R. Beddington, I. R. Crute et al., “Food security: the challenge of feeding 9 billion people,” Science, vol. 327, no. 5967, pp. 812–818, 2010.
- J. D. Glover, J. P. Reganold, L. W. Bell et al., “Increased food and ecosystem security via perennial grains,” Science, vol. 328, no. 5986, pp. 1638–1639, 2010.
- Anonomous, “Investing in soils,” Nature Geoscience, vol. 3, p. 295, 2010.
- Anonomous, “How to feed a hungry world,” Nature, vol. 466, pp. 531–532, 2010.
- S. S. Andrews, D. L. Karlen, and C. A. Cambardella, “The soil management assessment framework: a quantitative soil quality evaluation method,” Soil Science Society of America Journal, vol. 68, no. 6, pp. 1945–1962, 2004.
- F. Bastida, A. Zsolnay, T. Hernández, and C. García, “Past, present and future of soil quality indices: a biological perspective,” Geoderma, vol. 147, no. 3-4, pp. 159–171, 2008.
- J. F. Parr, R. I. Papendick, S. B. Hornick, and R. E. Meyer, “Soil quality: attributes and relationship to alternative and sustainable agriculture,” American Journal of Alternative Agriculture, vol. 7, no. 1-2, pp. 5–11, 1992.
- J. W. Doran and T. B. Parkin, “Defining and assessing soil quality,” in Defining Soil Quality for a Sustainable Environment, J. W. Doran, D. C. Coleman, D. F. Bezdicek, and B. A. Stewart, Eds., pp. 3–21, Soil Science Society of America, Madison, Wis, USA, 1994.
- J. L. Smith, J. L. Halvorson, and R. I. Papendick, “Using multiple-variable indicator kriging for evaluating soil quality,” Soil Science Society of America Journal, vol. 57, no. 3, pp. 743–749, 1993.
- F. Han, J. Zheng, W. Hu, F. Du, and X. C. Zhang, “Spatial variability and distribution of soil nutrients in a catchment of the Loess Plateau in China,” Acta Agriculturae Scandinavica Section B, vol. 60, no. 1, pp. 48–56, 2010.
- H. J. Chu, Y. P. Lin, C. S. Jang, and T. K. Chang, “Delineating the hazard zone of multiple soil pollutants by multivariate indicator kriging and conditioned Latin hypercube sampling,” Geoderma, vol. 158, no. 3-4, pp. 242–251, 2010.
- J. J. Lee, C. S. Jang, S. W. Wang, and C. W. Liu, “Evaluation of potential health risk of arsenic-affected groundwater using indicator kriging and dose response model,” Science of the Total Environment, vol. 384, no. 1–3, pp. 151–162, 2007.
- Y. Qi, J. L. Darilek, B. Huang, Y. Zhao, W. Sun, and Z. Gu, “Evaluating soil quality indices in an agricultural region of Jiangsu Province, China,” Geoderma, vol. 149, no. 3-4, pp. 325–334, 2009.
- N. Diodato and M. Ceccarelli, “Multivariate indicator Kriging approach using a GIS to classify soil degradation for Mediterranean agricultural lands,” Ecological Indicators, vol. 4, no. 3, pp. 177–187, 2004.
- C. S. Jang, S. K. Chen, and L. Ching-Chieh, “Using multiple-variable indicator kriging to assess groundwater quality for irrigation in the aquifers of the Choushui River alluvial fan,” Hydrological Processes, vol. 22, no. 22, pp. 4477–4489, 2008.
- J. P. E. Anderson and K. H. Domsch, “A physiological method for the quantitative measurement of microbial biomass in soils,” Soil Biology and Biochemistry, vol. 10, no. 3, pp. 215–221, 1978.
- H. Bolton Jr., L. F. Elliott, R. I. Papendick, and D. F. Bezdicek, “Soil microbial biomass and selected soil enzyme activities: effect of fertilization and cropping practices,” Soil Biology and Biochemistry, vol. 17, no. 3, pp. 297–302, 1985.
- J. L. Smith, “Cycling of nitrogen through microbial activity,” in Effects on Soil Quality, J.L. Hatfield and B.A. Stewart, Eds., Soil Biology, pp. 91–120, CRC Press, Boca Raton, Fla, USA, 1994.
- J. D. Snyder and J. A. Trofymow, “A rapid accurate wet oxidation diffusion procedure for determining organic and inorganic carbon in plant and soil samples,” Communications in Soil Science & Plant Analysis, vol. 15, no. 5, pp. 587–597, 1984.
- SPSS, SYSTAT 11.0 Statistics, SPSS, Chicago, Ill, USA, 2004.
- J. J. Halvorson, J. L. Smith, and R. I. Papendick, “Issues of scale for evaluating soil quality,” Journal of Soil and Water Conservation, vol. 52, no. 1, pp. 26–30, 1997.
- J. L. Smith and E. A. Paul, “The significance of soil microbial biomass estimations,” in Soil Biochemistry, J. M. Bollag and G. Stotzkey, Eds., vol. 6, pp. 357–396, Marcel Dekker, New York, NY, USA, 1990.
- R. E. Rossi, D. J. Mulla, A. G. Journel, and E. H. Franz, “Geostatistical tools for modeling and interpreting ecological spatial dependence,” Ecological Monographs, vol. 62, no. 2, pp. 277–314, 1992.
- P. Rochette, R. L. Desjardins, and E. Pattey, “Spatial and temporal variability of soil respiration in agricultural fields,” Canadian Journal of Soil Science, vol. 71, no. 2, pp. 189–196, 1991.
- M. Bonmati, B. Ceccanti, and P. Nanniperi, “Spatial variability of phosphatase, urease, protease, organic carbon and total nitrogen in soil,” Soil Biology and Biochemistry, vol. 23, no. 4, pp. 391–396, 1991.
- R. García-Ruiz, V. Ochoa, M. B. Hinojosa, and J. A. Carreira, “Suitability of enzyme activities for the monitoring of soil quality improvement in organic agricultural systems,” Soil Biology and Biochemistry, vol. 40, no. 9, pp. 2137–2145, 2008.
- F. Zvomuya, H. H. Janzen, F. J. Larney, and B. M. Olson, “A long-term field bioassay of soil quality indicators in a semiarid environment,” Soil Science Society of America Journal, vol. 72, no. 3, pp. 683–692, 2008.
- T. M. Nissen and M. M. Wander, “Management and soil-quality effects on fertilizer-use efficiency and leaching,” Soil Science Society of America Journal, vol. 67, no. 5, pp. 1524–1532, 2003.
- D. L. Karlen, E. G. Hurley, S. S. Andrews et al., “Crop rotation effects on soil quality at three northern corn/soybean belt locations,” Agronomy Journal, vol. 98, no. 3, pp. 484–495, 2006.
- T. M. Zobeck, A. D. Halvorson, B. Wienhold, V. Acosta-Martinez, and D. L. Karlen, “Comparison of two soil quality indexes to evaluate cropping systems in northern Colorado,” Journal of Soil and Water Conservation, vol. 63, no. 5, pp. 329–338, 2008.
- B. N. Moebius-Clune, H. M. Van Es, O. J. Idowu et al., “Long-term effects of harvesting maize stover and tillage on soil quality,” Soil Science Society of America Journal, vol. 72, no. 4, pp. 960–969, 2008.
- E. H. Isaaks and R. M. Srivastava, An Introduction to Applied Geostatistics, Oxford University Press, New York, NY, USA, 1989.