2. 四川省烟草公司凉山州公司, 四川西昌 615000;
3. 四川省烟草公司攀枝花市公司, 四川攀枝花 617000
氮素是烤烟生长发育过程中重要的营养元素,尤其是后期土壤氮含量和矿化量对烟叶的品质有直接的显著影响[1-2];当后期土壤供氮较多时会增加烤烟上部叶的烟碱含量[3]。此外,烤烟是喜硝态氮的作物,土壤硝态氮的含量对烤烟的生长发育有间接的影响[4-5]。烤烟生长过程中所需的氮素主要来源于当季施入的肥料氮和土壤有机氮的矿化,其中烤烟累积的氮素主要来源于土壤氮,占总量的50%[6]。谷海红等[7]研究表明,无论是水稻土还是红壤,烤烟生长后期土壤的供氮能力是生产优质烟叶的关键。因此,植烟土壤的供氮特性和氮素矿化能力对优质烟叶的生长有重要作用[8]。
土壤水分条件是土壤氮矿化的关键因素之一,湿润的土壤环境会促进土壤氮矿化,但含水量过高或过低会降低土壤氮矿化[9]。烤烟生长期间高温高湿的气候条件,会影响微生物数量和活性,进而影响植烟土壤的氮矿化程度[10],最终影响优质烤烟的生长发育状况。高真真等[8]对典型香型植烟土壤氮矿化特征研究表明,不同烟区土壤的含水量对土壤氮矿化影响各异,50% 田间持水量处理有利于质地相对黏重的土壤氮素矿化,而减少了质地相对较轻的土壤氮素矿化量和矿化速率。Stanford和Epstein[11]研究表明土壤含水量为饱和含水量50% ~ 60% 时,有利于土壤氮的硝化作用。赵长盛等[12]研究指出,土壤矿化速率在土壤含水量为25% 时达到最大值,且土壤含水量与土壤氨化速率呈正相关关系,而与硝化速率呈负相关关系。土壤酶活性是土壤肥力重要指标之一,其活性大小会影响土壤氮矿化[13]。不同水分条件和不同土壤类型会造成土壤酶活性的变化,最终影响土壤氮矿化。但目前对植烟土壤氮矿化及酶活性对土壤水分条件的响应特征研究较少,凉攀烟区(凉山州和攀枝花)是四川地区主要烤烟生产地,区域内地形多样,立体气候特征明显,雨热同期,干湿分明。本研究以凉攀烟区的土壤为研究对象,采用室内培养法,研究了不同水分条件下红壤、紫色土氮矿化特征及其相关的酶活性变化,以期为凉攀烟区合理调控土壤氮素和生产优质烤烟提供理论依据。
1 材料与方法 1.1 研究区概况研究区位于四川省凉山州和攀枝花市的主要烤烟生产地。凉山州(100°03′ ~ 103°52′E,26°03′ ~ 29°18′N)地处四川省西南部横断山区的东北部,青藏高原南缘,气候类型为亚热带季风气候,干湿分明,冬半年日照充足,少雨干暖;夏半年云雨较多,降水量集中在4—9月,气候凉爽,年均气温6 ~ 17 ℃。攀枝花(101°08′ ~ 102°15′E,26°05′ ~ 27°2l′N)地处四川西南川滇结合部,是以南亚热带为基带的立体气候小岛,干雨季分明,昼夜温差大,日照长,降雨量集中在6—10月。凉山和攀枝花植烟土壤类型主要为红壤和紫色土。
1.2 土样采集和试验设计 1.2.1 土样采集在凉山州和攀枝花市的主要烤烟生产区,采集红壤和紫色土样品。分别在烟苗移栽前,采集0 ~ 20 cm耕层土壤,在阴凉处风干后去除可见的植物残体和根系等。每个样品采集3 kg,取部分用于土壤基本指标测定,其余土壤样品风干过筛(2 mm)用于培养试验,土壤样品基本理化性质见表 1。
取风干后的植烟土壤样品均匀地铺装于容积为500 ml的塑料瓶中,瓶口覆膜,留3个通气孔,防止水分蒸发和通气。为活化土壤,在培养试验前,加无氮水至土壤质量含水量的60%,在25 ℃恒温预培养一周。基于采样区烤烟生长期气候条件,模拟干湿两季环境,本试验共设置两种水分条件:60%田间持水量(60% WHC)和淹水状态(FS)培养。60%WHC培养下,整个培养试验过程中保持土壤含水量为田间持水量的60%;FS培养下,整个培养试验过程中始终保持土壤在淹水状态,土壤表面有1 ~ 2 cm淹水。预培养结束后,根据不同水分条件调整土壤含水量,土壤氮矿化培养试验在光照培养室内进行,温度保持在25 ℃,空气湿度为60%,每2天用重量法保持土壤含水量,添加散失水分。整个试验共培养64 d,分别在培养试验的第0、2、8、16、32、64天取土样,试验为不破坏性取样,每个处理设置3个重复,共计72个试验处理。
1.3 样品测定及数据处理 1.3.1 土壤无机氮培养试验采集的土壤样品,用于测定土壤NH4+-N和NO3–-N含量。土壤NH4+-N和NO3–-N采用2 mol/L KCl提取,取滤液分别用靛酚蓝比色法和紫外分光光度计法测定,具体测定详细步骤参见文献[14]。本研究中的土壤无机氮为土壤NH4+-N和NO3–-N含量之和。
$ 土壤氮矿化速率(\mathrm{mg}/(\mathrm{kg·d}))=[\mathrm{N}]_{T}/(ΔT) $ |
式中:[N]T为培养第T天时土壤NH4+-N或NO3–-N的含量(mg/kg);ΔT为两次采样之间的天数间隔(d)。
土壤矿化的无机氮累积过程采用一级动力学方程进行非线性拟合[15]:
$ N_{t} = N_{0}(1-\mathrm{e}^{-kt}) $ |
式中:Nt为第t天的土壤无机氮累积矿化量;N0为氮矿化势;k为氮矿化速率常数。
1.3.2 土壤酶活性脲酶活性采用苯酚钠-次氯酸钠比色法测定,以24 h后1 g土壤中NH3-N的毫克数(mg)表示;蛋白酶活性采用福林法测定,以2 h后1 g土壤中酪氨酸的微克数(μg)表示。亚硝酸还原酶测定采用碳酸钙-亚硝酸钠比色法,以24 h后土壤中NO2-N的浓度(mg/L)表示。具体方法详见参考文献[16]。
1.3.3 数据处理采用SPSS 20.0、Origin 8.5和Excel 2010软件进行数据处理、绘图制表。所有的结果均用3次测定结果的平均值表示。不同处理之间的多重比较采用LSD法(P < 0.05)。
2 结果与分析 2.1 土壤铵态氮和硝态氮含量红壤和紫色土在FS和60%WHC条件下土壤NH4+-N和NO3–-N含量随培养时间的变化存在显著差异。如图 1A所示,无论是红壤还是紫色土,FS培养的土壤NH4+-N含量均高于60%WHC。FS培养下,红壤和紫色土的NH4+-N含量均在第2天时达到最大值,分别为12.91和12.44 mg/kg;随着培养时间的增加,土壤NH4+-N含量均呈下降趋势,在培养第64天,紫色土NH4+-N含量较红壤明显增加了2.49 mg/kg (P < 0.05)。60%WHC培养下,土壤NH4+-N含量随着培养时间的增加呈先增加后减少并趋于稳定的变化趋势;在培养的前16天,紫色土NH4+-N含量较红壤显著提高了1.5% ~ 63.0%(P < 0.05)。总体上,FS培养的土壤NH4+-N含量明显高于60%WHC,且培养期内紫色土NH4+-N平均含量高于红壤。
培养期内的土壤NO3–-N含量变化与NH4+-N含量存在差异,如图 1B所示,60%WHC培养的土壤NO3–-N含量明显高于FS。FS培养下,各培养阶段的红壤和紫色土NO3–-N含量存在显著差异,除培养第2天外,红壤的NO3–-N含量较紫色土增加了0.89 ~ 10.52 mg/kg。60%WHC培养下,红壤的NO3–-N含量在培养第8天达到最高值(91.69 mg/kg),而紫色土在培养第32天达到最高值(80.38 mg/kg)。总体上,60%WHC培养的NO3–-N含量明显高于FS,且培养期内红壤NO3–-N平均含量高于紫色土。
2.2 土壤氨化速率和硝化速率如图 2所示,培养期内,土壤氨化和硝化速率均随培养时间呈先下降后平稳的变化趋势。对土壤NH4+-N而言(图 2A),FS培养的氨化速率(0.07 ~ 6.46 mg/(kg·d))明显高于60%WHC(0.02 ~ 2.03 mg/(kg·d));FS培养下,红壤氨化速率除在第8天和第64天低于紫色土外,其他培养时间内均高于紫色土;60%WHC培养下,除第16天外,紫色土氨化速率略高于红壤;总体上,整个培养时期内,紫色土平均氨化速率(0.52和1.75 mg/(kg·d))均高于红壤(0.50和1.73 mg/(kg·d))。对土壤NO3–-N而言(图 2B),FS培养的硝化速率均低于60%WHC,平均减少了27.2% ~ 79.8%;FS培养下,除第2天外,红壤硝化速率(1.59 ~ 7.13 mg/(kg·d))均高于紫色土(1.26 ~ 6.60 mg/(kg·d));60%WHC培养下,红壤硝化速率在培养前期(第16天之前)明显高于紫色土;总体上,60%WHC培养的土壤硝化速率高于FS,且红壤在60%WHC培养下的硝化速率明显高于紫色土。
如图 3所示,培养期内,60%WHC培养的土壤无机氮累积矿化量明显高于FS,其中红壤60%WHC培养下的土壤无机氮累积矿化量最大,为42.12 ~ 394.06 mg/kg,其次是同一水分条件培养下的紫色土。通过一级动力学方程进行非线性拟合(表 2)发现,土壤氮矿化符合动力学方程,且相关性均达到极显著水平。60%WHC培养的土壤氮矿化势(N0)高于FS,且红壤的N0较紫色土提高了4.8%(60%WHC)和7.4%(FS);紫色土FS培养的土壤矿化率(k)高于60%WHC,而红壤60%WHC培养下的k略高于FS;施用综合指数(N0×k)表现为60%WHC-红壤 > FS-紫色土 > FS-红壤 > 60%WHC-紫色土,其中红壤60%WHC培养的供氮指数N0×k较同一水分条件培养下的紫色土增加了12.3 mg/(kg·d)。
从表 3可以看出,FS和60%WHC培养下,红壤的脲酶活性均高于紫色土。FS培养下,红壤和紫色土的脲酶活性均在第32天时达到最高值,分别为NH3-N 1.09和0.82 mg/(g·24 h);60%WHC培养下,随着培养时间的增加,红壤和紫色土脲酶活性呈先增加后减少的趋势,仍在第32天时达到最高值;整个培养期内,红壤60%WHC培养的土壤脲酶活性平均值最高,为NH3-N 0.86 mg/(g·24 h)。对土壤蛋白酶而言,FS培养下,红壤蛋白酶活性在第16天前均高于紫色土,而培养结束后(第64天),其蛋白酶活性较紫色土明显降低了53.6%;60%WHC培养下,土壤蛋白酶活性随着培养时间的增加呈下降趋势,且红壤的蛋白酶活性高于紫色土,2种土壤的蛋白酶均在培养初期(第2天)表现出较高的活性。从表 3还可以看出,紫色土的亚硝酸还原酶的活性高于红壤,培养期内土壤亚硝酸还原酶活性平均值表现为:FS-紫色土 > FS-红壤 > 60%WHC-紫色土 > 60%WHC-红壤。综上,60%WHC培养的土壤脲酶活性高于FS,尤其在红壤中效果较佳;而60%WHC培养的土壤蛋白酶活性和亚硝酸还原酶活性则低于FS,尤其是在紫色土培养中。
通过相关性分析发现(表 4),紫色土NO3–-N含量与土壤蛋白酶、亚硝酸还原酶活性存在显著的负相关关系,而红壤NO3–-N含量与亚硝酸还原酶活性存在显著的负相关关系;紫色土氨化速率与脲酶活性显著负相关,红壤和紫色土硝化速率与脲酶活性显著负相关;除紫色土无机氮累积矿化量与蛋白酶活性不存在显著相关外,2种土壤的无机氮累积矿化量均与脲酶、蛋白酶、亚硝酸还原酶活性存在显著相关。不同水分条件培养下,60%WHC培养的土壤NO3–-N含量与亚硝酸还原酶活性之间存在显著负相关;无论FS还是60%WHC培养,土壤氨化速率均与脲酶活性存在显著相关,除FS培养的土壤无机氮累积矿化量与蛋白酶活性不存在显著相关外,2种水分条件培养下的土壤无机氮累积矿化量均与脲酶、蛋白酶、亚硝酸还原酶活性间存在显著相关。
本研究中,水分条件对2种植烟土壤NH4+-N、NO3–-N含量和无机氮矿化存在显著影响,这与高真真等[8]、王常慧等[17]研究结果一致。一般认为,50% ~ 60% 田间持水量条件下土壤硝化作用最为活跃,因为该土壤水分含量既能满足微生物所需的水分条件,又能满足氧分分压[18]。本试验培养期间,土壤NH4+-N和NO3–-N含量存在“此消彼长”的关系,即淹水(FS)培养的土壤NH4+-N含量明显高于60% 田间持水量(60%WHC)培养,而FS培养的土壤NO3–-N含量明显低于60%WHC,可见FS培养不利于土壤NO3–-N累积,而增加了NH4+-N累积。一方面是因为FS培养相较于60%WHC培养土壤通透性较差,氧分分压较低,抑制了土壤氨氧化作用细菌的活性,从而促进NH4+-N累积[19];另一方面,FS培养的条件为反硝化作用提供了厌氧环境,促进了土壤NO3–-N还原,从而减少土壤NO3–-N累积,但可能会导致更多的NO3–-N以气态形式或淋溶损失,减少土壤氮的利用率[19-21]。
不同土壤类型的土壤氮矿化特征各异,本研究中,培养期内红壤NH4+-N平均含量和平均矿化速率均低于紫色土,而其NO3–-N平均含量和平均矿化速率、无机氮累积矿化量均高于紫色土,可见四川凉攀烟区红壤的矿化能力强于紫色土,这可能是因为不同类型的土壤性质,如有机质含量,导致了土壤的氮矿化潜力的差异[8, 22-23]。马兴华等[23]研究表明土壤有机氮的矿化率与土壤有机质含量显著相关,有机质含量高有利于土壤氮素矿化,提高了土壤的供氮能力。本试验条件下,红壤有机质含量(23.14 g/kg)高于紫色土(16.54 g/kg),有机氮源丰富,提高了微生物活性,因此在同一培养条件下红壤的矿化能力较强。但王树会和刘青丽[22]对云南主要植烟土壤的矿化特征研究表明,紫色土表层矿化能力和矿化量均高于红壤,这可能是因为云南主要烟区的紫色土潜在肥力高,阳离子交换量大,透气透水能力良好;此外该培养条件为间歇淋洗好气培养法,为微生物的生长发育提供了良好的微生态环境,加速了土壤矿化。土壤氮矿化势(N0)和矿化率(k)可以分别用来反映土壤的供氮容量和供氮速度;且施用综合指标(N0×k)可以更好地反映土壤真实的供氮情况[15]。一级动力学进行非线性拟合结果表明(表 2),无论FS或60%WHC培养下红壤矿化势均高于紫色土,而红壤的氮矿化率和供氮指标仅在60%WHC培养下高于紫色土,这也进一步说明了该研究区的红壤氮矿化潜力要强于紫色土,尤其是在60%WHC条件下,而含水量过高不利于土壤氮素矿化。
凉攀地区降水主要集中的4—9月,尤其是6—8月,如果烤烟生长后期仍有大量的氮素供应,会促进烟叶合成更多的烟碱,影响烤烟品质。本研究中,结合凉攀地区的降水特征,在降雨较少,土壤含水量较低的情况下,具有较强氮矿化潜力和较高硝化速率的红壤在烤烟前期可以释放大量的氮素供烟株吸收,但可能会引起烤烟后期仍有大量的氮素吸收,增加烟叶烟碱合成,导致烟碱含量过高。在降水量较多,土壤含水量过高的情况下,相较于红壤而言,氨化速率较高的紫色土可能在烤烟中期会提供较多的NH4+-N,由于NH4+-N可以促进芳香族挥发油,增加烤烟香味,但是过高的NH4+-N不利于烟株生长后期碳氮代谢的转化,最终影响烟叶品质[21]。此外,在降水量相对较少的烤烟后期,紫色土的供氮水平可能低于红壤,更有利于烟叶品质。因此,烤烟实际生产中,在推荐氮肥施用量时,不但要考虑不同土壤类型,还应该注意土壤含水量的调控,防止烤烟生长后期氮素供应能力过强,不利于烟叶成熟落黄和香气物质形成,降低了烟叶品质[24]。
土壤脲酶是土壤氮循环过程中重要的一种酰胺酶,可催化尿素生成易被植物吸收利用的无机氮和二氧化碳,其活性在一定程度上可以反映土壤供氮能力[25]。本研究中,紫色土的脲酶活性低于红壤,这可能与紫色土的有机质含量较低有关。有机质可以为微生物提供碳源,促进土壤微生物繁殖,促进土壤脲酶的分泌[13],从而增加了土壤无机氮的含量,因此无论是60%WHC还是FS培养,土壤脲酶活性均与土壤无机氮矿化累积量存在显著的正相关关系(r = 0.599 ~ 0.812,表 4)。土壤蛋白酶是一种胞外酶,可将蛋白质和肽类等分解形成氨基酸[26],FS培养的土壤蛋白酶活性高于60%WHC,可见水分条件是影响土壤蛋白酶活性的重要因素之一,尤其是在60%WHC培养条件下土壤蛋白酶活性显著影响了土壤无机氮矿化累积量。亚硝酸还原酶是反硝化作用中的关键酶,催化NO3–-N形成NO气体,其活性不仅受到自身成分的影响,还受到外界环境的影响,例如通气情况、水分等均会影响到酶活性[27]。淹水培养较60%WHC为反硝化作用提供了厌氧环境,提高了土壤反硝化微生物的活性,进而增加了土壤亚硝酸还原酶的活性。红壤和紫色土NO3–-N含量均与亚硝酸还原酶存在显著的负相关关系,因此,较高亚硝酸还原酶活性的紫色土催化了NO3–-N形成NO气体,从而减少了紫色土NO3–-N含量。
4 结论1) 淹水培养提高了土壤NH4+-N含量和氨化速率,但降低了土壤NO3–-N含量和硝化速率,且60% 田间持水量培养下的红壤硝化速率高于紫色土;该烟区的红壤氮矿化潜力要强于紫色土,尤其是在60% 田间持水量培养下,而含水量过高不利于土壤氮素矿化。
2) 不同土壤类型在不同水分条件下土壤酶活性各异。红壤的脲酶活性强于紫色土,尤其是在60% 田间持水量培养的效果更为显著;紫色土的蛋白酶活性和亚硝酸还原酶活性高于红壤,尤其是在淹水培养下的效果显著。土壤脲酶、蛋白酶和亚硝酸还原酶的活性能显著影响土壤无机氮矿化累积量。
[1] |
刘青丽, 任天志, 李志宏, 等. 植烟黄壤氮素矿化动态模拟研究[J]. 植物营养与肥料学报, 2010, 16(2): 400-406 (0) |
[2] |
李春俭, 张福锁, 李文卿, 等. 我国烤烟生产中的氮素管理及其与烟叶品质的关系[J]. 植物营养与肥料学报, 2007, 13(2): 331-337 DOI:10.3321/j.issn:1008-505X.2007.02.024 (0) |
[3] |
刘卫群, 郭群召, 张福锁, 晁逢春. 氮素在土壤中的转化及其对烤烟上部叶烟碱含量的影响[J]. 烟草科技, 2004, 37(5): 36-39 DOI:10.3969/j.issn.1002-0861.2004.05.014 (0) |
[4] |
孙榅淑, 杨军杰, 周骏, 等. 不同氮素形态对烟草硝态氮含量和TSNA形成的影响[J]. 中国烟草学报, 2015, 21(4): 78-84 (0) |
[5] |
李孝刚, 彭曙光, 靳志丽, 等. 有机物料对植烟土壤氮素矿化及微生物性质的影响[J]. 土壤学报, 2021, 58(1): 225-234 (0) |
[6] |
马兴华, 张忠锋, 荣凡番, 等. 高低土壤肥力条件下烤烟对氮素吸收、分配和利用的研究[J]. 中国烟草科学, 2009, 30(1): 1–4, 9 DOI:10.3969/j.issn.1007-5119.2009.01.001 (0) |
[7] |
谷海红, 李志宏, 李天福, 等. 不同来源氮素在烤烟体内的累积分配及对烟叶品质的影响[J]. 植物营养与肥料学报, 2009, 15(1): 183-190 DOI:10.3321/j.issn:1008-505X.2009.01.027 (0) |
[8] |
高真真, 段卫东, 胡坤, 等. 温度和水分对典型香型烟区植烟土壤氮素矿化的影响[J]. 土壤, 2019, 51(3): 442-450 (0) |
[9] |
于淑华, 张丽霞, 谢雪迎, 等. 不同水分模式对山东茶园土壤氮素动态的影响[J]. 水土保持学报, 2021, 35(4): 289-298 (0) |
[10] |
袁磊, 李文周, 陈文伟, 等. 戴云山自然保护区森林土壤氮转化特点研究[J]. 土壤, 2017, 49(2): 240-247 (0) |
[11] |
Stanford G, Epstein E. Nitrogen mineralization-water relations in soils[J]. Soil Science Society of America Journal, 1974, 38(1): 103-107 DOI:10.2136/sssaj1974.03615995003800010032x (0) |
[12] |
赵长盛, 胡承孝, 孙学成, 等. 温度和水分对华中地区菜地土壤氮素矿化的影响[J]. 中国生态农业学报, 2012, 20(7): 861-866 (0) |
[13] |
张璐, 任天宝, 阎海涛, 等. 不同有机物料对烤烟根际土壤碳库、酶活性及根系活力的影响[J]. 中国烟草科学, 2018, 39(2): 39-45 (0) |
[14] |
鲍士旦. 土壤农化分析(3版)[M].
中国农业出版社, 北京, 2000
(0) |
[15] |
潘金华, 吕成文, 吴红宝, 等. 无机改良剂对皖南植烟红壤氮矿化的影响[J]. 水土保持学报, 2016, 30(2): 193–199, 207 (0) |
[16] |
关松荫. 土壤酶及其研究法[M].
农业出版社, 北京, 1986
(0) |
[17] |
王常慧, 邢雪荣, 韩兴国. 温度和湿度对我国内蒙古羊草草原土壤净氮矿化的影响[J]. 生态学报, 2004, 24(11): 2472-2476 (0) |
[18] |
刘若萱, 张丽梅, 白刃, 等. 模拟条件下土壤硝化作用及硝化微生物对不同水分梯度的响应[J]. 土壤学报, 2015, 52(2): 415-422 (0) |
[19] |
Huang R, Wang Y Y, Liu J, et al. Variation in N2O emission and N2O related microbial functional genes in straw- and biochar-amended and non-amended soils[J]. Applied Soil Ecology, 2019, 137: 57-68 (0) |
[20] |
Lognoul M, Theodorakopoulos N, Hiel M P, et al. Impact of tillage on greenhouse gas emissions by an agricultural crop and dynamics of N2O fluxes: Insights from automated closed chamber measurements[J]. Soil and Tillage Research, 2017, 167: 80-89 (0) |
[21] |
谷海红. 水旱轮作植烟土壤供氮特征及对烤烟氮素积累分配的影响[D]. 北京: 中国农业科学院, 2008.
(0) |
[22] |
王树会, 刘青丽. 云南主要植烟土壤不同土层氮矿化研究[J]. 中国土壤与肥料, 2013(1): 14-19 (0) |
[23] |
马兴华, 荣凡番, 苑举民, 等. 云南曲靖烟区典型植烟土壤氮素矿化规律研究[J]. 中国烟草学报, 2011, 17(6): 66-69 (0) |
[24] |
龙世平, 李宏光, 曾维爱, 等. 湖南省主要植烟区域土壤有机氮矿化特性研究[J]. 中国烟草科学, 2013, 34(3): 6-9 (0) |
[25] |
龚丝雨, 聂亚平, 张启明, 等. 增施生物炭对烤烟成熟期根际土壤酶活性的影响[J]. 江西农业学报, 2017, 29(10): 54-57 (0) |
[26] |
李琰琰, 刘国顺, 向金友, 等. 氮营养水平对植烟土壤养分含量及土壤酶活性的影响[J]. 中国烟草学报, 2011, 17(2): 39-43 (0) |
[27] |
边雪廉, 赵文磊, 岳中辉, 等. 土壤酶在农业生态系统碳、氮循环中的作用研究进展[J]. 中国农学通报, 2016, 32(4): 171-178 (0) |
2. Liangshan Branch of Sichuan Tobacco Company, Xichang, Sichuan 615000, China;
3. Panzhihua Branch of Sichuan Tobacco Company, Panzhihua, Sichuan 617000, China