|
|
||||||||
a USDA-ARS Southeast Watershed Research Laboratory, Tifton, GA 31794
b Dep. of Agronomy and Soils, Auburn University, Auburn, AL 36849
c Global Hydrology and Climate Center, Huntsville, AL 35805
d Physics Dep., Auburn University Auburn, AL 36849
* Corresponding author (dgs{at}tifton.usda.gov)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: ATLAS, Airborne Terrestrial Applications Sensor CAI, cellulose absorption index CV, coefficient of variation NIR, near infrared RS, remote sensing SOC, Soil organic carbon TC, total carbon TIR, thermal infrared TOC, total organic carbon VIS, visible
| INTRODUCTION |
|---|
|
|
|---|
Problems with field-scale residue coverage assessments arise because obtaining spatially representative estimates of residue cover in a timely and cost efficient manner is difficult. Cover estimates are increasingly important due to eligibility and compliance with government cost-share programs, such as the Environmental Quality Incentives Program (EQIP). According to the 1985 Food Security Act, lands considered highly erodible must implement an acceptable conservation program to remain eligible for farm benefits. Furthermore, cost-share recipients for reduced tillage systems must maintain a minimum of 30 to 50% crop residue cover to receive program reimbursements. Current line-transect techniques are labor intensive and accuracy is often a function of line length and the number of data points collected. Remote sensing techniques using high spatial and spectral resolution sensors may facilitate field-scale and regional crop residue cover assessment.
Unlike growing vegetation, there is a general lack of information regarding spectral signatures associated with decaying crop residues. However, a fundamental understanding of molecular functional groups in growing and senescent vegetation provides a foundation for addressing residue spectra. Functional groups present in plant material, such as CH3, OH, and H2O, significantly affect spectral response properties via the presence of absorption bands within the 700- to 2600-nm range (Murray and Williams, 1988). During the initial stages of tissue chlorophyll loss, spectral response is greatest from 400 to 800 nm, as senescent plant tissues absorb incoming blue (300400 nm) and red (500600 nm) spectra while reflecting green (400500 nm). Presence of water at this stage masks absorbance features in the near infrared (NIR) associated with lignin and cellulose (Elvidge, 1990). As decay progresses, the relative abundance of lignin and cellulose present is evidenced by broad absorption bands throughout the 400- to 900-nm spectral region (Elvidge, 1990).
Studies conflict regarding the use of remote sensing data to reliably differentiate between residue and soil. Early attempts to differentiate between soil and residue spectra showed differences in spectral reflectance were greatest in the NIR (Gausman et al., 1973; Aase and Tanaka, 1991). These results are in congruence with a similar study conducted by McMurtrey et al. (1993). McMurtrey et al. (1993) developed a vegetation index using spectrophotometer data of five different crop residues and four different soil types. McMurtrey et al. (1993) found reflection at 450, 660, and 830 nm captured most differences between soil and crop residue spectra in a laboratory setting. In another study, Daughtry et al. (1995) utilized reflectance data to distinguish between a variety of crop residues and soils representing 14 suborders. Results showed visible (VIS) and NIR energy could not reliably distinguish soil from residue due to variability in soil properties, water content, and residue age. Nagler et al. (2000) applied the cellulose absorption index (CAI) developed by Daughtry et al. (1996) to differentiate among residue samples arranged on a black 45 by 45 by 2.5 cm plate using a controlled illumination source. As residue decomposed, CAI values decreased. Findings were confounded by water, litter type, and decomposition stage. More recently, Daughtry (2001) used the CAI to differentiate between spectral response patterns of corn (Zea mays L.), soybean (Glycine max L.), and wheat residue and five soil types at different soil water contents. Samples were arranged in 45-cm square trays consisting of soil or residue, mixed scenes of soil and residue were then simulated. Daughtry (2001) was able to discern the relative amount of residue present, via a positive CAI value. Although moist conditions yielded lower CAI values for residue, the CAI for soil remained constant.
Laboratory and field studies have had some success differentiating among residue coverages based on spectral response patterns. Under controlled laboratory conditions, properly calibrated red and NIR spectra may differentiate among degrees of residue cover. Thermal infrared spectra also show promise as a new method for assessing field scale variability in crop residue coverage. However, due to the expense associated with high-resolution TIR imagery, little has been done to investigate TIR as an alternate method of crop residue assessment. Rapid assessment of residue cover is particularly important in the southeastern USA where residue management may significantly impact soil quality and sustainability of these highly weathered soil systems. Thus, the goals of this study are two-fold: (i) evaluate handheld radiometer data and atmospherically corrected airborne imagery as a means to assess ground cover at a large plot scale, and (ii) evaluate high spectral resolution TIR spectra as a method for depicting crop residue coverage at a large plot scale.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Wheat straw residue applications were designed to mimic a conventional double-cropped system (wheatcotton or wheatsoybean] to evaluate residue spectra during periods of minimal warm season crop cover. Plots (15 by 15 m) were established on weed-free, fresh-tilled surfaces in March 2000 and repeated in April 2001. Pre- and postemergent herbicides were used to control weeds and grasses. Residue cover was calculated on a mass basis as a percentage of the amount of residue necessary for complete ground cover. Treatments consisted of five residue cover rates (0, 10, 20, 50, and 80%) arranged in a completely randomized design. A digital camera was used to acquire images of each plot at inception and classified in 2000 to ensure treatment coverages were met. Monthly digital images were taken in 2001 to monitor changes in residue cover over the year, and average estimates of cover per treatment were used in regression analyses
Laboratory
Composite soil samples were collected within each plot at the onset of the study (01 cm) before residue application to determine near-surface soil properties. Soils were air-dried and sieved to pass a 2-mm sieve. Analyses included total C via dry combustion on pulverized samples, citrate-dithionite extractable Fe (Jackson, 1975), and particle-size distribution on the <2-mm fraction (Kilmer and Alexander, 1949). Near surface soil attributes were similar across sites differing primarily by silt and sand content, with Appalachian Plateau soils having greater silt content and lesser amounts of sand (Table 1). Near-surface samples (01 cm) coincident with each remotely sensed data acquisition were also collected for gravimetric water content (
g).
|
Sensors
GER 1500 Spectroradiometer
Reflectance measurements were collected monthly April through June and October through December, on clear days using a hand-held GER 1500 spectroradiometer (GER Corp., New York). The GER 1500 uses a diffraction grating of silicon photo diodes with 512 individual detectors, and collects data between 350 to 1050 nm in 1.5-nm increments. Wavelengths utilized in this study encompassed the 520- to 900-nm spectrum to coincide with the spectral band passes of the ATLAS. Plot data were collected as close to solar noon as possible under clear conditions. Measurements were taken at nadir, within an 8° field of view, from a distance of 2.4 m above ground to approximate a spatial resolution of 0.30 m2. Data collection consisted of five measurements from within each plot. Measurements were converted to percent reflectance based on the reflectance properties of a spectralon reference plate. The spectralon reflectance plate was placed horizontally on the ground outside of each plot. Reference data were collected from a distance of approximately 25 cm from the spectralon plate.
Airborne Terrestrial Applications Sensor (ATLAS)
The ATLAS multispectral scanner acquired data onboard a Lear jet flown at approximately 1400 m. Airborne Terrestrial Applications Sensor collects data in 15 nominal bands ranging from 400 to 12500 nm, with an approximate spatial resolution of 2.5 m at nadir, and a 72° field of view (Birk, 1992) (Table 2). Simultaneous with acquisition the system records 6° of geometric data: latitude, longitude, pitch, roll, altitude, and heading. The onboard radiometric calibration subsystem consisted of three active sources: integrating sphere, hot black body, and cold black body. These are referenced on each revolution of the scan mirror. Data from the active calibration sources within the sensor were used to develop a system transfer function on a per scan line basis. The specific technique used accounts for sensor drift while eliminating high frequency noise. The result of this process converts each of the original airborne measurements recorded in eight bits per pixel to units of Watts cm2 sr1 of known irradiance at the sensor.
|
![]() | [1] |
Airborne Terrestrial Applications Sensor data were acquired for each site close to solar noon, under clear conditions on 2 June 2000 and 30 July 2001. Observations from the Coastal Plain site were limited to the 2000 data acquisition. Pixels lying completely within each plot were extracted, with each plot consisting of 16 pixels. Surface features and atmospheric attributes were assumed to be equal within each plot, thus the distribution of pixel values about mean plot values and percentage of coefficients of variation (CV) were used to assess sensor noise. Based on this analysis, during the 2000 acquisition <5% of the pixels sampled in most bands were greater than two standard deviations from mean plot values. However, in 2001, the percentage of pixels falling greater than two standard deviations from mean plot values was greater than in 2000. Coefficients of variation within plots were generally <10% in 2000, but as high as 36% in 2001. Furthermore, % CV showed that ATLAS bands 1 (450520 nm) and 8 (20802350 nm) exhibited the most variability in 2000 and bands 1 to 3 within the 450- to 630-nm range were most variable in 2001. In 2001, ATLAS Bands 3 (600630 nm) and 10 (82008600 nm) were faulty and excluded from analyses.
Statistical Analysis
Multivariate spectral data were first subject to band selection using a principal components analysis (PCA), since data encompassed multiple bands within the red, NIR, and TIR spectrum. Next, Duncan's least significant difference routine was used to delineate significant spectral differences (
= 0.10) and determine the magnitude of spectral differences between treatments. Based on these results, stepwise linear regression analysis was used to determine the degree of variability in wheat straw residue that could be explained via remotely sensed data.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
g < 2%) being more than twice that observed under wetter conditions (
g > 10%). Residue degradation also impacted spectral response patterns. Throughout the growing season, the total C content (TC) of residue gradually declined from 50 to <25% TC, suggesting decomposition of residue (Fig. 2)
. Residue cover differences were best observed during late-spring and fall, and related to decomposition of residue as treatment differences were best when residue TC
25%.
|
|
Residue Coverage
Differentiation between residue treatments was best observed using a combination of bands in the 600- to 760-nm range (redNIR). These results are consistent with previous studies, indicating red and NIR spectra best differentiate between residue cover differences (Biard and Baret, 1997; Nagler et al., 2000). Results showed that spectroradiometer data, at the spectral and spatial resolution used, could differentiate between plots receiving 20, 50, and 80% residue cover (Fig. 3)
. The relative magnitude of difference between bare soil plots and treatment plots varied during the collection period, with the greatest differences among treatments occurring during the October collection period, when TC content of the residue was <25%.
|
|
Spectral Response Curves
Visible and NIR spectral response patterns were similar to the bare soil line, differing only in magnitude of spectral response. Reflectance in the VIS portion of the spectrum slowly increases to a peak reflectance at approximately 760 nm, and then declines rapidly into the middle infrared regions (Fig. 4)
. Unlike VIS and NIR spectra, TIR spectral response curves differed in slope and magnitude of response when compared with the bare soil line. With <20% residue cover, soil spectral response dominated the shape of the spectral response curve with plots receiving residue distinguishable only by magnitude of spectral response (Fig. 4). However, as the amount of residue cover increased, the slope of the residue line for plots receiving 50 or 80% cover becomes positive in the 8200- to 8600-nm region and levels off in the 8600- to 9200-nm region. Beyond this point, emittance decreases rapidly and residue cover treatments are mostly indistinguishable.
|
g = 9.9%) compared with relatively drier (
g = 1.2%) surface conditions in 2000. Near-surface water absorbs a greater proportion of energy, thereby reducing the magnitude of difference in reflected energy for each treatment (Capehart and Carlson, 1997). At the Coastal Plain site, ATLAS Band 6 (red) best distinguished between 20, 50, and 80% cover, with no treatment differences between 0, 10, and 20% cover.
|
Regression analyses confirmed that a highly significant linear relationship (P < 0.0001) existed between emittance and residue cover at the Appalachian Plateau site. Thermal infrared spectra resulted in r2 peaking at 0.98 and 0.83 in 2000 and 2001, respectively (Table 4). Visible-red spectra were also useful, accounting for as much as 98% of the residue variability in 2000 and 74% of the residue variability in 2001. Under relatively dry conditions in 2000, TIR data explained 95% of the residue variability at the Coastal Plain site while, VIS-red spectra accounted for 7781% of the variability in residue cover (Table 4).
|
| CONCLUSION |
|---|
|
|
|---|
Comparison of ATLAS imagery to handheld spectroradiometer data demonstrate the role atmospherically corrected, high spectral resolution, airborne imagery can play in rapid assessment of residue cover. First, VIS and NIR ATLAS datasets more clearly differentiated among treatments compared with handheld spectroradiometer data. Second, ATLAS results showed that high spectral resolution TIR imagery more accurately delineated treatment differences compared with VIS and NIR spectra. Emitted energy is a function of the emissivity of an object and its temperature; hence emittance is intrinsic to the object of interest, whereas reflected energy is an indirect measure of the state of a given ground feature. Thus, plots receiving greater amounts of residue cover are distinguishable based on the differing heat capacities of organic (residue) and mineral (soil) surfaces. All other conditions being equal, our data suggest TIR is a more stable assessment residue cover compared with VIS and NIR. Another benefit of airborne imagery is a function of spatial representation. ATLAS data captured "whole-plots" compared with handheld datasets that consisted of a limited number of samples within each plot. Thus, the ATLAS sensor likely captured a more representative sample of the variability in residue distribution within the plot.
Overall, residue cover and spectral response exhibited a significant (p
0.05) linear relationship. Data acquired using ATLAS provided the most accurate results throughout the VIS, NIR, and TIR (r2 = 0.770.98). Results were best when surfaces were dry (
g < 2%). Monthly spectroradiometer data were useful in estimating residue cover by combining two bands in the red and NIR region using linear regression (r2 = 0.650.86).
Received for publication May 25, 2004.
| 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] |
||||
![]() |
D.G. Sullivan, D. Lee, J. Beasley, S. Brown, and E.J. Williams Evaluating a crop residue cover index for determining tillage regime in a cotton-corn-peanut rotation Journal of Soil and Water Conservation, January 1, 2008; 63(1): 28 - 36. [Abstract] [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] |
||||
![]() |
D. G. Sullivan, C. C. Truman, H. H. Schomberg, D. M. Endale, and T. C. Strickland Evaluating Techniques for Determining Tillage Regime in the Southeastern Coastal Plain and Piedmont Agron. J., August 3, 2006; 98(5): 1236 - 1246. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Vadose Zone Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||