|
|
||||||||
a USDA-ARS Southeast Watershed Research Lab., P.O. Box 748, Tifton, GA, 31794
b USDA-ARS South J. Phil Campbell Sr. Natural Resource Conservation Center, 1420 Experiment Station Rd., Watkinsville, GA 30677
* Corresponding author (dgs{at}tifton.usda.gov)
Received for publication October 21, 2005.
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: ATLAS, airborne terrestrial applications sensor CAI, cellulose absorption index CRC, crop residue cover CT, conventional tillage NIR, near infrared NT, no tillage RS, remote sensing ST, strip tillage TIR, thermal infrared TM, Thematic Mapper VIS, visible
| INTRODUCTION |
|---|
|
|
|---|
Estimates of residue cover are typically made via roadside surveys (Conservation Technology Information Center, 2004) or in-field line-transect measurements (Shelton et al., 1993). Thoma et al. (2004) compared three methods of residue classification including the roadside survey, in-field line transect, and satellite-derived estimates of residue coverage via Landsat Thematic Mapper (TM) data. In their study, the roadside survey failed to identify residue coverages between 25 and 35% cover nearly 58% of the time compared to in-field line-transect estimates. Limited accuracy of the roadside survey was attributed to difficulties in observing crop residue at oblique viewing angles, and generally resulted in an overestimation of percent residue cover. Thoma et al. (2004) reported that Landsat TM data were more efficient, providing rapid, unbiased estimates of residue cover into two classes (<30% and >30% crop residue cover).
Laboratory and field-scale RS of crop residues have yielded mixed results (McMurtrey et al., 1993; Chen and McKyes, 1993; Sullivan et al., 2004; Daughtry et al., 2005). Unlike growing vegetation, crop residue lacks a unique spectral signature in much of the visible (VIS) and NIR spectrum (McMurtrey et al., 1993; Streck et al., 2002). Instead, crop residues have spectral response features similar to soil spectra, increasing without inflection throughout the VIS and NIR, and differing only in magnitude of spectral response (Baumgardner et al., 1985; Aase and Tanaka, 1991; Daughtry et al., 1995; Sullivan et al., 2004). Difficulties in estimating crop residue cover are a function of soil physical properties, soil water content, crop residue type, crop residue water content, and surrounding green vegetation (Chen and McKyes, 1993; Daughtry et al., 1995; Nagler et al., 2000). In particular, soil background reflectance may be greater or less than crop residue reflectance depending on soil color and water content (Aase and Tanaka, 1991). This manifests a significant challenge in remote residue cover determinations based on differences in the magnitude of spectral response alone.
In much the same way as vegetative indices are designed to reduce soil background effects, researchers have begun investigating tillage indices designed to capture the spectral response of crop residue (McNairn and Protz, 1993; Daughtry et al., 1996; van Deventer et al., 1997; Gowda et al., 2001). Indices capitalize on differences in spectral response between residue and soil spectra within the NIR (Gausman and Allen, 1973; Aase and Tanaka, 1991). van Deventer et al. (1997) evaluated several Landsat 5 TM indices as a means to differentiate between CT and conservation tillage on 27 farms in Ohio. Results indicated the normalized difference tillage index using TM bands 5 (15501750 nm) and 7 (20802350 nm) best discriminated between CT and conservation tillage practices with 89% accuracy. Later, Gowda et al. (2001) applied logistic regression models developed by van Deventer et al. (1997) to Minnesota fields using 1997 Landsat TM imagery. Using logistic regression the percentage of conservation tillage fields classified correctly ranged from 42 to 77%, with indices containing TM band 5 having accuracies between 70 and 77%. Classification errors were attributed to field-scale variability in soil organic C, soil water content, and soil color.
More recently, Daughtry et al. (2005) evaluated RS indices, including the cellulose absorption index (CAI), to more specifically classify tillage practices by the extent of crop residue coverage. The CAI is designed to take advantage of absorption bands centered on 2100 nm, which are highly correlated with the presence of cellulose and lignin in organic materials (Elvidge, 1990; Daughtry et al., 1996). Results from Daughtry et al. (2005) indicated that Landsat TM vegetation and tillage indices were not well correlated with small changes in crop residue coverage. However, the CAI was linearly related to increasing amounts of crop residue coverage having an r2 = 0.88 when the vegetative cover fraction was <0.30. In other studies, the CAI has been shown to be effective even in the presence of little crop residue coverage (Nagler et al., 2003). Earlier techniques used to estimate crop residue cover include fluorescence and the "soil-line" approach (Daughtry et al., 1995; Biard and Baret, 1997). However, data were collected in the laboratory or under artificial field conditions.
Thermal infrared (TIR) spectra also show promise as a new method for assessing field-scale variability in crop residue coverage. In an early study, Aase and Tanaka (1991) used infrared thermometer data to quantify varying degrees of residue cover under wet and dry conditions in the Great Plains. Results showed that under moist conditions, TIR data more accurately quantified residue cover compared to VIS and NIR spectra. Sullivan et al. (2004), evaluated the high spatial and spectral resolution airborne terrestrial applications sensor (ATLAS) to differentiate among wheat residue covers (0, 10, 20, 50, and 80%) in 15 by 15 m plots in Alabama. Results demonstrated that although red and NIR spectra could be used to assess crop residue coverage, TIR data more accurately differentiated between plots receiving 10, 20, 50, and 80% wheat residue cover. Moreover, spectral response curves indicate unique spectral signatures associated with increasing wheat residue cover were present in the 8200 to 9200 nm spectral regions. Authors attribute differences in TIR emittance to contrasting heat capacities of mineral vs. organic materials.
Few studies have evaluated the potential for RS data to depict residue cover in the southeastern USA, where conservation tillage practices are becoming increasingly common. Water quality, conservation effects assessment, and eligibility in federal conservation programs necessitates an accurate and rapid method to measure crop residue distributions. Our study was designed to (i) assess the impact of surface conditions on our ability to remotely discriminate between tillage regimes in two intensively row cropped physiographic regions, (ii) evaluate new RS indices to assess residue cover following cover crop kill and bed preparation in two distinct soils, and iii) compare line-transect crop residue cover estimates with RS residue cover estimates.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The second site was located in the Piedmont region of Georgia, at the USDA, ARS, J. Phil Campbell, Sr., Natural Resource Conservation Center near Watkinsville, GA (33°54' N, 83°24' W). The soil was a Cecil sandy loam (clayey, kaolinitic, thermic Typic Kanhapludult). Treatments consisted of CT and conservation tillage in the form of NT. Plots (10 by 30 m) were arranged in a completely randomized design and replicated three times. A winter rye cover crop was planted following corn harvest on 25 Oct. 2004. Corn stalks were mowed and CT plots were disk harrowed before planting rye. On 14 Mar. 2005, a contact herbicide was used to kill the winter rye before planting corn on 18 Apr. 2005. All plots were planted using a 76-cm row spacing. Corn was planted directly into the killed cover crop, thus rye residues were distributed evenly across each treatment. Conventional tillage plots were mowed and disked to turn rye and prepare beds for planting.
Ground Truth
Ground truth and RS data were collected three times at each site over a 4-wk period. This time frame corresponded to 24 May to 16 June 2004 at the Coastal Plain site, and 19 Apr. to 9 May 2005 at the Piedmont site. Sampling times were chosen to minimize crop canopy interferences based on planting dates and crop growth patterns. Ground truth consisted of digital images, soil water content (05 cm), and soil texture.
Two digital images were taken at nadir from random locations within each plot to quantify the extent of residue cover. Digital images were acquired without a flash, using a 5-mega pixel Olympus C-505 Zoom (Olympus, London, UK). Images were acquired from a height of 1.5 m, centered directly over the row, and represent an area of 1.4 m2 on the ground. Images were classified into four classes: shadow, soil, residue and vegetation using ERDAS Imagine 8.4 (Leica Geosystems, Heerbrugg, Switzerland). Percent residue cover was calculated by dividing pixels classified as residue by the total pixel count in each image (5 million). To determine the validity of the classified images, an accuracy assessment was conducted in ERDAS Imagine using a random selection of four images (conservation tillage plots only) from each RS acquisition date. For each classified image, the assessment randomly chooses 40 points. Each point was then referenced as shadow, soil, crop residue, or vegetation based on a visual interpretation of the unclassified digital image. Results of the accuracy assessment were used to calculate the percentage of each class that was accurately identified. Based on these results, 80% of the classified points were accurately identified as shadow, soil, crop residue, or vegetation at each study site and RS acquisition. Average estimates of residue cover per RS acquisition were used in statistical analyses (Table 1).
|
v, 05 cm) was obtained coincident with each RS acquisition using a Wet Sensor probe (Dynamax, Houston, TX). The Wet Sensor probe uses a measure of the dielectric constant of the soil matrix to estimate volumetric water content (Topp et al., 1980; Whalley, 1993). The general equation can be solved to estimate volumetric water content:
![]() | [1] |

