Effects of mowing modes on biomass, crude protein and crude fat contents of Broussonetia papyrifera in karst area
-
摘要:目的 构树是重要的生态和经济树种,在喀斯特地区大面积种植,但是目前刈割模式对其产量及品质的影响研究比较缺乏,导致构树经济林资源利用率低,产量和质量不高。本研究通过对不同刈割频度和留茬高度下构树生物量和营养品质的研究,明确最佳的刈割模式,为构树林采收利用提供依据。方法 本研究以贵州省贞丰县2年生构树饲用林为研究对象,分3个刈割频度(一年刈割3次、4次、5次)和3个留茬高度(齐地割0 cm、留茬15 cm、留茬30 cm)进行刈割处理,研究不同处理间发枝数,收获的枝、叶生物量,及粗蛋白、粗脂肪含量的差异。结果 一年刈割4次和5次构树的发枝数明显高于3次。随着刈割次数和留茬高度的增加,茎叶比逐渐降低。累积生物量、粗蛋白含量均不受刈割频次的影响,但在不同留茬高度间差异显著。其中,留茬30 cm时生物量最高,为2.63 t/hm2,齐地割时粗脂肪含量最高。结论 年刈割3 ~ 5次,构树产量和品质差异不大;留茬30 cm,构树生物量最好,但粗脂肪含量略有下降。因此,综合考虑成本等因素,推荐喀斯特地区构树饲用林采用一年刈割3 ~ 4次、留茬30 cm的模式进行采收利用。Abstract:Objective Broussonetia papyrifera is an important ecological and economic species, which is widely planted in karst area. However, at present, research on the impact of mowing patterns on its yield and quality is seldom, resulting in low resource utilization, low yield and quality of B. papyrifera stand. Therefore, this research studied the effects of different mowing frequencies and stubble heights on the yield and nutrient contents of B. papyrifera, in order to determine the optimal mowing modes and guide the harvesting and utilization of B. papyrifera.Method Research objects were two-year-old B. papyrifera in Zhenfeng County, Guizhou Province of southwestern China. Three mowing frequencies (cutting 3 times, 4 times and 5 times a year) and three stubble heights (stubble 0 cm, 15 cm and 30 cm) were randomly set up. The number of sprouting, biomass of harvested branches and leaves, crude protein and crude fat contents among different treatments were studied.Result The number of resprouting of B. papyrifera under 4 and 5 time mowing a year was significantly higher than 3 times. With the increase of mowing times and stubble heights, stem to leaf ratio decreased gradually. The biomass, crude protein were not affected by mowing frequency. However, stubble heights significantly affected biomass, and showed that the biomass of stubble 30 cm was the highest (2.63 t/ha). The crude fat content of neat mowing was significantly higher than that of stubble 15 cm and 30 cm.Conclusion There is little difference in yield and quality of B. papyrifera among 3−5 times mowing a year. Yield with stubble of 30 cm is the highest, but the crude fat content decreases slightly. Therefore, it is recommended that mowing 3−4 times a year and stubble height of 30 cm for harvesting and utilization in B. papyrifera forage forest.
-
Keywords:
- Broussonetia papyrifera /
- mowing frequency /
- stubble height /
- biomass /
- crude protein /
- crude fat
-
森林的天然更新是一个持久而复杂的生态学过程[1-2],是森林生态系统自我繁衍及生态恢复的重要方式[3]。幼苗是森林生态系统的重要组成部分[4],在森林的天然更新过程中发挥着重要作用[5]。因此,深入了解木本植物幼苗的组成、结构、存活、生长规律及更新的限制因子对森林生态系统管理和生物多样性保护意义重大[6]。
幼苗更新的组成与结构特征是天然更新研究的基础。刘帅[7]对长白山针阔混交林乔木幼苗进行动态研究,发现群落乔木幼苗年际间物种组成变化较小,幼苗个体的新生和死亡导致群落乔木幼苗数量有较大的年际和月际变化;Comita等[8]研究表明,在大树多度长期保持不变的情况下,幼苗多度应与周围大树多度呈线性关系。在天然林中,更新的幼苗幼树多呈现聚集分布格局。李艳丽等[9]通过长白山云冷杉针阔混交林天然更新空间分布格局的研究发现,样地中更新幼苗幼树呈现斑块聚集分布;李雪云等[10]对闽楠天然次生林幼树幼苗研究表明更新层树种空间分布格局总体上呈现聚集分布的状态,且聚集强度较大。幼苗的组成与分布受多种因素影响,适宜的环境条件是幼苗更新的关键。幼苗存活及空间分布与光照、土壤湿度、凋落物等因子息息相关[11-12]。同时,幼苗的存活率在不同的地形中存在很大的差异[13]。地形因子如海拔、坡向和坡度通过影响光照和水分的再分配,从而影响幼苗的更新[14-15]。而食草动物和昆虫等通过对幼苗的取食和损伤限制幼苗活力,引起幼苗死亡[16]。另外,植物病原体(如细菌、真菌、病毒等)以及资源竞争引起的负密度制约效应[17]也是幼苗死亡的重要驱动因素。当然,充足的种源是保证森林更新的物质前提和基础[18]。在森林群落中,无论是种群水平还是群落尺度,森林下层幼苗的组成、数量和存活等动态特征随时间和空间发生较大变异[19-21],因此,研究长时间尺度上幼苗动态格局,对揭示森林群落稳定机制意义重大。茂兰喀斯特常绿落叶阔叶混交林是我国亚热带地区主要的天然林生态系统和喀斯特区域的非地带性顶极植被,蕴藏着丰富的生物多样性,具备重要的生态功能[22]。但其生态系统十分脆弱,成土较慢且土层浅薄不连续,植物生长缓慢[23-24],研究该系统幼苗组成和更新特征对深入揭示其生态系统稳定机制具有十分重要的科学价值和现实意义[25]。目前,国内学者对茂兰喀斯特常绿落叶阔叶混交林的研究多集中在物种组成[26-27]、退化森林的修复[28]、种子繁殖[29]、种群的结构动态[22]、物种分布与环境的关系[30]等方面,但有关茂兰喀斯特森林木本植物幼苗的组成及年际变化研究报道较少[31]。
本研究以茂兰喀斯特常绿落叶阔叶混交林为研究对象,基于2017—2019年连续3次的幼苗调查数据,分析该群落木本植物幼苗的组成、分布及其年际动态变化,探讨幼苗密度在空间上的分布特征。研究结果有助于揭示和预测该森林群落组成的动态变化,为喀斯特森林的合理管理和保护提供科学依据。
1. 研究地概况与研究方法
1.1 研究地概况
茂兰国家级自然保护区位于贵州省南部黔、桂交界处,总面积约212.85 km2,地理坐标为25°09′20″N、107°52′10″E,海拔为430 ~ 1 078 m,与广西木伦国家级自然保护区相邻。茂兰保护区属于中亚热带石灰岩区常绿落叶阔叶混交林生态系统,是全球喀斯特地貌区幸存的连片面积最大、保存完整且原生性强的喀斯特森林。保护区主要是由纯质石灰岩及白云岩构成的典型喀斯特地貌,局部地点覆盖有少量砂页岩,岩石裸露率高达90%以上。该区具有中亚热带山地湿润气候的特点,年平均气温为15.3 ℃,年均降雨量为1 320.5 mm,太阳辐射年总量为63 289.8 kW/m2。茂兰保护区内区系成分复杂,生物种类丰富。保护区内有维管束植物154科514属1 203种,种子植物143科494属1 172种,裸子植物6科,蕨类植物11科20属31种,森林覆盖率约为87%。主要树种有齿叶黄皮(Clausena dunniana)、粗糠柴(Mallotus philippensis)、贵州鹅耳枥(Carpinus kweichowensis)、香叶树(Lindera communis)、黄梨木(Boniodendron minus)、南方桂樱(Laurocerasus australis)、朴树(Celtis sinensis)、青冈(Cyclobalanopsis glauca)、狭叶润楠(Machilus rehderi)、杨梅叶蚊母树(Distylium myricoides)、西南卫矛(Euonymus hamiltonianus)、裂果卫矛(Euonymus dielsianus)、榆树(Ulmus pumila)等。
1.2 研究方法
1.2.1 样地建立与调查
按照CTFS(the Center for Tropical Forest Science)样地建设标准在茂兰国家级自然保护区常绿落叶阔叶混交林中,选择具有代表性地段设置面积为1.28 hm2(160 m × 80 m)的固定样地。将该样地等面积划分为128个10 m × 10 m样方,采用相邻网格法调查每个样方内所有胸径 ≥ 1 cm的活立木胸径、树高、树种、坐标等信息。在距每个10 m × 10 m样方起始点1 m处分别建立1个1 m × 1 m的小样方(共128个),调查每个小样方内所有胸径小于1 cm的木本植物,包括乔木、灌木和木质藤本幼苗的地径、高度和树种名。调查工作分别在2017年、2018年和2019年的8月进行,同时记录样地的海拔和坡度等地形因子。
根据幼苗的高度,将其划分为不同的高度等级:0 ~ 5 cm、5 ~ 10 cm、10 ~ 20 cm、20 ~ 30 cm、30 ~ 40 cm、40 ~ 50 cm、50 ~ 60 cm、60 ~ 70 cm、70 ~ 80 cm、80 ~ 90 cm、90 ~ 100 cm、≥ 100 cm,分析幼苗的个体数及物种数在不同高度级上的分布。根据样地的海拔和坡度范围,将海拔分为4级(Ⅰ级:790 ~ 810 m、Ⅱ级:810 ~ 830 m、Ⅲ级:830 ~ 850 m、Ⅳ级:850 ~ 870 m);坡度分为平缓坡(0 ~ 16°)、斜坡(16° ~ 26°)、陡坡(26° ~ 36°)、急险坡(> 36°)。
1.2.2 植物种类的确定
参照《中国植物志》《贵州植物志》及相关文献,对调查的木本植物幼苗进行鉴定,建立植物信息数据库。
1.2.3 数据分析
采用相对多度和相对频度来计算各物种的重要值:幼苗重要值 = (相对多度 + 相对频度)/2。其中,相对多度 = (某种的多度/所有种的多度和) × 100%,相对频度 = (某种的频度/所有种的频度和) × 100%。分别统计每个幼苗样方中的木本植物幼苗个体数和物种数,并计算每个样方的幼苗密度:幼苗密度 = 幼苗个体数/样方面积。计算幼苗新生率和死亡率,新生率 = (调查年份当年生个体数/调查年份总个体数) × 100%,死亡率 = (经过1年或者几年死亡的个体数/调查年份总个体数) × 100%[32]。
地统计学中的半方差函数(也称为变异函数)经常用于描述空间异质性[33],其计算公式如下:
λ(h)=12N(h)N(h)∑i=1[Z(xi)−Z(xi+h)]2 (1) 式中:λ(h)为变异函数,N(h)为
以h为间距的所有观测点对数量,Z(xi)为系统某属性Z在空间位置xi处的值,Z(xi + h)为Z在(xi + h)处的值,h为空间距离。通过半方差函数及曲线图可得到以下重要参数:基台值(sill,C0 + C)、块金值(nugget,C0)和变程(range,A)。基台值(C0 + C)是系统或系统属性中最大变异程度,块金值(C0)代表不能被参数所解释的随机变化,表示随机因素产生的空间异质性;变程反映空间变异的尺度,在变程内,变量具有空间自相关,反之不存在[34]。空间结构比(C0/(C0 + C))常用来反映随机部分引起的空间变异在总空间异质性中的比例[33]。 利用单因素方差分析(one-way ANOVA)和Ducan多重比较分析幼苗个体数和种类在不同海拔和坡度上的差异。利用地统计学软件GS+(Gramma design software)9.0对幼苗密度进行克里格空间插值分析。采用Microsoft Excel 2016和SPSS 24软件进行数据统计分析。
2. 结果与分析
2.1 幼苗数量与物种组成
2017—2019年的3次调查共计录到161种木本植物幼苗,分属于63科124属(表1)。树种较多的是蔷薇科(Rosaceae)、豆科(Fabaceae)、紫金牛科(Myrsinaceae)、芸香科(Rutaceae)、大戟科(Euphorbiaceae)、樟科(Lauraceae)等。2017—2019年记录到的乔木幼苗、灌木幼苗和木质藤本幼苗个体数共3 080株(表1)。2019年记录到的幼苗个体数最多,为1 342株(10 株/m2),其次是2017年,为971株(8 株/m2),2018年最少,为767株(6 株/m2)。调查期间,物种数和个体数在年际间存在一定的差异,2018年调查到的灌木幼苗和木质藤本幼苗是3年调查中种类最少的。
表 1 3年调查到的幼苗科、属、种Table 1. Seedling families, genera and species investigated in the past three years年份
Year所有幼苗 All seedlings 乔木幼苗 Tree seedling 灌木幼苗 Shrub seedling 木质藤本幼苗
Seedling of woody vine幼苗数量 Seedling number 2017 55科96属121种55 families,
96 genera, 121 species29科47属55种 29 families, 47 genera, 55 species 31科42属49种 31 families, 42 genera, 49 species 15科16属17种 15 families, 16 genera, 17 species 971 2018 40科66属78种 40 families,
66 genera, 78 species25科38属44种 25 families, 38 genera, 44 species 15科21属24种 15 families, 21 genera, 24 species 8科10属10种 8 families,
10 genera, 10 species767 2019 39科70属87种 39 families,
70 genera, 87 species22科34属42种 22 families, 34 genera, 42 species 19科27属30种 19 families, 27 genera, 30 species 13科15属15种 13 families, 15 genera, 15 species 1 342 2017—2019 63科124属161种 63 families, 124 genera, 161 species 34科61属74种 34 families, 61 genera, 74 species 35科50属62种 35 families, 50 genera, 62 species 20科25属26种 20 families, 25 genera, 26 species 3 080 2017—2019年,共有17个幼苗样方没有记录到任何木本植物幼苗,27个幼苗样方仅调查到1个树种,203个幼苗样方调查到2 ~ 5个树种,6个树种以上的幼苗样方有127个(图1)。2017、2018和2019年分别有1个、10个和6个样方没有调查到木本植物幼苗。3次调查的幼苗样方中,物种数主要集中在1 ~ 7种。10个物种以上的幼苗样方在2017年、2018年和2019年分别有12个、0个和12个。幼苗个体数主要集中在2 ~ 7株,仅在2019年调查到1个幼苗样方的个体数为42株。2018年调查到2株幼苗的样方数最多,为16个,2017年5株的样方数最多,为14个,2019年5株和6株最多,都为12个。2018年调查的幼苗个体数大于15株的样方数明显少于2017年和2019年。
从表2中可知,茂兰喀斯特常绿落叶阔叶混交林林下幼苗物种丰富,络石、珍珠莲、齿叶黄皮、杨梅叶蚊母树、狭叶润楠和南方桂樱等树种在3次调查中重复出现,且数量较多。络石和珍珠莲为木质藤本。络石在3年调查中的数量最多,且重要值排名第一,分别为6.80、10.24和8.35。珍珠莲3次调查的重要值也较大。乔木幼苗中,齿叶黄皮的重要值在2017—2019年都排名前5,分别为4.30、4.98和4.78。杨梅叶蚊母树、狭叶润楠和南方桂樱等树种的重要值也较大。
表 2 茂兰喀斯特森林木本植物幼苗的种类及数量特征Table 2. Species and quantity characteristics of woody plant seedlings in karst forest of Maolan树种 Species 科 Family 属 Genus 幼苗数量
Number of seedling出现样方数
Number of quadrat presented重要值
Important value2017 2018 2019 2017 2018 2019 2017 2018 2019 齿叶黄皮 Clausena dunniana 芸香科 Rutaceae 黄皮属 Clausena 38 36 67 30 24 31 4.30 4.98 4.78 杨梅叶蚊母树 Distylium myricoides 金缕梅科 Hamamelidaceae 蚊母树属 Distylium 45 30 27 27 16 14 4.43 3.71 2.04 针齿铁仔 Myrsine semiserrata 紫金牛科 Myrsinaceae 铁仔属 Myrsine 35 6 12 25 4 7 3.76 0.83 0.96 朴树 Celtis sinensis 大麻科 Cannabaceae 朴属 Celtis 30 16 17 25 13 15 3.50 2.47 1.74 狭叶润楠 Machilus rehderi 樟科 Lauraceae 润楠属 Machilus 38 24 34 19 15 13 3.44 3.21 2.23 南方桂樱 Laurocerasus australis 蔷薇科 Rosaceae 桂樱属 Laurocerasus 28 26 29 19 18 15 2.93 3.67 2.19 粗糠柴 Mallotus philippensis 大戟科 Euphorbiaceae 野桐属 Mallotus 13 11 22 13 9 14 1.69 1.70 1.85 广西密花树 Rapanea kwangsiensis 紫金牛科 Myrsinaceae 密花树属 Rapanea 18 4 20 9 2 12 1.63 0.48 1.63 榆树 Ulmus pumila 榆科 Ulmaceae 榆属 Ulmus 14 7 7 10 7 6 1.50 1.22 0.70 贵州鹅耳枥 Carpinus kweichowensis 桦木科 Betulaceae 鹅耳枥属 Carpinus 11 10 14 6 6 11 1.04 1.31 1.33 野柿 Diospyros kaki var. sylvestris 柿科 Ebenaceae 柿属 Diospyros 8 3 6 7 3 4 0.96 0.52 0.52 小叶柿 Diospyros mollis 柿科 Ebenaceae 柿属 Diospyros 6 15 35 6 12 24 0.78 2.29 3.07 西南卫矛 Euonymus hamiltonianus 卫矛科 Celastraceae 卫矛属 Euonymus 5 19 3 5 9 3 0.65 2.23 0.33 海桐 Pittosporum tobira 海桐科 Pittosporaceae 海桐属 Pittosporum 6 17 8 4 3 6 0.62 1.44 0.74 伞花木 Eurycorymbus cavaleriei 无患子科 Sapindaceae 伞花木属 Eurycorymbus 3 3 2 3 3 1 0.39 0.52 0.15 香叶树 Lindera communis 樟科 Lauraceae 山胡椒属 Lindera 2 13 25 2 7 14 0.26 1.62 1.96 圆叶乌桕 Triadica rotundifolia 大戟科 Euphorbiaceae 乌桕属 Triadica 2 2 9 2 1 5 0.26 0.24 0.70 梾木 Swida macrophylla 山茱萸科 Cornaceae 山茱萸属 Cornus 1 3 4 1 3 6 0.13 0.52 0.59 黄梨木 Boniodendron minus 无患子科 Sapindaceae 黄梨木属 Boniodendron 9 17 25 9 17 19 1.17 2.97 2.33 角叶槭 Acer sycopseoides 槭树科 Aceraceae 槭属 Acer 16 1 3 4 1 20 1.14 0.17 1.59 九里香 Murraya exotica 芸香科 Rutaceae 九里香属 Murraya 9 10 35 8 9 3 1.09 1.64 1.53 长叶冻绿 Rhamnus crenata 鼠李科 Rhamnaceae 鼠李属 Rhamnus 1 2 5 1 1 5 0.13 0.24 0.56 异叶花椒 Zanthoxylum dimorphophyllum 芸香科 Rutaceae 花椒属 Zanthoxylum 8 5 4 8 3 2 1.04 0.65 0.30 小叶女贞 Ligustrum quihoui 木犀科 Oleaceae 女贞属 Ligustrum 5 10 18 4 8 14 0.57 1.53 1.70 十大功劳 Mahonia fortunei 小檗科 Berberidaceae 十大功劳属 Mahonia 5 3 3 3 3 3 0.49 0.52 0.33 南天竹 Nandina domestica 小檗科 Berberidaceae 南天竹属 Nandina 45 18 17 24 10 10 4.20 2.27 1.37 长梗罗伞 Brassaiopsis glomerulata 五加科 Araliaceae 罗伞属 Brassaiopsis 9 3 5 5 2 4 0.85 0.41 0.48 棕竹 Rhapis excelsa 棕榈科 Arecaceae 棕竹属 Rhapis 13 6 7 7 4 3 1.22 0.83 0.48 九管血 Ardisia brevicaulis 紫金牛科 Myrsinaceae 紫金牛属 Ardisia 3 2 34 2 2 14 0.31 0.35 2.30 箬竹 Indocalamus tessellatus 禾本科 Poaceae 箬竹属 Indocalamus 3 12 21 2 2 3 0.31 1.00 1.00 硃砂根 Ardisia crenata 紫金牛科 Myrsinaceae 紫金牛属 Ardisia 5 7 1 3 4 1 0.49 0.89 0.11 苎麻 Boehmeria nivea 荨麻科 Urticaceae 苎麻属 Boehmeria 8 5 7 4 3 3 0.72 0.65 0.48 雀梅藤 Sageretia thea 鼠李科 Rhamnaceae 雀梅藤属 Sageretia 5 6 3 5 6 3 0.65 1.05 0.33 寒莓 Rubus buergeri 蔷薇科 Rosaceae 悬钩子属 Rubus 5 10 9 2 6 5 0.41 1.31 0.70 络石 Trachelospermum jasminoides 夹竹桃科 Apocynaceae 络石属 Trachelospermum 85 105 151 31 31 37 6.80 10.24 8.35 球兰 Hoya carnosa 夹竹桃科 Apocynaceae 球兰属 Hoya 16 5 12 6 5 5 1.29 0.87 0.82 菝葜 Smilax china 菝葜科 Smilacaceae 菝葜属 Smilax 23 17 40 13 13 25 2.20 2.53 3.33 珍珠莲 Ficus sarmentosa
var. henryi桑科 Moraceae 榕属 Ficus 19 24 106 13 15 27 2.00 3.21 5.94 常春藤 Hedera nepalensis var. sinensis 五加科 Araliaceae 常春藤属 Hedera 28 5 11 20 5 6 3.01 0.87 0.85 链珠藤 Alyxia sinensis 夹竹桃科 Apocynaceae 链珠藤属 Alyxia 1 19 40 1 12 15 0.13 2.55 2.60 小计 Sub-total 624 537 898 388 317 428 62.49 69.71 64.99 其他 Others 347 230 444 251 139 250 37.51 30.29 35.01 2.2 幼苗高度级分布
3年调查的木本植物幼苗个体数和物种数都呈现出偏峰现象,在10 ~ 20 cm的高度级达到高峰,在 ≥ 100 cm高度级又出现一个小高峰(图2)。10 ~ 20 cm高度级的幼苗个体数最多的是2019年(404株),其次是2017年(221株),2018年为207株;物种数最多的是2017年(68种),2018年和2019年分别为52种和60种。2017年在 ≥ 100 cm高度级的幼苗个体数和物种数都是最多的。幼苗个体数和物种数主要集中分布在5 ~ 30 cm高度级间。
2.3 主要树种幼苗新生和死亡数量组成
对齿叶黄皮、粗糠柴、贵州鹅耳枥、黄梨木、南方桂樱、朴树、西南卫矛、狭叶润楠、杨梅叶蚊母树和榆树10个主要树种幼苗的新生和死亡情况进行统计(表3)。2017—2019年间,齿叶黄皮、粗糠柴、黄梨木、南方桂樱、朴树、狭叶润楠和杨梅叶蚊母树新生幼苗数量均大于10,新生率均达到80%以上;齿叶黄皮、南方桂樱、狭叶润楠和杨梅叶蚊母树死亡幼苗数量均大于10,其中,南方桂樱的死亡率高达80%以上。2018年和2019年新生幼苗总数变化不大,但死亡幼苗总数变化较大。2018年,齿叶黄皮、粗糠柴、贵州鹅耳枥、南方桂樱、朴树、狭叶润楠、杨梅叶蚊母树和榆树的死亡数量都大于新生数量。齿叶黄皮、粗糠柴、南方桂樱、朴树和狭叶润楠的幼苗数量在2017—2019年间先减少后增加。新生幼苗数量较多时,死亡幼苗数量也较多,如齿叶黄皮、南方桂樱和杨梅叶蚊母树等树种。
表 3 10个主要树种的幼苗物种新生和死亡数量组成Table 3. Composition of seedling species regeneration and death of 10 main tree species树种 Species 幼苗数量
Number of seedling出现样方数
Number of quadrat presented新生幼苗数量(新生率) Number of newborn seedling (newborn rate) 死亡幼苗数量(死亡率) Number of dead seedling
(mortality rate)2017 2018 2019 2017 2018 2019 2017—2018 2018—2019 2017—2018 2018—2019 齿叶黄皮 Clausena dunniana 38 36 67 30 24 31 29(80.6%) 48(71.6%) 31(81.6%) 17(47.2%) 粗糠柴 Mallotus philippensis 13 11 22 13 9 14 10(90.9%) 18(81.8%) 12(92.3%) 7(63.6%) 贵州鹅耳枥 Carpinus kweichowensis 11 10 14 6 6 11 10(100.0%) 9(64.3%) 11(100.0%) 5(50.0%) 黄梨木 Boniodendron minus 9 17 25 9 17 19 14(82.4%) 22(88.0%) 6(66.7%) 14(82.4%) 南方桂樱 Laurocerasus australis 28 26 29 19 18 15 22(84.6%) 26(89.7%) 24(85.7%) 23(88.5%) 朴树 Celtis sinensis 30 16 17 25 13 15 10(62.5%) 10(58.8%) 24(80.0%) 9(56.3%) 西南卫矛 Euonymus hamiltonianus 5 19 3 5 9 3 19(100.0%) 3(100.0%) 5(100.0%) 19(100.0%) 狭叶润楠 Machilus rehderi 38 24 34 19 15 13 11(45.8%) 21(61.8%) 25(65.8%) 11(45.8%) 杨梅叶蚊母树 Distylium myricoides 45 30 27 27 16 14 24(80.0%) 11(40.7%) 39(86.7%) 14(46.7%) 榆树 Ulmus pumila 14 7 7 10 7 6 5(71.4%) 6(85.7%) 12(85.7%) 6(85.7%) 小计 Sub-total 231 196 245 163 134 141 154(78.6%) 174(71.0%) 189(81.8%) 125(63.8%) 2.4 幼苗的分布特征
2.4.1 幼苗的地理分布
幼苗个体数和物种数均随着海拔和坡度的升高呈现增加的趋势(图3)。Ⅰ级海拔、Ⅱ级海拔和Ⅳ级海拔之间幼苗个体数分布存在显著差异,Ⅰ级海拔和Ⅳ级海拔之间幼苗物种数分布存在显著差异。在急险坡地段上分布的幼苗个体数和物种数均高于平缓坡、斜坡和陡坡的。
2.4.2 幼苗密度的空间分布
2017年幼苗密度的最优拟合模型为高斯模型,2018年和2019年幼苗密度的最优拟合模型都为线性模型,R2值在0.50 ~ 0.90之间(表4),模型拟合效果较好。3年幼苗密度的块金值(C0)存在明显差异,如2017年幼苗密度的块金值为3.200,远小于2018年的21.050和2019年的63.482。基台值(C0 + C)2019年最大,为84.584,其次为2017年的33.630,2018年最小,为25.126。
表 4 茂兰喀斯特森林木本植物幼苗密度半方差函数的模型类型及参数Table 4. Fitted model types and parameters for semivariograms of woody plant seedling density from karst forest of Maolan年份 Year 最优模型
Optimal model块金值
Nugget (C0)基台值
Sill (C0 + C)空间结构比
Spatial structure ratio C0/(C0 + C)变程
Range (A)/mR2 2017 高斯模型 Gaussian model 3.200 33.630 0.095 14.896 0.654 2018 线性模型 Linear model 21.050 25.126 0.838 76.850 0.521 2019 线性模型 Linear model 63.482 84.584 0.751 76.850 0.863 注:C 表示拱高。Note: C means arch height. 幼苗密度的空间分布具有一定的斑块性(图4)。2017年幼苗密度主要分布在5 ~ 10 株/m2、20 ~ 25 株/m2和25 ~ 30 株/m2的幼苗零星分布在样地中。2018年幼苗密度分布在0 ~ 5 株/m2和5 ~ 10 株/m2,密度为20 ~ 25 株/m2分布在样地的西北和东南方向;2019年幼苗密度分布在20 株/m2以上的数量较2017年和2018年都大,主要分布在样地的西北、北和东南方向,密度为0 ~ 5 株/m2的分布在样地边缘。
3. 讨 论
3.1 幼苗组成及动态变化
2017—2019年间,对喀斯特常绿落叶阔叶混交林内128个1 m2的小样方进行3次调查共记录到木本植物幼苗63科124属161种3 080株(8株/m2)。各年调查到的物种均较多,表明喀斯特常绿落叶阔叶混交林幼苗的物种组成较为丰富,与刘检明[35]在八大公山亚热带山地常绿落叶阔叶混交林1.2 hm2样地共调查到的幼苗28科43属68种3 101株相比数量上相似,物种组成存在差异。齿叶黄皮、南方桂樱、狭叶润楠和杨梅叶蚊母树等树种在3次调查中重复出现且重要值在前10位,表明林下幼苗的主要树种成分与该森林群落内主要树种组成保持着一定的相似性,主要树种的幼苗在林下同样保持着优势地位[6, 36-37]。2017—2019年3次调查共记录到木质藤本幼苗20科25属26种814株,其中络石和珍珠莲的重要值较高,说明木质藤本幼苗是茂兰喀斯特常绿落叶阔叶混交林林下幼苗群落中的重要组成部分。在森林群落中,木质藤本具有木质化的茎,但是不能直立生长,通过攀附或缠绕其他植物(主要是乔木)的茎到达森林冠层[38],与树木在地上竞争空间和光照,在地下争夺水分和养分[39]。同时木质藤本也可以通过机械损伤和资源竞争来降低幼苗的存活与生长[40]。杨梅叶蚊母树幼苗的数量在逐年递减,可能受到林下木质藤本的影响。
2017年、2018年和2019年的幼苗密度分别为8、6和10 株/m2。不同年份幼苗数量的波动现象主要是由幼苗新生和死亡所驱动[41-43]。齿叶黄皮、南方桂樱和杨梅叶蚊母树等幼苗大量产生后,其对应的死亡率也出现相对较高的现象,这与张健等[6]、姚杰等[36]对林下幼苗更新研究得到的“幼苗大量发生、大量死亡”的研究结果一致。种子是幼苗更新的来源,幼苗的数量和分布很大程度上依赖于种子的可获得性[44]。2018年调查到的幼苗物种数较多,但是幼苗个体数仅有767株,可能和这一年树种的种子产量较少或者是种子发芽能力低有关。种子萌发需要的外界基本条件是水份、温度和氧气[45]。适宜的温度能促进种子的萌发和幼苗的生长[46],温度过高或过低均会影响种子活力造成发芽和出苗不良[47],特别是在温度不断降低的趋势之下,种子的发芽态势、发芽率、活力指数、发芽指数以及发芽速度指数都会随之降低[48]。经查阅相关气象资料,发现喀斯特常绿落叶阔叶混交林2018年1月和2月的最低温度都小于0 ℃,且低温天气持续时间较2017年和2019年长。因此,温度变化是导致喀斯特常绿落叶阔叶混交林幼苗数量年际差异的主要原因。
3.2 幼苗高度级结构
幼苗个体数和物种数都是随着高度级的增加呈先增后减趋势,10 ~ 20 cm高度级达到峰值。在早期阶段,种子在母树周围容易形成较高的幼苗密度,并且能成功定居[41],所以数量不断增加。当密度达到一定值时,资源竞争加剧,幼苗个体间的种间竞争和种内的相互作用导致幼苗死亡率升高[49],所以在高度级大于10 ~ 20 cm后,幼苗数量出现下降的现象。经过环境筛选后,幼苗的适应能力增强[50],当其高度 ≥ 100 cm时,渡过对环境反映最敏感的时期[17],死亡率下降,个体数量增多,这在一定程度上证实了负密度制约效应出现在植物生活史早期阶段[51]。
3.3 幼苗的分布
不同树种幼苗出现的样方数表现出一定差异。张建等[6]研究表明,种子良好的扩散能力可扩大种子和幼苗与母树的距离。Harms等[52]认为幼苗间出现样方数的差异可能是由于扩散限制的存在,种子不能扩散到适合其生长的地区,限制其种群数量和物种多度。如榆树的种子主要通过风力传播,西南卫矛的种子主要通过重力传播,扩散能力有限。同时,同种成年大树和其他物种的成年大树的分布[53]以及动物对种子和幼苗的取食[54]都影响幼苗的分布。
地形的变化会导致环境因子如光照、水分、温度、土壤的养分含量等发生空间变异[55]。研究发现,茂兰喀斯特常绿落叶阔叶混交林木本植物幼苗个体数和物种数都表现为随着海拔、坡度的升高呈现增加的趋势。结果与李洁琼等[4]对玉龙雪山样地研究得到幼苗物种数量随海拔的升高逐渐减少,以及黄萍等[56]对北京山油松林相同林型的研究得到更新数量与坡度呈负相关关系不同。张忠华等[57]对茂兰喀斯特常绿落叶阔叶混交林树种的空间分布格局及其分形特征研究结果表明,林内乔木树种的种群空间分布格局多呈现聚集分布,有些树种集中分布于干旱瘠薄、海拔较高的生境中。本研究得到的幼苗在高海拔上分布较多的结果可能是林内一些成年大树树种聚集分布在海拔较高的地段,从而使林下幼苗较多。草本是比木本更进化的性状,对环境的适应能力更强[58]。在坡度较缓处,小生境条件较坡度大的地段好,因此会存在较多的草本植物,其会抑制幼苗高度生长[59],木本植物幼苗竞争不过草本植物,从而导致幼苗难以存活。
2017—2019年幼苗密度在空间上呈现复杂的斑块分布,且不同年份间的空间异质性程度存在差异,这与研究区地形的变化有关。研究区属于喀斯特地貌,地形较复杂,土层浅薄不连续,岩石裸露率高[60],部分区域无法为幼苗萌发提供一个适宜的生存环境。同时,幼苗更新过程中存在负密度制约效应[49],使幼苗分布不均。林隙也会影响幼苗的分布[61],它为新个体提供存活与生长空间[62]。研究表明林隙下的幼苗密度一般高于非林隙[63]。另外,幼苗的生长还受其他因素的制约,如幼苗与成树之间的关系[36],种子产量[64]、腐殖质和枯落物[65]等,有待深入探索。
4. 结 论
本研究以茂兰喀斯特常绿落叶阔叶混交林为研究对象,基于2017—2019年连续3次的幼苗调查数据,对茂兰喀斯特常绿落叶阔叶混交林木本植物幼苗的组成、分布及其年际动态变化进行分析,得出以下结论:茂兰喀斯特常绿落叶阔叶混交林木本植物幼苗的数量和物种组成在年际之间都表现出明显的波动,木质藤本幼苗是林下幼苗的重要组成部分。幼苗新生和死亡的数量在年际间变化较大,林下幼苗更新存在一定的负密度制约效应。幼苗个体数和物种数随着海拔和坡度的增大呈现增加的趋势,且幼苗密度在空间上呈现复杂的斑块分布。
-
图 1 刈割频次和留茬高度对构树发枝数的影响
大写字母表示同一刈割频次不同留茬高度之间差异显著,小写字母表示同一留茬高度不同频次之间差异显著(P < 0.05)。下同。Capital letters indicate significant difference among varied stubble heights under the same mowing frequency. Lowercase letters indicate significant difference among varied mowing frequencies under the same stubble height (P < 0.05). The same below.
Figure 1. Effects of mowing frequency and stubble height on the resprouting number of Broussonetia papyrifera
图 4 刈割频次和高度对粗蛋白和粗脂肪含量的影响
大写字母表示同一刈割频次不同留茬高度之间差异显著,小写字母表示同一留茬高度不同频次之间差异显著(P < 0.05)。Capital letters indicate significant difference among varied stubble heights under the same mowing frequency. Lowercase letters indicate significant difference among varied mowing frequencies under the same stubble height (P < 0.05).
Figure 4. Influence of mowing frequency and stubble height on crude protein and crude fat contents
表 1 不同频次刈割时间表
Table 1 Mowing time of different cutting frequencies
刈割次数
Mowing time刈割日期
Mowing date间隔期
Interval time/d3 2017−05−22,2017−08−07,2017−10−22 50 4 2017−05−22,2017−07−12,2017−09−02,2017−10−22 42 5 2017−05−22,2017−06−29,2017−08−07,2017−09−14,2017−10−22 38 表 2 不同刈割频次下各茬次生物量 t/hm2
Table 2 Biomass of each stubble under different mowing frequencies
t/ha 刈割次数
Mowing time第1茬
1st stubble第2茬
2nd stubble第3茬
3rd stubble第4茬
4th stubble第5茬
5th stubble合计
Total3 0.60 ± 0.08A 0.75 ± 0.24A 0.41 ± 0.09A 1.75 ± 0.32A 4 0.50 ± 0.04B 0.54 ± 0.09B 0.64 ± 0.10B 0.18 ± 0.4A 1.86 ± 0.37A 5 0.60 ± 0.11C 0.33 ± 0.06AB 0.30 ± 0.08AB 0.46 ± 0.16BC 0.06 ± 0.01A 1.79 ± 0.19A 注:不同大写字母表示同一刈割频次下各茬次之间差异显著(P < 0.05)。Note: different capital letters indicate significant difference among stubbles under the same mowing frequency (P < 0.05). 表 3 不同刈割频次下各茬次粗蛋白和粗脂肪含量
Table 3 Crude protein and crude fat contents of each stubble under different mowing frequencies
养分含量
Nutrient content/%刈割次数
Mowing time第1茬
1st stubble第2茬
2nd stubble第3茬
3rd stubble第4茬
4th stubble第5茬
5th stubble粗蛋白
Crude protein3 24.40 ± 0.44B 23.65 ± 0.94B 19.98 ± 0.65A 4 24.69 ± 0.42C 19.12 ± 0.28A 25.47 ± 0.22C 22.70 ± 0.48B 5 24.31 ± 0.48BC 25.31 ± 0.42BC 26.24 ± 0.71C 21.10 ± 0.98A 23.41 ± 0.64B 粗脂肪
Crude fat3 11.75 ± 0.33B 8.01 ± 1.14A 9.53 ± 1.70AB 4 12.07 ± 1.53B 8.41 ± 0.82A 11.47 ± 0.78AB 10.19 ± 0.83AB 5 12.59 ± 0.93AB 11.54 ± 0.83AB 10.42 ± 0.98A 12.52 ± 0.88AB 13.89 ± 0.60B 注:不同大写字母表示同一刈割频次下各茬次之间差异显著(P < 0.05)。Note: different capital letters indicate significant difference among varied stubbles under the same mowing frequency (P < 0.05). -
[1] 杨凯. 四倍体刺槐刈割利用技术研究[D]. 杨凌: 西北农业科技大学, 2015. Yang K. Study on cutting and utilization techniques of tetraploid Robinia pseudoacacia[D]. Yangling: Northwest Agricultural University of Science and Technology, 2015.
[2] 张国君, 李云, 徐兆翮, 等. 栽培模式对四倍体刺槐生物量和叶片营养影响[J]. 北京林业大学学报, 2010, 32(5): 102−106. Zhang G J, Li Y, Xu Z H, et al. Effects of cultivation models on biomass and leaf nutritive value of tetraploid Robinia pseudoacacia[J]. Journal of Beijing Forestry University, 2010, 32(5): 102−106.
[3] Lin W, Wang J, Xu C, et al. Effects of Epichloë endophyte and repeated cutting on nutrition compositions of Festuca sinensis[J]. Plant, Soil and Environment, 2020, 66: 526−532. doi: 10.17221/444/2020-PSE
[4] 张德强, 周国逸, 温达志, 等. 刈割时间和次数对牧草产量和品种的影响[J]. 热带亚热带植物学报, 2000(增刊Ⅲ): 43−51. Zhang D Q, Zhou G Y, Wen D Z, et al. Effect of cutting frequency on yield and quality of forage grasses[J]. Journal of Tropical and Subtropical Botany, 2000(Suppl.Ⅲ): 43−51.
[5] 郭霞, 颜淑云, 周志宇, 等. 刈割次数对紫穗槐生长和氮磷钾含量的影响[J]. 草业科学, 2011, 28(3): 431−435. doi: 10.3969/j.issn.1001-0629.2011.03.016 Guo X, Yan S Y, Zhou Z Y, et al. Effects of cutting frequency on growth and NPK content of Amorpha fruticose[J]. Pratacultural Science, 2011, 28(3): 431−435. doi: 10.3969/j.issn.1001-0629.2011.03.016
[6] 常玉卉, 赵威军, 张阳, 等. 不同刈割次数对甜高粱产草量的影响[J]. 山西农业科学, 2017, 45(5): 753−755,768. doi: 10.3969/j.issn.1002-2481.2017.05.21 Chang Y H, Zhao W J, Zhang Y, et al. Effect of different cutting times on the grass yield of sweet sorghum[J]. Journal of Shanxi Agricultural Sciences, 2017, 45(5): 753−755,768. doi: 10.3969/j.issn.1002-2481.2017.05.21
[7] 游永亮, 赵海明, 李源, 等. 刈割制度对海河平原区紫花苜蓿产量和品质的影响[J]. 中国草地学报, 2018, 40(6): 47−55. You Y L, Zhao H M, Li Y, et al. Effect of cutting system on the forage yield and quality of alfalfa in Haihe Plain area[J]. Chinese Journal of Grassland, 2018, 40(6): 47−55.
[8] 杨秀芳, 梁庆伟, 娜日苏, 等. 年刈割次数对科尔沁沙地不同秋眠级紫花苜蓿品种产量、品质和越冬率的影响[J]. 草地学报, 2019, 27(3): 637−643. Yang X F, Liang Q W, Narisu, et al. The effects of cutting times on yield, quality and winter survival rate of the different fall dormancy rating alfalfa in Keerqin Sandy Land[J]. Acta Agrestia Sinica, 2019, 27(3): 637−643.
[9] Giambalvo D, Amato G, Stringi L. Effects of stubble height and cutting frequency on regrowth of berseem clover in a Mediterranean semiarid environment[J]. Crop Science, 2011, 51(4): 1808−1814. doi: 10.2135/cropsci2010.05.0271
[10] 李新媛, 王汝富, 王国生, 等. 留茬高度对人工草地苜蓿产量和质量的影响[J]. 草地与牧草, 2020, 40(6): 43−45, 54. Li X Y, Wang R F, Wang G S, et al. Effect of stubble height on yield and quality of alfalfa artificial grassland[J]. Grassland and Pasture, 2020, 40(6): 43−45, 54.
[11] 刘国辉. 刈割高度对高寒沼泽湿地植被再生影响的研究[J]. 青海草业, 2016, 25(4): 2−5. doi: 10.3969/j.issn.1008-1445.2016.04.001 Liu G H. Effects of mowing height on vegetation regrowth in alpine mire wetland[J]. Qinghai Prataculture, 2016, 25(4): 2−5. doi: 10.3969/j.issn.1008-1445.2016.04.001
[12] Mukangango M, Nduwamungu J, Naramabuye F X, et al. Biomass production and nutrient content of three agroforestry tree species growing on an acid anthropic ferralsol under recurrent harvesting at different cutting heights[J]. Agroforestry Systems, 2020, 94(3): 857−867. doi: 10.1007/s10457-019-00455-8
[13] Zheng Y X, Zhang Y P, Wu J C. Yield and quality of Moringa oleifera under different planting densities and cutting heights in southwest China[J]. Industrial Crops and Products, 2016, 91: 88−96. doi: 10.1016/j.indcrop.2016.06.032
[14] 付志慧, 格根图, 贾玉山, 等. 茬次与品种对多花黑麦草在北方地区种植的产量及营养价值的影响[J]. 草地学报, 2021, 29(7): 1571−1577. Fu Z H, Ge G T, Jia Y S, et al. Effects of cutting times and varieties on the yield and nutrient value of Italian ryegrass planted in northern regions[J]. Acta Agrestia Sinica, 2021, 29(7): 1571−1577.
[15] 包乌云, 赵萌莉, 红梅, 等. 刈割对人工草地产量和补偿性生长的影响[J]. 中国草地学报, 2015, 37(5): 46−51. doi: 10.3969/j.issn.1673-5021.2015.05.008 Bao W Y, Zhao M L, Hong M, et al. Effect of cutting on yield and compensatory growth of artificial pasture[J]. Chinese Journal of Grassland, 2015, 37(5): 46−51. doi: 10.3969/j.issn.1673-5021.2015.05.008
[16] 巩皓, 杨柳, 李丹丹, 等. 寒地黑土农区紫花苜蓿生产与品质对施肥和刈割频次的响应及效益分析[J]. 中国农业科学, 2020, 53(13): 2657−2667. doi: 10.3864/j.issn.0578-1752.2020.13.014 Gong H, Yang L, Li D D, et al. Response of alfalfa production and quality to fertilization and cutting frequency and benefit analysis in mollisol agricultural area in cold region[J]. Scientia Agricultura Sinica, 2020, 53(13): 2657−2667. doi: 10.3864/j.issn.0578-1752.2020.13.014
[17] Valencia J A G, Rivera R C, Aparicio Y V, et al. Alfalfa (Medicago sativa L.) biomass yield at different pasture ages and cutting frequencies[J]. Revista Mexicana de Ciencias Pecuarias, 2019, 10(2): 353−366. doi: 10.22319/rmcp.v10i2.4319
[18] 郭秀卿, 崔福柱, 杜天庆, 等. 4种甜秆饲草高粱对密度与刈割次数的响应[J]. 山西农业大学(自然科学版), 2015, 35(6): 596−602. Guo X Q, Cui F Z, Du T Q, et al. Response of 4 sweet forage sorghum for different density and mowing times[J]. Shanxi Agriculture University (Natural Science Edition), 2015, 35(6): 596−602.
[19] Partey S T. Effect of pruning frequency and pruning height on the biomass production of Tithonia diversifolia (Hemsl) A. Gray[J]. Agroforestry Systems, 2011, 83(2): 181−187. doi: 10.1007/s10457-010-9367-y
[20] Raj R M, Raj A K, Kunhamu T K. Fodder yield and nutritive value of mulberry (Morus indica L.) under varying plant density and pruning frequency in coconut garden[J]. Range Management and Agroforestry, 2019, 40(2): 255−261.
[21] Isiam M R, Carcia S C, Horadagoda A, et al. Management strategies for forage rape (Brassica napus L. cv. Goliath): impact on dry-matter yield, plant reserves, morphology and nutritive value[J]. Grass and Forage Science, 2020, 75(1): 96−110. doi: 10.1111/gfs.12462
[22] de Sousa G J, Alexandrino E, dos Santos A C, et al. Megathyrsus maximus cv. Massai at different cutting frequencies[J]. Semina-Ciencias Agrarias, 2019, 40(5): 1913−1923. doi: 10.5433/1679-0359.2019v40n5p1913
[23] 杨晓松, 高丽娟, 斯日古楞, 等. 科尔沁区天然打草场植被特征对刈割频次的响应[J]. 畜牧与饲料科学, 2016, 37(9): 32−34. doi: 10.3969/j.issn.1672-5190.2016.09.010 Yang X S, Gao L J, Siriguleng, et al. The responses of vegetation characteristic of natural meadow on mowing frequency in Keerqin District[J]. Animal Husbandry and Feed Science, 2016, 37(9): 32−34. doi: 10.3969/j.issn.1672-5190.2016.09.010
[24] 刘晓静, 张进霞, 李文卿, 等. 施肥及刈割对干旱地区紫花苜蓿产量和品质的影响[J]. 中国沙漠, 2014, 34(6): 1516−1526. Liu X J, Zhang J X, Li W Q, et al. Effects of nitrogen and phosphorus addition and cuttings on yield and quality of alfalfa in dry region of Gansu, China[J]. Journal of Desert Research, 2014, 34(6): 1516−1526.
[25] 闫士元, 李迎, 王晓宇, 等. 刈割次数对苜蓿相关表型性状与饲草产量的影响[J]. 草原与草业, 2021, 33(1): 6−11. doi: 10.3969/j.issn.2095-5952.2021.01.002 Yan S Y, Li Y, Wang X Y, et al. Effects of cutting times on phenotypic traits and forage yield of alfalfa[J]. Grassland and Prataculture, 2021, 33(1): 6−11. doi: 10.3969/j.issn.2095-5952.2021.01.002
[26] 张峰, 赵天启, 乔荠瑢, 等. 刈割留茬高度对大针茅草原生产力及可持续利用的影响[J]. 草地学报, 2021, 7(29): 1491−1498. Zhang F, Zhao T Q, Qiao J R, et al. Effects of stubble height on productivity and sustainable utilization in the Stipa grandis steppe[J]. Acta Agrestia Sinica, 2021, 7(29): 1491−1498.
[27] Elena M, Andreas K, John C B, et al. Influence of cutting height on biomass yield and quality of miscanthus genotypes[J]. Global Change Biology Bioenergy, 2021, 13(10): 1675−1689. doi: 10.1111/gcbb.12881
[28] Hu R, Zhang G, Li Y, et al. Effects of cultivation measures on biomass in tetraploid Robinia pseudoacacia[J]. Nonwood Forest Research, 2011, 29(4): 46−53.
[29] 胡建忠, 卢顺光. 我国主要野生沙棘叶片的饲用价值研究[J]. 现代畜牧兽医, 2021(9): 19−24. Hu J Z, Lu S G. Study on the feeding value of main wild Hipphphae leaves in China[J]. Modern Journal of Animal Husbandry and Veterinary Medicine, 2021(9): 19−24.
[30] 孙扣忠, 邹海洋, 李亚芳, 等. 刈割次数对江苏沿海地区多花黑麦草产量及品质的影响[J]. 安徽农业科学, 2020, 48(2): 106−108, 112. doi: 10.3969/j.issn.0517-6611.2020.02.028 Sun K Z, Zou H Y, Li Y F, et al. Effects of cutting times on the yield and quality of Lolium multiflorum in coastal areas of Jiangsu Province[J]. Journal of Anhui Agriculture Science, 2020, 48(2): 106−108, 112. doi: 10.3969/j.issn.0517-6611.2020.02.028
[31] 赵京东, 宋彦涛, 徐鑫磊, 等. 施氮和刈割对辽西北退化草地牧草产量和品质的影响[J]. 草业学报, 2021, 30(8): 36−48. doi: 10.11686/cyxb2020292 Zhao J D, Song Y T, Xu X L, et al. Effects of nitrogen application and mowing on yield and quality of forage in degraded grassland in northwest Liaoning Province[J]. Acta Prataculturae Sinica, 2021, 30(8): 36−48. doi: 10.11686/cyxb2020292
[32] Inyang U, Vendramini J M B, Sellers B, et al. Harvest frequency and stubble height affect herbage accumulation, nutritive value, and persistence of ‘Mulato II’ brachiariagrass[J]. Forage and Grazinglands, 2010, 8(1): 7.