留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

沙棘林退化对林木根系和土壤抗冲性的影响

舒乔生 谢立亚

舒乔生, 谢立亚. 沙棘林退化对林木根系和土壤抗冲性的影响[J]. 浙江农林大学学报, 2014, 31(3): 380-385. DOI: 10.11833/j.issn.2095-0756.2014.03.008
引用本文: 舒乔生, 谢立亚. 沙棘林退化对林木根系和土壤抗冲性的影响[J]. 浙江农林大学学报, 2014, 31(3): 380-385. DOI: 10.11833/j.issn.2095-0756.2014.03.008
HUANG Haiyan, TANG Huimin, JIN Lu, et al. Effects of tree species importance, diversity and soil physicochemical properties on soil amino sugars[J]. Journal of Zhejiang A&F University, 2024, 41(4): 778-786. DOI: 10.11833/j.issn.2095-0756.20230527
Citation: SHU Qiaosheng, XIE Liya. Root system distribution and soil anti-scouring capability with Hippophae rhamnoides degradation[J]. Journal of Zhejiang A&F University, 2014, 31(3): 380-385. DOI: 10.11833/j.issn.2095-0756.2014.03.008

沙棘林退化对林木根系和土壤抗冲性的影响

DOI: 10.11833/j.issn.2095-0756.2014.03.008
基金项目: 

辽宁省自然科学基金资助项目 20102109,20120210

重庆市水利局重大项目 渝水科201306

详细信息
    作者简介: 舒乔生,副研究员,博士,从事水土保持科学研究。E-mail:cysqs@163.com
    通信作者: 谢立亚,副研究员,从事水土保持科学研究。E-mail:xly2004@126.com
  • 中图分类号: S714.5;S157.1

