Abstract

Little work has been done to characterize and quantify the residue traits affecting decomposition of winter and spring canola (Brassica napus L.) residue in dryland farming systems of the Pacific Northwest United States. Traditional methods of characterizing residue fiber and nutrients are time-consuming and expensive and require large quantities of chemical reagents. The goal of this research was to determine whether near-infrared spectroscopy (NIRS) could accurately predict neutral detergent fiber (NDF), acid detergent fiber (ADF), acid detergent lignin (ADL), carbon (C), and nitrogen (N) of canola stems, litter, and roots and decomposition of canola stems. Canola residue varied in decomposition, fiber, and nutrients by year, location, and type. NIRS predictions were successful for NDF and ADF in 2011 (standard error of prediction ; ) and NDF, ADF, and N in 2012 (; ). Other predictions for residue fiber and nutrient characteristics were considered moderately successful. Prediction of canola residue decomposition with NIRS was useful for screening purposes. Near-infrared spectroscopy shows promise for rapidly and reproducibly predicting some canola residue fiber and nutrient traits and may be useful for estimating residue decomposition potential in dryland conservation cropping systems.

1. Introduction

Production of winter and spring canola (Brassica napus L.) is increasing in the Pacific Northwest (PNW) United States as marketing opportunities and profitability for growers surge upward. Worldwide, canola production has increased dramatically, and demand for canola oil and meal continues to rise [1]. In the United States in 2008-2009, domestic consumption of canola oil was more than 2.5 times domestic production [1]. Cultivars of winter and spring canola perform well in conventional and conservation farming systems, and higher yielding cultivars and herbicide resistant varieties have made these crops more economically feasible. A review by Johnston et al. [2] notes that canola is well-suited in conservation tillage systems, including no-till, and can result in increased wheat yields in rotation compared to continuous wheat crops.

In the dryland farming areas of the inland PNW, average annual rainfall varies from 150 mm to 500 mm. Soils in the eastern portion of the region with higher rainfall and steeper slopes are cropped annually and are susceptible to erosion by water. The traditional cropping system in the western segment is a winter wheat-summer fallow system with conventional, high disturbance tillage, which leaves soils prone to wind erosion. Growers throughout the region are working to implement conservation farming practices to curb erosion and enhance soil quality; however, residue management is of concern to growers who wish to employ management practices such as no-till farming. There is limited information available on the decomposition of oilseed crop residue in both conventional and no-till systems [3, 4], and little is known about the residue characteristics of canola cultivars that are currently grown or in development for production in dryland and irrigated cropping systems of the PNW.

A combination of factors must be considered in determining decomposition of crop residues. Among those characteristics are fiber [5] and nutrient content of the residue [6, 7]. Neutral detergent fiber (NDF) of residue includes hemicellulose, cellulose, and lignin, which are the insoluble cell wall components. Acid detergent fiber (ADF) consists of the cellulose and lignin portions, with hemicellulose removed; and acid detergent lignin (ADL) is the portion remaining after cellulose is removed. Along with fiber components of residue, nutrient content, especially the C/N of residue, is critical in determining potential of residue to decompose under field conditions. Traditional wet chemistry and combustion methods of fiber, carbon, and nitrogen analysis are time-consuming, expensive, and destructive to the sample. Near-infrared spectroscopy (NIRS) is a potential method for estimating decomposability, fiber, and nutrient content of crop residues that is reproducible, inexpensive, rapid, and nondestructive to the sample [8]. Previous research has shown NIRS to be an acceptable method to predict fiber components [911] and nitrogen [12, 13] of crop residue.

The objectives of this research were to develop calibration equations for winter and spring canola residue NDF, ADF, ADL, N, C, C/N, and decomposition and to determine the feasibility of using NIRS to predict those characteristics in canola residue from various cultivars produced in different years in distinct growing locations. Our goal is to develop NIRS as a rapid screening method for predicting residue composition and decomposition potential in dryland oilseed crops. As growers seek diversification in their wheat-based cropping systems and marketing opportunities for oilseed crops advance, information on canola residue decomposition will be useful to growers who wish to incorporate winter and spring canola into their rotations. Analysis of canola residue nutrient content can be used in determining nitrogen fertilizer efficiency. Additionally, canola residue may be managed for the greatest economic success and soil quality benefits in both conventional and conservation farming systems.