is the square root of the dielectric constant,
v is volumetric soil water content, a0 is the intercept, and a1 is the slope. Using default calibration parameters for a mineral soil, the Wet Sensor has an accuracy of ±3 to 5% volumetric water content. Because the probe was 7.6 cm in length, it was inserted at a 45° angle to ensure only the upper 5 cm of soil water content was measured. Wet Sensor measurements were made at four random locations and composited within each plot. Because soil water content can vary from 0 to 5 cm, precipitation data preceding RS data acquisitions have been provided (Fig. 1
). In addition, composite soil samples were collected within each plot at the onset of the study (020 cm) for soil texture on the <2 mm fraction (Kilmer and Alexander, 1949).
|
The line transect was a flat, plastic measuring stick marked with tape beginning at zero and continuing at 0.63-cm intervals to 0.75 m (120 tick marks). A tick mark was counted each time a piece of residue touched the outside, left edge of the tape (Shelton et al., 1993; Eck et al., 2001). Only crop residues having a width >0.25 cm were counted (Shelton et al., 1993). Percent cover was calculated by dividing the counted number of ticks by total ticks (n = 120) along the transect and multiplying by 100. To evaluate variability in the line transect approach, two transects were established for each image: (i) upper left corner to lower right corner, and (ii) lower left corner to upper right corner.
CropScan Multispectral Radiometer
Reflectance measurements were collected using a hand-held CropScan Multispectral Radiometer (CropScan, Rochester, MN). The CropScan uses narrow band interference filters to select discrete bands in the VIS and NIR regions of the electromagnetic spectrum. Eight bands were measured in this study within the 485 to 1650 nm range (Table 2). The CropScan is equipped with upward and downward looking sensors in each band, and simultaneously acquires irradiance as well as radiance over the target. It is assumed that irradiance over the sensor head is equal to irradiance over the target. Radiance and irradiance were measured in millivolts, adjusted for temperature of the Cropscan, and converted to an energy term. Percent reflectance was determined using the following equation:
![]() | [2] |
|
Because TIR data were acquired over approximately 1 h, it was necessary to adjust all output for changes in ambient air temperature. Ambient air temperatures were recorded using a HOBO Pro Temp/RH Weatherproof Recorder and radiation shield (Onset Computer Corp., Bourne, MA). Temperatures were recorded every 2 min throughout each RS acquisition and used to calibrate TIR data. Since surface temperatures were highly correlated (r = 0.91, P < 0.10) with ambient air temperature, each Ti30 measurement was adjusted using a simple difference approach (Sadler et al., 2002). Thus, each TIR measurement was adjusted by adding or subtracting the change in ambient air temperature from initial conditions.
Statistical Analysis
Using the Statistical Analysis System (SAS Inst., Cary, NC), an analysis of variance (
= 0.10) was used to evaluate differences in tillage regime using line-transect, visible, NIR, or TIR methods of estimation. Visible and NIR indices included the greenness normalized difference vegetation index (GNDVI) (Gitelson et al., 1996), which was calculated as:
![]() | [3] |
![]() | [4] |
Next, linear regression analyses were used to determine the degree of variability between tillage treatments that could be explained via the line transect, VIS/NIR indices or TIR methods. It should be noted that a significant linear relationship between RS data (VIS, NIR, and TIR) and increasing crop residue cover has been established (Biard and Baret, 1997; Nagler et al., 2003; Sullivan et al., 2004). Thus, extreme residue cover conditions (NT or ST vs. CT) were sufficient to establish a relationship between residue cover and RS data. Because tillage regimes differed between sites (NT vs. ST) statistical analyses were run individually for each site. Average crop residue cover estimates for NT and ST treatments during each RS acquisition were used in the regression analyses.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
In the TIR, no treatment differences were observed (P < 0.10). Direct measures of ground temperature showed no significant differences between treatments and suggest that heat capacities of mineral soil and crop residues had not yet been reached. Results contradict previous findings by Sullivan et al. (2004), which showed significant differences between bare soil emittance and varying degrees of wheat residue cover (1080% cover). Conflicting results were likely associated with differences in the time of RS acquisition. Sullivan et al. (2004) collected airborne imagery at 1430 h EST, compared to 1100 h EST in this study. Perhaps the earlier acquisition time in our study may not have been adequate to capture differences in surface emittance associated with contrasting heat capacities of soil and crop residue. Future work is necessary to evaluate the impact of image acquisition time for TIR assessments of cover.
Piedmont
At the Piedmont site, the overall shape of the spectral response curve was similar for NT and CT treatments, steadily increasing from 485 to 1650 nm (Fig. 3
). Differences in spectral response were greatest in the NIR compared to the VIS. As a result of higher clay content at this site, surface soil water contents were relatively higher in the Piedmont compared to the Coastal Plain site, ranging from 8 to 18 cm3 cm3 (Table 3). Thus, compared to the Coastal Plain site, more irradiant energy was absorbed at the soil surface and NT treatments were more reflective compared CT treatments. Data demonstrate the impact that surface soil properties can have on spectral response and our ability to accurately differentiate between conservation tillage and CT systems in two different physiographic regions.
|
|
In the TIR, no significant differences between treatments were observed. As previously mentioned, this may be related to pre-noon RS acquisitions. Additional research is necessary to identify timing for optimal TIR acquisitions.
Cover Estimates
Remotely sensed crop residue cover indices were compared to line-transect estimates of cover to evaluate the utility and accuracy of a rapid, RS residue cover assessment.
Line Transect
Pairs of transects were compared to evaluate variability in the line-transect approach. At both sites, no significant differences were observed between pairs of transects. However, differences in estimated residue cover over time were observed (P < 0.10). Since cover estimates were acquired over a short sampling period, differences in estimated cover were likely due to variability in sample location within a plot. Although treatment differences between NT or ST and CT treatments were significant, variability in cover estimates over time suggest that point-based sampling methodologies may not yield spatially representative estimates.
Using the line-transect technique in linear regression, the line-transect explained >95% of the variability in residue coverage at both sites despite differences in surface conditions at the time of data acquisition (Table 4).
|
![]() | [5] |
![]() | [6] |
![]() | [7;] |
At the Coastal Plain study site, significant differences between ST and CT treatments were observed using all RS indices, except the NDVI (Fig. 4 ). Differences in spectral response between CT and ST treatments were best observed before 25% canopy closure. Once the peanut canopy exceeded 25% cover (NDVI > 0.33), reflectance from the developing peanut canopy masked crop residue reflectance. Furthermore, increasing canopy cover was positively correlated with GNDVI and NDVI indices (r = 0.82, P < 0.10), which limits our ability to accurately assess crop residue cover in the presence of developing vegetation. Daughtry et al. (2005) also reported a linear relationship between vegetative indices and canopy cover, which limited the utility of vegetative indices in crop residue cover determination.
|
Compared with vegetative indices, CRC indices more consistently differentiated between tillage treatments over time. Treatment differences were best observed during the 24 May and 8 June RS acquisitions. Although CRC indices use a portion of the NIR spectrum, the correlation between CRC indices and canopy cover (r < 0.56, P < 0.10) was generally lower compared to vegetation indices. Before canopy closure, CRC index values varied as much as 10 to 38% as a function of soil water content (Table 3) with CRC1 and CRC3 exhibiting the greatest stability (Fig. 4). However, given the low range in soil water content studied here, future research is necessary to determine the effects of changing soil water content on our ability to distinguish between CT and ST using the CRC1 or CRC3.
Threshold values, based on separability of ST and CT plots, were established for each of the three crop residue indices for the 24 May and 8 June 2004 RS acquisition dates. During this time, tillage treatments were separable using a unique CRC threshold value. Specifically, ST treatments exhibited a CRC index value greater than the established threshold of 0.58, 0.38, or 0.43 for the CRC1, CRC2, and CRC3, respectively (Fig. 4).
At the Piedmont study site, significant treatment differences were a function of RS index and surface (residue and soil) conditions at the time of data acquisition. Treatment differences were best observed using the NDVI, CRC1, and CRC3 indices (Fig. 5 ). Keeping in mind that canopy interference was minimal at this site, the NDVI accurately and consistently differentiated between CT and NT treatments, despite differences in soil water content between RS acquisitions. No-tillage treatments typically exhibited a threshold NDVI > 0.16. Moreover, volumetric water content ranged from 8 to 18 cm3 cm3 with NDVI values fluctuating within 2% of the benchmark index calculated using RS data acquired on 19 Apr. 2005. Results suggest that the NDVI, if acquired proximate to planting, is relatively stable under the range in soil water content studied here, and may be used to differentiate between tillage regimes in the Southern Piedmont physiographic region.
|
Crop residue cover index values varied from 10 to 30% of the benchmark index calculated on 19 Apr. 2005 (Fig. 5). The CRC1 provided the most consistent results, exceeding benchmark CRC1 values by 10% under moist soil conditions (Table 3). Using the CRC1, NT treatments were separated using a threshold value <0.65. Significant treatment differences were also observed using CRC3 for both April RS acquisitions, however, the CRC3 was more sensitive to changes in soil water content. For NT treatments the maximum observed CRC3 value ranged from 0.46 under moist soil conditions to 0.53 under dry soil conditions, compared to 0.53 and 0.57 for CT treatments under dry and wet conditions, respectively (Fig. 5). Because the observed CRC3 values for CT and NT treatments overlap as soil conditions change, it would be difficult to differentiate between tillage systems using CRC3 without a priori knowledge of surface soil water content. Other studies confirm, that variability in surface conditions at the time of RS data acquisition significantly impact estimates of vegetative cover (Daughtry et al., 1995; Guerif and Duke, 2000; Nagler et al., 2003).
Linear regression was used to compare the amount of variability between tillage treatments explained using the line transect approach and RS indices (Table 4). Only RS indices that exhibited significant (P < 0.10) differences between CT and ST or NT treatments were used in the analysis. At the Coastal Plain site the GNDVI and CRC indices explained >84% of the variability in crop residue cover. In Piedmont, the NDVI, CRC1, and CRC3 were best during the April RS acquisitions, having coefficients of determination from 0.60 to 0.96 with CRC1 explaining the greatest degree of variability between tillage treatments (Table 4).
Crop residue cover indices, particularly CRC1 and CRC3, performed similarly to the line transect method of crop residue estimation (Table 4). Although, the line transect method was not sensitive to changes in soil water content and canopy cover, the line transect approach may not provide spatially representative estimates of crop residue cover distribution.
| CONCLUSIONS |
|---|
|
|
|---|
Surface soil property variability between sites, was perhaps the greatest single variable affecting the crop residue cover indices. In our study, crop residue reflectance generally ranged from 8 to 40% (4851650 nm), however, the magnitude of response was greater or less than bare soil reflectance as a function of surface soil texture and soil water content. Surface soil texture was the greatest determining factor of crop residue cover index threshold values. Variability in soil water content was secondary to differences in soil texture, however, additional information is necessary to determine the robustness of CRC Index 1 over a wider range of soil water content.
Considering that line-transect estimates are time and resource intensive, results are promising and suggest that threshold RS index values, when used in combination with commonly available soil survey data, may be used to more rapidly differentiate between CT and ST systems compared to the line-transect approach. Remotely-derived crop residue cover maps are necessary to determine eligibility for federal conservation program resources, assess changes in watershed hydrology, and better determine the impact that conservation tillage practices have on soil and water quality.
| NOTES |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D.G. Sullivan, T.C. Strickland, and M.H. Masters Satellite mapping of conservation tillage adoption in the Little River experimental watershed, Georgia Journal of Soil and Water Conservation, May 1, 2008; 63(3): 112 - 119. [Abstract] [PDF] |
||||
![]() |
T. L. Potter, C. C. Truman, T. C. Strickland, D. D. Bosch, and T. M. Webster Herbicide Incorporation by Irrigation and Tillage Impact on Runoff Loss J. Environ. Qual., May 1, 2008; 37(3): 839 - 847. [Abstract] [Full Text] [PDF] |
||||
![]() |
D.G. Sullivan, J.N. Shaw, A. Price, and E. van Santen Spectral Reflectance Properties of Winter Cover Crops in the Southeastern Coastal Plain Agron. J., November 6, 2007; 99(6): 1587 - 1596. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Crop Science | Vadose Zone Journal | |||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||