Root system distribution and soil anti-scouring capability with Hippophae rhamnoides degradation

  • 摘要: 近年来,辽宁建平县亚洲最大的人工沙棘Hippophae rhamnoides林呈现大面积死亡的现象。通过标准样地法调查沙棘的死亡程度,测定与对比不同沙棘死亡率下林地的根系密度和土壤抗冲性,研究沙棘林退化对水土流失的影响。结果表明:0~30 cm范围内,直径d<1 mm的根系密度随深度增大急剧减小,72%~80%根量分布在0~20 cm土层;0~30 cm范围内,随着沙棘死亡率的提高,d<1 mm根密度显著降低(P<0.05);随着深度增大,抗冲性也逐步降低,沙棘的死亡率对土壤抗冲性的影响亦达到显著水平(P<0.05),这主要是由根系密度的相应变化引起的;经曲线拟合,土壤抗冲性与d<1 mm根密度的关系可用幂函数表示。目前该区主要应通过采取有效措施降低沙棘林退化的生态影响。
  • 竹材和竹制品的环保特性已成为竹产业的标志属性。毛竹Phyllostachys edulis的微观结构是由维管束纤维与基体组成的两相复合材料,维管束纤维是影响毛竹宏观力学性能的关键。现有竹纤维分离方法可分为高温蒸煮和机械冲击摩擦2类,如蒸汽爆破、碾压捶打、机械梳解、化学分离等[18]。以上分离方法关注的焦点均在于所得纤维产品的性能优劣,而忽视了维管束与薄壁组织的分离效果评价。在分离过程中,常存在纤维热损伤、比强度降低、生产效率低、纤维尺寸一致性差等问题。

    数控加工中心可通过数控程序精确控制刀具路径和切削方向,并通过合理设置刀具路径和铣削条件精确控制纤维形状。杨永福等[9]分析了不同刀具前角、切削速度以及进给量等切削参数对竹片平面直角自由切削过程中切削力的影响;郭莹洁等[1011]研究了铣削加工参数的改变对竹片表面质量及超前劈裂的影响。OGAWA等[1213]通过设定合适的加工中心切削参数,切削提取得到了长度均匀、无热损伤的竹纤维。上述研究主要围绕竹制品的切削特性和切削提取原竹纤维的加工条件展开,并未深入探讨天然竹维管束纤维制取过程中切削力及切削特性等变化。

    为揭示天然毛竹维管束纤维提取过程中铣削参数对切削力及维管束纤维形态的影响,本研究以切削加工三要素为变量,使用双刃直槽硬质合金木工雕刻刀对竹板开展单向顺铣正交切削试验,运用极差、方差分析方法,验证分析铣削参数对切削力的影响,并建立切削力经验公式;基于瞬态切削几何模型和单因素试验,探究切削速度(Vc)、切削深度(Ap)和每齿进给量(fz)对维管束纤维提取质量的影响,以期为合理选取切削参数,高效获取优质竹纤维提供理论以及技术指导。

    选取浙江湖州5年生毛竹,并截取直径约110 mm,竹壁厚度约10 mm的竹段作为试材。去除竹节后放入烘箱中,在70~80 ℃下烘烤8 h。然后,用破竹器将每段试材分成8根竹条,并去除0.8 mm厚度的竹青表皮,最后对竹条进行切割,得到高度分别为36、40和44 mm的竹板毛坯。加工机床采用宁波其锐达机械有限公司生产的KMD-80120型三轴高速雕铣机,主轴最高转速达到12 000 r·min−1;试验刀具为直径10 mm的双刃直槽硬质合金木工雕刻刀。在试验过程中,刀刃径向跳动调整在6 μm以下。切削力测量系统为宁波灵元测控工程有限公司生产的压电晶体测力仪。

    对于纤维增强复合材料,具有大长径比的纤维是理想的[1417],因此本研究以纤维长度、纤维直径和纵横比评判维管束纤维质量。设计正交切削试验并获取各参数下的切削力数据,再基于切削力数据分析切削参数对切削力及维管束纤维质量的影响规律。

    参考以往研究[1820],本研究将竹壁沿径向划分为图1所示的内、中、外3层,其中外部维管束纤维呈半开放型,纤维长轴平均长度为0.39 mm,短轴平均长度为0.23 mm。经测量,竹壁外部维管束纤维在水平方向的平均间隙为0.2~0.4 mm。切削过程中,为得到尽可能完整的维管束纤维,切削步进距(Ae)设置为0.4 mm。在以上加工参数的基础上,将切削速度、每齿进给量及切削深度分别定义为因素ABC,制定正交试验方案(表1)。试验采用刀刃平行于维管束纤维,沿X轴单向顺铣,干式切削方式进行,对每个试件连续切削3个步距,每组参数重复切削5个试件。通过以上正交试验,可得到135组切削力数据。

    图 1  竹壁分层示意图
    Figure 1  Schematic diagram of bamboo wall layering
    表 1  正交试验及切削力
    Table 1  Forces of orthogonal experimental cutting
    序号 切削速度(A)/
    (m·min−1)
    每齿进给量(B)/
    mm
    切削深度(C)/
    mm
    切削力/N 序号 切削速度(A)/
    (m·min−1)
    每齿进给量(B)/
    mm
    切削深度(C)/
    mm
    切削力/N
    Fx Fy Fx Fy
    1 62.8 0.2 4 23.70 −8.89 6 125.7 0.4 4 34.78 −10.86
    2 62.8 0.3 12 81.04 −25.87 7 188.5 0.2 8 37.22 −24.62
    3 62.8 0.4 8 69.56 −23.88 8 188.5 0.3 4 27.01 −15.82
    4 125.7 0.2 12 55.84 −19.59 9 188.5 0.4 12 85.16 −42.63
    5 125.7 0.3 8 51.97 −16.12
    下载: 导出CSV 
    | 显示表格

    在切削速度为62.8 m·min−1,每齿进给量为0.2 mm,切削深度为4 mm,切削步进距为0.4 mm条件下,对高度为36 mm的竹块毛坯进行切削试验,每个切削刃的部分切削分力(FxFyFz)的动态分布曲线如图2所示。切削力呈周期性变化。由于竹壁维管束纤维呈错列排布,铣削过程中维管束纤维的切削厚度不尽相同,因此加工过程中切削力的最大值变化显著。1个切削周期内,Fx先增后减,且为正值,表明在切削过程中,切削刃始终沿进给方向对竹板产生拉力,且随着切削厚度增至最大,切削力也增至最大;Fy向切削刃始终对竹板产生挤压作用,且顺铣时切削刃切出工件侧面过程中,切削厚度逐渐减小为0,故Fy始终为负值。此外,不同于螺旋铣刀,双刃直槽硬质合金木工雕刻刀在切削过程中对材板无Z轴方向作用力,但机床系统振动会导致Z轴方向作用力产生如图2C所示的微小波动[21],因此本研究忽略Z轴方向作用力的影响[2223]。为保证试验数据的一致性,均在切削力曲线中选取重复频率较高、数值较大的切削力数值进行记录。每组试验参数记录15组切削力数据,并计算其平均值作为该组正交试验参数的切削力,正交试验方案及切削力如表1所示。

    图 2  切削力动态分布
    Figure 2  Dynamic distribution of cutting force

    图3为维管束纤维切削瞬态几何模型[24]。定义纤维纵横比为α=L/amax,用以标定维管束纤维的长度和窄度,其中L为切割弧长,amax为最大纤维切削厚度。amaxL可根据下式得到:

    图 3  维管束纤维切削瞬态几何关系
    Figure 3  Transient geometry of vascular bundle fiber cutting
    $$ \mathit{a} _{ \mathrm{\max }} \mathrm= \mathit{f} _{ \mathrm{z}} \mathrm{sin} \mathit{\varphi } _{ \mathrm{0}} \mathrm{\text{;}} $$ (1)
    $$ \mathit{L} \mathrm{=2{\text{π}} } \mathit{R\varphi } _{ \mathrm{0}} \mathrm{/360。} $$ (2)

    式(1)~(2)中:fz为每齿进给量,φ0为最大接触角,R为刀具半径。定义维管束纤维理论直径为De=(L+amax)/2。从每组试验中随机选取的100根维管束纤维,使用游标卡尺分别测量纤维厚度、切割弧长和纤维长度(L0),并分别以各自平均值作为对应参数的最终数据。

    在每齿进给量为0.4 mm,切削深度为12 mm,切削步进距为0.4 mm条件下,切削速度与切削力平均值的关系如图4所示。切削过程中,Fx为正值,为主要切削力。Fy为负值,切削刃对工件产生挤压作用。随着切削速度的增加,刀具单次切削竹维管束的体积减小,单次切削功率也随之减小,因此Fx随切削速度的增加而减小;Fy在−40 N附近轻微波动。由此可见:切削速度的变化对Fy方向的切削力影响不大,其轻微波动主要由切削系统共振造成。

    图 4  切削速度对切削力影响
    Figure 4  Influence of cutting speed on cutting force

    在每齿进给量为0.4 mm,切削深度为12 mm, 切削步进距为0.4 mm条件下,切削速度对维管束纤维形态的影响如图5所示。随着切削速度的增加,维管束纤维直径(D0)先减小后增大,且数值只有维管束纤维理论直径的0.5~0.7倍(图5A)。这是因为理论切割弧长大于维管束纤维直径,即使维管束纤维被完整切割也无法达到理论值。各参数铣削获得的维管束纤维长度比较稳定,且与理论值的误差极小(图5B)。维管束纤维纵横比随着切削速度的增大而增大,且当切削速度大于125.7 m·min−1时,纤维纵横比有较明显的提高(图5C)。由纤维直径、纤维长度及纤维纵横比的变异系数(CV)可知(图5D):随着切削速度的变化,纤维长度并未产生明显的变化,纤维直径和纤维纵横比均在切削速度为160.0~188.5 m·min−1时出现最小值且变化幅度较小。综上,提高切削速度有助于获得尺寸稳定且有较大纵横比和直径的维管束纤维。

    图 5  切削速度对竹纤维形态的影响
    Figure 5  Effect of cutting speed on the shape of bamboo fiber

    在切削速度为125.7 m·min−1,切削深度为12 mm,切削步进距为0.4 mm条件下,每齿进给量与切削力平均值的关系如图6所示。随着每齿进给量增加,单位时间内切下切屑的体积增大,切削功率也随之增大,因此随着每齿进给量的增加,Fx呈近似线性增大,Fy呈近似线性减小。

    图 6  每齿进给量对切削力的影响
    Figure 6  Influence of feed rate per tooth on cutting force

    在切削速度为125.7 m·min−1,切削深度为12 mm, 切削步进距为0.4 mm条件下,每齿进给量对维管束纤维形态的影响如图7所示。不同每齿进给量下获得维管束纤维直径先减小后增大,且当每齿进给量为0.2 mm时出现最大值(图7A);各参数铣削获得的维管束纤维长度比较稳定,且与理论值接近(图7B);纤维纵横比随着每齿进给量增加,呈先下降后上升再下降的趋势,且当每齿进给量为0.2 mm时出现最大值(图7C);纤维长度、纤维直径和纤维纵横比随每齿进给量变化而产生明显的变化,但均在每齿进给量为0.2 mm时获得最优维管束纤维(图7D)。综上,减小每齿进给量有利于获得较高质量的维管束纤维。

    图 7  每齿进给量对竹纤维形态的影响
    Figure 7  Effect of feed rate per tooth on bamboo fiber morphology

    在切削速度为188.5 m·min−1,每齿进给量为0.3 mm,切削步进距为0.4 mm条件下,切削深度与切削力平均值的关系如图8所示。随着切削深度增加,单位时间内切下切屑的体积增大,切削功率也随之增大,因此Fx随切削深度的增加而增大,而Fy随切削深度的增加而减小。

    图 8  切削深度对切削力的影响
    Figure 8  Influence of cutting depth on cutting force

    在切削速度为188.5 m·min−1,每齿进给量为0.3 mm, 切削步进距为0.4 mm条件下,切削深度对维管束纤维形态的影响如图9所示。不同切削深度下获得维管束纤维直径变化不大,且为理论值的一半左右(图9A);维管束纤维长度随切削深度增加而增加,且曲线斜率为1,表明纤维长度与理论值保持一致(图9B);纤维纵横比随着切削深度增加,呈先下降后上升的趋势(图9C);图9D为切削深度对纤维形态变异系数的影响,其中纤维长度和纤维直径的变异系数基本保持恒定,纤维纵横比的变异系数先降低后升高且在切削深度大于10 mm后基本保持恒定。综上,增加切削深度有利于获得较高质量的维管束纤维。

    图 9  切削深度对竹纤维形态的影响
    Figure 9  Effect of cutting depth on the shape of bamboo fiber

    对切削力正交试验结果进行极差分析,结果如表2所示。由表2可知:切削深度(C)对FxFy的影响最大,切削速度(A)和每齿进给量(B)分别对FxFy的影响最小。因此,可以得到对FxFy影响最小的最优加工参数方案为A2B1C1

    表 2  切削力正交试验极差
    Table 2  Range of cutting force orthogonal experiment
    切削力 参数 切削参数 切削力 参数 切削参数
    切削速度(A)/
    (m·min−1)
    每齿进给量(B)/
    mm
    切削深度(C)/
    mm
    切削速度(A)/
    (m·min−1)
    每齿进给量(B)/
    mm
    切削深度(C)/
    mm
    Fx K1 174.301 116.764 85.497 Fy K1 −58.638 −53.102 −35.573
    K2 142.582 160.022 158.749 K2 −46.577 −57.819 −64.624
    K3 149.399 189.497 222.037 K3 −83.078 −77.373 −88.096
    k1 58.100 38.921 28.499 k1 −19.546 −17.701 −11.858
    k2 47.527 53.341 52.916 k2 −15.526 −19.273 −21.541
    k3 49.800 63.166 74.012 k3 −27.693 −25.791 −29.365
    R 10.573 24.244 45.514 R 12.167 8.090 17.507
    主次因素 ApfzVc 主次因素 ApVcfz
    最优方案 A2B1C1 最优方案 A2B1C1
      说明:Ki. 满足要求的单元格求和;ki. Ki的算术平均值;R. 极差。
    下载: 导出CSV 
    | 显示表格

    对切削力数据进行方差分析,结果如表3所示。在单向顺铣过程中,切削深度对Fx影响极显著(P<0.01),对Fy影响显著(P<0.05)。每齿进给量对Fx影响显著(P<0.05),对Fy影响不显著。切削速度对FxFy切削力影响均不显著。因此,单向顺铣过程中,切削深度是影响切削力的最主要因素,每齿进给量对切削力的影响大于切削速度。方差分析结果与极差分析结果一致,表明正交试验设计的有效性。

    表 3  切削力方差分析
    Table 3  Variance analysis of cutting force
    方差来源离差平方和均方F显著性方差来源离差平方和均方F显著性
    FxVc245.647122.82314.0320.067FyVc230.558115.27918.2910.052
    fz860.307430.15349.1440.020*fz110.41255.2068.7600.102
    Ap3 125.7551 562.877178.5560.006**Ap461.497230.74836.6130.027*
    误差17.5068.753误差12.6056.302
      说明:*. 差异显著(P<0.05),**. 差异极显著(P<0.01)。
    下载: 导出CSV 
    | 显示表格

    对不同切削状态下的切削数据进行多元非线性回归分析[2526]Fx非线性回归方程决定系数(R2)为0.956,说明对应模型能解释95.6%的变异,模型拟合效果很好。Fy非线性回归方程R2为0.697,说明对应模型仅能解释69.7%的变异。

    $$ {F}_{x}=37.818\;3{{A}_{\mathrm{p}}}^{0.865}{{f}_{\mathrm{z}}}^{0.688}{{V}_{\mathrm{c}}}^{-0.135}{\text{,}}{R}^{2}=0.956{\text{;}} $$ (3)
    $$ {F}_{y}=1.832\;4{{A}_{\mathrm{p}}}^{0.842}{{f}_{\mathrm{z}}}^{0.65}{{V}_{\mathrm{c}}}^{0.312}{\text{,}}{R}^{2}=0.697\mathrm{。} $$ (4)

    式(3)~(4)中:FxFy为切削分力;Ap为切削深度度;fz为每齿进给量;Vc为切削速度;R2为决定系数。正交试验与回归方程预测的对比结果如表4所示。Fx的预测值与试验值的误差为±6%左右,说明回归方程具有较高的可靠性。Fy的预测值与试验值的偏差较大,是因为切削刃对竹板反复产生挤压作用,导致切削系统振动。但对比Fy的试验值与预测值可知,预测值与理论值之间的差值在±5 N以内,且Fy的切削力数值相对较小,因此Fy的回归模型能够一定程度上反映其实际切削状态。

    表 4  正交试验参数回归方程预测值检验
    Table 4  Test of predicted value of regression equation of orthogonal experimental parameters
    序号Fx/NFy/N序号Fx/NFy/N
    试验值预测值误差/%试验值预测值误差/%试验值预测值误差/%试验值预测值误差/%
    124.3123.70−2.57−8.89−7.53−18.06633.1634.784.66−10.86−14.6625.92
    277.8181.043.99−25.87−24.71−4.69737.2537.22−0.08−24.62−19.01−29.51
    373.8169.56−6.11−23.88−21.17−12.80825.4127.015.92−15.82−13.80−14.64
    457.6355.84−3.21−19.59−23.5716.89984.0085.161.36−42.63−41.97−1.57
    549.0051.975.71−16.12−21.8026.06
    下载: 导出CSV 
    | 显示表格

    为进一步验证回归模型的可靠性,设计各切削要素条件下的非正交试验,试验结果如表5所示。回归方程对各参数的切削力预测与上述非正交试验结果基本一致,说明回归方程具有较高的可靠性。

    表 5  非正交试验参数回归方程预测值检验
    Table 5  Test of predicted values of regression equations of non-orthogonal experimental parameters
    序号切削速度/(m·min−1)每齿进给量/mm切削深度/mmFx/NFy/N
    试验值预测值误差/%试验值预测值误差/%
    162.80.3431.3332.724.25−9.80−13.7328.62
    262.80.41298.78104.345.33−29.79−35.1315.20
    3125.70.2421.5922.242.92−9.35−8.35−11.98
    4125.70.2839.3239.29−0.08−16.75−14.65−14.33
    5125.70.31273.8073.09−0.97−30.67−29.17−5.14
    6125.70.4863.3462.86−0.76−26.29−20.00−31.45
    7125.70.41289.9590.530.64−36.98−37.371.04
    8188.50.21252.8656.426.31−26.75−37.7829.20
    9188.50.3849.2046.38−6.08−24.74−22.04−12.25
    下载: 导出CSV 
    | 显示表格

    切削深度对维管束纤维切削力的影响最大,平行于刀具进给方向受力主要为切削力,每齿进给量相对于切削速度对切削力的影响更明显,垂直于刀具进给方向主要受挤压力作用,切削速度比每齿进给量对切削力的影响大。对不同切削状态下的切削数据进行多元非线性回归分析,得到切削力非线性回归模型。试验值与模型预测值的对比结果显示:该模型可以较为准确地计算各方向的切削力。

    采用铣削方法可提取长度一致性好的维管束纤维,且可通过控制铣刀的切削深度,得到不同长度的维管束纤维。采用较高切削速度、较大切削深度和较小每齿进给量有助于获得较大长径比和直径的维管束纤维。

  • 图  1  各样地不同深度根径d<1mm 的根密度

    Figure  1  Distribution tend of root number diameter less than 1 mm in different soil depth

    图  2  各样地土壤冲刷量随时间的变化

    Figure  2  Relationship between soil scoured mass and time

    表  1  各标准地活立木数量

    Table  1.   Numbers of living trees in sample plots

    样地沙棘/株榆树/株沙棘死亡率/%
    F11723
    F21432216.86
    F3834151.74
    F4495971.51
    F587695.35
    下载: 导出CSV

    表  2  各样地针式抗冲仪测定的抗冲性

    Table  2.   Anti-scouring abilities with needle-type instrument

    土层深度/cm各样地针式抗冲仪测定的抗冲性/(g·cm-2mm-1)
    F1F2F3F4F5
    0〜106.25±0.913 33±0.482.68±0.292 38±0.241 79±0.19
    10〜204.00±0.752 98±0.352 14±0.271.89±0.211.56±0.20
    20〜302.78±0.592 24±0.361 79±0.281.39±0.191.59±0.18
    30〜402.04±0.341 76±0.271.59±0.161 21±0.131.11±0.19
    平均3.77±1.84 aA2 58±0.71 bB2.04±0.48 cBC1 72±0.53 cdBC1.51±0.29 dC
    下载: 导出CSV

    表  3  各样地槽式抗冲仪测定的抗冲指数

    Table  3.   Anti-scouring abilities with trough-type instrument

    深度/cm槽式抗冲仪测定的抗冲指数/(s·g-1)
    F1F2F3F4F5
    0〜1014.42±1.0612.39±1.1511.56±1.268.44±1.276.52±0.68
    10-2014.13±0.9812.15±1.0910.90±1.648.04±1.236.40±0.72
    20-3013.86±1.2012.01±1.3510.64±1.347.77±1.126.27±0.59
    30-4013.59±2.0111.77±1.2710.04±1.287.40±1.036.15±0.67
    平均14.00±0.36 A12.08±0.26 B10.78±0.63 C7.91±0.44 D6.34±1.62 E
    下载: 导出CSV

    表  4  土壤抗冲性(%)与根径d<1mm根密度(x)的幂函数(y=axb)参数

    Table  4.   Parameters for power functions between anti-scouring abilities and root density which diameter less than 1 mm

    样地针式抗冲仪槽式抗冲仪
    abR2PabR2P
    F10.530.400.990.00312.680.0220.960.008
    F20.730.270.920.00111.190.0200.960.003
    F31.000.150.980.0058.680.0500.950.009
    F40.390.310.980.0026.240.0500.960.004
    F50.840.140.920.0105.820.0200.960.020
    下载: 导出CSV
  • [1] 查小春,贺秀斌. 土壤物理力学性质与侵蚀关系研究进展[J]. 水土保持研究,1999,6(2):98-104.

    ZHA Xiaochun,HE Xiubin. Progress on the study of soil erosion from soil physical property and erosive force[J]. Res Soil Water Convers,1999,6(2):98-104.
    [2] 朱显谟. 黄土高原植被因素对于水土流失的影响[J]. 土壤学报,1960,8(2):110-120.

    ZHU Xianmo. Effect of vegetation on soil and water loss in Loess Plateau[J]. Acta Pedol Sin,1960,8(2):110-120.
    [3] 李勇,徐晓琴,朱显谟,等. 植物根系与土壤抗冲性[J]. 水土保持学报,1993,7(3):11-18.

    LI Yong,XU Xiaoqin,ZHU Xianmo,et al. Plant roots and anti-scouring ability[J]. J Soil Water Conserv,1993,7(3):11-18.
    [4] 宫海志,王军. 建平县沙棘林大面积死亡成因分析及经营对策[J]. 中国水土保持,2003(4):17-18.

    GONG Haizhi,WANG Jun. Analysis on reasons of large-acreage death of Seabuckthorn forest and its operation strategies in Jianping County[J]. Chin Soil Water Conserv,2003(4):17-18.
    [5] 李树彬,党福江. 建平县沙棘林大面积死亡原因调查分析[J]. 水土保持科技情报,2001(6):4-6.

    LI Shubin,DANG Fujiang. Investigation and analysis on reasons of large-acreage death of Seabuckthorn forest in Jianping County[J]. Sci Technol Inform Soil Water Conserv,2001(6):4-6.
    [6] 惠兴学,张连翔,孔繁轼,等. 建平县沙棘林大面积死亡成因调查分析及对策[J]. 防护林科技,2002(2):53-55.

    HUI Xingxue, ZHANG Lianxiang,KONG Fanshi,et al. Investigation and analysis on reasons of large-acreage death of Seabuckthorn forest and its control in Jianping County[J]. Prot For Sci Technol,2002(2):53-55.
    [7] 张连翔,惠兴学,黄立华,等. 建平县沙棘林大面积死亡原因及其治理对策[J]. 沙棘,2002,15(3):26-29.

    ZHANG Lianxiang,HUI Xingxue,HUANG Lihua,et al. A large-acreage death of Seabuckthorn forest and its control in Jianping County[J]. Hippophae,2002,15(3):26-29.
    [8] 闫芊,何文珊,陆健健. 崇明东滩湿地植被演替过程中生物量与氮含量的时空变化[J]. 生态学杂志,2006,25(9):1019-1023.

    YAN Qian,HE Wenshan,LU Jianjian. Temporal and spatial changes of plant biomass and its N content during vegetation succession at Chongming Dongtan wetland[J]. Chin J Ecol,2006,25(9):1019-1023.
    [9] 宋洪涛,张劲峰,田昆,等. 滇西北亚高山地区黄背栎林植被演替过程中的林地土壤化学响应[J]. 西部林业科学,2007,36(2):65-70.

    SONG Hongtao,ZHANG Jinfeng,TIAN Kun,et al. Response of soil chemical properties to vegetation succession of Quercus pannosa forest in subalpine area of northwest Yunnan[J]. J West China For Sci,2007,36(2):65-70.
    [10] 王韵,王克林,邹冬生,等. 广西喀斯特地区植被演替对土壤质量的影响[J]. 水土保持学报,2007,21(6):130-134.

    WANG Yun,WANG Kelin,ZOU Dongsheng,et al. Effects of vegetation succession on soil quality in Karst region of Guangxi,China[J]. J Soil Water Conserv,2007,21(6):130-134.
    [11] 周正朝,上官周平. 子午岭次生林植被演替过程的土壤抗冲性[J]. 生态学报,2006,26(10):3270-3275.

    ZHOU Zhengchao,SHANGGUAN Zhouping. Soil anti-scouring ability during vegetation succession of Ziwuling secondary forest[J]. Acta Ecol Sin,2006,26(10):3270-3275.
    [12] 李勇,徐晓琴,朱显谟,等. 草类根系对土壤抗冲性的强化效应[J]. 土壤学报,1992,29(3):302-309.

    LI Yong,XU Xiaoqin,ZHU Xianmo,et al. Strengthening effects of grass root on soil anti-scouring ability[J]. Acta Pedol Sin,1992,29(3):302-309.
    [13] DILLAHA T A. Vegetative filter strips for agricultural nonpoint-source pollution control[J]. Tran ASAE,1989,32(3):513-519.
    [14] CHAUBEY I,EDWARDS D R,DANIEL T C,et al. Effectiveness of vegetative filter strips in controlling losses of surface-applied poultry litters constituents[J]. Tran ASAE,1995,38(6):1687-1692.
    [15] GONZALEZ H J C,RAVENTOS J,ECHEVARRIA M T. Comparison of sediment ratio curves for plants with different architectures[J]. Catena,1997,29(3/4):333-340.
  • [1] 杨馥羽, 陈奇伯, 黎建强, 龚顺清, 付钇珊.  计划烧除对云南松林土壤抗蚀和抗冲性的影响 . 浙江农林大学学报, 2023, 40(1): 188-197. doi: 10.11833/j.issn.2095-0756.20220321
    [2] 彭健健, 王增, 张勇, 刘海英, 顾光同, 彭欣怡, 吴家森, 叶子豪, 张申, 尚世宇.  杨梅人工林相容性单株生物量模型构建 . 浙江农林大学学报, 2022, 39(2): 272-279. doi: 10.11833/j.issn.2095-0756.20210272
    [3] 魏继华, 李佳益, 刘宏, 张建国, 罗红梅, 何彩云.  沙棘根瘤内生菌株库构建与微生物多样性分析 . 浙江农林大学学报, 2022, 39(2): 356-363. doi: 10.11833/j.issn.2095-0756.20210246
    [4] 方伟, 余晓, 王晶, 徐秋芳, 梁辰飞, 秦华, 陈俊辉.  施加石灰石粉和微生物肥料对发病山核桃林土壤化学性质和微生物群落的影响 . 浙江农林大学学报, 2020, 37(2): 273-283. doi: 10.11833/j.issn.2095-0756.2020.02.011
    [5] 郭益昌, 庄舜尧, 胡昱彦, 桂仁意.  埋管通气对雷竹林土壤氧气体积分数的影响 . 浙江农林大学学报, 2020, 37(1): 69-75. doi: 10.11833/j.issn.2095-0756.2020.01.009
    [6] 姜仕昆, 周运超, 谭伟, 陈竹, 黄剑峰.  马尾松林近自然不同经营管理措施下土壤肥力 . 浙江农林大学学报, 2020, 37(5): 876-882. doi: 10.11833/j.issn.2095-0756.20190549
    [7] 李智超, 张勇强, 厚凌宇, 宋立国, 孙启武.  杉木人工林土壤微生物对林分密度的响应 . 浙江农林大学学报, 2020, 37(1): 76-84. doi: 10.11833/j.issn.2095-0756.2020.01.010
    [8] 何姗, 刘娟, 姜培坤, 周国模, 王会来, 李永夫, 吴家森.  经营管理对森林土壤有机碳库影响的研究进展 . 浙江农林大学学报, 2019, 36(4): 818-827. doi: 10.11833/j.issn.2095-0756.2019.04.023
    [9] 许俊丽, 张桂莲, 张希金, 高志文, 仲启铖, 张亚萍, 宋坤, 达良俊.  上海市人工林土壤理化性质与群落特征的相关性 . 浙江农林大学学报, 2018, 35(6): 1017-1026. doi: 10.11833/j.issn.2095-0756.2018.06.004
    [10] 李晨晨, 周再知, 梁坤南, 黄桂华, 杨光.  不同林药复合经营模式对杉木生态公益林土壤理化性质的改良效果 . 浙江农林大学学报, 2018, 35(1): 51-59. doi: 10.11833/j.issn.2095-0756.2018.01.007
    [11] 竹万宽, 陈少雄, RogerARNOLD, 王志超, 许宇星, 杜阿朋.  不同种桉树人工林土壤呼吸速率时空动态及其影响要素 . 浙江农林大学学报, 2018, 35(3): 412-421. doi: 10.11833/j.issn.2095-0756.2018.03.004
    [12] 辜忠春, 李光荣, 李军章, 杜业云, 王宵.  正交试验优化分光光度法测定森林土壤有机质 . 浙江农林大学学报, 2017, 34(2): 239-243. doi: 10.11833/j.issn.2095-0756.2017.02.006
    [13] 高国日, 刘娟娟, 陈道国, 董礼隆, 何彩云.  土壤干旱对2个沙棘品种叶片水势和光合特性的影响 . 浙江农林大学学报, 2017, 34(6): 999-1007. doi: 10.11833/j.issn.2095-0756.2017.06.005
    [14] 孙鹏跃, 徐福利, 王渭玲, 王玲玲, 牛瑞龙, 高星, 白小芳.  华北落叶松人工林地土壤养分与土壤酶的季节变化及关系 . 浙江农林大学学报, 2016, 33(6): 944-952. doi: 10.11833/j.issn.2095-0756.2016.06.004
    [15] 朱仁欢, 李玮, 郑子成, 李廷轩, 洪月, 何秋佳, 田宗渠.  退耕植茶地土壤碳氮磷生态化学计量学特征 . 浙江农林大学学报, 2016, 33(4): 612-619. doi: 10.11833/j.issn.2095-0756.2016.04.009
    [16] 徐莉, 周伟, 俞元春, 王玫, 侯文军, 徐长柏, 戴成.  沼液施用对杨树林地土壤性质及林分生长的影响 . 浙江农林大学学报, 2015, 32(2): 204-207. doi: 10.11833/j.issn.2095-0756.2015.02.006
    [17] 郭帅, 徐秋芳, 沈振明, 李松昊, 秦华, 李永春.  雷竹林土壤氨氧化微生物对不同肥料的响应 . 浙江农林大学学报, 2014, 31(3): 343-351. doi: 10.11833/j.issn.2095-0756.2014.03.003
    [18] 叶玲燕, 傅伟军, 姜培坤, 李永夫, 张国江, 杜群.  浙江省森林表层土壤基本化学性质和有机碳储量的空间变异 . 浙江农林大学学报, 2012, 29(6): 803-810. doi: 10.11833/j.issn.2095-0756.2012.06.001
    [19] 叶耿平, 刘娟, 姜培坤, 周国模, 吴家森.  集约经营措施对毛竹林生长季土壤呼吸的影响 . 浙江农林大学学报, 2011, 28(1): 18-25. doi: 10.11833/j.issn.2095-0756.2011.01.004
    [20] 李正才, 傅懋毅, 杨校生.  经营干扰对森林土壤有机碳的影响研究概述 . 浙江农林大学学报, 2005, 22(4): 469-474.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2014.03.008

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2014/3/380