2. Materials and Methods

2.1. Residue Collection

Canola residue was collected after harvest from University of Idaho Brassica Breeding and Research Program canola variety trials in 2011 and 2012. Winter canola residue was collected from trials at Odessa, WA, Moscow, ID, and Genesee, ID, and spring canola residue was sampled from trials located at Davenport, WA, Colfax, WA, and Moscow, ID [14]. Shoot residue from seven cultivars and four replications was sampled at each site, and root residue was sampled to a depth of 30 cm from select locations. Canola litter (leaves, pods, and small stems) was collected from each site. Of the seven spring canola cultivars sampled in 2011, only three were included in the 2012 trials, and so four new spring canola cultivars were substituted in the 2012 sampling. The sampling locations were selected to represent varying amounts of average annual precipitation in eastern Washington and northern Idaho. Winter and spring canola was produced under rainfed conditions at all sites except the irrigated winter canola site at Odessa, WA. Each cultivar was grown under the same management (tillage, planting, seeding rate, fertilization, and weed control) as the grower/cooperator. The trial design at each location was a randomized complete block with four replications. The experimental design was a completely randomized split split-plot, where year was the main plot and location and cultivar were the subplots. Stems and roots were separated, and roots were washed with water and dried. For fiber, nutrient, and NIRS analyses, residue samples were ground to pass a one mm sieve.

2.2. Laboratory Reference Data

Ground canola residue was enclosed in filter bags (ANKOM Technology Corp., Fairport, NY) for determination of neutral detergent fiber (NDF), acid detergent fiber (ADF), and acid detergent lignin (ADL) using a modification of the Van Soest et al. [15] procedure. Neutral detergent fiber and acid detergent fiber were determined sequentially following processing with an ANKOM 200 Fiber Analyzer (ANKOM Technology Corp., Fairport, NY). Acid detergent lignin was determined following digestion in 72% H2SO4. Ground residue samples were analyzed by dry combustion using a LECO TruSpec Analyzer (St. Joseph, MI) to determine total carbon (C) and nitrogen (N) and calculate C/N ratios. For consistency, we chose to analyze litter (leaves, pods, and small stems), stem, and root residues separately, and only stem residue was used in the decomposition study.

Decomposition of canola stem residue was determined in a laboratory incubation study using a procedure similar to Summerell and Burgess [16]. Canola stem residue was cut into 5 cm lengths and oven-dried. A layer of Ritzville silt loam soil (oven-dried and brought to field capacity) was placed in a one-pint glass canning jar, followed by approximately 2.0 g of oven-dry canola stem residue, and covered with a second layer of soil to simulate buried residue. Another 2.0 g of stem residue was placed on the soil surface to simulate surface residue. Residue was dried and weighed prior to placing in soil. Jars were sealed with parafilm to allow air exchange and incubated at 22°C. Residue from 2012 was used for the study and included the seven cultivars each of winter and spring canola from each site. Each cultivar was replicated three times. After an 11-week incubation period, canola residue was removed from soil, washed, dried, and reweighed to determine mass of residue lost during incubation.

2.3. Near-Infrared Spectroscopy

