第一作者:杨天辉(1989-),男,宁夏同心人,助理研究员,硕士,研究方向为草地栽培。E-mali:[email protected]
对2012-2013年黄土高原种植的13个牧草品种、780份干草样品的营养成分建立了近红外光谱(near infrared reflectance spectroscopy,NIRS)的检测模型。豆科牧草的粗脂肪(EE)、酸性洗涤纤维(ADF)和粗灰分(Ash)建模结果最好,其定标决定系数(RSQ)>0.94,交叉验证相关系数(1-VR)>0.7最高,定标标准分析误差(SEC)在0.071~0.713,交叉校验定标标准分析误差(SECV)在0.160~2.751。禾本科牧草的EE和可溶性糖(WSC)建模结果最好,RSQ分别达0.916和0.859,1-VR分别为0.609和0.810,SEC和SECV分别是0.250、1.488和0.505、3.172。菊科和车前科牧草的模型,除ADF外,其它指标预测的稳定性和准确性较为理想,RSQ在0.85以上,1-VR在0.70以上,SEC和SECV分别在0.361~3.557和0.495~4.602。NIRS对豆科粗蛋白(CP)和WSC的数值预测较差,RSQ仅>0.55,对禾本科CP、ADF、中性洗涤纤维(NDF)、Ash及菊科和车前科的ADF的预测稍差,RSQ>0.7。
A visible/near-infrared reflectance spectroscopy (visible/NIRS) model was developed to determine the forage cultivar nutritional composition of 780 hay samples under simulated rotational grazing. Hay samples ( n=780) from 13 forage crop varieties under simulated grazing in Loess Plateau during the 2012 to 2013 growing season were evaluated using calibration methods for prediction of nutrient contents using NIRS. The following results were obtained. The optimal calibrations in Leguminosae were ether extract (EE), acid detergent fiber (ADF), and crude ash(Ash). The multiple correlation coefficients(RSQ) and 1-variance ratio (1-VR) were > 0.94 and > 0.7, and standard error of calibration (SEC) and standard error of cross validation (SECV) were 0.071~0.713 and 0.160~2.751, respectively. Optimal calibrations in Gramineae were EE and water-soluble carbohydrate content(WSC). RSQ were 0.916 and 0.859 and 1-VR were 0.609 and 0.810 for EE and WSC, respectively, and SEC was 0.250 and 1.488 and SECV was 0.505 and 3.172, respectively. For the other species, the results for nutrient predication were reasonably good, with the exception of ADF. RSQ and 1-VR were >0.85 and >0.70, and SEC and SECV were 0.361~3.557 and 0.495~4.602, respectively. These results indicate that the accuracy of prediction using NIRS was acceptable for 13 forage crop nutrients, although the crude protein(CP), ADF, neutral detergent fibre(NDF), Ash of Gramineae and the ADF of others species (RSQ > 0.7) may require further calibration. The accuracies of the predictions for CP and WSC in Leguminosae (RSQ > 0.55) were not acceptable and thus more samples and greater precision during measurement will be required in further investigations.
随着计算机技术和化学计量学的进步, 近红外检测技术(near infrared reflectance spectroscopy, NIRS)因对样品无损伤、用量少、速度快、精度高而应用范围不断扩大[1], 目前已成功应用于动物、植物、微生物、土壤等生物和非生物样品的成分检测[2]。近红外光谱数据分析的主旨是定量光谱数据与样品成分含量之间的关系, 核心是模型的建立与评价[3]。最小二乘法(PLS)是NIRS数学建模中应用最多的方法[4]; 模型建立后, 其校正和评价的主要参数包括定标决定系数、交叉验证决定系数、外部验证决定系数、交叉验证标准差等, 而且决定系数越大, 标准差越小, 模型的定标及预测效果越好[3]。
NIRS 技术在牧草常规营养成分分析的应用虽然起步较晚, 但是已经成为一个重要分支[5]。1976年, Norris等[6]利用NIRS技术测定了牧草的粗蛋白(CP)。此后, 该技术相继用于分析紫花苜蓿(Medicago sativa)[7]、红三叶(Trifolium pretense)[8]、多年生黑麦草(Lolium perenne)[9]和紫羊茅(Festuca rubra)[10]等牧草的营养成分检测上, NIRS技术分析样品粗蛋白含量的相关系数均在0.90以上。而且, 青贮饲料样品连续4年营养成分建模的决定系数均大于0.90[11]。国外学者[3]对牧草、秸秆类粗饲料中CP、酸性洗涤纤维(ADF)、中性洗涤纤维(NDF)和其它成分均成功建立了NIRS 预测模型。20 世纪末, 我国先后对紫花苜蓿[12]、燕麦(Avena sativa)[13]、黑麦草[14]、羊草(Leymus chinense)[15]等牧草的常规营养指标, 建立了效果良好的NIRS定标模型。以往研究多对天然草原牧草混合样品的营养成分进行实验室检测, 结合近红外光谱技术建立定标模型, 或者对单一的栽培牧草品种进行测定和定标[2]。不同生态区域和生长季节, 各种牧草的养分含量可能大不相同, 而在不同生育期, 对多种利用方式的栽培牧草和混播草地的样品建立近红外定标检查模型的研究鲜有报道。为此, 本研究选取13种优质牧草, 通过定期刈割模拟轮牧, 采集牧草单一样品和混合样品参与定标, 研究混合定标的可行性, 以期为NIRS技术的改进、应用范围的扩展提供基础数据。
2012年4月29日, 12种多年生牧草品种和燕麦在兰州大学榆中草地农业实验站播种(表1), 完全随机区组设计, 每小区面积4 m× 10 m, 每个品种4次重复。每组小区间距1.5 m, 组内小区间隔0.5 m。播种前翻耕, 耕深30~40 cm, 施二胺150 kg· hm-2、尿素300 kg· hm-2作为底肥。每小区平均分为两个裂区, 设置两个处理, 用定期刈割来模拟轮牧, 生长季末收获一次干草。2012年6月29日开始首次轮牧, 牧草株高30 cm左右, 此后每20 d刈割一次, 留茬高度8 cm, 末次齐地面刈割, 全年共计刈割6次。每次模拟轮牧后, 追肥7.5 g· m-2尿素, 并灌水, 灌溉量约为42.24 m3· hm-2。收获干草的裂区刈割时间与模拟轮牧裂区末次刈割时间相同, 齐地面刈割, 施肥、灌溉的时间和量同模拟刈割。2013年, 燕麦品种Makura复播, 方法均与2012年相同; 根据其它他牧草长势, 模拟轮牧的首次刈割时间为6月1日, 末次为10月1日, 共7次, 其它管理均与2012年相同。所有草样经过60 ℃烘干, 粉碎后过0.37 mm筛, 用于测定营养成分和扫描近红外光。
1.2.1 化学成分测定方法 CP含量采用PROXIMA流动分析仪测定, 粗脂肪(EE)采用ANKOMAXT15i型全自动脂肪分析仪滤袋法测定, NDF和ADF采用ANKOMA220型半自动纤维素分析仪进行测定, 粗灰分(Ash)采用TM-O91OP型马弗炉进行测定[16, 17, 18]; 可溶性糖(WSC)采用FOSS流动分析仪测定[19]。
1.2.2 样品光谱测定 用Foss公司XDS-RCA可见/近红外光谱仪扫描牧草样品[2], 光谱波长范围400-2 498 nm, 光谱分辨率2 nm, 波长准确度< 0.05 nm, 数据采样频率2 f· s-1。每样品重复装样扫描3次, 计算其平均值, 用软件自动合并各处理重复和相似数值以增加精确度。双检测器系统:硅检测器(400-1 100 nm), 硫化铅检测器(110-2 500 nm)。WinISⅢ 光谱处理软件建立牧草成分的标准曲线并验证。
1.2.3 定标和校验 把全部样品分成两份, 其中585份作为定标集, 195份作为预测集, 在WinISⅢ 定标软件中对585份样品先进行聚类分析, 筛选出具有代表性的样品514份, 测定营养成分后输入。采用软件自动合并各处理重复值, 将牧草干草分为豆科、禾本科和其它科, 先用定标样品集建立预测模型, 再做交叉检验, 根据交叉验证相关系数(1-VR)、交叉校验定标标准分析误差(SECV)和定标标准分析误差(SEC)等指标对模型进行优化, 并最终确定最佳定标方式和模型。然后, 利用195份预测样品集对模型进行外部验证, 检测其预测效果, 评价模型的外部验证能力。
对定标样品各营养成分NIRS模型进行随机样品验证(表2), 豆科牧草、禾本科和其它科牧草验证集的预测值与化学测定值间的平均偏差(Bias)分别为-0.10~0.21、-0.16~0.24和-0.06~0.05, 牧草定标模型预测效果较为理想。
用WinISⅢ 定标软件对光谱数据进行预处理, 选择全谱范围(400―2 500 nm), 确定主因子数。采用NIRS处理软件自带的共12种回归、光谱和数学处理方法组合对定标集样品光谱进行定标[2](表3)。
豆科牧草定标集样品进行NIRS扫描后, 分别采用改进偏最小二乘法、最小二乘法建立模型, 选择最优模型和光谱处理模式(表4); 其中EE、ADF和Ash的定标决定系数(RSQ)均在0.9以上, 交叉验证相关系数(1-VR)在0.7以上, 定标标准分析误差(SEC)在0.071~0.713, 交叉校验定标标准分析误差(SECV)在0.160~2.751, 说明豆科牧草所建模型在EE、ADF和Ash数值的预测具有较高的稳定性和准确性; 所建定标模型在NDF数值预测方面稍差, RSQ达0.831, 1-VR为0.631, SEC和SECV分别是3.331和4.624, 需对定标集进一步扩充和完善以进一步提高预测准确度; 模型对CP和WSC的预测能力最差, RSQ分别仅为0.557和0.559, 1-VR分别为0.543和0.404, SEC和SECV分别为1.889、2.017和2.461、2.224。
禾本科牧草定标集样品处理参数及模型选择方法与豆科相同(表4); 其中EE和WSC的RSQ分别达到0.916和0.859, 1-VR分别为0.609和0.810, SEC和SECV分别是0.250、1.488和0.505、3.172; 说明所建模型在EE和WSC数值的预测较为稳定和准确; 定标模型在CP、ADF、NDF和Ash的数值预测方面稍差, RSQ在0.75~0.78, 1-VR在0.62~0.65, SEC和SECV分别在0.744~3.472和0.909~5.378。
其它科牧草定标集样品所建模型的处理参数和模型选择方式与上述两类相同(表4); 除ADF数值预测方面稍差RSQ为0.708, 1-VR为0.801, SEC和SECV分别是2.799和4.709; 其它指标数值预测稳定性和准确性较为理想, RSQ在0.85以上, 1-VR在0.70以上, SEC和SECV分别在0.361~3.557和0.495~4.602。
牧草和饲料的常规养分分析是NIRS技术的传统应用领域[3], 国际上建立了系统的NIRS模型库[20]。NIRS技术测定反刍家畜粗饲料、多年生黑麦草和紫羊茅(F. rubra)的CP含量, 模型的相关系数在0.97~0.98[8, 21]。NIRS技术还可以测定多花黑麦草的CF、NDF、ADF含量[22], 玉米秸秆的NDF和ADF含量[23], 紫花苜蓿的Ash等[7, 24]。但是, 由于牧草种类构成、营养成分含量以及环境等诸多因素的差异, 目前还没有一个较为完善的NIRS模型可以适应于各种牧草的生产和多种草地类型[25]。本研究建立了2012-2013年模拟轮牧条件下, 13种优质牧草品种干草营养成分的NIRS模型; 各品种牧草营养成分的NIRS模型, 除豆科CP和WSC外, 定标决定系数均在0.7以上, 交叉验证相关系数均在0.63以上(表4), 均达实用标准, 这与Marten等的结果较为一致[10], 他们虽然建立了4种豆科牧草常规营养成分的NIRS 模型, 并证实NIRS 可分析混合牧草营养成分, 但未涉及不同科牧草NIRS模型的比较。本研究发现, 将NIRS模型运用于不同科牧草品种营养成分的预测效果较为理想, 也表明定标模型对样品集的变异描述能力较强, 定标模型建立较为成功。
然而, 单一牧草样品建模较混合样品更准确[26]。本研究中, 豆科牧草NIRS模型, CP和WSC的决定系数较其它模型小, RSQ仅> 0.55, 禾本科CP、ADF、NDF和Ash及其它牧草的ADF的预测稍差, RSQ> 0.7(表4), 既可能与混合样品有关, 也可能因为实验室中样品测定值存在误差, 需要扩充样本量、降低实验室化学分析的误差等, 以校正定标方程。