图(2) / 表(4)
计量
  • 文章访问数:  3781
  • HTML全文浏览量:  684
  • PDF下载量:  523
  • 被引次数: 0
出版历程
  • 收稿日期:  2013-09-21
  • 修回日期:  2013-10-15
  • 刊出日期:  2014-06-20

沙棘林退化对林木根系和土壤抗冲性的影响

doi: 10.11833/j.issn.2095-0756.2014.03.008
    基金项目:

    辽宁省自然科学基金资助项目 20102109,20120210

    重庆市水利局重大项目 渝水科201306

    作者简介:

    舒乔生,副研究员,博士,从事水土保持科学研究。E-mail:cysqs@163.com

    通信作者: 谢立亚,副研究员,从事水土保持科学研究。E-mail:xly2004@126.com
  • 中图分类号: S714.5;S157.1

摘要: 近年来,辽宁建平县亚洲最大的人工沙棘Hippophae rhamnoides林呈现大面积死亡的现象。通过标准样地法调查沙棘的死亡程度,测定与对比不同沙棘死亡率下林地的根系密度和土壤抗冲性,研究沙棘林退化对水土流失的影响。结果表明:0~30 cm范围内,直径d<1 mm的根系密度随深度增大急剧减小,72%~80%根量分布在0~20 cm土层;0~30 cm范围内,随着沙棘死亡率的提高,d<1 mm根密度显著降低(P<0.05);随着深度增大,抗冲性也逐步降低,沙棘的死亡率对土壤抗冲性的影响亦达到显著水平(P<0.05),这主要是由根系密度的相应变化引起的;经曲线拟合,土壤抗冲性与d<1 mm根密度的关系可用幂函数表示。目前该区主要应通过采取有效措施降低沙棘林退化的生态影响。