Near-infrared spectroscopy (NIRS) is a rapid, low cost, nondestructive secondary method that was used to predict fiber and nutrient content of canola root and shoot residue. Calibration and prediction equations were developed for canola residue enclosed in stationary metal ring cups (36 mm inside diameter) and scanned with a FOSS XDS Rapid Content Analyzer (Foss NIRSystems, Laurel, MD) using ISIscan software, version 3.10 (Infrasoft International, State College, PA). Samples were scanned twice using the wavelength range 400–2498 nm at 2 nm intervals. The mean of the two scans was used for data analysis. Canola stem, root, and litter data for the two years of the study were analyzed separately and each set was randomly divided into two parts for developing calibration equations (Table 1; 2011, ; 2012, ) and for validation of equations (Table 2; 2011, ; 2012, ). Separate calibration equations and validation of equations were performed for the canola residue decomposition study (Table 3). For both the calibration and validation sets, buried residue samples and surface-placed residue samples. Reference analysis and NIRS analysis were performed on each of the samples in both sets. Random selection of sample sets and calibration and validation statistics were completed using WinISI software, version 4.0 (Infrasoft International, State College, PA). Calibration equations were derived using modified partial least-squares (MPLS) and cross-validation techniques. The scatter correction of standard normal variant and detrend (SNV-D) was applied, along with several different math treatments for derivative order number, gap, and first smoothing. The second smoothing was set at 1 to indicate no second smoothing. Principal component analysis was used to identify and remove spectral outliers. Samples having spectra with Mahalanobis distance values greater than 3.0 were considered outliers and were removed from the file. The appropriate calibration equation for each component was determined by selecting the one with the lowest standard error of cross-validation (SECV) and the 1 − variance ratio () closest to 1 [8]. The ratio of standard deviation (SD)/SECV was calculated [17] and used to determine calibration equations that were acceptable for quantitative prediction of fiber characteristics (Tables 1 and 3). Correlations between NIRS-predicted values and reference values (fiber wet chemistry, LECO TruSpec, or residue decomposition) were determined using Pearson correlation coefficients [18] with JMP software [19].

3. Results and Discussion

3.1. Calibration Results

In the present study, the residue of commonly grown cultivars of winter and spring canola produced at multiple locations were analyzed for neutral detergent fiber (NDF), acid detergent fiber (ADF), acid detergent lignin (ADL), carbon (C), and nitrogen (N) contents and decomposition in lab incubations. Traditional methods to quantify fiber and nutrient content of plant residues are laborious; destroy the plant sample; and are unfriendly to the environment. Near-infrared spectroscopy (NIRS) represents a secondary method to estimate fiber and nutrient content which is of lower cost and nondestructive to the crop sample. Near-infrared spectroscopy is reliable and accurate and requires very little sample preparation [8], while allowing for rapid analysis of multiple properties at one time [20]. Sample characteristics are quantitatively predicted by NIRS using the near-infrared absorbance of a sample to measure organic functional groups. In the present study, stem, root, and litter residue from winter and spring canola cultivars grown at multiple locations over two years were scanned with a FOSS XDS Rapid Content Analyzer, and those results were compared to laboratory reference data from fiber, nutrient, and decomposition experiments to develop calibration equations for prediction of those characteristics. Predictions of canola shoot, root, and litter residue characteristics, as well as canola stem residue decomposition, using NIRS were conducted as part of the present study.

Calibration statistics for prediction equations of winter and spring canola residue traits are listed in Table 1. Separate calibrations were done for each year of the study, with Table 1(a) showing results for residue from 2011 and Table 1(b) for residue from 2012. The math treatments which resulted in most optimal calibration statistics were used for the traits of NDF, ADF, ADL, N, C, and C/N. Using equations developed with specific math treatments improves the predictability for individual data sets [21]. Equations that yielded the values closest to 1 and the lowest SECV values were chosen for predictions. For 2011 winter and spring canola residue, NDF, ADF, and ADL, values were 0.95, 0.97, and 0.88, and SECV values were 2.24, 1.61, and 0.83. Nitrogen, C, and C/N and SECV values for 2011 were 0.97, 0.94, and 0.93 and 0.09, 0.74, and 9.93, respectively. The NDF, ADF, and ADL and SECV values for 2012 canola residue were 0.95, 0.94, and 0.83 and 1.59, 1.40, and 0.83. The 2012 and SECV for N, C, and C/N were 0.97, 0.89, and 0.92 and 0.09, 0.69, and 12.99. Deaville et al. [17] developed guidelines for the predictive ability of calibration equations using SD/SECV values. Under their protocol, calibrations with SD/SECV ratios >3.0 are acceptable for prediction of characteristics using NIRS. Ratios >2.5 but <3.0 indicate calibrations that would be useful for screening, and equations with SD/SECV <2.5 are not useful for making predictions. In this study, all SD/SECV values for 2011 were greater than 3.0 except ADL (2.72; Table 1), indicating that successful calibrations were developed for the traits of NDF, ADF, N, C, and C/N (SD/SECV ranging from 3.71 to 5.55). The same was true for 2012, with SD/SECV values ranging from 3.03 to 5.41 for all traits but ADL (2.41; Table 1). In 2011, the SD/SECV of 2.72 for ADL fell within the range that would be acceptable for screening of traits; however, in 2012 the SD/SECV for ADL was only 2.41, which is outside the range considered acceptable for a screening prediction.

