一种精确计算多年生禾本科牧草越冬率的方法
【技术领域】
[0001]本发明涉及农业科技领域，特别是涉及一种精确计算多年生禾本科牧草越冬率的方法。
【背景技术】
[0002]禾本科牧草是畜牧养殖业生产中必不可少的主要饲草料之一，在我国分布广泛、资源量大，且营养丰富、易加工，便于储藏运输，也可以与其它成分混合制作不同用途饲料，为畜禽喜食。禾本科牧草抗逆性强、耐刈割、耐践踏，也是我国北方地区庭园绿化、防风固沙及水土保持的优良植物，在草地农业与草原生态系统中具有重要作用。禾本科牧草包括野生和栽培两类，依据寿命长短可以分为一年生、二年生和多年生，其中，二年生和多年生禾本科牧草在人工草地建植与牧草栽培生产中占有绝对优势地位。然而我国北方大部分地区气候寒冷，冬季少雪，尤其是西北干旱区、河北坝上高原及东北地区，极端气温可达零下20-40 V，给普通禾本科牧草及国外引进优质多年生禾本科牧草越冬带来极大挑战，由此，如何确保栽培牧草在寒冷的冬春季安全越冬存活并返青，便成为生产实践的关键问题。
[0003]由牧草栽培学或草地生态学可知，衡量禾本科牧草越冬存活效率的主要指标是越冬率，g卩，单位面积上禾本科牧草越冬返青后的存活数量占越冬前数量的比例，理论上讲，该定义是越冬率的本质含义，也是生产实际中迫切求知的指标，但该指标在实际操作中却难以准确获得，其主要原因与牧草生长特性有关，多年生禾本科牧草类型多样、分蘖与新生枝条繁殖方式复杂，且整个生长阶段植株数量始终处于动态变化状态，即，在进行植株数量统计时，哪一时间点记录的数据才能成为最具代表性的基础数据，将成为越冬率准确计算关键点和难点，由此，如果有一种方法能够精准确定越冬前后牧草植株数量统计的时间点，或者多年生禾本科牧草越冬率计算过程中，植株数量统计均以统一方法进行标准计数，那么，不同种类的多年生禾本科牧草越冬率数据之间才可能存在有效的可比性。
[0004]2006年10月颁布实施的中华人民共和国农业行业标准“草品种审定技术规程”(NY/T 1091-2006)中，对禾本科牧草越冬率的描述是“在小区中选择有代表性的样段两处，每段长lm，在越冬前及第二年返青后分别计算样段中植株总数及返青数，便可统计越冬率”，其中，越冬前及第二年返青后的具体计数时间指定不明，实际测定中难以操作，其原因是多年生禾本科牧草越冬前及第二年返青后物候时间较长，尤其是第二年返青期有的甚至达到1-2个月，在此期间，植株数量一直处于变化状态，不同时间点样段内植株数量不同，则选择不同的时间点计数时，计算结果差异很大，则不同物种或处理间得出的越冬率缺乏可比性，另外，每个试验处理仅仅取两样段，则不能进行统计学的方面的差异显著性分析，再次降低了该试验数据的有效性和实用性。
[0005]2014年7月颁布实施的中华人民共和国国家标准“草品种审定技术规程”(GB/T30395-2013)中的描述是“在同一区组的小区中随机选择有代表性的样段3处，每段长lm，在越冬前后分别计数样段中植株数量，计算越冬率”，其中，将2006版中“在小区中选择有代表性的样段两处”修改为“在同一区组的小区中随机选择有代表性的样段3处”，避免了差异显著性分析时的统计数据不足，但将“越冬前及第二年返青后分别计算样段中植株总数及返青数”修改为“越冬前后分别计数样段中植株数量”，使实际测定时植株计数时间点的选择更加宽泛和模糊，进一步降低了不同物种或处理间计算结果的可比性，由此，需要结合多年生禾本科牧草生长发育等生物学特性和物候特征，选择一种准确性高、可比性强、执行中容易操作的测定方法，才是解决越冬率计算问题的关键所在。
【发明内容】
[0006]为了解决目前多年生禾本科牧草越冬率计算过程中，植株计数时间点难以准确选择、越冬率计算结果不唯一、测定数据缺乏有效性和可比性等主要问题，进一步提高该指标测定中的可操作性及其在生产实践中的指导意义，本发明的目的是提供一种精确计算多年生禾本科牧草越冬率的方法。
[0007]本发明提供一种精确计算多年生禾本科牧草越冬率的方法，包括以下步骤:
[0008]I)选定区域，在植株冬前休眠期时统计越冬前植株数量，越冬后在植株返青数量达到10 %时，统计首次越冬后植株数量，之后每10天计数一次，共统计M次越冬后植株计数；
[0009]2)以步骤I)统计的数据拟合如下方程:
[0010]N=K/[l+e(a_rt)]
[0011]其中，t为测定日期，首次测定日期在计算公式中可以计数为10，其余测定日期均为10的倍数，N为测定日计数的植株数量，然后根据连续若干次的测定日期和植株计数数据，求得对应的函数关系式N=f(t);
[0012]3)以步骤2)的结果，求函数N的二阶导数等于零时的t值和N值，S卩，t是植株计数的测定日期，N是曲线拐点时样段内越冬植株的拐点数量；
[0013 ] 4)以步骤3)所求N值，按如下公式计算越冬率:
[0014]越冬率=越冬后样段内植株拐点数量/越冬前休眠期样段植株数量X100 %。
[0015]其中，所述M次为5次及以上。
[0016]其中，所述选定区域为选Im长的样条。样条宽度以播种后的自然生长单垄为条宽。
[0017]其中，所述选定区域，同一处理设3个重复。
[0018]本发明还提供所述精确计算多年生禾本科牧草越冬率的方法在牧草种植中的应用。
[0019]本发明的关键点如下:
[0020]1、越冬前后植株计数时间点的确定:越冬前后植株计数是获得越冬率准确计算结果的基础，其中，越冬前的休眠期，植株形态与生理变化仅限于正常的生理代谢，数量基本不变，该时间段内测定的数据相对稳定;另外，越冬返青后，由于植株生长较快，不同时间点测定的数量差异较大，所以需要多次连续测定，才能最终确定合理的计算结果。
[0021]2、拟合方程式的选择:栽培后的多年生禾本科牧草草地，群落物种比较单一，基本上是由栽培种构建的单一种群，其种群数量的发展受到环境资源影响，变化趋势符合生态学常用的连续增长模型一一逻辑斯谛方程，所以，以该方程为模型，求得方程中的主要参数，即得植株数量变化公式，这也是栽培牧草返青后的植株数量实际发展轨迹。
[0022]3、函数方程式二阶求导:严格地讲，越冬后计数的植株与越冬前应该有一定的对应性，是越冬率计算的理论基础，但是，实际环境中，植株越冬返青后数量的变化一般分为2个阶段，第一个阶段是越冬休眠芽在春季返青后生长发育为正常成株，而这一批成株也正是越冬后应该计入统计数据的植株，第二个阶段是生长发育后的成株产生新一级分蘖，且新一级分蘖又生长发育，最终成为又一批新的成株，同时，这两个阶段形成的正常成株在形态上难以区分，从而为植株的正确选择计数带来很大困难。不过，从高等数学知识可知，对函数求二阶导数，并使二阶导数为零，即可获得函数曲线的拐点值，也就是曲线斜率变化的极值，该点对应的自变量是测定日期，对应函数则是该测定日的植株拐点数量，正常情况下，牧草植株第一阶段与第二阶段的生长发育存在一定时间间隔，并在整体数量变化规律上表现为连续性的暂时打破，即为一个转折点，这也符合函数曲线拐点特征，由此，选择方程式进行二阶求导，并令导数为零，即可获得曲线变化中的拐点值，这也是越冬植株返青后完成第一阶段生长发育而成为正常成株关键点，同时也是第二阶段植株分蘖开始生长发育的起点，因此，由拐点处植株数量值除以越冬前休眠期植株数量值即为准确的植株越冬率计算结果。
[0023]本发明具有以下有益效果:
[0024]根据多年生禾本科牧草自身生长特性和物候期特征，在植株越冬前休眠期和越冬返青后多次定期计数，符合植株生长发育基本规律，再通过生态学常用的种群增长模型逻辑斯谛方程，寻找返青植株数量增长拐点，计算