English Abstract

舒乔生, 谢立亚. 沙棘林退化对林木根系和土壤抗冲性的影响[J]. 浙江农林大学学报, 2014, 31(3): 380-385. DOI: 10.11833/j.issn.2095-0756.2014.03.008
引用本文: 舒乔生, 谢立亚. 沙棘林退化对林木根系和土壤抗冲性的影响[J]. 浙江农林大学学报, 2014, 31(3): 380-385. DOI: 10.11833/j.issn.2095-0756.2014.03.008
HUANG Haiyan, TANG Huimin, JIN Lu, et al. Effects of tree species importance, diversity and soil physicochemical properties on soil amino sugars[J]. Journal of Zhejiang A&F University, 2024, 41(4): 778-786. DOI: 10.11833/j.issn.2095-0756.20230527
Citation: SHU Qiaosheng, XIE Liya. Root system distribution and soil anti-scouring capability with Hippophae rhamnoides degradation[J]. Journal of Zhejiang A&F University, 2014, 31(3): 380-385. DOI: 10.11833/j.issn.2095-0756.2014.03.008
  • 对于降雨侵蚀力相当的同一地区来说,水力侵蚀的强弱主要取决于土壤本身特性及地形、植被等影响因子。土壤抗冲性是表征土壤抗蚀能力的重要指标,它指土壤抵抗径流的机械破坏和推动下移的能力,抗力来源于土粒间、微结构间的胶结力及其抵抗离散力[1]。朱显谟[2]研究表明: 土壤抗冲性取决于根系的缠绕和固结作用,这种作用使土体有较高的水稳结构和抗蚀强度,从而不易被径流带走。近年来,由于人口剧增,陡坡毁林毁草、垦荒耕种,破坏了植物根系对土体的固结作用,导致土壤抗冲性降低,是加剧全球性土壤侵蚀、生态环境恶化的一个重要内在原因[3]。辽宁西部的建平县人工沙棘Hippophae rhamnoides林面积为6.67万hm2,是全国乃至亚洲最大的人工沙棘林[4],对保持水土、改善生态环境、促进经济发展等方面起到了积极的作用。近年来,由于自然老化,过度干旱,病虫危害,经营管理不善等原因[5-6],该区沙棘死亡率达58.5%[4, 7],生态环境遭到严重破坏。然而关于沙棘林死亡对水土保持影响的定量评价及机制仍没有进行系统的研究。本研究通过样地测试,分析沙棘林林分退化对林地根系和土壤抗冲性的影响。

    • 辽宁省建平县地处辽西低山丘陵区内,属北温带大陆季风气候区,多年平均气温为8.4 ℃,年平均降水量 450.0 mm,其中6-9月占70%,是辽宁省最严重的水土流失地区之一。

      试点位于建平县建平镇胡家店村红金台小流域,土壤为黄土母质发育的褐土,水土保持树种主要有沙棘Hippophae rhamnoides,榆树Ulmus pumila,油松Pinus tabulaeformis,柠条Caragana korshinskii等,属辽西低丘区的典型小流域。该流域纯沙棘林有20年生以上的植株,出现逐步死亡的现象,取而代之是榆树的入侵,表现出近似演替的特征。

    • 根据植被演替研究的通用方法——时空替代法[8-10],在空间上选择纯沙棘林(未死亡,F1),榆树入侵期(沙棘仍为优势种,F2),榆树发展期(沙榆势均力敌,F3),榆树优势期(有少量沙棘,F4),纯榆树林(F5)5个标准地,3次重复,分别作为沙棘死亡及近似演替的5个阶段,标准地尺寸为20 m × 20 m,各标准地平均坡度为7°,死亡的沙棘均被农民搬运移出。各标准地活立木数量见表 1

      表 1  各标准地活立木数量

      Table 1.  Numbers of living trees in sample plots

      样地沙棘/株榆树/株沙棘死亡率/%
      F11723
      F21432216.86
      F3834151.74
      F4495971.51
      F587695.35
    • 在每个标准地内,挖长2.0 m,宽2.0 m,深1.0 m的土壤剖面,对根系分布进行调查,通过土壤取样并对根系进行冲洗,利用游标卡尺将根系按直径大小分为d<1 mm,1≤d<2 mm,2≤d<5 mm,5≤d<10 mm和d≥l0 mm共5个径级,统计各径级根系的条数,以1 dm2中根的条数表示根系密度。

    • ①针式抗冲仪:采用NLKC-90 针式土壤抗冲仪,在0~40 cm深度范围内按照10 cm间隔,在101.325 kPa压力(1个标准大气压)下,以0.7 mm直径的水柱冲击土层1 min,使其产生水蚀穴,每10个水蚀穴深(cm)与宽(cm)乘积均值的倒数,即为该土层的抗冲指数。②槽式抗冲仪:在0~40 cm深度范围内以10 cm间隔,将每个土壤剖面用15 cm × 30 cm × 10 cm的取样器在土壤剖面取土样,带回室内置于底部有少量水的盆内浸泡12 h(水面不能高于土表面),使土柱自然吸水达到饱和状态。利用NLKC-101槽式抗冲仪将土样坡度调整为10°,在101.325 kPa水压力下对土样进行冲刷,对泥沙过滤、烘干称量,将冲刷时间除以冲刷的土量,即冲刷单位质量的土所需时间,作为抗冲性度量指标(单位为s·g-1)。

    • 从各径级的根系分布情况来看,径级为1≤d<2 mm,2≤d<5 mm,5 ≤d<10 mm和d≥l0 mm的根密度分布较小,因此将根系分为直径d≥1 mm和d<1 mm等2种径级进行分析。

    • 经统计(详细结果未列出),直径d≥1 mm根密度随土层深度增加逐步减少,0~40 cm土层内递减较快,在40 cm以下土层递减较慢;不同沙棘死亡率下直径d≥1 mm的根量分布差异不显著(P>0.05)。

    • 直径d<1 mm根密度的分布见图 1,在0~30 cm范围内,根密度随深度增大急剧减少,30~100 cm范围内逐渐减少并维持在极低水平,总根数的72%~80%分布在0~20 cm土层,20~40 cm土层内根系仅占总根数的10%~20%。这主要是由于该区缺水非常严重,且有愈演愈烈之势;加上该区为黄土状母质上发育的淋溶褐土,表土质地粗孔隙大,而表层以下为质地较为坚硬的黄土,结构不良,渗水性能差,极易形成超渗产流;同时该区地下水埋藏极深,沙棘根系无法吸收深层地下水,其生长主要靠天然降水。因此沙棘的根系主要位于浅层土壤中,随着深度增加,根系急剧减少,这也是沙棘出现大面积死亡的原因之一。分析表明,在地表下0~30 cm深度范围内,随着沙棘死亡率的提高,根径d<1 mm的根密度呈显著的降低趋势(P<0.05),而大于30 cm的区域,差异不显著(P>0.05),表明沙棘死亡会导致表层土壤的根系分布变得稀疏。若按照植被演替的规律,沙棘死亡和榆树生长意味着植被的正向演替,根密度会随之增大[11]。然而本试验却出现了相反的现象,这主要是由于沙棘林对浅层根系分布起决定作用,而榆树数量较少且根系分布相对较深,对浅层根系起次要作用。

      图  1  各样地不同深度根径d<1mm 的根密度

      Figure 1.  Distribution tend of root number diameter less than 1 mm in different soil depth

    • 针式土壤抗冲仪测定结果见表 2。结果表明:随着沙棘死亡率的提高,土壤抗冲性呈逐步降低趋势,纯沙棘林与其他样地差异为极显著(P<0.01),除榆树优势期外,其他3个样地之间的差异均达到显著水平(P<0.05),说明沙棘死亡对土壤的抗冲性具有显著的影响;并且随着土层深度的增大,抗冲性也逐步降低。这主要是由于表层根系分布密度大和枯枝落叶层较厚,虽然土壤较疏松,但根系对土壤的固持作用较强;随着深度增加,根径d<1 mm的根密度迅速减少,致使土壤缺乏根系对集中水流的缓冲作用,冲刷作用较强[12]

      表 2  各样地针式抗冲仪测定的抗冲性

      Table 2.  Anti-scouring abilities with needle-type instrument

      土层深度/cm各样地针式抗冲仪测定的抗冲性/(g·cm-2mm-1)
      F1F2F3F4F5
      0〜106.25±0.913 33±0.482.68±0.292 38±0.241 79±0.19
      10〜204.00±0.752 98±0.352 14±0.271.89±0.211.56±0.20
      20〜302.78±0.592 24±0.361 79±0.281.39±0.191.59±0.18
      30〜402.04±0.341 76±0.271.59±0.161 21±0.131.11±0.19
      平均3.77±1.84 aA2 58±0.71 bB2.04±0.48 cBC1 72±0.53 cdBC1.51±0.29 dC

      槽式抗冲仪测定的土壤冲刷量见图 2。从图 2可以看出:冲刷量随着时间的延续呈增加趋势,0~120 s之间冲刷量增加较快,随后增速变缓,至300 s时达到较平稳状态。同时随着沙棘死亡率的升高,土壤冲刷量逐步呈增加趋势,纯榆树林的冲刷量远高于其他样地,差异为极显著(P<0.01),榆树发展期与榆树入侵期之间差异不显著。通过计算单位冲刷量所需的时间,得到抗冲指数见表 3

      图  2  各样地土壤冲刷量随时间的变化

      Figure 2.  Relationship between soil scoured mass and time

      表 3  各样地槽式抗冲仪测定的抗冲指数

      Table 3.  Anti-scouring abilities with trough-type instrument

      深度/cm槽式抗冲仪测定的抗冲指数/(s·g-1)
      F1F2F3F4F5
      0〜1014.42±1.0612.39±1.1511.56±1.268.44±1.276.52±0.68
      10-2014.13±0.9812.15±1.0910.90±1.648.04±1.236.40±0.72
      20-3013.86±1.2012.01±1.3510.64±1.347.77±1.126.27±0.59
      30-4013.59±2.0111.77±1.2710.04±1.287.40±1.036.15±0.67
      平均14.00±0.36 A12.08±0.26 B10.78±0.63 C7.91±0.44 D6.34±1.62 E

      结果表明:无论是用水穴尺寸还是用冲刷量表示,随着沙棘的死亡和榆树生长,土壤抗冲性均呈逐步降低的趋势,这与周正朝等[11]植被演替增强土壤抗冲性的结论不一致。这主要由于土壤抗蚀性主要通过沙棘的根系起作用,随着沙棘的死亡,土壤中的根系数量减少,加上农民及时将死亡的沙棘移出,死亡的沙棘又无法变成腐殖质进一步保护土壤。虽然沙棘的死亡伴随着榆树的生长,单株乔木根系固结土壤和阻止径流的能力较强[13-14],然而水土流失的整体防治效果与植被结构密切相关,层次结构较好的植被对土壤的保护作用要优于单层植被[15];因此,密度较小的榆树林木对土壤的固结作用仅限于其根系周边,无法起到沙棘全面保护土壤的作用。另外,由于榆树根系较沙棘深,它更侧重于保护深层次的土壤,而土壤抗冲性主要体现在表层或近表层土壤,因此在这种情况下,榆树对土壤抗冲性的提高要逊于沙棘。

      通过土壤抗冲性与根密度进行曲线拟合,表明抗冲性与根径d<1 mm根密度之间可用幂函数(y=axb)表示,与李勇等[12]的研究结论相类似。各样地参数a和b表 4

      表 4  土壤抗冲性(%)与根径d<1mm根密度(x)的幂函数(y=axb)参数

      Table 4.  Parameters for power functions between anti-scouring abilities and root density which diameter less than 1 mm

      样地针式抗冲仪槽式抗冲仪
      abR2PabR2P
      F10.530.400.990.00312.680.0220.960.008
      F20.730.270.920.00111.190.0200.960.003
      F31.000.150.980.0058.680.0500.950.009
      F40.390.310.980.0026.240.0500.960.004
      F50.840.140.920.0105.820.0200.960.020
    • 通过对辽西建平县沙棘林死亡率的调查,并对不同死亡程度的林地进行土壤根系和土壤抗冲性的测定,研究沙棘林退化对水土流失的影响。结果表明:随着沙棘死亡率的提高,根径d<1 mm根密度逐步降低,沙棘死亡对较粗根系分布影响不显著;沙棘林具有较好的防治水土流失作用,随着沙棘活立木的减少,土壤抗蚀性呈逐步降低趋势,纯沙棘与其他时期差异为极显著;随着土层深度的增大,抗冲性逐步降低,土壤抗冲性与根径d<1 mm根系密度之间的关系可用幂函数表示。林地根系分布的减少和抗蚀性的减弱,水土流失使土壤结构产生破坏,因此土壤的水肥条件必然进一步恶化,从而加快沙棘的死亡,形成恶性循环。

      近年来,当地林业和水土保持部门对人工沙棘林的退化现象非常重视,但一直仅就其林业方面的损失进行评估,对其引起的生态问题未进行系统的观测。本项研究从根系分布和土壤抗冲性入手,探讨沙棘退化对水土流失的影响,可以作为生态损失评估的基础。由于该区自然条件的恶劣性,加上沙棘林本身的老化现象,沙棘退化的趋势不可避免,需对其产生的影响进行全面评估,同时针对当地不良水分条件、土壤特性,以及稀疏乔木无法快速发挥水土保持效益的现状,采取措施抑制其退化速度或补植其他水土保持效果较好的灌木,将沙棘退化的生态影响降到最低。

参考文献 (15)

目录

/

返回文章
返回