3.2. Validation Results

Reference measurements and validation results are shown in Table 2(a) for 2011 winter and spring canola residue and Table 2(b) for 2012 canola residue samples. In both years, acceptable predictions were attained for each of the residue traits. With 2011 residue, values ranged from 0.83 to 0.96, and values for 2012 residue ranged from 0.75 to 0.92. For all traits in both years, slope values were near 1, with the lowest slope = 0.87 for ADL in 2012. Comparisons of laboratory reference values for each trait as compared to NIRS-predicted values are shown in Figures 16. Because similar results were obtained for the two years of the study, only results from 2011 are shown. In determining whether NIRS is an acceptable method for prediction of residue traits, Mathison et al. [22] used an additional calculation, ratio of prediction to deviation (RPD), which is the measured standard deviation (SD) divided by the standard error of prediction (SEP) of the validation results. When RPD values are >2.5, the predictions are considered to be acceptable, while RPD values of 10 are excellent [22]. In this study, all RPD values indicated acceptable predictions except for ADL in 2012 () and C in both years of the study (2011 RPD = 2.40; 2012 RPD = 2.40). The highest RPD values were NDF (4.33) and ADF (5.20) in 2011. Malley et al. [23] include a separate calculation, RER, in determining suitable predictions using NIRS. The RER is equal to the range of measured reference values divided by the SEP with RER values of >20 considered to be excellent, RER 15–20 successful, RER 10–15 moderately successful, and RER 8–10 moderately useful [23]. In this study, we calculated RER values ranging from 12.73 to 37.21 (data not shown); however, these values gave conflicting results when compared to the guidelines using and RPD values in determining success of predictions. Roggo et al. [24] suggested that RPD values of >3 and RER values >10 may have use for screening characteristics in their study of sugar beet (Beta vulgaris) quality. In our study, only NDF and ADF consistently met that level for RPD, while the RER values for each characteristic would indicate acceptable predictions for screening fiber and nutrient traits of canola residue.

3.3. Prediction of Canola Residue Decomposition

Canola stem residue decomposition was measured in a laboratory incubation study [14], and the percent decomposition was compared to NIRS-predicted values. The calibration statistics and validation results are shown in Table 3. The SD/SECV values were 2.07 for buried residue, as in a conventionally tilled farming system, and 1.89 for residue that was surface-placed, as in no-till farming. Using the guidelines proposed by Deaville et al. [17], these calibration equations would be unacceptable for predicting canola residue decomposition using NIRS. Validation results yielded values of 0.63 for buried residue and 0.67 for surface-placed residue. The RPD values were 1.63 for buried canola stems and 1.71 for surface-placed residue, and both values are considered unsuitable for screening; however, Pearson correlation coefficients were all >0.77 when comparing laboratory canola residue decomposition with NIRS-predicted values (; data not shown). Difficulties in cleaning soil from decomposing residue may have led to variability in the reference measurements [16]. Shepherd et al. [25] found that NIRS has the potential for accurately determining decomposition of organic residues. Additional studies with a larger number of canola residue samples may be needed to develop stronger calibrations in order to use NIRS for predicting canola residue decomposition in soil.

3.4. Correlation between Laboratory Reference Measurements and NIRS-Predicted Measurements

NIRS-predicted values were significantly correlated () with measured reference values for each of the residue fiber and nutrient characteristics tested when calculated across the two years, six study locations, two crop types, and all plant components (Table 4). The six traits of NDF, ADF, ADL, N, C, and C/N were all well-correlated with one another, with correlation coefficients ≥0.92 (level of significance <0.0001). Nitrogen was negatively correlated with each of the other characteristics. Correlations between the fiber fractions (NDF, ADF, and ADL) were higher than those with the nutrient contents (N, C, and C/N). We previously found that NIRS was a successful method for predicting NDF, ADF, and ADL in cereal crop residue, but not for C and N [26]. Velasco and Mollers [12] successfully used NIRS to accurately predict N concentration in Brassica napus L. tissues. Additional work has shown the usefulness of NIRS in predicting features of canola seed [27] and/or meal [28]. Long et al. [29] found that NIRS could measure the oil content of canola seed prior to crushing. Wittkop et al. [30] used NIRS to predict NDF, ADF, and ADL in intact seed of various genotypes of oilseed rape (Brassica napus L.). They found good values for ADL and adequate predictions for NDF and ADF when looking for low fiber genotypes. Others have studied properties of canola residue for decomposition [4, 31], biofuel production [32], and pulp or paper production [33]; however, we are unaware of attempts to develop NIRS calibrations for prediction of canola residue properties for those purposes.

