Effects of fresh seaweed bio-enzymatic organic fertilizer on the growth and soil physicochemical properties of Atraphaxis bracteata
-
摘要:目的为提高荒漠沙地植被的存活率和覆盖度，本文以毛乌素沙地西南缘的沙木蓼为研究对象，通过施用鲜海带生物酶解有机液肥，研究其对沙地土壤理化性质和沙木蓼生长及群落物种多样性的影响。方法采用不同梯度的鲜海带生物酶解有机液肥对沙木蓼进行1年的试验，利用土壤水分测定仪测量土壤含水率，凯氏定氮法测定土壤全氮含量，重铬酸钾滴定法测定土壤有机质含量，公式法计算群落物种多样性。结果(1) 该有机液肥有效地增加了土壤0~10 cm层的含水率，且施用冲施肥研究区的土壤含水率较叶面肥研究区的高10.6%。(2)土壤0~20 cm层的全氮和有机质的含量随着施肥量的增加而增加，并在施肥量为15 mL/m2时，达到最高，与对照间差异显著(P < 0.05)，且叶面肥处理的有机质含量远高于冲施肥。(3)该有机液肥有效地促进了沙木蓼的生长，沙木蓼地径、株高和新梢增长量均随着施肥量增加而呈先增加后减小的趋势，在施肥量为6 mL/m2时达到最大，与对照差异显著(P < 0.05)，且冲施肥的效果要比叶面肥更优。(4)该有机液肥增加了沙木蓼群落的物种组成、Margalef物种丰富度指数、Simpson优势度指数和Shannon-Wiener多样性指数。结论鲜海带生物酶解有机液肥能够增加毛乌素沙地土壤水分与养分的含量，促进沙木蓼的生长，有助于毛乌素沙地植被的恢复，更好地发挥其生态效益。Abstract:ObjectiveIn order to improve the survival rate and coverage of desert vegetation, a field experiment by adding fresh seaweed bio-enzymatic liquid organic fertilizer on the Atraphaxis bracteata located in the southwest of Mu Us Desert, Ningxia of northwestern China, was carried out to study its effect on soil physicochemical properties, growth and community species diversity of Atraphaxis bracteata in desert ecosystems.MethodUsing a different gradient of fresh seaweed bio-enzymatic organic fertilizer to conduct a one-year test on Atraphaxis bracteata, the soil moisture was measured by soil moisture meter, soil total nitrogen was measured by Kjeldahl method, soil organic matter was measured by potassium dichromate titration, and species diversity was calculated by formula method.Result(1) Fresh seaweed biological enzyme liquid organic fertilizer can increase soil moisture content in 0-10 cm layer effectively, and the soil moisture content of flush fertilization was 10.6% higher than that of foliar fertilizer application area. (2) The content of soil total nitrogen and organic matter in 0-20 cm layer increased with the increase of fertilizing amount and reached peak value when the fertilizing amount was 15 mL/m2, and there was significant difference compared with control (P < 0.05). The content of organic matter treated with foliar fertilizer was much higher than that of flush fertilizer. (3) The fresh liquid organic fertilizer promoted the growth of Atraphaxis bracteata effectively, the growth of its ground diameter, plant height and new shoot all increased first and then decreased with the increase of fertilizing amount, and reached peak value when the fertilizing amount was 6 mL/m2, and there was significant difference compared with control (P < 0.05). The effect of flush fertilization was better than foliar fertilizer. (4) The fresh seaweed bio-enzymatic liquid organic fertilizer increased the species composition, Margalef species richness index, Simpson dominance index and Shannon-Wiener diversity index of Atraphaxis bracteata.ConclusionFresh seaweed bio-enzymatic liquid organic fertilizer can increase the content of soil moisture and nutrients in Mu Us Desert, promote the growth of Atraphaxis bracteata, contributes to the restoration of vegetation in Mu Us Desert, and then the vegetation can give better play to the ecological benefits.
-
松栎柱锈菌(Cronartium orientale)是导致樟子松(Pinus sylvestris var. mongolica)、兴凯湖松(Pinus takahasii)和马尾松(Pinus massoniana)等2针松瘤锈病的森林有害生物，是一种转主寄生菌，其春孢子阶段寄生在松树枝干部[1-2]。当松树被松栎柱锈菌侵染致病后，其病原菌就会长期寄生在松树枝干部，并且病部肿瘤会逐年增大，给林分健康和林木生长造成极大危害[3]。黄酮类化合物具有较高的生物活性，是一种强抗氧化剂，能有效清除自由基从而延缓衰老和阻止细胞衰老、癌变，对人类心脑血管疾病也有很好的治疗和预防效果，具有较高的医疗和保健作用[4-8]。已有较多研究成果表明松科植物树皮与松针中黄酮类化合物含量丰富，具有较好的开发利用潜力[9-10]。国内学者如李石清等对南方红豆杉(Taxus chinensis var. mairei)枝叶总黄酮的研究、冯慧英等对巴山冷杉(Abies fargesii)总黄酮的研究、王飞等对青扦(Picea wilsonii)针叶总黄酮的研究等均表明其易提取且具有较强的抗氧化能力[11-13]。
担子菌门(Basidiomycota)的许多菌类代谢产物已经被开发利用，关于担子菌的研究已经成为当代生命科学、食品科学和医药领域的热点[14-15]。锈菌目(Uredinales)隶属于担子菌门，目前却尚未发现任何关于锈菌目真菌开发利用的研究，这也表明我们对于锈菌目真菌的潜在价值知之甚少。本研究选择了锈菌目中寄生在松科植物上的易于采集的松栎柱锈菌春孢子为研究对象，进行了松栎柱锈菌春孢子黄酮类化合物(以下简称CQF)的提取与体外抗氧化活性的研究，旨在为锈菌目菌物开发利用、森林有害生物开发利用和提升与拓宽森林有害生物防治理念提供科学依据。
2. 材料与方法
2.1 研究材料与预处理
松栎柱锈菌春孢子于2015年5月中旬采自大兴安岭地区，过80目筛后于45℃烘干至恒质量，密封冷藏。所用试剂盒均购自南京建成生物研究所；芦丁标准品，购自上海源叶生物科技有限公司；其他试剂均为国产分析纯；KQ-600B型超声波仪，昆山市超声波仪器有限公司；分析天平(感应量0.0001g)；Gene Quant 1300型分光光度计，美国GE公司。
2.2 测定波长的选择与标准曲线的建立
精密称取芦丁标准品10.0mg，使用体积分数为75%的乙醇溶液配制成50mL 0.2mg/mL的标准溶液。取上述标准溶液稀释为质量浓度为0.04，0.06，0.08，0.10，0.12，0.14，0.16，0.18，0.20mg/mL的芦丁溶液，分别吸取5mL上述各质量浓度芦丁溶液，置于1，2，3，4，5，6，7，8，9号共9支25mL容量瓶中，分别加入5%亚硝酸钠(NaNO2)溶液1.0mL，摇匀后室温下静置6min，再分别加入10%硝酸铝(Al(NO3)3)溶液1.0mL，摇匀后室温下静置6min，再加入1mol/L的氢氧化钠(NaOH)溶液8mL，用75%乙醇溶液定容至25mL。显色15min后取5号容量瓶进行紫外全波长扫描，以确定吸收峰。用蒸馏水调零后在300~900nm范围内的扫描，截取其在400~650nm范围内的吸光度变化示意图，结果如图 1所示。可知经过本方法处理过的黄酮类化合物在500~520nm处出现了明显的吸收峰，因此判断本方法能够用于CQF的测定。为减小红移和蓝移对检测结果的影响，选择510nm作为测定波长。
不同质量浓度芦丁标准品溶液按照本方法操作后，在510nm处测定各质量浓度对应的吸光度A，以芦丁标准品质量浓度(mg/mL)为横坐标，吸光度A值为纵坐标，绘制标准曲线，结果如图 2所示。标准曲线的回归方程为y=2.5617x-0.001，回归方程的相关系数R2=0.9995，说明该检测方法可行，所得回归方程能够用来测定CQF质量浓度。
2.3 单因素试验
设定松栎柱锈菌春孢子黄酮类化合物的初始提取条件为：提取温度A=40℃、料液比B=1:40、提取时间C=1h、乙醇体积分数D=50%(v/v)；每次试验所用春孢子粉200mg；超声频率40kHz，超声功率600W。将所得提取液5000r/min离心8min，吸取上层清液按照2.2的方法进行测定样品吸光值A。根据2.2中所得芦丁标准曲线：y=2.5617x-0.001推算CQF提取液中CQF的质量浓度，并通过式(1)计算CQF的得率。采用控制变量法依次考察该提取条件下上述4个单因子对CQF得率的影响。
T=NnVm×100% (1) 式中：T为CQF得率，%；N为根据标准曲线计算出的CQF的质量浓度，mg/mL；n为稀释倍数；V为提取液的体积，mL；m为供试春孢子质量，mg。
2.4 得率建模试验
在单因子考察试验基础上，根据中心组合实验(central composite design)原理使用软件Design-Expert 8.0.6对提取温度(A)，料液比(B)，提取时间(C)，乙醇体积分数(D)这4个因子较优的水平进行了CQF得率的响应面建模试验。完成建模试验后使用建模试验数据建立CQF得率关于以上4个影响因素的二次多项式回归方程数学模型，对所得模型及各个影响因子进行方差分析。通过软件Design-Expert 8.0.6与所得数学模型得出最优CQF提取条件，并验证优化结果。建模试验各因素水平与实际值对应关系见表 1。
表 1 建模试验因子水平表Table 1. Factors and levels of flavonoid extraction experiment水平
Level影响因子Influencing factor 温度Temperature(A)/℃ 料液比Solid-liquid ratio(B)/(g·mL-1) 时间Time(C)/min 乙醇体积分数
Ethanol contentration(D)/%-2 35 1:100 20 50 -1 45 1:120 30 60 0 55 1:140 40 70 1 65 1:160 50 80 2 75 1:180 60 90 2.5 CQF抗氧化能力测定
除CQF的还原力与DPPH·清除能力之外，清除超氧阴离子(O2-·)活性、清除羟基自由基(·OH)活性和总抗氧化活性均采用从南京建成生物研究所购买的对应的试剂盒测定，并严格按照试剂盒使用说明执行。本文所有CQF体外抗氧化活性测定实验均以抗坏血酸(Vc)做阳性对照，用以评价CQF的体外抗氧化活性。
2.5.1 CQF还原力测定
用体积分数为50%的乙醇溶液将在建模试验中所收集的CQF溶液稀释为原来质量浓度的0.2、0.4、0.6、0.8和1.0倍，分别向装有1.0mL各质量浓度CQF溶液的试管中加入0.2mol/L的pH 6.6的磷酸缓冲溶液2.5mL，快速混匀后加入质量分数1%的铁氰化钾溶液2.5mL并在50℃下水浴20min。迅速冷却反应液并加入质量分数10%的三氯乙酸溶液2.5mL，充分混匀，3000r/min离心10min后取上清液3mL，依次加入3mL蒸馏水和1.0mL的0.1%的三氯化铁溶液，充分混匀后静置10min，于700nm波长下测定溶液吸光度。吸光度越大还原能力也越大，计算公式如(2)。
A700=A样品−A空白 (2) 2.5.2 DPPH·清除能力测定
用95%乙醇溶液配置DPPH·浓度为0.2mmol/L的乙醇溶液，置于冰箱中避光冷藏备用。测定方法：测定管用2mL DPPH·乙醇溶液中加入2mL各质量浓度的CQF样品溶液，CQF质量浓度设置梯度与2.5.1相同。对照管用95%乙醇代替DPPH·溶液；空白管用蒸馏水代替CQF溶液。混匀器混匀后避光反应30min。波长调至517nm处用0.5mL蒸馏水和2.5mL 95%的乙醇调零，然后测定反应液吸光值A，每个质量浓度设置3个平行实验，结果取平均值。DPPH·自由基的清除率通过公式(3)计算得出。
DPPH⋅清除率=A0−(A1−A2)A0×100% (3) 式中：A0为2mL DPPH·乙醇溶液+2mL超纯水的吸光值；A1为2mL DPPH·乙醇溶液+2mL样品溶液后的吸光值；A2为2mL样品溶液+2mL超纯水的吸光值。
3. 结果与分析
3.1 单因子试验
3.1.1 乙醇体积分数对CQF得率的影响
由图 3可知，在其他影响因子不变的条件下，乙醇体积分数在50%~90%范围内变化导致CQF的得率先逐渐增大，当乙醇体积分数增至70%的时候，得率最大，但是仍然不足0.5%。当乙醇体积分数超过70%以后，随着体积分数的增加导致CQF得率骤降并逐渐趋于稳定。黄酮类化合物难溶于水，易溶于有机溶剂，水和乙醇可以以任意比例互溶。在水中加入不同量的乙醇，溶液的极性也随之不同，因而对某种或某类物质的溶解能力也不相同。根据所得实验结果绘制的折线图可以看出乙醇体积分数为70%是较优水平。
![]()
图 3 乙醇体积分数对黄酮类化合物(CQF)得率的影响Figure 3. Effects of ethanol concentration on extraction rate of flavonoids3.1.2 料液比对CQF得率的影响
根据图 4可以看出在其他条件不变的情况下，随着料液比的减小，CQF得率在不断升高，并逐渐趋于平稳，在料液比达到1:140以后，CQF得率增长幅度变小，得率增幅减小，因此选择1:140为较优水平。这一结果与CQF在50%乙醇溶液中的溶解度相关，当料液比较大时，CQF只有少部分溶于了乙醇溶剂，此时的提取液处于一种近饱和状态；随着提取液的增多，也就是料液比减小，未溶出的CQF部分也逐渐溶出，导致CQF的得率上升。当料液比减小到一定程度时，能溶于50%乙醇溶液的部分溶解殆尽，所以这时的得率逐渐趋于稳定，靠近极值。
![]()
图 4 料液比对CQF得率的影响Figure 4. Effects of solid-liquid ratio on extraction rate of flavonoids3.1.3 提取温度对CQF得率的影响图
在其他试验条件不变的条件下，考察了温度在15~75℃之间变化对CQF得率的影响，通过考察发现CQF得率先随着温度升高而增大，随着温度升高至65℃以后又逐渐降低。当提取液温度较低时，分子获得的能量相对较少，分子运动缓慢，有效提取CQF耗费的时间增长，这一点可以通过分子扩散现象很好地解释。但是，温度的变化同样会影响溶剂分子之间的作用力，从而影响溶剂溶解溶质的能力，导致CQF得率也随之发生变化。也就是说温度的变化同时影响着CQF分子运动速率与溶剂，因而CQF得率受此二者的交互作用影响，CQF得率随提取温度的变化如图 5所示。
![]()
图 5 提取温度对CQF提取率的影响Figure 5. Effects of extraction temperature on extraction rate of flavonoids3.1.4 提取时间对CQF得率的影响
在其他试验因子不变的条件下，考察提取时间对CQF得率的影响，结果显示提取时间对CQF得率的影响微弱，时间在20~60min变化时，CQF得率的变化始终处在0.46%~0.468%之间。CQF分子溶解到提取液中需要一定的时间，所以在提取时间40min以内时，CQF得率随着时间增加而增高，随着提取时间的延长，CQF的得率反而逐渐开始降低。这或许是因为部分原生质体内的CQF分子进入了细胞壁，因为细胞壁本身就可以容纳黄酮类化合物，另一方面，松栎柱锈菌春孢子表面具有油球，对黄酮类化合物具有一定的溶解能力。除此之外，随着溶剂分子运动时间的增加，乙醇分子扩散增多引起溶剂中乙醇浓度降低，从而导致CQF提取率也降低。CQF得率随提取时间的变化如图 6所示。
![]()
图 6 提取时间对CQF提取率的影响Figure 6. Effects of extraction time on extraction rate of flavonoids3.2 响应面法分析
3.2.1 建模试验设计及结果
使用软件Design-Expert 8.0.6对提取温度(A)，料液比(B)，提取时间(C)，乙醇体积分数(D)这4个因子的较优水平进行CQF得率的响应面建模试验。试验因子水平设计如表 1所示，中心组合实验设计与结果见表 2。
表 2 中心组合实验设计与结果Table 2. Design and results of central composite design试验号
Run No.温度
Temperature(A)/℃料液比
Solid-liquid ratio(B)/(g·mL-1)时间
Time(C)/min乙醇体积分数
Ethanol concentration(D)/%得率
Yield(T)/%1 0 -2 0 0 1.6329 2 1 -1 -1 1 1.0071 3 -1 1 1 1 0.9619 4 1 1 1 -1 1.4116 5 2 0 0 0 1.2952 6 0 0 0 0 1.1149 7 -1 1 1 -1 1.1555 8 0 0 0 0 1.1586 9 1 -1 1 1 1.1196 10 0 0 0 0 1.0930 11 1 1 -1 -1 1.3928 12 0 0 0 2 0.7542 13 0 0 0 0 1.1477 14 -1 -1 1 1 0.9697 15 1 -1 -1 -1 1.2507 16 -1 -1 -1 -1 1.1524 17 -1 -1 -1 1 0.8713 18 1 1 1 1 1.1805 19 -1 1 -1 1 0.9993 20 0 2 0 0 1.1453 21 -1 -1 1 -1 1.1055 22 0 0 -2 0 1.2133 23 1 -1 1 -1 1.3959 24 0 0 2 0 0.1969 25 1 1 -1 1 1.0555 26 0 0 0 -2 1.1040 27 0 0 0 0 1.1258 28 -2 0 0 0 1.3675 29 0 0 0 0 1.1477 30 -1 1 -1 -1 1.2304 3.2.2 数学模型的建立与分析
以CQF的得率为响应值，通过使用软件Design-Expert 8.0.6对得率T与4个影响因子回归拟合之后，得到CQF关于4个影响因子的数学模型方程：
T=1.13+0.051A−0.019B−0.071C−0.11D+1.171e−3AB+0.029AC−0.015AD−0.017BC−3.513e−3BD+0.016CD+0.059A2+0.073B2−0.098C2−0.042D2 对该数学模型进行方差分析(表 3)可知该模型回归显著(P＜0.05)，失拟项不显著(P=0.149)，说明未知因素对CQF得率影响较小，模型拟合效果良好。在我们设置的4因子水平梯度范围内，乙醇体积分数D对于CQF得率的影响显著(P＜0.05)，这与黄酮类化合物易溶于有机溶剂难溶于水的化学性质密不可分。由F值可得出各个影响因子对CQF得率的影响程度大小顺序为：乙醇体积分数(D)＞提取时间(C)＞提取温度(A)＞料液比(B)。
表 3 数学模型方差分析Table 3. Parameter estimate of regression model and variance analysis方差来源
Variance source平方和Sum
f squares自由度
df均方
Mean squareF值
F valueP值
P value显著性
Significance模型Model 1.150 14 0.082 7.469 0.043 1 显著Significant 温度Temperature(A) 0.062 1 0.062 1.520 0.237 1 料液比Solid-liquid ratio(B) 0.009 1 0.009 0.210 0.650 1 时间Time(C) 0.120 1 0.120 2.900 0.109 1 乙醇体积分数Ethanol concentration(D) 0.290 1 0.290 7.010 0.018 3 AB 0.000 1 0.000 0.001 0.981 9 AC 0.013 1 0.013 0.320 0.577 2 AD 0.004 1 0.004 0.093 0.765 2 BC 0.005 1 0.005 0.120 0.736 6 BD 0.000 1 0.000 0.005 0.945 7 CD 0.004 1 0.004 0.100 0.756 6 A2 0.095 1 0.095 2.300 0.149 8 B2 0.150 1 0.150 3.580 0.078 1 C2 0.260 1 0.260 6.380 0.023 3 D2 0.048 1 0.048 1.160 0.297 5 残差Residual 0.620 15 0.041 失拟项Lack of fit 0.610 10 0.061 0.017 0.149 不显著Not significant 纯误差Pure error 0.003 5 0.001 总误差Total error 1.770 29 3.2.3 交互项对CQF得率的影响
响应面是以响应值与影响响应值的因子构建起了空间三维图形，可以让我们更直观的判断出两两不同影响因子交互作用对响应值的影响趋势以及响应值的变化范围。响应面图中颜色的变化表示响应值(CQF得率)的变化，颜色变化快，响应值变化也快，一般较低的响应值使用蓝色表示，较高的响应值用红色表示，其他颜色均为中间过渡颜色。同时，在一定方向和范围内响应面3D图形中坡度越大表示响应值变化的越快，即对CQF得率的影响更为显著。这一点在等高线图形中也有很好的体现，当某条线段在垂直其中一条轴的方向上移动时，在一定范围内与它所相交的等高线越多，则说明在该范围内某一因子对响应值的影响越显著。等高线表示在此条线上所有试验方案都能够得到相同的响应值，即虽然等高线所对应的试验条件不同，但是能够产生相同的实验结果。4个影响因子的两两之间交互作用对CQF得率的影响示意图见图 7。
![]()
图 7 各因素交互作用对CQF得率影响的响应面与等高线图Figure 7. Response surface and contour plots for the interaction effects of various factors on extraction rate of CQF3.2.4 最佳提取工艺及验证
根据在3.2.2中对模型进行的方差分析，已经得知我们拟合出的回归方程模型显著，失拟项不显著，说明该方程是适合用于CQF得率与提取工艺的数学模型，所以我们可以借助于该方程确定最佳的CQF提取工艺。经软件Design-Expert 8.0.6分析，得知模型预测的最佳提取工艺为：提取温度65℃，料液比1:120，提取时间38min，提取液乙醇体积分数60%。在此条件下有对该优化结果进行了3次验证实验，3次验证实验所得出的CQF得率为1.382%，实际得率与在该条件下模型所预测的得率1.416%相差2.46%，这表明我们优化出来的CQF提取条件真实可靠，能够应用于实践。
3.3 抗氧化活性分析
3.3.1 CQF的还原力测定
还原力是评价物质抗氧化能力的重要指标之一，其原理是具有较强还原力的物质将Fe3+还原为Fe2+，Fe2+在酸性环境中能与FeCl3产生有颜色的普鲁士蓝。在该检测体系下物质的还原力越强，对应的吸光值也越大。根据图 8可以看出，CQF与Vc的还原力与质量浓度均呈现很好的线性关系，拟合度均超过0.998。在图示浓度范围内，CQF和Vc的质量浓度与还原力均表现出较好的量效关系，相同质量浓度条件下CQF的还原力更强，而且CQF质量浓度与还原力的线性方程斜率更大，量效关系更为明显。在图示范围内，不论是量效关系还是还原力，CQF都强于Vc，说明CQF具备很好的还原力。
3.3.2 CQF对DPPH·的清除能力
在图 9中，根据受试物质量浓度与DPPH·清除率的散点图，我们可以看出Vc对DPPH·的清除能力较CQF强，并迅速逼近清除率最大值。在散点图的基础上，可以得出CQF与Vc清除DPPH·的IC50分别为0.0492mg/mL和0.0335mg/mL，也就是说对CQF对DPPH·清除能力较Vc差。但是随着CQF质量浓度的逐步增大，CQF对DPPH·清除能力也逐步增强，量效关系越来越明显。对于该种人工合成自由基，CQF的反应不够灵敏。
![]()
图 9 受试物清除DPPH·能力的比较Figure 9. Comparison of DPPH· radical scavenging capacity of CQF and Vc3.3.3 CQF总抗氧化能力
生物活性物质总抗氧化能力是指活性物质清除各种自由基的有效总和。虽然已有多种测定抗氧化活性的方法，但这些方法只是从某个角度来衡量天然抗氧化物的活性。事实上，抗氧化作用涉及到多种作用机制，它涉及一系列的问题，目前尚未形成标准方法[16]。本实验测定总抗氧化能力采用的南京建成生物研究所试剂盒原理为FRAP法，在该种检测方法下，由图 10可以看出CQF与Vc的总抗氧化能力相当。根据散点图经计算得知，CQF的半抑制率IC50为0.0303mg/mL，而Vc的半抑制率为0.0240mg/mL，可以看出Vc的总抗氧化能力较CQF强，但是CQF与Vc趋势线在0.0654mg/mL处相交。整体来看，两种物质达到最大清除率都不需要较大浓度，量效关系相近，其总抗氧化能力并不具有很大的差异。所以说，CQF的总抗氧化能力还是非常理想的。
3.3.4 CQF对·OH的清除能力
羟基自由基(·OH)属于一种强氧化剂(E0=2.80V，3.06V)，氧化能力远大于一般氧化剂，是机体中最活跃的一种活性氧分子，具有很强的破坏性。它几乎可以和活细胞中的任何生物分子发生反应而造成损伤，而且反应速度极快[17]。因此，清除机体内过多的·OH具有非常重要的生物学意义。通过与Vc的对比测试，我们可以看出Vc在质量浓度0.016 ~0.080mg/mL范围内呈现线性变化规律，但对·OH较强的清除能力远低于CQF。由图 11可以非常明显的看出Vc与CQF在·OH清除能力上的差距。虽然Vc具有很强的还原性，不仅能够直接清除·OH，还能够阻断Fenton反应，从而减少体系中的·OH，但是CQF对·OH的敏感程度与清除能力依然远胜Vc。经测定，CQF对·OH清除能力的IC50仅为0.0129mg/mL，因为Vc对·OH的清除能力远小于CQF，故未予测定与计算其对·OH的半抑制率。
3.3.5 CQF对O2-·的清除能力
超氧阴离子O2-·经过一系列反应可以生成其他氧自由基，是大多数氧自由基的母体，能够攻击细胞DNA，具有较强的毒性，所以清除生物机体内过多的O2-·具有重要意义。理论上来讲，CQF能够提供电子从而有效清除自由基，但本课题在测定CQF清除O2-·能力时却发现了截然相反的实验结果。由图 12可知，两种受试物对清除O2-·都具有一定的量效关系，但是测定结果不在同一象限内。这也就是说在该中检测体系下受试物CQF非但没有清除O2-·，反而促进了O2-·的产生。所使用试剂盒检测原理为模拟机体中黄嘌呤与黄嘌呤氧化酶反应系统，黄嘌呤在黄嘌呤酶的催化下转化为尿酸，在此过程中能够产生O2-·，然后通过衡量O2-·的量来推算加入反应体系的物质时清除还是产生O2-·。黄酮类化合物作为电子提供物质，不会产生O2-·，但是CQF提取液是未经过除杂的，一些会对该反应体系产生影响的物质可能一并被提取了，比如说嘌呤或黄嘌呤氧化酶以及一些自氧化物质，这些物质生成O2-·的能力远超黄酮类物质清除O2-·的能力，所以才产生了如图 12所示的测定结果[18-20]。这一问题没有得到有效解决，所以尚不清楚导致该检测结果的具体原因，需要在后续工作中进一步研究导致该结果的具体原因与解决方法。黄嘌呤酶氧化法目前是测定O2-·最完善最普遍的方法，但是这种方法应该对待测样品设置一定的要求或合适的处理条件以排除干扰因子，需要我们进一步完善这种测定方法。
4. 结论与讨论
从微生物孢子中提取黄酮类化合物的例子尚不多见，在此背景下，本研究采用控制变量法逐个考察了能够影响CQF得率的4个因子。然后使用软件Design-Expert 8.0.6设计了建模试验，完成建模试验之后利用可靠的CQF得率数学模型得到了适用的最优CQF提取条件：提取温度65℃，料液比1:120，提取时间38min，提取液乙醇体积分数60%，在该提取条件下CQF得率为1.382%。虽然松栎柱锈菌春孢子具有厚且坚硬的几丁质细胞壁，能够有效阻止细胞内含物的流失以抵抗不良环境因素，但是我们使用一般植物总黄酮提取的方法并辅以超声振荡也提取到了CQF，说明醇提取法同样能够用于其他孢子总黄酮的提取并且达到了预期的CQF得率。总黄酮得率既取决于提取方法，也取决于实验材料总黄酮含量。一般来说针叶树实验材料总黄酮得率高于菌物但不足5%。闫琳娜通过对落叶松(Larix gmelinii)松针黄酮的研究发现落叶松松针总黄酮含量较高且在不同提取条件下得率都能达到6%以上[21]；孙琼优化了在菌物中黄酮类化合物含量相对较高的杏鲍菇(Pleurotus eryngii)黄酮类化合物的提取条件，但是得率仍然不足3%[22]。同时我们也看到了黄嘌呤酶氧化法在实际应用中暴露出来的问题，没能有效测定出CQF清除O2-·的能力，需要进一步明确并排除干扰因子或者改进O2-·检测方法。
黄晓冬等对阔叶树种白桦(Betula platyphylla)活体寄生真菌桦褐孔菌(Inonotus obliquus)总黄酮的研究表明其具有较强的体内外抗氧化活性[23]，Boonsong等研究了5种食用菌总酚与类黄酮的提取与抗氧化活性，发现乙醇溶液更适合菌物总酚的提取，其提取物均表现出较强的抗氧化活性，然而很多的食用菌不是活体寄生或不含有黄酮类物质[24]。整体来看，对于针叶树总黄酮的研究成果与热情远胜担子菌，而活体寄生的担子菌总黄酮的研究很鲜见。本研究所选研究对象将针叶树与担子菌门菌物总黄酮的研究很好地有机结合，旨在为锈菌目菌物开发利用、森林有害生物开发利用和转变森林有害生物防治理念提供科学依据。迄今为止，人类对锈菌的认知还很不完善，尚未看到国内外任何关于锈菌开发利用的报道。本研究首次完成了CQF的提取与体外抗氧化活性测定工作，发现CQF具有很强的还原力，还原力较Vc强；其次CQF的总抗氧化能力和DPPH·清除能力与Vc相当；CQF对·OH具有非常好的清除作用，这一能力远胜Vc，与Treml等对黄酮类化合物清除·OH能力的研究现状综述相一致[26]。体外试验结果为体内试验提供了一定的科学依据，这表明CQF具有开发利用的潜在价值。
长久以来，人类对于有害生物的认识是不全面的，往往只是单纯的从治理的角度去对待，总是尝试使用各种各样的手段降低或控制有害生物的种群数量。森林有害生物确实为林业生产活动和森林健康产生了很多负面的影响，如果说发现并证实某种有害生物有其独特的利用价值，那么它的某种功能或价值就能在一定程度上抵偿它所造成的损失，从而引导我们转变传统的有害生物防治理念。
-
![]()
图 1 试验小区布设图
CK.对照组；Y.叶面肥；C.冲施肥。数字表示肥料用量，单位为mL/m2。
Figure 1. Layout of experimental plots
CK, control; Y, foliar fertilizer; C, flush fertilizer. The numbers represent the amount of fertilizer, the unit is mL/m2.
![]()
图 2 2016年7—8月降雨量与土壤蒸发量
Figure 2. Rainfall and soil evaporation during July to August in 2016
![]()
图 3 冲施肥对沙木蓼土壤含水率的影响
不同字母表示处理间差异显著(P＜0.05)。下同。
Figure 3. Effects of flush fertilizer on the soil water content of Atraphaxis bracteata
Different letters indicate significant difference among treatments at P < 0.05 level. The same below.
![]()
图 4 叶面肥对沙木蓼土壤含水率的影响
Figure 4. Effects of foliar fertilizer on soil water content of Atraphaxis bracteata
![]()
图 7 施肥对沙木蓼地径增长量的影响
Figure 7. Effects of fertilization on the growth of ground diameter of Atraphaxis bracteata
![]()
图 8 施肥对沙木蓼株高增长量的影响
Figure 8. Effects of fertilization on the tree height increment of Atraphaxis bracteata
![]()
图 9 施肥对沙木蓼新梢生长的影响
Figure 9. Effects of fertilization on the growth of Atraphaxis bracteata new shoots
![]()
9 施肥对沙木蓼群落物种多样性的影响
不同字母表示处理间差异显著(P＜0.05)。
9. Effects of fertilization on the species diversity of Atraphaxis bracteata
Different letters indicate significant difference among treatments at P < 0.05 level.
表 1 沙木蓼地径与株高的原始值分析
Table 1 Analysis on original values of ground diameter and tree height of Atraphaxis bracteata
cm 施肥处理
Fertilization treatment植物指标
Plant index肥料用量Fertilizing amount 0 3 mL/m2 6 mL/m2 9 mL/m2 12 mL/m2 15 mL/m2 冲施肥Flush fertilizer 地径Ground diameter 0.50a 0.49a 0.53a 0.54a 0.57a 0.53a 株高Plant height 104.87a 103.87a 119.67a 104.93a 106.40a 104.93a 叶面肥Foliar fertilizer 地径Ground diameter 0.50a 0.53a 0.49a 0.54a 0.56a 0.55a 株高Plant height 104.87a 103.33a 100.33a 108.73a 108.60a 104.67a 注：不同小写字母表示不同处理之间差异显著(P>0.05)。Note: different lowercases mean significant difference among varied treatments at P>0.05 level. 
下载: 导出CSV
表 2 鲜海带生物酶解有机液肥主要成分
Table 2 Main ingredients of fresh seaweed bio-enzymatic liquid organic fertilizer
液肥
Liquid fertilizer主要成分
Main component叶面肥
Foliar fertilizer海带酶解液、海藻多糖、海藻蛋白、甘露醇、甜菜碱、赤霉素、吲哚乙酸、脱落酸、氨基酸等Kelp hydrolysate, seaweed polysaccharide, seaweed protein, mannitol, betaine, gibberellin, indoleacetic acid, abscisic acid, amino acid, etc 冲施肥
Flush fertilizer海带纤维、有机质、海藻多糖、吲哚乙酸、海藻蛋白、甘露醇、甜菜碱等Kelp fiber, organic matter, seaweed polysaccharide, indoleacetic acid, seaweed protein, mannitol, betaine, etc
下载: 导出CSV
表 3 小区施肥用量
Table 3 Fertilizing amount in 3 plots
肥料用量
Fertilizing amount/(mL·m-2)小区肥料用量
Community fertilizing amount/mL小区用水量
Community water consumption/mL小区补水量
Community water recharge/mL0 0 27 000 0 3 27 5 400 21 600 6 54 10 800 16 200 9 81 16 200 10 800 12 108 21 600 5 400 15 135 27 000 0
下载: 导出CSV
表 4 不同梯度冲施肥处理的物种相对丰度
Table 4 Relative abundance of species in different gradient treatments of flush fertilizer
物种名称Species name CK 3 mL/m2 6 mL/m2 9 mL/m2 12 mL/m2 15 mL/m2 杨柴Hedysarum fruticosum 44.84 23.88 22.09 19.06 16.08 14.53 白蒿Artimisiae sphaerocephala 25.82 21.70 0 6.36 22.99 39.40 油蒿Artemisia ordosica 20.92 9.20 0 10.00 0 0 冰草Agropyron cristatum 8.42 0 22.42 24.57 0 9.53 芦苇Phragmites australis 0 14.48 6.30 0 0 12.15 狗尾草Setaria viridis 0 17.70 19.76 18.66 15.54 12.37 砂蓝刺头Echinops gmelini 0 13.49 0 0 0 6.36 沙蓬Agriophyllum squarrosum 0 14.56 5.39 0 0 0 角蒿Incarvillea sinensis 0 0 0 0 4.53 0
下载: 导出CSV
表 5 不同梯度叶面肥处理的物种相对丰度
Table 5 Relative abundance of species in different gradient treatments of foliar fertilizer
物种名称Species name CK 2 mL/m2 4 mL/m2 6 mL/m2 8 mL/m2 10 mL/m2 杨柴Hedysarum fruticosum 58.33 29.67 15.06 18.62 13.47 12.75 白蒿Artimisiae sphaerocephala 0 34.71 7.10 31.36 27.99 19.52 油蒿Artemisia ordosica 41.67 6.56 8.87 23.11 11.92 17.27 冰草Agropyron cristatum 0 24.56 21.89 21.07 22.57 16.53 芦苇Phragmites australias 0 0 12.32 0 6.57 12.52 狗尾草Setaria viridis 0 9.56 11.42 14.18 15.94 7.28 针茅Stipa capillata 0 0 10.59 0 0 6.40 花棒Hedysarum scoparium 0 0 0 10.61 0 0
下载: 导出CSV
-
[1] 孙保平.荒漠化防治工程学[M].北京:中国林业出版社, 2000:1. Sun B P. Desertification engineering [M]. Beijing: China Forestry Publishing House, 2000:1.
[2] 韩东锋, 孙德祥, 周广阔, 等.半荒漠风沙区5种优良沙生灌木造林效果比较[J].西北农业学报, 2009, 18(5): 312-315. doi: 10.3969/j.issn.1004-1389.2009.05.068 Han D F, Sun D X, Zhou G K, et al. Comparison among the afforest effects of 5 excellent desert shrubs in semi-desert area [J]. Acta Agriculturae Boreali-Occidentalis Sinica, 2009, 18 (5):312-315. doi: 10.3969/j.issn.1004-1389.2009.05.068
[3] 赵岩, 周文渊, 孙保平, 等.毛乌素沙地三种荒漠灌木根系分布特征与土壤水分研究[J].水土保持研究, 2010, 17(4):129-133. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=stbcyj201004028 Zhao Y, Zhou W Y, Sun B P, et al. Root distribution of three desert shrubs and soil moisture in Mu Us Sand Land[J]. Research of Soil and Water Conservation, 2010, 17(4):129-133. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=stbcyj201004028
[4] Harpole W S, Potts D L, Suding K N. Ecosystem responses to water and nitrogen amendment in a California grassland[J]. Global Change Biology, 2007, 13(11): 2341-2348. doi: 10.1111/gcb.2007.13.issue-11
[5] 周国英, 陈桂琛, 赵以莲, 等.施肥和围栏封育对青海湖地区高寒草原影响的比较研究Ⅱ地上生物量季节动态[J].草业科学, 2005, 22(1): 59-63. doi: 10.3969/j.issn.1001-0629.2005.01.019 Zhou G Y, Chen G C, Zhao Y L, et al. Comparative studies on influence of chemical fertilizer application and enclosure on alpine steppes in Qinghai Lake (Ⅱ): seasonal and annual biomass dynamics[J]. Pratacultural Science, 2005, 22(1): 59-63. doi: 10.3969/j.issn.1001-0629.2005.01.019
[6] 王大明.试论高寒人工草地施氮的增产效应[J].中国草地, 1994, 16(2): 76-80. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400569880 Wang D M. A discussion on incremental response of to alpine sown pasture N application[J]. Grassland of China, 1994, 16(2): 76-80. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400569880
[7] 孙斌.三种改良措施对高寒退化草地植被的影响[J].甘肃农业大学报, 2005, 40(6): 797-801. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gsnydxxb200506017 Sun B. Three grassland improvement measures on alpine degraded grassland vegetation[J]. Journal of Gansu Agricultural University, 2005, 40(6): 797-801. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=gsnydxxb200506017
[8] 纪亚君.青海高寒天然草地施氮肥试验[J].中国草地, 2005, 27(6): 69-71. http://d.old.wanfangdata.com.cn/Periodical/zgcd200506014 Ji Y J. Effects of nitrogen on grass yield in Qinghai alpine meadow[J]. Grassland of China, 2005, 27(6): 69-71. http://d.old.wanfangdata.com.cn/Periodical/zgcd200506014
[9] 仁青吉, 罗燕江, 王海洋, 等.青藏高原典型高寒草甸退化草地的恢复:施肥刈割对草地质量的影响[J].草业学报, 2004, 13(2):43-49. doi: 10.3321/j.issn:1004-5759.2004.02.007 Ren Q J, Luo Y J, Wang H Y, et al. Restoration of degraded typical meadowland on the Qinghai-Tibetan Plateau: effect of fertilizing and cutting on grassland quality[J]. Acta Prataculturae Sinica, 2004, 13(2):43-49. doi: 10.3321/j.issn:1004-5759.2004.02.007
[10] 李燕婷, 肖艳, 李秀英, 等.作物叶面施肥技术与应用[M].北京:科学出版社, 2009:11. Li Y T, Xiao Y, Li X Y, et al. The foliar fertilization technology and application[M]. Beijing: Science Press, 2009:11.
[11] 王周琼, 李述刚, 程心俊.有机无机复合与荒漠化防治[M].北京:科学出版社, 2003:92-212. Wang Z Q, Li S G, Cheng X J. Organic inorganic compound and prevention of desertification [M]. Beijing: Science Press, 2003:92-212.
[12] 马红梅, 王锁民, 张建全, 等.专用有机肥对荒漠植物梭梭的生长及耐旱性的影响[J].草业科学, 2016, 33(2):240-248. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=caoyekx201602009 Ma H M, Wang S M, Zhang J Q, et al. The special organic fertilizer in improving the growth and drought tolerance of desert plant Haloxylon ammodendron[J]. Pratacultural Science, 2016, 33(2):240-248. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=caoyekx201602009
[13] Goldberg D E, Miller T E. Effects of different resource additions of species diversity in an annual plant community [J]. Ecology, 1990, 71(1): 213-225. doi: 10.2307/1940261
[14] 汪金松, 赵秀海, 张春雨, 等.模拟氮沉降对油松林土壤有机碳和全氮的影响[J].北京林业大学学报, 2016, 38(10): 88-94. http://d.old.wanfangdata.com.cn/Periodical/bjlydxxb201610011 Wang J S, Zhao X H, Zhang C Y, et al. Effects of simulated nitrogen deposition on soil organic carbon and total nitrogen content in plantation and natural forests of Pinus tabuliformis[J]. Journal of Beijing Forestry University, 2016, 38(10):88-94. http://d.old.wanfangdata.com.cn/Periodical/bjlydxxb201610011
[15] Xu W W, Liu H Y, Tan X S, et al. Effects of seaweed bio-organic fertilizer on growth and yield of winter wheat[J]. Agricultural Science & Technology, 2016, 17(11):2555-2559. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hnnykjtx-e201611027
[16] Liu P J, Wang F, Zhang S Q. Effect of seaweed bio-organic fertilizer on seed germination and seedling growth of vegetables[J]. Agricultural Science & Technology, 2014, 15(3):407-411.
[17] 张军红.毛乌素沙地油蒿群落土壤水分分布与动态[D].北京: 中国林业科学研究院, 2013. Zhang J H. Distribution and dynamics of soil moisture in Artemisia ordosica communities in Mu Us Sandy Land [D]. Beijing: Chinese Academy of Forestry, 2013.
[18] 李博.生态学[M].北京:高等教育出版社, 2004:118-119. Li B. Ecology [M]. Beijing: Higher Education Press, 2004:118-119.
[19] 王莉, 秦树高, 张宇清, 等.生物土壤结皮对毛乌素沙地油蒿群落土壤水分的影响[J].北京林业大学学报, 2017, 39(3):48-56. doi: 10.13332/j.1000-1522.20160383 Wang L, Qin S G, Zhang Y Q, et al. Influence of biological soil crusts on soil moisture in Artemisia ordosica community in Mu Us Desert, northwestern China[J]. Journal of Beijing Forestry University, 2017, 39(3):48-56. doi: 10.13332/j.1000-1522.20160383
[20] 穆家伟, 查天山, 贾昕, 等.毛乌素沙地典型沙生灌木对土壤蒸发的影响[J].北京林业大学学报, 2016, 38(12): 39-45. doi: 10.13332/j.1000-1522.20160192 Mu J W, Zha T S, Jia X, et al. Influence of typical sandy shrubs on soil evaporation in Mu Us Sandland, northwestern China[J].Journal of Beijing Forestry University, 2016, 38(12):39-45. doi: 10.13332/j.1000-1522.20160192
[21] 王改兰, 段建南, 贾宁凤, 等.长期施肥对黄土丘陵区土壤理化性质的影响[J].水土保持学报, 2006, 20(4): 82-85, 89. doi: 10.3321/j.issn:1009-2242.2006.04.019 Wang G L, Duan J N, Jia N F, et al. Effects of long-term fertilization on soil physical and chemical property in loess hilly area[J]. Journal of Soil and Water Conservation, 2006, 20(4):82-85, 89. doi: 10.3321/j.issn:1009-2242.2006.04.019
[22] 苏洁琼, 李新荣, 鲍婧婷.施氮对荒漠化草原土壤理化性质及酶活性的影响[J].应用生态学报, 2014, 25(3): 664-670. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yystxb201403006 Su J Q, Li X R, Bao J T. Effects of nitrogen addition on soil physico-chemical properties and enzyme activities in desertified steppe[J]. Chinese Journal of Applied Ecology, 2014, 25(3):664-670. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yystxb201403006
[23] 胡静.绿色纤维的应用及发展[J].科技创新导报, 2013 (11):66. doi: 10.3969/j.issn.1674-098X.2013.11.053 Hu J. The application and development of green fiber[J]. Science and Technology Innovation Herald, 2013(11):66. doi: 10.3969/j.issn.1674-098X.2013.11.053
[24] Vourlitis G L, Pasquini S, Zorba G. Plant and soil N response of southern Californian semi-arid shrub lands after 1 year of experimental N deposition[J].Ecosystems, 2007, 10:263-279. doi: 10.1007/s10021-007-9030-2
[25] 涂利华, 胡庭兴, 张健, 等.华西雨屏区苦竹林土壤酶活性对模拟氮沉降的响应[J].应用生态学报, 2009, 20(12): 2943-2948. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yystxb200912015 Tu L H, Hu T X, Zhang J, et al. Soil enzyme activities in a Pleioblastus amurus plantation in rainy area of west China under simulated nitrogen deposition [J]. Chinese Journal of Applied Ecology, 2009, 20(12):2943-2948. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yystxb200912015
[26] Patterson T B, Guy R D, Dang Q L.W. Hole-plant nitrogen and water-relations trais, and their associated tradeoffs, in adjacent muskeg and upland boreal spruce species[J]. Oecologia, 1997, 110(2):160-168. doi: 10.1007/s004420050145
[27] 刘彤, 王琛, 徐璐, 等.叶面施氮对东北红豆杉幼苗生长及光合特性的影响[J].北京林业大学学报, 2016, 38(8): 64-70. doi: 10.13332/j.1000-1522.20160099 Liu T, Wang C, Xu L, et al. Effects of foliar nitrogen application on growth and photosynthetic characteristics of Taxus cuspidata seedings[J]. Journal of Beijing Forestry University, 2016, 38(8):64-70. doi: 10.13332/j.1000-1522.20160099
[28] 李宜伟, 蒋进, 宋春武, 等.发酵肥料对沙地水分与植物生长的影响[J].干旱区资源与环境, 2012, 26(2): 161-166. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ghqzyyhj201202028 Li Y W, Jiang J, Song C W, et al. Effects of fermented manure on sandy land soil moisture and growth of plants[J]. Journal of Arid Land Resources and Environment, 2012, 26(2):161-166. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=ghqzyyhj201202028
[29] 陈亚明, 李自珍, 杜国祯.施肥对高寒草甸植物多样性和经济类群的影响[J].西北植物学报, 2004, 24(3): 424-429. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xbzwxb200403009 Chen Y M, Li Z Z, Du G Z. Effects of fertilization on plant diversity and economic herbage groups in alpine meadow [J].Acta Botanica Boreali-Occidentalia Sinica, 2012, 24(3):424-429. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=xbzwxb200403009
[30] 郑华平, 陈子萱, 王生荣, 等.施肥对玛曲高寒沙化草地植物多样性和生产力的影响[J].草业学报, 2007, 16(5): 34-39. doi: 10.3321/j.issn:1004-5759.2007.05.004 Zheng H P, Chen Z X, Wang S R, et al. Effects of fertilizer on plant diversity and productivity of desertified alpine grassland at Maqu, Gansu [J]. Acta Prataculturae Sinica, 2007, 16(5):34-39. doi: 10.3321/j.issn:1004-5759.2007.05.004
[31] 王堃, 吕进英, 邵新庆.沙漠化草地的恢复与重建途径[J].草地学报, 2004, 12(3):240-245. http://d.old.wanfangdata.com.cn/Periodical/cdxb200403016 Wang K, Lü J Y, Shao X Q. Measures of restoring and rebuilding desertified grassland[J]. Acta Agrestia Sinica, 2004, 12(3):240-245. http://d.old.wanfangdata.com.cn/Periodical/cdxb200403016
[32] 王长庭, 龙瑞军, 王启基, 等.不同类型高寒草地群落物种特征和均匀度的重要性[J].草地学报, 2005, 13(4): 320-323. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=cdxb200504013 Wang C T, Long R J, Wang Q J, et al. Importance of species characteristics and evenness in different type grassland[J]. Acta Agrestia Sinica, 2005, 13(4):320-323. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=cdxb200504013
计量
- 文章访问数: 2163
- HTML全文浏览量: 479
- PDF下载量: 49
