Next Article in Journal
Distillery Anaerobic Digestion Residues as Fertilizers for Field Vegetable Crops: Performance and Efficiency in Mid-term Successions
Previous Article in Journal
Allelopathic Activity of Spearmint (Mentha spicata L.) and Peppermint (Mentha × piperita L.) Reduces Yield, Growth, and Photosynthetic Rate in a Succeeding Crop of Maize (Zea mays L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rapid Estimation of Wheat Straw Decomposition Constituents Using Near-Infrared Spectroscopy

by
Nathan S. Nielsen
1,
Tami L. Stubbs
2,
Kimberly A. Garland-Campbell
3 and
Arron H. Carter
1,*
1
Department of Crop and Soil Sciences, Washington State University, Pullman, WA 99164-6420, USA
2
Department of Plant Pathology, Washington State University, Pullman, WA 99164-6430, USA
3
Wheat Health, Genetics and Quality Research Unit, United States Department of Agriculture, Agricultural Research Service, Pullman, WA 99164-6420, USA
*
Author to whom correspondence should be addressed.
Agronomy 2019, 9(8), 462; https://doi.org/10.3390/agronomy9080462
Submission received: 29 June 2019 / Revised: 9 August 2019 / Accepted: 13 August 2019 / Published: 16 August 2019
(This article belongs to the Section Farming Sustainability)

Abstract

:
Adoption of no-till systems in Eastern Washington has been slow due to the difficulty of managing wheat (Triticum aestivum L.) straw residue and the unknown decomposition potential of cultivars. We hypothesize that by analyzing wheat straw with near-infrared spectroscopy (NIRS), calibration models can be developed to accurately predict fiber and chemical constituents of wheat, determining straw decomposition potential. Straw from a panel of 480 soft winter wheat cultivars adapted to the Pacific Northwest are analyzed for neutral detergent fiber (NDF), acid detergent fiber (ADF), acid detergent lignin (ADL), cellulose, hemicellulose, carbon (C), and nitrogen (N). Using modified partial least squares regression and cross validation techniques, specific environment and broad-based NIRS models are calibrated and predictive ability is validated. R2cal values from broad models are better than the specific models, and are 0.85 (NDF), 0.86 (ADF), 0.65 (ADL), 0.88 (cellulose), 0.42 (hemicellulose), 0.67 (C), and 0.73 (N). The corresponding SEP values are 1.68% (NDF), 1.54% (ADF), 0.62% (ADL), 1.14% (cellulose), 1.11% (hemicellulose), 1.23% (C), and 0.06% (N). A Finch × Eltan breeding population is used to further validate models and prediction accuracies for variety selection within a breeding program scenario. The broad NIRS models prove useful for estimating high and low ranges of NDF, ADF, and cellulose in wheat cultivars which translate into characteristics of slow and fast decomposition potential.

1. Introduction

The rainfed regions of Eastern Washington are recognized for their highly successful dryland wheat (Triticum aestivum L.) production systems. However, soil erosion threatens the future productive stability of these cropping systems [1,2]. In regions receiving high precipitation (>300 mm), annual cropping is generally practiced with a three-year rotation of winter wheat/spring cereal/spring legume [3]. These regions are distinguished by slopes as steep as 45% [4] and soils high in clay content [3]. Conventional tillage practices are common and can accelerate erosion issues [5]. The combination of steep slopes, clay soils, and tillage makes these cropping systems prone to serious water erosion [3,6]. Water erosion is especially serious on newly planted winter wheat following a low-residue spring legume [3].
In contrast, the production regions of central Washington receive less than 300 mm of annual precipitation and the grain yield potential is low [3]. Due to the low amount of annual precipitation, a regularly practiced winter wheat/summer fallow cropping system is used, intended to conserve moisture [7]. This rotation involves a 13-month period of barren soil between the harvest and planting of wheat [8] that allows the soil to store water. Tillage is conducted during summer months to break the capillary pores in the soil and it creates a soft dust mulch, which minimizes water losses due to evaporation [9]. Unfortunately, this practice produces fine particulate matter ≤10 µm in diameter (PM10) which are prone to intense wind erosion [10]. This is the dominant rotation in the low precipitation region because no other system has proven able to compete agronomically or economically [11].
Leaving straw residue on top of the soil through reduced tillage or through planting directly into the stubble of the previous crop (colloquially termed ‘no-till’) is an excellent solution to minimize soil erosion in both high and low rainfall regions of Eastern Washington. Accumulated runoff and soil erosion have been observed to be up to 54 times lower in a no-till system than in a conventional tillage system [5]. Long term benefits in soil quality have been observed in no-till systems [4]. Excessive surface residue can, however, create difficulties in planting and can delay seed germination [12]. Seeding into a no-till field is often only possible for farmers who can afford to invest in specialized planting equipment designed for this purpose [13]. Many growers resort to alternative management practices for straw which can be detrimental to air and soil quality, such as burning or baling [14,15]. Managing straw residue in sustainable ways will be important for continued production of wheat in the Pacific Northwest while minimizing soil erosion.
Depending on the cropping system each grower is using, the needs for straw residue management will differ. Within no-till production systems in high rainfall regions, straw will need to decompose rapidly over the winter months to avoid the planting complications that excessive straw residue can create in the spring. On the other hand, straw residue in the low rainfall regions needs to decompose slowly to cover the soil during the entirety of the fallow season. Wheat fiber and chemical constituents have previously been associated with straw decomposition [16,17]. The standard method for determining decomposition potential in winter wheat is based upon a fiber and nutrient analysis [18] that quantifies hemicellulose, cellulose, and acid detergent lignin (ADL) through a wet chemistry process that measures neutral detergent fiber (NDF), acid detergent fiber (ADF), and ADL. Nitrogen (N), carbon (C), and C:N ratios are analyzed in straw by dry combustion and are also used to determine decomposition [19]. High NDF, ADF, ADL, and low N are all indicators of slow decomposition [16]. However, these fiber and chemical analysis methods are slow, expensive, and destructive to the samples.
Near-infrared spectroscopy (NIRS) is potentially an inexpensive, rapid, and nondestructive method for predicting fiber and nutrient content of wheat residue. Prediction models can be made by developing a partial least squares regression to predict phenotypic data based on NIRS- generated spectra. In wheat, NIRS calibration models have been successful for determining protein and gluten content [20,21]. NIRS has also been evaluated to predict the degradability and ash content of wheat straw for use in bioethanol production [22]. Furthermore, wheat and rice straw has been evaluated for use in energy conversion, with NIRS being tested as a predictive measure [23]. The ruminal animal industry has successfully used NIRS to determine chemical composition and nutritive value of forages and straw [24,25]. NIRS also successfully predicted NDF, ADF, and N in canola [26], NDF and ADF in dryland cereal cultivars [14], and ADF and ADL in rice [27]. The use of NIRS is routine for prediction of various traits in many other agricultural crops such as postharvest quality of apples [28], postharvest ripeness and quality of mango [29,30], and carbohydrate content in zucchini fruit [31].
The objective of this experiment was to analyze a diverse winter wheat population with NIRS to develop prediction models that would facilitate the determination of residue decomposition rates. By standardizing NIRS as a primary method of predicting residue decomposition rates, breeders will be able to select for this trait in their breeding programs and farmers can be provided with information regarding the decomposition potential of released winter wheat cultivars. The risk involved with transitioning to no-till systems will be reduced because farmers will be able to make informed decisions when selecting a cultivar that meets their specific needs for production.

2. Materials and Methods

2.1. Sample Collection and Preparation

Two separate populations were used in this study. The first was a panel of 480 advanced soft white winter wheat cultivars from breeding programs in the Pacific Northwest (Washington State University, University of Idaho, Oregon State University, USDA-ARS). This panel was harvested from Pullman, WA in 2016 and from Pullman, Central Ferry, and Mansfield, WA in 2017. Pullman is located at 46.73° N, 117.18° W, on the eastern border of Washington state at 717 m elevation above sea level. The average annual precipitation exceeds 500 mm, soil types range from silt loams to silty clay loams [32], and winter wheat grain yields average 7900 kg/ha [33]. Central Ferry is located approximately 95 km southwest of Pullman at 46.62° N, 117.79° W and 195 m elevation. Central Ferry averages approximately 6700 kg/ha of grain yield under irrigation (600 mm) on silt loam soils [32]. Mansfield, located 320 km northwest of Pullman at 47.81° N, 119.64° W, has an annual precipitation of less than 300 mm and grain yields average around 3550 kg/ha [33]. The soils are classified as ashy fine sandy loam [32] and the elevation is 692 m. At each location, the population and checks were planted in an augmented design with repeating checks every 20 entries and one rep per location. The second population included 167 recombinant inbred lines (RIL) developed through single seed descent after crossing Finch (PI 628640) and Eltan (PI 536994) winter wheat cultivars [34,35,36]. These two cultivars were selected for crossing because they are cultivars released for production in the low rainfall areas of Washington and stark contrasts in decomposition potential exist between the two [16]. This population was harvested from Pullman, Mansfield, and Waterville, WA in 2015, whereas samples were harvested in 2017 from Pullman and Mansfield. Waterville is located approximately 50 km southwest of Mansfield, has mostly silty loam and sandy loam soils [32], and has grain yields comparable to Mansfield. The coordinates of Waterville are 47.65° N, 120.07° W and the elevation is 800 m. At all locations, the Finch/Eltan RIL population was planted as a randomized complete block with two reps per location and a repeating check every 20 entries. In all experiments 0.5 m rows of straw from each plot was cut just above ground level at harvest maturity (Stage 11.4 on Feekes’ scale) [37,38] and placed in brown paper sacks with the grain removed. For consistency, the leaves and the nodes were removed, leaving only the internode portion of the straw. The internode residue was cut into 1–2 cm pieces and then ground to pass through a one mm sieve using a FOSS Cyclotec Sample Mill (FOSS North America, Eden Prairie, MN, USA). Samples were stored in a low humidity chamber prior to analysis to ensure similar moisture content.

2.2. Fiber and Nutrient Analysis

Neutral detergent fiber, ADF, and ADL were determined by analyzing 0.5–0.55 gram of ground winter wheat straw using the vanSoest et al. [18] procedure modified slightly by using an ANKOM automated system with specialized filter bags (ANKOM Technology, Macedon, NY, USA). The NDF procedure removed starches, sugars, free amino acids, and other water-soluble components, leaving hemicellulose, cellulose, and ADL. The ADF procedure removed the hemicelluloses, leaving only the cellulose and ADL. The ADL procedure removed the cellulose from the straw. The NDF and ADF procedures were performed sequentially using an ANKOM 200 Fiber Analyzer (ANKOM Technology, Macedon, NY, USA). The straw samples were then digested in 72% H2SO4 to determine the ADL. After each procedure, the samples were dried overnight in a fume hood and then were oven dried in a hybridization incubator for a minimum of 6 h at 64 °C. The samples were removed from the incubator, placed in desiccator pouches to cool, and individually weighed. The hemicellulose value was calculated as the difference between NDF and ADF whereas the cellulose value was determined by calculating the difference between ADF and ADL. Carbon and N were determined using dry combustion with a LECO TruSpec Analyzer (LECO Corp., St. Joseph, MI, USA) as described by Gazulla et al. [39].

2.3. Near-Infrared Spectroscopy

Finely ground winter wheat residue from the internode samples described above was enclosed in metal ring cups (36 mm inside diameter) and scanned with a FOSS XDS Rapid Content Analyzer (FOSS North America, Eden Prairie, MN, USA) using ISIscan software, version 3.10 (Infrasoft International, State College, PA, USA). At the start of each day, check samples were run to calibrate the machine. Each sample was then scanned twice, with the cup rotated 90° between the first and second scan, using the wavelength range 400–2498 nm at 2 nm intervals. The two resulting spectra from each sample were averaged. Background spectra was recorded with each sample analyzed.

2.4. Statistics

The data from each of the four environments from the first population of 480 cultivars were combined into one set to develop broad range prediction models for each trait, with the exception of one cultivar which was removed from analysis due to lack of emergence under field conditions. After laboratory reference data was combined with the corresponding spectra, each data set was randomly divided into two parts to develop calibration models (n = 1566) and to validate the models (n = 350) using WinISI software, version 4.0 (Infrasoft International, State College, PA, USA). Roughly 18% of the samples were used in the validation models using a random selection procedure within the WinISI software. Principle component analysis was used to eliminate spectral outliers, which were defined as spectra with Mahalanobis distance values greater than 3.0. Modified partial least-squares (MPLS) and fold-cross validation methods using four iterations, as implemented in the WinISI software, were used to develop prediction models [40]. The modification to PLS scales the reference method data and the reflectance data at each wavelength to ensure a standard deviation (SD) of 1.0 before each PLS term [40]. Standard normal variant and detrend (SNV-D) as a scatter correction, first and second derivative spectra, as well as mathematical treatments for derivative order number, gap, first smoothing, and second smoothing were applied. The derivative order improved spectral resolution by enhancing spectral differences [41,42]. The gap (listed as the second value in the math treatment) was the length in nm between points over which the derivative was calculated and the smoothing treatments were used to reduce random noise in the spectral data [42]. The degree of first smoothing is expressed by the segment of wavelength points used as a moving average to ‘smooth’ the spectral output, and the number of wavelength points used is listed as the third value in the math treatment. For all models, the second smoothing (fourth value in the math treatment) was set at 1 to indicate that no second smoothing was used. Numerous combinations of math treatments were tested to develop prediction models. The best prediction equation was determined by identifying that which displayed the highest 1-variance ratio (1-VR) and lowest standard error of cross validation (SECV). The ratio of the SD of phenotypic data to the SECV was calculated, as done by Deaville et al. [43], to determine if an equation was acceptable for quantitative prediction, screening only, or was not useful. The standard error of the calibration (SEC) and the coefficient of determination (R2) were also computed. The unit of measurement for SEC, SECV, and SEP are all recorded in percentages. The second population of RIL was used to confirm the best prediction equation performance in a breeding population closely related to the first population, and analyses were repeated as described above.

3. Results and Discussion

3.1. Calibration Models

Summary statistics for fiber and chemical constituents from the four locations are presented in Table S1. Calibration models were developed for each trait based on the highest 1-VR and lowest SECV, after all environments were combined into one dataset. The mathematical treatments accounting for derivative, gap, and first smoothing differed for each trait when developing calibration models (Table 1). The 1-VR values were highest for NDF, ADF, and cellulose models and lowest for ADL and hemicellulose models (Table 1). SECV values ranged from 0.05% for N to 1.52% for NDF. The SD/SECV is a guideline for determining whether an equation can be validly used for making predictions, as utilized by Deaville et al. [43]. Calibrations with SD/SECV ratios >3.0 are useful for quantitative prediction, models with ratios >2.5 and <3.0 are acceptable for screening purposes, and any equation with a ratio <2.5 is not considered to be useful. A high SD indicates high variation within a population and is desirable because increased data variation will improve the robustness of a calibration [28] and result in an increased SD/SECV ratio. A low SECV indicates that the NIRS calibration is predicting similarly during each iteration of the cross-validation process [41]. According to these guidelines, our NDF and ADF models are useful for sample screening whereas the cellulose equation is acceptable for quantitative prediction (Table 1).
Initially, individual calibration models were developed for each separate environment (Table S2). The 1-VR for NDF was best in Pullman 2016, best in Central Ferry 2017 for ADF, and moderately high across all environments for cellulose. However, the SD/SECV ratio threshold for acceptable performance was not met for the individual environment calibrations. Due to the decreased number of data points and variability within the populations, the individual environment calibration models were much poorer in comparison to the combined environment models (Table 1; Table S2). A broad-based NIRS calibration has advantages over specific calibrations within environments, with increased reliability and robustness [44].

3.2. Validation of Models

The best combined environment calibration models were used to predict a randomly selected validation set (n = 350) that was removed from the population prior to equation development (Table 2). Cellulose was predicted with the highest accuracy and the predictions for ADF and NDF were also acceptable. Nitrogen was predicted with reasonable accuracy, although the SD/SECV ratio of the N equation was below 2.5 which indicates that the usefulness of the equation may be questionable. The R2pred values for ADL and hemicellulose were lower than all other traits included in this study. Additionally, the slope of each trait neared 1 (Table 2; Figure 1).
In previous literature concerning NIRS predictions of fiber characteristics in other crops, there are findings that are both similar and dissimilar to ours. Excellent NIRS prediction models were developed by Jin and Chen [45] for cellulose of rice straw and NDF and ADF of barley straw were successfully predicted by Mathison et al. [22]. Campo et al. [46] reported high R2pred values for NDF (0.91) and ADF (0.91) in maize whereas Wittkop et al. [47] reported an R2pred value for NDF of oilseed rape as 0.62. Several other studies that used NIRS to predict ADL in numerous grain crops reported R2pred values higher than ours [27,48]. Stubbs and Kennedy [26] developed excellent calibration models for NDF, ADF, ADL, C, and N in canola and all traits, with the exception of ADL, which had models that surpassed the SD/SECV threshold of 3.0.
Correlations between reference data and NIRS predicted values for all traits were statistically significant (p < 0.05). Laboratory reference values and predicted values had high Pearson correlation coefficients (above 0.90) for NDF, ADF, and cellulose (Table 3). Reference and NIRS predicted values for ADL, C, and N were also moderately high, ranging from 0.81 to 0.86. The Pearson correlation coefficient was lowest for hemicellulose.
Other guidelines for acceptability of NIRS predictions have been proposed besides the SD/SECV ratio. Williams [49] recommended the residual prediction deviation (RPD) guideline, which is calculated as the ratio of the standard deviation of laboratory reference values and the standard error of prediction (SEP) of the validation set. An RPD value greater than or equal to 2.5 was considered acceptable for NIRS predictions [24], whereas an RPD of 10 or greater was considered excellent [24,49]. In our study, the RPD of NDF (2.61), ADF (2.66), and cellulose (2.90) predictions would be considered successful. The range error ratio (RER), or the ratio of the range of reference data to the SEP, was another NIRS guideline originally developed by Starr et al. [50]. Malley et al. [51] developed guidelines to use the RER alongside the R2pred and RPD for additional insight into the acceptability of NIRS predictions on environmental samples. They proposed that excellent predictions have R2pred > 0.95, RPD > 4, and RER > 20, successful predictions have R2pred between 0.90 and 0.95, RPD = 3–4, and RER = 15–20, and moderately successful predictions have R2pred between 0.8 and 0.9, RPD = 2.25–3, and RER = 10–15. Moderately useful predictions have R2pred from 0.7 to 0.8, RPD = 1.75–2.25, and RER = 8–10. Our predictions of NDF (R2pred = 0.85, RPD = 2.61, RER = 12.24), ADF (R2pred = 0.86, RPD = 2.66, RER = 12.73), and cellulose (R2pred = 0.88, RPD = 2.90, RER = 13.97) would be considered moderately successful whereas the prediction for N (R2pred = 0.73, RPD = 1.83, RER = 9.33) would be considered only moderately useful. Model indices for C, hemicellulose, and ADL are not reported as they did not meet the 2.5 RPD threshold for acceptable NIRS predictions.
Calibration models developed for individual environments were also used to predict validation sets that had been removed from each population before the models were created (Table S3).
Models developed for Central Ferry 2017 best predicted NDF, ADF, cellulose, and C whereas predictions for N, ADL, and hemicellulose were less accurate. For Pullman 2016, the best prediction accuracies were again observed for NDF, ADF, cellulose, and C. Hemicellulose and ADL were predicted with less accuracy, whereas the R2pred value for N was the lowest for all traits in this environment. In Pullman 2017, N had a higher R2pred than all other traits besides cellulose. This is contrary to what was predicted in other environments as prediction accuracies for N were generally lower than most traits. The predictions in Mansfield 2017 were similar to those of Central Ferry 2017 and Pullman 2016, with NDF, ADF, and cellulose being predicted with the highest accuracies. When compared to the R2pred values from the broad calibration model predictions, each trait was predicted with less accuracy in every environment. The R2pred values for C in Central Ferry 2017 and Pullman 2016 were similar to the R2pred for C from the broad model, but were still slightly lower. These results demonstrated that a narrow, specific calibration will generally predict with lower accuracy than a broad, robust calibration.

3.3. Predictions within Breeding Populations

Our calibration models were developed with the purpose of predicting decomposition constituents of winter wheat cultivars from a broad range of environments and germplasm. For this reason, the models were further validated by predicting NDF, ADF, ADL, cellulose, hemicellulose, C, and N of straw samples in a breeding population derived from a cross of Finch and Eltan, two winter wheat cultivars from the Pacific Northwest that were also included in the first population. Across all environments, the highest R2pred values were generally found in NDF, ADF, and cellulose (Table 4). Hemicellulose was consistently difficult to predict and NIRS was generally unsuccessful in predicting C and N. Based upon R2pred values, our equation best predicted NDF in Waterville 2015, ADL, hemicellulose, and C in Pullman 2015, and ADF, cellulose and N in Pullman 2017. The lowest R2pred values for each trait were observed in Mansfield 2017, with the exception of N which was lowest in Waterville 2015. Like the R2pred values, Pearson correlation coefficients were generally lower in this population than they were in the first validation set (Table 5). Even though these values are lower than in the original validation sets, they are expected since we are moving from a diversity panel to a breeding population. Within a breeding context, these values are still acceptable for making selections within segregating populations and identifying lines with either fast or slow decomposition potential.
In order to derive accurate predictions, the spectral data of a validation population must be relatively similar to the spectra of the calibration population. As a general rule, spectra with a Mahalanobis distance above 3.0 are considered outliers and may be difficult for an equation to predict [41]. The very low prediction accuracy of the Finch × Eltan population in Mansfield 2017 for all traits (R2pred range = 0.00–0.44) was surprising due to similarities between the spectral data of this population and the population used for the calibration equation. Approximately 85% of the spectra from Mansfield 2017 samples had Mahalanobis distances below 3.0 when compared with the calibration set. The Mansfield 2017 location did have some variability across the field in early spring growth due to snowdrifts that laid in the field, resulting in some lines being under snow cover for a longer period, and could have affected the straw composition. The SEP was quite high for most traits, indicating that there may be a lack of agreement between spectral information and the reference method [39]. However, it is still possible to get accurate predictions for samples that have a Mahalanobis distance above 3.0 [41]. The spectral data from Waterville 2015 was quite dissimilar to that of the calibration population but was predicted with the best accuracy for NDF (Table 4). Every spectra from samples within this environment had a Mahalanobis distance above 3.0. Nevertheless, to improve the predictive ability and increase the reliability of our models on breeding populations, it is likely that the number of environments included in the calibration will need to be increased.

4. Conclusions

The predictive ability of a NIRS equation is more reliable and accurate when a broad and robust range of data is used in the calibration. Our broad NIRS models were successful in predicting NDF, ADF, and cellulose of the first validation set with high accuracy whereas hemicellulose and ADL were predicted with lower accuracy. The overall predictive ability of NIRS decreased when used to predict the same traits in the Finch × Eltan breeding population but was still moderately high for NDF, ADF, and cellulose and suitable for our purposes of estimating decomposition potential. The C and N prediction accuracies were too low across every environment to be trusted. Using NIRS for screening, rather than prediction, will identify whether a sample fits into a range of high or low NDF, ADF, and cellulose values, which will be sufficient for recommendations and breeding purposes. Grouping NDF, ADF, and cellulose values into a “high” or “low” category, accompanied with C and N values obtained through TruSpec analysis, will provide an estimate of fast or slow decomposition for individual cultivars. Additional samples from a wider range of environments would be beneficial for increasing the predictive accuracy and reliability of the NIRS models and would ultimately assist the transition to conservation farming.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4395/9/8/462/s1, Table S1: Summary statistics for neutral detergent fiber (NDF), acid detergent fiber (ADF), acid detergent lignin (ADL), cellulose (CELL), hemicellulose (HEMI), carbon (C), and nitrogen (N) in Pullman 2016, Central Ferry 2017, Pullman 2017, and Mansfield 2017, Table S2: Near-infrared spectroscopy (NIRS) calibration statistics of the specific equations developed from the population of 480 in individual environments. Traits include neutral detergent fiber (NDF), acid detergent fiber (ADF), acid detergent lignin (ADL), cellulose (CELL), hemicellulose (HEMI), carbon (C), and nitrogen (N). Table includes trait, math treatment, number of samples (n), mean, standard deviation (SD), standard error of calibration (SEC), coefficient of determination (R2), standard error of cross-validation (SECV), 1 minus the variance ratio (1-VR), standard deviation to standard error of cross-validation ratio (SD/SECV), Table S3. Near-infrared spectroscopy (NIRS) validation statistics displaying laboratory reference measurements vs. NIRS predicted results of the specific equations developed from the population of 480 in individual environments. Measurements for neutral detergent fiber (NDF), acid detergent fiber (ADF), acid detergent lignin (ADL), cellulose (CELL), hemicellulose (HEMI), C, and N. Includes trait, number of samples (n), range of reference data, lab mean, lab standard deviation (Lab SD), NIRS predicted mean, bias, coefficient of determination (R2), standard error of prediction (SEP), and slope.

Author Contributions

Conceptualization, A.H.C.; methodology, N.S.N., T.L.S.; formal analysis, N.S.N., K.A.G.-C.; writing—original draft preparation, N.S.N.; writing—review and editing, A.H.C., K.A.G.-C., T.L.S.; supervision, A.H.C.; funding acquisition, A.H.C., N.S.N.

Funding

This material is based upon work that is supported by the National Institute of Food and Agriculture, U.S. Department of Agriculture, under award number 2017-38640-26913 through the Western Sustainable Agriculture Research and Education program under subaward number 201207-595. USDA is an equal opportunity employer and service provider. This research was also funded by Biologically Intensive Agriculture and Organic Farming (BIOAg) grant program from the Center for Sustaining Agriculture and Natural Resources, the National Institute of Food and Agriculture, U.S. Department of Agriculture, under award number 2016-68004-24770 and hatch project 1014919.

Acknowledgments

We thank the technical assistance of Gary Shelton, Kyall Hagemeyer, and Adrienne Burke for maintenance of field locations, assistance with sample preparation, and laboratory technical assistance.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Schillinger, W.F.; Papendick, R.I. Then and now: 125 years of dryland wheat farming in the Inland Pacific Northwest. Agron. J. 2008, 100, S166–S182. [Google Scholar] [CrossRef]
  2. Greer, R.C.; Wu, J.Q.; Singh, P.; McCool, D.K. WEPP simulation of observed winter runoff and erosion in the U.S. Pacific Northwest. Vadose Zone J. 2006, 5, 261–272. [Google Scholar] [CrossRef]
  3. Schillinger, W.F.; Papendick, R.I.; McCool, D.K. Soil and water challenges for Pacific Northwest agriculture. In Soil and Water Conservation Advances in the United States, SSSA Special Publication 60; Zobeck, T.M., Schillinger, W.F., Eds.; SSSA: Madison, WI, USA, 2010; pp. 47–79. [Google Scholar]
  4. Kennedy, A.C.; Schillinger, W.F. Soil quality and water intake in traditional-till vs. no-till paired farms in Washington’s Palouse region. Soil Sci. Soc. Am. J. 2006, 70, 940–949. [Google Scholar] [CrossRef]
  5. Williams, J.D.; Wuest, S.B.; Long, D.S. Soil and water conservation in the Pacific Northwest through no-tillage and intensified crop rotations. J. Soil Water Conserv. 2014, 69, 495–504. [Google Scholar] [CrossRef] [Green Version]
  6. Frazier, B.E.; McCool, D.K.; Engle, C.F. Soil erosion in the Palouse: An aerial perspective. J. Soil Water Conserv. 1983, 38, 70–74. [Google Scholar] [CrossRef]
  7. Schillinger, W.F.; Young, D.L. Best management practices for summer fallow in the world’s driest rainfed wheat region. Soil Sci. Soc. Am. J. 2014, 78, 1707–1715. [Google Scholar] [CrossRef]
  8. Feng, G.; Sharratt, B.; Young, F. Soil properties governing soil erosion affected by cropping systems in the U.S. Pacific Northwest. Soil Till. Res. 2011, 111, 168–174. [Google Scholar] [CrossRef]
  9. Papendick, R.I.; Lindstrom, M.J.; Cochran, V.L. Soil mulch effects on seedbed temperature and water during fallow in Eastern Washington. Soil Sci. Soc. Am. Proc. 1973, 37, 307–314. [Google Scholar] [CrossRef]
  10. Sharratt, B.; Feng, G.; Wendling, L. Loss of soil and PM10 from agricultural fields associated with high winds on the Columbia Plateau. Earth Surf. Process. Landf. 2007, 32, 621–630. [Google Scholar] [CrossRef]
  11. Schillinger, W.F.; Papendick, R.I.; Guy, S.O.; Rasmussen, P.E.; van Kessel, C. Dryland Cropping in the Western United States. In Dryland Agriculture, Agronomy Monograph No. 23, 2nd ed.; Peterson, G.A., Unger, P.W., Payne, W.A., Eds.; American Society of Agronomy, Crop Science Society of America, Soil Science Society of America: Madison, WI, USA, 2006; pp. 365–393. [Google Scholar]
  12. Sullivan, P. Conservation tillage. ATTRA Sustain. Agric. 2003, 105, 1–5. [Google Scholar]
  13. Kok, H.; Papendick, R.I.; Saxton, K.E. STEEP: Impact of long-term conservation farming research and education in Pacific Northwest wheatlands. J. Soil Water Conserv. 2009, 64, 253–264. [Google Scholar] [CrossRef] [Green Version]
  14. Stubbs, T.L.; Kennedy, A.C.; Fortuna, A. Using NIRS to predict fiber and nutrient content of dryland cereal cultivars. J. Agric. Food Chem. 2010, 58, 398–403. [Google Scholar] [CrossRef] [PubMed]
  15. Lafond, G.P.; Stumborg, M.; Lemke, R.; May, W.E.; Holzapfel, C.B.; Campbell, C.A. Quantifying straw removal through baling and measuring the long-term impact on soil quality and wheat production. Agron. J. 2009, 101, 529–537. [Google Scholar] [CrossRef]
  16. Stubbs, T.L.; Kennedy, A.C.; Reisenauer, P.E.; Burns, J.W. Chemical composition of residue from cereal crops and cultivars in dryland ecosystems. Agron. J. 2009, 101, 538–545. [Google Scholar] [CrossRef]
  17. Douglas, C.L., Jr.; Allmaras, R.R.; Rasmussen, P.E.; Ramig, R.E.; Roager, N.C., Jr. Wheat straw composition and placement effects on decomposition in dryland agriculture of the Pacific Northwest. Soil Sci. Soc. Am. J. 1980, 44, 833–837. [Google Scholar] [CrossRef]
  18. van Soest, P.J.; Robertson, J.B.; Lewis, B.A. Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J. Dairy Sci. 1991, 74, 3583–3597. [Google Scholar] [CrossRef]
  19. Johnson, J.M.-F.; Barbour, N.W.; Meyers, S.L. Chemical composition of crop biomass impacts its decomposition. Soil Sci. Soc. Am. J. 2007, 71, 155–162. [Google Scholar] [CrossRef]
  20. Başlar, M.; Ertugay, M.F. Determination of protein and gluten quality-related paramenters of wheat flour using near-infrared reflectance spectroscopy (NIRS). Turk. J. Agric. For. 2009, 35, 139–144. [Google Scholar] [CrossRef]
  21. Albanell, E.; Miñarro, B.; Carrasco, N. Detection of low-level gluten content in flour and batter by near infrared reflectance spectroscopy (NIRS). J. Cereal Sci. 2012, 56, 490–495. [Google Scholar] [CrossRef]
  22. Bruun, A.; Jensen, J.W.; Magid, J.; Lindedam, J.; Engelsen, S.B. Prediction of the degradability and ash content of wheat straw from different cultivars using near infrared spectroscopy. Ind. Crop. Prod. 2010, 31, 321–326. [Google Scholar] [CrossRef]
  23. Huang, C.; Han, L.; Yang, Z.; Liu, X. Ultimate analysis and heating vlaue prediction of straw by near infrared spectroscopy. Waste Manag. 2008, 29, 1793–1797. [Google Scholar] [CrossRef] [PubMed]
  24. Mathison, G.W.; Hsu, H.; Soofi-Siawash, R.; Recino-Diaz, G.; Okine, E.K.; Helm, J.; Juskiw, P. Prediction of composition and ruminal degradability characteristics of barley straw by near infrared reflectance spectroscopy. Can. J. Anim. Sci. 1999, 79, 519–523. [Google Scholar] [CrossRef]
  25. Redshaw, E.S.; Mathison, G.W.; Milligan, L.P.; Weisenburger, R.D. Near infrared reflectance spectroscopy for predicting forage composition and voluntary consumption and digestibility in cattle and sheep. Can. J. Anim. Sci. 1986, 66, 103–115. [Google Scholar] [CrossRef]
  26. Stubbs, T.L.; Kennedy, A.C. Prediction of canola residue characteristic using near-infrared spectroscopy. Int. J. Agron. 2017, 2017, 1–9. [Google Scholar] [CrossRef]
  27. Kong, X.; Xie, J.; Wu, X.; Huang, Y.; Bao, J. Rapid prediction of acid detergent fiber, neutral detergent fiber, and acid detergent lignin of rice materials by near-infrared spectroscopy. J. Agric. Food Chem. 2005, 53, 2843–2848. [Google Scholar] [CrossRef] [PubMed]
  28. Bobelyn, E.; Serban, A.; Nicu, M.; Lammertyn, J.; Nicolai, B.M.; Saeys, W. Postharvest quality of apple predicted by NIR-spectroscopy: Study of the effect of biological variability on spectra and model performance. Postharvest Biol. Tec. 2010, 55, 133–143. [Google Scholar] [CrossRef]
  29. Rungpichayapichet, P.; Mahayothee, B.; Nagle, M.; Khuwijitjaru, P.; Müller, J. Robust NIRS models for non-destructive prediction of postharvest fruit ripeness and quality in mango. Postharvest Biol. Tec. 2016, 111, 31–40. [Google Scholar] [CrossRef]
  30. Nordey, T.; Joas, J.; Davrieux, F.; Chillet, M.; Léchaudel, M. Robust NIRS models for non-destructive prediction of mango internal quality. Sci. Hortic. 2017, 216, 51–57. [Google Scholar] [CrossRef]
  31. Pomares-Viciana, T.; Martínez-Valdivieso, D.; Font, R.; Gómez, P.; del Río-Celestino, M. Characterisation and prediction of carbohydrate content in zucchini fruit using near infrared spectroscopy. J. Sci. Food Agric. 2018, 98, 1703–1711. [Google Scholar] [CrossRef]
  32. Web Soil Survey. Available online: https://websoilsurvey.sc.egov.usda.gov/App/HomePage.htm (accessed on 16 January 2019).
  33. Wheat & Small Grains. Available online: http://smallgrains.wsu.edu/variety/variety-2018-data/ (accessed on 27 February 2019).
  34. Peterson, C.J.; Allan, R.E.; Rubenthaler, G.L.; Line, R.F. Registration of ‘Eltan’ wheat. J. Plant Regist. 1991, 31, 1704. [Google Scholar] [CrossRef]
  35. Garland Campbell, K.A.; Allan, R.E.; Anderson, J.; Pritchett, J.A.; Little, L.M.; Morris, C.F.; Line, R.F.; Chen, X.; Walker-Simmons, M.K.; Carter, B.P.; et al. Registration of ‘Finch’ wheat. J. Plant Regist. 2005, 45, 1656–1657. [Google Scholar] [CrossRef]
  36. Balow, K.; Shelton, G.; Burke, A.; Hagemeyer, K.; Klarquist, E.F.; Froese, P.S.; Carle, S.W.; Roa, A.; Nielsen, N.; Carter, A.H.; et al. Registration of the Finch-Eltan winter wheat recombinant inbred mapping population. J. Plant Regist. 2019, 13, 287–293. [Google Scholar]
  37. Large, E.C. Growth stages in cereals illustrations of the Feekes scale. Plant Pathol. 1954, 3, 128–129. [Google Scholar] [CrossRef]
  38. Feekes, W. De Tarwe en haar milieu. Tech. Tarwe Comm. 1941, XVII, 560–561. [Google Scholar]
  39. Gazulla, M.F.; Rodrigo, M.; Orduña, M.; Gómez, C.M. Determination of carbon, hydrogen, nitrogen, and sulfur in geological materials using elemental analysers. Geostand. Geoanal. Res. 2012, 36, 1397–1402. [Google Scholar] [CrossRef]
  40. Shenk, J.S.; Westerhaus, M.O. Population structuring of near infrared spectra and modified partial least squares regression. Crop Sci. 1991, 31, 1548–1555. [Google Scholar] [CrossRef]
  41. Foss North America. ISIscan & WinISI Software Training Class; Foss North America: Eden Prairie, MN, USA, 2017. [Google Scholar]
  42. Li, J.; Li, Q.; Li, D.; Chen, Y.; Wang, X.; Yang, W.; Zhang, L. Use of near-infrared reflectance spectroscopy for the rapid determination of the digestible energy and metabolizable energy content of corn fed to growing pigs. J. Anim. Sci. Biotechnol. 2016, 7, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Deaville, E.R.; Humphries, D.J.; Givens, D.I. Whole crop cereals 2. Prediction of apparent digestibility and energy value from in vitro digestion techniques and near infrared reflectance spectroscopy and of chemical composition by near infrared reflectance spectroscopy. Anim. Feed Sci. Tech. 2009, 149, 114–124. [Google Scholar] [CrossRef]
  44. Filho, P.A.C.; Volery, P. Broad-based versus specific NIRS calibration: Determination of total solids in fresh cheese. Anal. Chim. Acta 2005, 554, 82–88. [Google Scholar] [CrossRef]
  45. Jin, S.; Chen, H. Near-infrared analysis of the chemical composition of rice straw. Ind. Crop. Prod. 2007, 26, 207–211. [Google Scholar] [CrossRef]
  46. Campo, L.; Manteagudo, A.B.; Salleres, B.; Castro, P.; Moreno-Gonzales, J. NIRS determination of non- structural carbohydrates, water soluble carbohydrates and other nutritive quality traits in whole plant maize with wide range variability. Span. J. Agric. Res. 2013, 11, 463–471. [Google Scholar] [CrossRef]
  47. Wittkop, B.; Snowdon, R.J.; Friedt, W. New NIRS Calibrations for fiber fractions reveal broad genetic variation in Brassica napus seed quality. J. Agric. Food Chem. 2012, 60, 2248–2256. [Google Scholar] [CrossRef] [PubMed]
  48. Bruno-Soares, A.M.; Murray, I.; Paterson, R.M.; Abreu, J.M.F. Use of near infrared reflectance spectroscopy (NIRS) for the prediction of the chemical composition and nutritional attributes of green crop cereals. Anim. Feed Sci. Tech. 1998, 75, 15–25. [Google Scholar] [CrossRef]
  49. Williams, P. Variables Affecting Near-Infrared Reflectance Spectroscopic Analysis. In Near-Infrared Technology in the Agricultural and Food Industries; Williams, P., Norris, K., Eds.; American Association of Cereal Chemists, Inc.: St. Paul, MN, USA, 1987; pp. 143–167. [Google Scholar]
  50. Starr, C.; Morgan, A.G.; Smith, D.B. An evaluation of near infra-red reflectance analysis in some plant breeding programmes. J. Agric. Sci. 1981, 97, 107–118. [Google Scholar] [CrossRef]
  51. Malley, D.F.; Ben-Dor, E.; Martin, P.D. Application in analysis of soils. In Near Infrared Spectroscopy in Agriculture, Agronomy Monograph 44; Roberts, C., Workman, J., Reeves, J.B., III, Eds.; American Society of Agronomy, Crop Science Society of America, Soil Science Society of America: Madison, WI, USA, 2004; Volume 44, pp. 729–784. [Google Scholar]
Figure 1. Linear regression relationship (includes R2pred, slope, and SEP) of laboratory reference data and NIRS predicted values from the first validation set for (a) neutral detergent fiber (NDF), (b) acid detergent fiber (ADF), (c) acid detergent lignin (ADL), (d) cellulose, (e) hemicellulose, (f) carbon, and (g) nitrogen using the broad NIRS models developed when combining a population of 480 cultivars in four environments.
Figure 1. Linear regression relationship (includes R2pred, slope, and SEP) of laboratory reference data and NIRS predicted values from the first validation set for (a) neutral detergent fiber (NDF), (b) acid detergent fiber (ADF), (c) acid detergent lignin (ADL), (d) cellulose, (e) hemicellulose, (f) carbon, and (g) nitrogen using the broad NIRS models developed when combining a population of 480 cultivars in four environments.
Agronomy 09 00462 g001aAgronomy 09 00462 g001bAgronomy 09 00462 g001c
Table 1. Near-infrared spectroscopy (NIRS) calibration statistics of the broad models developed when combining a population of 480 cultivars in four environments.
Table 1. Near-infrared spectroscopy (NIRS) calibration statistics of the broad models developed when combining a population of 480 cultivars in four environments.
TraitMath TreatmentnMeanSDSEC (%)R2calSECV (%)1-VRSD/SECV
NDF2,6,4,1151779.484.121.490.871.520.862.71
ADF2,6,4,1151151.793.921.320.891.380.882.84
ADL1,10,10,115046.420.970.550.680.560.671.73
CELL2,6,4,1151345.343.170.970.911.000.903.17
HEMI2,10,10,1151327.631.431.060.451.070.441.34
C3,5,5,1150843.322.191.080.761.120.741.96
N2,4,4,115110.160.100.050.750.050.732.00
SNV and detrend was used for scatter correction. For traits, NDF is neutral detergent fiber, ADF is acid detergent fiber, ADL is acid detergent lignin, CELL is cellulose, HEMI is hemicellulose, C is carbon, N is nitrogen. The four numbers listed for each trait under Math Treatment are sequentially: derivative number, gap, first smoothing, second smoothing. n is number of samples in calibration set; SD is standard deviation; SEC is standard error of calibration; R2cal is coefficient of determination; SECV is standard error of cross validation; 1-VR is 1 minus the variance ratio; SD/SECV is ratio of standard deviation to standard error of cross validation.
Table 2. Near-infrared spectroscopy (NIRS) validation statistics displaying laboratory reference measurements vs. NIRS predicted results for the broad models and validation set.
Table 2. Near-infrared spectroscopy (NIRS) validation statistics displaying laboratory reference measurements vs. NIRS predicted results for the broad models and validation set.
Laboratory MeasurementsValidation Results
TraitnRangeLab MeanLab SDNIRS MeanBiasR2predSEP (%)Slope
NDF35066.64–87.2178.884.3879.07−0.1960.851.681.02
ADF34940.54–60.1451.474.1051.52−0.0530.861.541.03
ADL3503.91–9.216.371.056.370.0050.650.621.08
CELL34935.20–51.1345.093.3145.12−0.0300.881.141.02
HEMI34922.84–31.3527.451.4527.53−0.0740.421.110.94
C34939.39–47.5643.352.1243.41−0.0610.671.230.94
N3470.00–0.560.170.110.18−0.0030.730.060.95
NDF is neutral detergent fiber, ADF is acid detergent fiber, ADL is acid detergent lignin, CELL is cellulose, HEMI is hemicellulose, C is carbon, N is nitrogen. n is number of samples in calibration set; SD is standard deviation; Bias is the difference between the mean of reference data and the mean of NIRS predicted values; R2pred is coefficient of determination; SEP is standard error of prediction.
Table 3. Coefficients of Pearson Correlation between laboratory reference data and NIRS predicted values for NDF, ADF, ADL, cellulose, hemicellulose, C, and N in the first validation set.
Table 3. Coefficients of Pearson Correlation between laboratory reference data and NIRS predicted values for NDF, ADF, ADL, cellulose, hemicellulose, C, and N in the first validation set.
NDFADFADLCELLHEMICN
Correlation0.920.930.810.940.650.820.86
p-value<0.001<0.001<0.001<0.001<0.001<0.001<0.001
NDF is neutral detergent fiber, ADF is acid detergent fiber, ADL is acid detergent lignin, CELL is cellulose, HEMI is hemicellulose, C is carbon, N is nitrogen. α = 0.05 for the probability level.
Table 4. NIRS validation statistics for all fiber and chemical constituents in five environments from the Finch × Eltan Breeding Population.
Table 4. NIRS validation statistics for all fiber and chemical constituents in five environments from the Finch × Eltan Breeding Population.
Lab MeasurementsValidation Results
TraitnLab MeanLab SDNIRS MeanBiasR2predSEP (%)Slope
F × E Pullman 2015
NDF14286.970.8385.381.5870.531.710.68
ADF14259.301.1856.452.8540.562.960.89
ADL14210.990.637.433.5520.523.580.98
CELL14248.340.8048.83−0.4970.500.770.75
HEMI14227.650.7529.16−1.5080.501.600.93
C15446.140.5746.32−0.1780.290.770.35
N1580.1550.0620.206−0.0520.130.080.72
F × E Waterville 2015
NDF14783.951.8982.061.8940.772.150.77
ADF14755.851.7752.303.5500.643.720.84
ADL14710.370.636.713.6570.413.690.91
CELL14745.481.3245.210.2740.640.900.76
HEMI14728.161.0129.41−1.2460.381.480.84
C14546.087.2946.10−0.0150.037.47−1.37
N1470.1450.0570.151−0.006<0.010.07−0.05
F × E Mansfield 2015
NDF17080.041.8776.303.7410.544.040.62
ADF17050.601.7646.943.6550.583.890.66
ADL1708.300.755.612.6860.322.760.86
CELL17042.281.2541.580.7010.551.250.60
HEMI17029.320.7228.590.7300.251.000.55
C13242.771.6443.67−0.9060.031.940.32
N1530.3810.1690.2140.1660.140.231.22
F × E Pullman 2017
NDF34281.663.2180.930.7250.751.771.08
ADF34253.992.7353.270.7230.691.691.10
ADL3427.020.847.09−0.0640.240.740.83
CELL34246.972.2946.460.5140.791.171.06
HEMI34227.661.6527.420.2450.431.270.96
C33645.410.4243.122.2920.172.390.24
N3410.2330.0810.1730.0590.420.090.73
F × E Mansfield 2017
NDF34375.963.9774.621.3450.333.530.91
ADF34347.313.6345.921.3930.403.140.96
ADL3435.961.005.700.2560.140.980.64
CELL34341.352.9340.331.0290.442.430.99
HEMI34328.651.8128.490.1590.101.740.67
C33645.310.5142.602.7050.002.960.001
N3451.0055.8240.1240.8810.065.8728.32
NDF is neutral detergent fiber, ADF is acid detergent fiber, ADL is acid detergent lignin, CELL is cellulose, HEMI is hemicellulose, C is carbon, N is nitrogen. n is number of samples in calibration set; SD is standard deviation; Bias is the difference between the mean of reference data and the mean of NIRS predicted values; R2pred is coefficient of determination; SEP is standard error of prediction.
Table 5. NIRS Coefficients of Pearson Correlation (and P values) of laboratory reference data and NIRS predicted values for NDF, ADF, ADL, cellulose, hemicellulose, C, and N in five environments from the Finch × Eltan breeding population.
Table 5. NIRS Coefficients of Pearson Correlation (and P values) of laboratory reference data and NIRS predicted values for NDF, ADF, ADL, cellulose, hemicellulose, C, and N in five environments from the Finch × Eltan breeding population.
F × E Pullman 2015
TraitNDFADFADLCelluloseHemicellCarbonNitrogen
Correlation0.730.750.720.710.710.540.36
p value<0.001<0.001<0.001<0.001<0.001<0.001<0.001
F × E Waterville 2015
TraitNDFADFADLCelluloseHemicellCarbonNitrogen
Correlation0.880.800.640.800.61−0.17−0.03
p value<0.001<0.001<0.001<0.001<0.0010.040.73
F × E Mansfield 2015
TraitNDFADFADLCelluloseHemicellCarbonNitrogen
Correlation0.730.760.570.740.500.170.38
p value<0.001<0.001<0.001<0.001<0.0010.05<0.001
F × E Pullman 2017
TraitNDFADFADLCelluloseHemicellCarbonNitrogen
Correlation0.870.830.490.890.660.410.65
p value<0.001<0.001<0.001<0.001<0.001<0.001<0.001
F × E Mansfield 2017
TraitNDFADFADLCelluloseHemicellCarbonNitrogen
Correlation0.570.630.370.660.32<0.010.24
p value<0.001<0.001<0.001<0.001<0.0010.94<0.001
NDF is neutral detergent fiber, ADF is acid detergent fiber, ADL is acid detergent lignin, CELL is cellulose, HEMI is hemicellulose, C is carbon, N is nitrogen. α = 0.05 for the probability level.

Share and Cite

MDPI and ACS Style

Nielsen, N.S.; Stubbs, T.L.; Garland-Campbell, K.A.; Carter, A.H. Rapid Estimation of Wheat Straw Decomposition Constituents Using Near-Infrared Spectroscopy. Agronomy 2019, 9, 462. https://doi.org/10.3390/agronomy9080462

AMA Style

Nielsen NS, Stubbs TL, Garland-Campbell KA, Carter AH. Rapid Estimation of Wheat Straw Decomposition Constituents Using Near-Infrared Spectroscopy. Agronomy. 2019; 9(8):462. https://doi.org/10.3390/agronomy9080462

Chicago/Turabian Style

Nielsen, Nathan S., Tami L. Stubbs, Kimberly A. Garland-Campbell, and Arron H. Carter. 2019. "Rapid Estimation of Wheat Straw Decomposition Constituents Using Near-Infrared Spectroscopy" Agronomy 9, no. 8: 462. https://doi.org/10.3390/agronomy9080462

APA Style

Nielsen, N. S., Stubbs, T. L., Garland-Campbell, K. A., & Carter, A. H. (2019). Rapid Estimation of Wheat Straw Decomposition Constituents Using Near-Infrared Spectroscopy. Agronomy, 9(8), 462. https://doi.org/10.3390/agronomy9080462

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

Article Metrics

Back to TopTop