Like other crop plants [34, 35], various components of the canola postharvest residue have differing nutrient and fiber contents [14] and have varying rates of residue decomposition [31]. We found that NIRS-predicted values for fiber and nutrient content of canola stems, litter, and roots were well-correlated when compared to laboratory reference methods (Table 5). The highest correlation coefficient for stem residue was for ADF (0.97) and the lowest was for ADL (0.90). Canola litter, which included leaves, pods, and very small stems, had the highest correlation for ADF as well (0.92) and the lowest for C (0.62). The highest Pearson correlation coefficient for canola roots was with N (0.98), and the lowest was with ADL (0.80).

Crop residues, including canola, vary in fiber and nutrient content among years [3, 36], crop types [3, 37], and growing locations [36, 38]. There was significant correlation between laboratory reference measurements and NIRS predictions for all traits in both years, both crop types, and the six field experiment locations (Table 6; ; level of significance < 0.0001 for all comparisons). The highest correlations were found for samples when separated by year, with the highest correlation for N in 2011 (0.98) and the lowest for ADL in 2011 (0.83). When samples were separated by crop type, the highest correlation was found for N in winter canola (0.95) and the lowest for NDF in spring canola (0.75). For the spring canola locations, the correlation coefficient was highest for C from the Colfax, WA, location (0.96) and lowest for NDF from Moscow, ID, spring canola (0.67). In winter canola, the highest correlation coefficient was for N from Genesee, ID, residue (0.97), and the lowest coefficient was also from the Genesee location for ADL (0.55). Bruun et al. [39] were able to distinguish differences in ash content and degradability among cultivars of wheat straw grown at two sites in Denmark using NIRS predictions.

In previous work, we found that there were no clear differences in fiber or nutrient characteristics among residue samples from various cultivars of winter and spring canola grown in eastern Washington and northern Idaho [14]. In the present study, values were too low to develop calibration equations with adequate sets of validation data for individual winter and spring canola cultivars (data not shown). Redaelli and Berardo [11] found that NIRS could accurately predict fiber components of oat (Avena sativa L.) hulls grown in Italy and that fiber components varied with cultivar, but location played a lesser role in that variation. Oilseed rape genotypes showed variability in ADL with NIRS analysis according to Wittkop et al. [30]. Future work in this area might focus on obtaining a greater number of canola residue samples over more years or experimental locations to determine whether NIRS predictions are feasible for defining differences in residue traits among cultivars, establishing usefulness of NIRS for prediction of fiber and nutrient content of individual plant components, and verifying fertilizer uptake and nutrient efficiency in Brassica crops.

4. Conclusions

NIRS was successful in predicting the characteristics of NDF, ADF, N, and C/N of winter and spring canola residue from two years of field studies. NIRS-predicted values were less consistent for ADL, where the prediction was successful for only one of the two years of the study. In this work, we found that NIRS was useful only for screening canola residue C and decomposition in the laboratory incubation studies, not for prediction of these two traits. NIRS predictions show promise for rapidly discerning differences in certain fiber and nutrient characteristics among canola residue samples while using smaller quantities of chemicals and performing less costly analyses. Ultimately, the rapid screening of crop residues for nutrient and fiber characteristics will provide information to crop producers and scientists for designing conservation cropping systems that make efficient use of fertilizers, increase soil organic matter, prevent soil erosion, and lead to soils with greater productivity.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the funding received from the Pacific Northwest Regional Canola Research Program and the assistance from Jim Davis, Jack Brown, and Megan Wingerson of the University of Idaho Brassica Breeding and Research Program. They are grateful to Jeremy Hansen, Brianna Kemp, and Mackenzie Owen for excellent technical assistance.