Remediation Effect of Biological Agents on Rhizosphere Soil on Land of Continuously Cropped Panax notoginseng
-
摘要: 本研究探讨有益微生物复合菌剂（Effective Micoorganisms，EM）与木霉菌剂（Trichoderma，TRI）对三七连作障碍的消减作用及其作用机理。结果表明：EM和TRI处理的重茬三七出苗率、保苗率、产量均显著提高；EM处理的重茬土壤中的各营养指标均有所提高，TRI处理的重茬土壤中的速效钾等营养指标显著提高；HPLC分析发现，重茬三七根际土壤的酚酸含量未出现累积效应，可见三七连作障碍的发生并非是由化感物质的直接作用造成的，可能与化感物质对根际环境的间接作用有关；变性梯度凝胶电泳（DGGE）、荧光定量PCR（qRT-PCR）原位分析验证结果发现两种菌剂处理均能使连作土壤中的有益菌含量增加，有害菌含量减少，其中TRI的处理效果更为显著。说明这两种菌剂处理均可有效修复三七根际土壤，缓解连作障碍，同时显著提高三七产量，达到了预期效果。Abstract: Two biological agents, EM and TRI, were applied in soil on a continuous cropping field of Panax notoginseng to study the remediating effect on the soil and the mechanism associated with the application. Treatments with the agents were found significantly increased the seedling emergence rate and survival rate as well as the root yield. All nutrients in the soil treated with EM were increased, while the available potassium in the soil treated with TRI significantly improved. The HPLC analysis showed no accumulation of phenolic acids in the rhizosphere soil by these treatments indicating that allelochemicals were not directly, but might indirectly, involved in the detrimental effect on P. notoginseng caused by continuous cropping. The DGGE and qRT-PCR test results showed that both EM and TRI, especially TRI, could increase the population of beneficial, while depressed the harmful, bacteria in the soil. Thereby, the rhizosphere soil was effectively repaired by alleviating the deficiency brought about by the continuous cropping practice, and consequently, the growth and yield of P. notoginseng was significantly recovered or improved.
-
耕地资源是保障粮食安全的重要战略资源，是人类赖以生存的根基。由于我国人口众多、人均耕地少，耕地生态系统在人类长期活动下，日渐演变成一个社会-经济-生态的高度耦合的复合系统[1-2]。在城镇化进程加快、工业化水平提高的压力背景下，人类与耕地资源间的矛盾日益尖锐，建设用地占用耕地、工业废弃物污染土壤等问题使得耕地生态环境问题日益凸显。因此，开展耕地生态安全评价研究，对于改善耕地生态环境，协调人地关系和保障粮食安全具有重要的现实意义。
在生态文明的大前提下，耕地生态安全是指基于耕地资源可持续利用，既能满足当代人对耕地资源的生活所需及合理利用，又能使子孙后代从中受益的一种状态[3-4]。国外的研究范围主要是将耕地生态安全与区域生态安全相结合，单独分析耕地生态安全则较少[5-6]。目前，国内学者的研究范围主要涵盖耕地生态安全评价、影响因素、时空演变特征以及调控政策等方面[7-11]。郭荣中等运用改进的PSR模型结合组合赋权方法构建评价体系，定量评价了长株潭地区耕地生态安全状态[12]。杨俊等以长江中下游粮食主产区为评价对象，基于PSR模型选取耕地压力指数等36个评价指标，采用综合评价值计量等方法探索研究区域耕地生态安全现状及影响因素的内在机理[13]。张锐等在明确耕地生态安全内涵之后，针对传统综合评价研究中存在评价信息缺失和评价主观性的弊端，采用改进的熵值法结合物元模型评价我国耕地生态安全[14]。王军通过分析石家庄市近10年来的耕地数量、质量和生态环境动态变化，运用人工神经网络法评价石家庄市1998年、2002年和2007年的耕地生态安全状况[15]。赵文晖等采用可拓优度评价模型分析评价了河北省保定市2003-2012年的耕地生态安全状况[16]。综上，耕地生态安全评价研究方法主要有：人工神经网络法、综合评价法、物元分析法等，在研究尺度上主要以国家、省、市级尺度较多，而从山区县域层面展开耕地生态安全评价的研究甚少。
2014年福建省被确立为全国首个生态文明先行示范区，可见福建生态省的建设已然上升为国家战略[17]。但福建省近九成陆地面积为丘陵山地，地形复杂，耕地碎片化，素有“八山一水一分田”之称，清流县作为福建省山区农业县的典型代表，所面临的耕地生态环境问题具有显著普遍性。本文在前人的研究基础上，从耕地现状、社会经济和生态环境三方面综合提取评价指标，建立基于熵权法确权和物元模型的耕地生态安全评价模型，评价福建省清流县在2010-2016年的耕地生态安全状况，以期为山区县域耕地生态安全评价和县域耕地保护政策提供参考依据。
1. 研究区概况与数据来源
1.1 研究区概况
清流县位于福建省西部，北纬25°46′～26°22′、东经116°38′～117°10′，属亚热带季风气候，兼具大陆性和海洋性气候特点，秋冬季多为西北风，春夏季则为东南风，全年气候温暖而湿润。全县土地总面积为1 806.3 km2，下辖5个镇、8个乡，境内低山丘陵分布广泛，地势从西向东逐级提高，平均海拔为350 m，是福建省典型的山区农业县。
截至2016年末，全县总人口达155 994人，地区生产总值为89亿元，固定资产投资额达到98.8亿元。耕地面积为1.32万hm2，人均耕地面积0.085 1 hm2，坡度大于25度耕地面积比重为8.2%。年内增加耕地57.12 hm2，减少耕地33.54 hm2。
1.2 数据来源
耕地数据主要来自清流县国土局，人口密度和人口自然增长率来源于清流县卫计局，GDP年增长率等社会经济数据以及化肥施用量、森林覆盖率等生态环境数据来源于《清流县统计年鉴》(2011-2017)和《三明市统计年鉴》(2011-2017)。
2. 生态安全评价指标体系构建
2.1 指标选取
明确耕地生态安全概念前提下，参照国家环保总局制定的《生态县、生态市、生态省建设指标(试行)》关于生态县建设指标[18]，借鉴国内已有的相关成果，同时依据各评价指标数据的获取难易程度，从耕地现状、社会经济、生态环境3个方面综合选取了人均耕地面积、农民人均纯收入年增长率、GDP年增长率、污水处理率等16个评价指标。
2.2 权重确定
权重反映了评价指标对被评价对象的重要程度，评价结果的科学与否离不开确权的客观程度。确权的方法包括主观和客观两个方面，涵盖包含层次分析法、德尔菲法、熵权法在内的10多种方法，相较于德尔菲专家评分法，熵权法避免了由于主观因素对赋权的干扰，结果较为客观，因此采用熵权法对指标赋权[19]，通过熵权法公式计算，得到赋权结果见表 1。
表 1 耕地生态安全评价指标体系及权重Table 1. Weighted indicators in evaluating ecological security of cultivated land in Qingliu county因素 指标 公式 权重 影响 耕地现状 X1人均耕地面积/(hm2·人-1) 耕地面积/总人口 0.0837 (+) X2耕地减少面积与耕地面积比例/% 耕地减少面积/耕地面积 0.0425 (-) X3耕地增加面积与耕地面积比例/% 耕地增加面积/耕地面积 0.2017 (+) X4大于25度耕地面积比重/% - 0.0705 (-) 社会经济 X5人口密度/(人·km-2) 总人口/土地总面积 0.0176 (-) X6人口自然增长率/‰ - 0.0778 (-) X7农民人均纯收入年增长率/% - 0.0941 (+) X8 GDP年增长率/% - 0.0520 (+) X9单位面积耕地粮食产量/(kg·hm-2) 粮食总产量/粮食作物播种面积 0.0246 (+) X10固定资产投资总额年增长率/% - 0.0806 (+) 生态环境 X11森林覆盖率/% - 0.0603 (+) X12单位耕地面积化肥施用量/(kg·hm-2) 化肥施用量/耕地总面积 0.0725 (-) X13单位耕地面积农药施用量/(kg·hm-2) 农药施用量/耕地总面积 0.0766 (-) X14工业固体废物综合利用率/% - 0.0213 (+) X15污水处理率/% - 0.0042 (+) X16生活垃圾处理率/% - 0.0200 (+) 3. 评价模型及确定经典域
20世纪80年代初，为研究事物拓展的可能性和解决不相容的复杂问题，蔡文教授提出了物元模型，同时这种模型也适用于多指标评价过程中。耕地生态系统是由社会、经济、生态间相互作用而成的多因素复杂系统，涉及评价指标较多且指标间的内在关系不明确，因此可将其看作一个不相容的问题，采用物元分析方法构建评价模型[20]。
3.1 构建物元模型
3.1.1 确定耕地生态安全物元
依据物元模型的相关概念，描述耕地生态安全的基本元由耕地生态安全Z，特征值X以及量值v以有序的三元组N=(Z, X, v)构成。其中，Y=X(N)反映事物间质和量的关系。若耕地生态安全Z有n个特征X1, X2，…，Xn和对应的模糊量值为v1，v2，…，vn，则称作n维模糊物元。相应矩阵为：
N=|ZX1v1X2v2⋯⋯Xnvn|=|N1N2⋯Nn| (1) 3.1.2 确定经典域和节域物元矩阵
以Z0表示所划分耕地生态安全的评价等级，X表示评价指标，各评价指标对应等级的量值范围称谓经典域，即v0=(u0i，v0i)(i=1, 2, …，n)。经典域的物元矩阵可表示为：
N0=(Z0,v0)=|Z0X1(u01,v01)X2(u02,v02)⋯⋯Xn(u0n,v0n)| (2) 以Zp表示评价节域物元，vp表示节域物元关于评价指标X的全部量值范围，vp=(upi，vpi)(i=1, 2, …，n)。耕地生态安全Z节域的物元矩阵可表示为：
Np=(Zp,X,vp)=|ZpX1(up1,vp1)X2(up2,vp2)⋯⋯Xn(upn,vpn)| (3) 3.1.3 确定关联度函数及关联度
K(Ci)={−ρ(C,C0)|C0|C∈C0ρ(C,C0)ρ(C,Cp)−ρ(C,C0)C∉C0 (4) ρ(C,C0)=|C−12(u0+v0)|−12(v0−u0) (5) ρ(C,Cp)=|C−12(up+vp)|−12(v0−u0) (6) 其中，ρ(C, C0)表示点C与有限区间C0=[u0, v0]的距离；ρ(C, Cp)表示表示点C与有限区间Cp=[up, vp]的距离；|C0|=|v0-u0|；C、C0、Cp分别表示待评价物元的数值、经典域物元的取值范围和节域物元的取值范围。当Kj(Ci)=max[Kj(Ci)](j=1, 2，…，n), 得出待评价物元Zx的第i个指标属于等级j。
3.1.4 计算综合关联度并确定评价等级
Kj(Zc)=n∑i=1wiKj(Ci) (7) 上述公式中，Kj(Zc)表示待评价物元Zc关于等级j的综合关联度；Kj(Cj)表示Zc关于等级j的单指标关联度；wi为评价指标的权重。当Kj(Ci)=max[Kj(Zi)](j=1, 2，…，n), 得出待评价物元Zx属于等级j。
关联度数值表明评价指标或评价对象在相应等级内的稳定程度，数值越大越稳定，反之，数值越小，表明该评价指标或待评价对象有往相邻次级转化的趋势，若相邻两个等级的关联度值越接近则转化可能性越大。当K(C)∈(-∞，-1)时，表示被评价对象不符合该级评价标准，且没有转化为该等级标准的条件；当K(C)∈[-1, 0)时，表示被评价对象不完全符合该级别评价标准，但具备转化条件，数值越大，越容易转化；当K(C)∈[0, 1)时，表示被评价对象符合该级别评价标准，数值越大，越接近该标准要求；当K(C)∈[1, +∞)时，评价对象超过该级别评价标准，数值越大，符合程度越高。
3.2 确定耕地生态安全评价的经典域和节域
基于耕地生态安全的可拓性，将评价等级划分为五个级别: Ⅰ表示理想安全、Ⅱ表示较安全、Ⅲ表示临界安全、Ⅳ表示较不安全、Ⅴ表示极不安全，从而确定经典域。经典域的确定主要依据《生态县、生态市、生态省建设指标(试行)》[18]关于生态县建设的标准值、全国平均水平值、福建省平均水平值以及三明市平均水平值。若有的评价指标已获得国际公认标准，则采用国际通行标准。所确定的经典域和节域见表 2。
表 2 基于占补平衡的耕地生态安全评价经典域和节域取值范围Table 2. Intervals of data collection on classical and segment fields of area cultivated land for ecological security evaluation based on occupation and compensation balance评价指标 经典域 节域 Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ X1 [0.107，0.137] [0.077, 0.107) [0.047，0.077) [0.023, 0.047) [0, 0.023) [0, 0.137] X2 [0，0.17) [0.17, 0.34) [0.34，1.15) [1.15, 1.96) [1.96，2.77] [0, 2.77] X3 [2.68, 3.22] [2.14, 2.68) [1.60, 2.14) [1.06, 1.60) [0, 1.06) [0, 3.22] X4 [0, 2) [2, 8) [8, 10) [10, 24) [24, 35] [0, 35] X5 [0, 100) [100, 300) [300, 1000) [1000, 2000) [2000, 3000] [0, 3000] X6 [0, 5) [5, 10) [10, 15) [15, 20) [20, 25] [0, 25] X7 [16, 20]> [12, 16) [8, 12) [4, 8) [0, 4) [0, 20] X8 [16, 20] [12, 16) [8, 12) [4, 8) [0, 4) [0, 20] X9 [6500,7500] [5920, 6500) [4970，5920) [2490, 4970) [0，2490) [0, 7500] X10 [45,50] [35, 45) [15, 35) [5, 15) [0, 5) [0, 50] X11 [80,85] [75, 80) [60, 75) [20, 60) [0, 20) [0, 85] X12 [50, 255) [255, 465) [465, 685) [685, 800) [800, 1000] [50,1000] X13 [0，10) [10，20) [20，25) [25，40) [40，80] [0, 80] X14 [97,100] [85, 97) [80, 85) [70, 80) [0, 70) [0, 100] X15 [97,100] [85, 97) [80, 85) [70, 80) [0, 70) [0, 100] X 16 [97,100] [85, 97) [80, 85) [70, 80) [0, 70) [0, 100] 4. 结果与分析
4.1 物元模型总体评价分析
将待评价物元依次输入公式(1)～公式(7)中，以2010年的人均耕地面积(x1)为例，得到该指标在对应5个等级下的关联度为K1(x1)=-0.2857，K2(x2)=0.3333，K3(x3)=-0.1667，K4(x4)=-0.4444，K5(x5)=-0.5614，可以判定2010年该指标属于级别Ⅱ(较安全)。同理可得其他年份相应指标的关联度如表 3。
表 3 清流县2010-2016年耕地生态安全指标关联度等级Table 3. Degree of association of indicators for ecological security evaluation on cultivated land in Qingliu county from 2010 to 2016关联度 2010年 2011等级 2012等级 2013等级 2014等级 2015等级 2016等级 Ⅰ级 Ⅱ级 Ⅲ级 Ⅳ级 Ⅴ级 等级 Kj(x1) -0.2857 0.3333 -0.1667 -0.4444 -0.5614 Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Kj(x2) -0.4220 -0.3152 0.3580 -0.4522 -0.6786 Ⅲ Ⅴ Ⅴ Ⅱ Ⅱ Ⅲ Ⅱ Kj(x3) -0.4963 -0.3692 -0.1563 0.4630 -0.1768 Ⅳ Ⅲ Ⅰ Ⅴ Ⅴ Ⅴ Ⅴ Kj(x4) -0.4323 -0.0455 0.2000 -0.1600 -0.6500 Ⅲ Ⅲ Ⅲ Ⅲ Ⅲ Ⅲ Ⅲ Kj(x5) 0.1400 -0.1400 -0.7133 -0.9140 -0.9570 Ⅰ Ⅰ Ⅰ Ⅰ Ⅰ Ⅰ Ⅰ Kj(x6) 0.0900 -0.0900 -0.5450 -0.6967 -0.7725 Ⅰ Ⅰ Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Kj(x7) -0.3540 -0.0876 0.2125 -0.2625 -0.4469 Ⅲ Ⅰ Ⅱ Ⅲ Ⅲ Ⅲ Ⅳ Kj(x8) -0.2368 0.4500 -0.2750 -0.5167 -0.6375 Ⅱ Ⅱ Ⅱ Ⅱ Ⅲ Ⅲ Ⅲ Kj(x9) -0.3679 -0.2536 0.1443 -0.0542 -0.5224 Ⅱ Ⅲ Ⅲ Ⅲ Ⅲ Ⅲ Ⅲ Kj(x10) -0.4081 -0.0640 0.0550 -0.5400 -0.6422 Ⅲ Ⅲ Ⅲ Ⅲ Ⅲ Ⅲ Ⅳ Kj(x11) -0.5000 0.2000 -0.1000 -0.6400 -0.8615 Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Kj(x12) -0.1940 0.3095 -0.3494 -0.5748 -0.6400 Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Ⅱ Kj(x13) -0.3735 -0.1615 0.0460 -0.0092 -0.3808 Ⅲ Ⅲ Ⅲ Ⅳ Ⅳ Ⅳ Ⅳ Kj(x14) -0.4658 -0.2650 -0.1285 0.3460 -0.2180 Ⅳ Ⅳ Ⅰ Ⅰ Ⅰ Ⅱ Ⅱ Kj(x15) -0.4667 -0.2733 -0.1439 0.4040 -0.1987 Ⅳ Ⅱ Ⅳ Ⅳ Ⅲ Ⅴ Ⅳ Kj(x16) 0.0000 -1.0000 -1.0000 -1.0000 -1.0000 Ⅰ Ⅱ Ⅰ Ⅰ Ⅰ Ⅰ Ⅰ 将表 3得到各指标的关联度值结合相应权重各指标对应的关联度值以及各指标的权重代入公式(7)中计算，得出清流县2010年耕地生态安全所有指标对应5个评价等级的综合关联度为：K1(Z2010)=-0.3353，K2(Z2010)=-0.0535，K3(Z2010)=-0.0886，K4(Z2010)=-0.2130，K5(Z2010)=-0.5071，可以判定清流县2010年耕地生态安全级别处在Ⅱ级别，但由于-1 < K2(Z2010) < 0，表示清流县2010年耕地生态安全不完全符合Ⅱ级别(较安全)的评价标准，但具备转化为该级别的条件，耕地生态安全正在向较安全状态转化。运用相同方法计算得出清流县2010-2016年耕地生态安全评价结果(表 4)。
表 4 清流县2010-2016年耕地生态安全评价结果Table 4. Ecological security evaluation on cultivated land in Qingliu county, 2010-2016关联度 Ⅰ级 Ⅱ级 Ⅲ级 Ⅳ级 Ⅴ级 安全级别 Kj(Z2010) -0.3353 -0.0535 -0.0886 -0.2130 -0.5071 向较安全转化 Kj(Z2011) -0.2590 -0.0731 -0.1128 -0.4020 -0.5632 向较安全转化 Kj(Z2012) -0.1542 -0.0702 -0.3100 -0.5067 -0.6455 向较安全转化 Kj(Z2013) -0.4097 -0.1086 -0.2124 -0.4417 -0.3931 向较安全转化 Kj(Z2014) -0.4324 -0.1393 -0.1745 -0.3971 -0.3666 向较安全转化 Kj(Z2015) -0.4714 -0.1667 -0.1815 -0.3286 -0.3200 向较安全转化 Kj(Z2016) -0.4993 -0.2029 -0.2595 -0.2876 -0.3293 向较安全转化 由表 4分析可得，从2010-2016年，清流县耕地生态安全均处于向较安全转化级别，但在Ⅱ级(较安全)状态内有小幅动态变化，变化幅度分别为：-0.019 6，0.002 9，-0.038 4，-0.030 7，-0.027 4，-0.036 2。清流县耕地生态安全水平在2011年到2012年间有所改善，综合关联度提升0.002 9。但之后几年均呈现下降趋势，总的来说，根据Kj(Z2010)>Kj(Z2012)>Kj(Z2011)>Kj(Z2013)>Kj(Z2014)>Kj(Z2015)>Kj(Z2016)，近几年清流县耕地生态安全水平关联度呈下降趋势，说明清流县耕地生态安全水平在Ⅱ级(较安全)状态内并没有朝改善的方向变化，反而呈现恶化趋势。
综合6年来的评价指标信息，工业固体废物综合利用率出现不同等级上升趋势，说明它是清流县耕地生态安全水平的贡献因素。而耕地增加面积与耕地面积比例、人口自然增长率、GDP年增长率、单位面积耕地粮食产量、固定资产投资总额年增长率、单位耕地面积农药施用量出现不同等级的下降趋势，说明以上评价指标是清流县耕地生态安全水平提升的限制因素。
4.2 物元模型阶段评价分析
4.2.1 指标分析信息角度
2011-2012年，综合关联度提升0.002 9，清流县耕地生态安全水平得到改善。此时期的耕地增加面积与耕地面积比例、工业固体废物综合利用率、生活垃圾处理率3项评价指标在关联度有安全等级的提升，为耕地生态安全的改善做出贡献。人口自然增长率、农民人均纯收入年增长率、污水处理率3项指标安全等级有所下降，成为该年份耕地生态安全水平的限制因素。
2012年后，清流县耕地生态安全水平连续两年呈恶化趋势。这一时期，耕地增加面积与耕地面积比例、农民人均纯收入年增长率、GDP年增长率、固定资产投资总额年增长率、工业固体废物综合利用率、单位耕地面积农药施用量6项指标的安全等级有所下降，成为此时期耕地生态安全水平的限制因素，其他指标趋于稳定状态。
4.2.2 指标权重角度
从评价指标权重来看，耕地增加面积与耕地面积比例(X3=0.2017)、农民人均纯收入年增长率(X7=0.0941)、人均耕地面积(X1=0.0837)、固定资产投资总额年增长率(X10=0.0806)、人口自然增长率(X6=0.0778)、单位耕地面积农药施用量(X13=0.0766)、单位耕地面积化肥施用量(X12=0.0725)、大于25度耕地面积比重(X4=0.0705)是清流县耕地生态安全的主要影响因子。
2012-2016年清流县耕地生态安全水平有所下降，主要原因还在于清流县城镇化进程加快的同时，耕地减少面积与耕地面积比例从Ⅰ级(理想安全)降到Ⅴ级(极不安全)，农民人均纯收入年增长率从Ⅱ级(较安全级别)下降到Ⅳ级(较不安全)。此外，固定资产投资总额年增长率、单位耕地面积农药施用量和工业固体废物综合利用率等因素安全级别均下降，使得清流县耕地生态安全水平在向较安全转化级别内呈现恶化趋势。
4.3 多指标综合评价分析
由于耕地的生态安全状况由多种因素共同决定，而多指标综合评价法旨在根据不同评价目的，选择多个指标，通过评价方法将多个评价指标转化为能反映待评价对象总体特征的信息[21]，是目前较为通用的评价方法。运用多指标综合评价方法对清流县耕地生态安全进行评价，用以验证物元模型的评价结果。计算公式如下：
Z=17∑i=1WiXi (8) Z为评价综合值；Wi为指标权重；Xi为标准化后的指标数值。
根据上述公式，得到清流县各年份耕地生态安全综合总分值(图 1)。
![]()
图 1 清流县耕地生态安全差异Figure 1. Differences on ecological security of cultivated lands in Qingliu county清流县耕地生态安全水平由高到低排序依次为：Z2012>Z2011>Z2010>Z2013>Z2014>Z2015>Z2016。两种方法的评价结果基本一致，从2012年后生态安全水平均呈现下降趋势，和物元分析方法的评价结果：Kj(Z2010)>Kj(Z2012)>Kj(Z2011)>Kj(Z2013)>Kj(Z2014)>Kj(Z2015)>Kj(Z2016)相比，均存在Z2012>Z2013>Z2014>Z2015>Z2016的变化趋势。但也存在一定差异：在物元模型分析中，2010-2011年的耕地生态安全水平呈现恶化趋势，在2012年存在微弱改善；而在多指标综合分析方法中，2010-2012年耕地生态安全水平连续改善，且在2012年的耕地生态安全水平高于其他年份。可能因为两种评价方法关注点不一致，物元分析方法更加关注单个指标的分析，而多指标综合评价分析方法注重对综合质量的评价。
5. 结论及建议
5.1 结论
以福建省典型丘陵地区清流县作为研究区域，在分析全县社会经济和耕地现状的基础上，构建清流县耕地生态安全评价指标体系，将耕地生态安全、特征值和量值作为基本元，结合《生态县、生态市、生态省建设指标(试行)》[18]关于生态县建设的标准值、全国平均水平值、福建省平均水平值以及三明市平均水平值构建评价经典域和节域，计算各指标间的关联度以及评价等级关联度，从而构建清流县耕地生态安全的物元模型。从熵权法确权到建立物元模型，此方法克服了评价过程中人为因素的干扰，量化处理评价指标，提高了评价结果的精度。并将评价结果与传统的多指标综合评价法进行比较，得出以下结论：
清流县耕地生态安全级别在评价年份内均为“向较安全转化”级别，但在Ⅱ级(较安全)状态内有小幅动态变化，总体而言，清流县耕地生态安全水平在Ⅱ级(较安全)状态内并没有朝改善的方向变化，反而呈现恶化趋势。其中工业固体废物综合利用率出现不同等级上升趋势，说明它是清流县耕地生态安全水平的贡献因素。而耕地增加面积与耕地面积比例、人口自然增长率、GDP年增长率、单位面积耕地粮食产量、固定资产投资总额年增长率、单位耕地面积农药施用量出现不同等级的下降趋势，说明以上评价指标是清流县耕地生态安全水平提升的限制因素。
2012年之后，随着清流城镇化进程加快，耕地减少面积与耕地面积比例、农民人均纯收入年增长率、固定资产投资总额年增长率、单位耕地面积农药施用量和工业固体废物综合利用率5项指标的安全等级有所下降，导致清流县耕地生态安全水平在向较安全转化级别内呈现恶化趋势。从指标权重来看，最重要的两个制约因素是耕地减少面积与耕地面积比例和农民人均纯收入年增长率。
将熵权物元评价与多指标综合评价相比，其评价结果基本一致，虽然个别年份出现差异，可能是因为两种评价方法关注点不一致，物元分析方法更加关注单个指标的分析，而多指标综合评价分析方法注重对综合质量的评价。
将物元模型应用到丘陵地区进行耕地生态安全评价，基本达到了预期目的。在评价指标体系方面，尝试在文献检索的基础上，综合选取指标，但由于生态系统受到多方面因素影响，有些指标例如水土流失、土壤重金属污染等数据获取度有限，尚未列入评价范围；另一方面，在经典域和节域确定上主要依据地方平均水平等，没有办法量化每一个指标的影响度，因此，经典域和节域的取值范围有待进一步探讨。
5.2 建议
研究结论显示，清流县耕地生态安全在“向较安全转化”级别内呈现恶化趋势，并未完全符合较安全级别的标准。耕地增加面积与耕地面积比例、人口自然增长率、GDP年增长率、单位面积耕地粮食产量、固定资产投资总额年增长率、单位耕地面积农药施用量安全等级在评价年份内均有下降。对耕地资源生态安全有所制约，有可能限制清流县未来的可持续发展。因此，应立足于清流县发展实际，在改善清流人民经济生活水平的同时，合理利用耕地资源，将经济发展与耕地生态安全结合起来，最终实现人与自然的和谐统一。提出以下几方面建议：
5.2.1 考核制度与耕地生态安全监测并行
在严格控制非农建设用地占用耕地的前提下，确保耕地资源数量安全的同时，加大高标准农田建设强度，通过合理的土地整治手段维持耕地资源的质量安全。将耕地资源安全纳入到政府部门绩效考核体系中，将耕地生态安全作为发展的战略高度，增强政府干部维护耕地生态安全的责任感和使命感。此外，还应完善环境保护的相关政策法规，严令禁止工业“三废”不合理排放行为，提高县域污水处理率，提倡绿色生态可持续发展。
5.2.2 创新补偿机制，协同多方力量共同守护耕地生态安全
由耕地生态安全的受益者向耕地生态安全提供者实施补偿行为[22]。耕地资源作为一种公共资源，多方利益既得者享受着耕地生态安全带来的好处，但在以往的补偿制度中，多以政府部门为主体，实际上，工厂企业的生产发展在一定程度上以耕地作为生产资料，公民也从耕地生态安全中收益。应将补偿机制的主体从原来的政府为主转化为政府-企业-公民参与的形式，避免产生外部不经济行为。通过惠农政策，农村合作社等方式鼓励农民耕种，种植生态价值较高的作物，避免耕地抛荒，在保障耕地资源数量的同时维持耕地生态安全。
5.2.3 加强对耕地生态保护的宣传
耕地生态安全与每一个人关系密切，农户兼有耕地生态安全的受益者和提供者双重角色，政府应当加强对耕地生态保护的宣传，通过线下与线上结合的创新宣传方式，找出有效的宣传突破点，充分利用乡贤在农村治理中的地位，广泛发动农户参与其中。此外，鼓励农户使用有机肥代替化肥，合理控制农药用量，倡导绿色生态种植，在保障耕地质量的同时提高农户保护耕地生态安全的意识。
-
![]()
图 1 混合酚酸与土样酚酸HPLC出峰
注：1为没食子酸；2为香豆酸；3为丁香酸；4为香兰素；5为阿魏酸；6为苯甲酸。
Figure 1. HPLC on mixed and soil phenolic acids samples
![]()
图 2 不同处理下的三七出苗率、保苗率及产量
Figure 2. Seedling emergence and survival rate and root yield of P. notoginseng under treatments
![]()
图 3 不同处理下连作三七根际土壤酚酸含量
Figure 3. Phenolic acids in rhizosphere soil of P. notoginseng under monoculture regime
![]()
图 4 不同处理三七根际土壤DGGE图谱
注：A为TRI；B为CK；C为NP-DP；D为NP；E为SP；F为EM。
Figure 4. DGGE profile of rhizosphere soil of P. notoginseng under treatment
![]()
图 5 DGGE数据主成分分析(PCA)及聚类分析
Figure 5. Principal component and clustering analyses on DGGE data
![]()
图 6 不同处理下连作三七根际土壤中部分特异菌群含量的原位qRT-PCR分析
注：①A为Anteholosticha gracilis，B为Scutellospora savannicola，C为Anteholosticha manca，D为Trichoderma harzianum，E为Thermococcus barophilus，F为Monacrosporium janus，G为Paraphaeosphaeria verruculosa，H为Geminibasidium donsium，I为Pleosperales，J为Sphaerobolus iowensis。②图中数据为每克鲜土中的特异菌含量。
Figure 6. In situ qRT-PCR analysis on population of specific microbial flora in rhizosphere soil of monocultured P. notoginseng under different treatments
表 1 不同处理对连作三七根际土壤中营养元素含量的影响
Table 1 Effect of treatments on soil nutrient content in monocultured P. notoginseng rhizosphere soil
处理 全氮/(g·kg-1) 全磷/(g·kg-1) 全钾/(g·kg-1) 碱解氮/(mg·kg-1) 速效磷/(mg·kg-1) 速效钾/(mg·kg-1) pH CK 1.29±0.80a 0.22±0.01c 0.75±0.06e 112.58±4.40a 29.54±0.32b 140.83±2.55b 6.49±0.03b NP 0.66±0.83c 0.22±0.01c 3.66±0.06c 42.58±11.65b 24.84±0.12d 155.87±2.55a 6.72±0.11a NP-DP 1.01±0.89b 0.25±0.02b 0.95±0.09d 109.67±9.95a 34.65±0.11a 125.25±1.67c 6.43±0.05b SP 0.22±0.84d 0.12±0.01d 3.57±0.07c 18.08±1.01d 22.70±0.21f 41.75±1.67f 6.03±0.05c EM 0.29±0.90d 0.14±0.01d 4.07±0.08a 20.42±2.02cd 23.72±0.18e 57.34±0.96e 5.95±0.02c TRI 0.58±0.35c 0.37±0.01a 3.83±0.10b 32.67±9.63bc 26.80±0.64c 109.11±0.96d 6.72±0.15a 注：①CK为空白对照；NP为头茬；NP-DP为头茬发病；SP为重茬未处理；EM为有益微生物处理；TRI为木霉菌剂处理。②同列数据后不同字母表示处理间差异显著(P < 0.05)。表 2同。 
下载: 导出CSV
表 2 不同根际土壤微生物多样性指数
Table 2 Microbial diversity index of rhizosphere soil under treatments
处理 SIMPSON(J) SHANNON(H) 均匀度 BRILLOUIN McIntosh(Dmc) CK 1.0540c 4.313c 0.992a 3.014a 1.131c NP 1.056b 4.407b 0.988b 2.615bc 1.138b NP-DP 1.060a 4.482a 0.990ab 2.740ab 1.149a SP 1.049d 4.221d 0.988b 2.896bc 1.118d EM 1.056bc 4.362b 0.993a 2.699bc 1.136bc TRI 1.042e 4.056e 0.992a 3.010a 1.100e
下载: 导出CSV
表 3 部分特异条带BLAST比对结果
Table 3 BLAST on specific bands
编号 真菌拉丁名 中文名 A1 Anteholosticha gracilis A2 Trichoderma harzianum 哈茨木霉菌 A3 Scutellospora savannicola 盾巨孢囊霉丛枝菌根 A4 Anteholosticha manca A6 Scutellospora savannicola 盾巨孢囊霉丛枝菌根 A9 Pleosperales 格孢腔菌目 B1 Bryum arcticum 极地真藓 B2 Pleosporales 格孢腔菌目 B3 Ignatius tetrasporus 伊格内修斯子囊菌 C1 Anteholosticha gracilis C2 Pleosperales 格孢腔菌目 C4 Bryum arcticum 极地真藓 D1 Bryum arcticum 极地真藓 D2 Scutellospora savannicola 盾巨孢囊霉丛枝菌根 D3 Entrophospora 稀有内养囊霉 E1 Pleosperales 格孢腔菌目 E2 Scutellospora savannicola 盾巨孢囊霉丛枝菌根 E4 Sphaerobolus iowensis 弹球菌属 E5 Ignatius tetrasporus 伊格内修斯子囊菌 F1 Thermococcus barophilus 嗜热球菌属 F2 Monacrosporium janus 两栖单顶孢菌 F3 Ignatius tetrasporus 伊格内修斯子囊菌 F4 Paraphaeosphaeria verruculosa F6 Geminibasidium donsium F7 Scutellospora savannicola 盾巨孢囊霉丛枝菌根 注：每个菌种比对结果同源性均在98%以上。
下载: 导出CSV
-
[1] 杨志刚, 陈阿琴, 俞颂东.三七药理研究新进展[J].上海中医药杂志, 2005, 39(4):59-62. DOI: 10.3969/j.issn.1007-1334.2005.04.025 [2] 王莹, 禇扬, 李伟, 等.三七中皂苷成分及其药理作用的研究进展[J].中草药, 2015, 46(9):1381-1392. http://d.old.wanfangdata.com.cn/Periodical/zcy201509023 [3] 孟祥霄, 黄林芳, 董林林, 等.三七全球产地生态适宜性及品质生态学研究[J].药学学报, 2016, 51(9):1483-1493. http://mall.cnki.net/magazine/Article/YXXB201609021.htm [4] 张子龙, 李凯明, 杨建忠, 等.轮作对三七连作障碍的消减效应研究[J].西南大学学报(自然科学版), 2015, 37(8):39-46. http://d.old.wanfangdata.com.cn/Periodical/xnnydxxb201508007 [5] 张子龙, 王文全, 缪作清, 等.主成分分析在三七连作土壤质量综合评价中的应用[J].生态学杂志, 2013, 32(6):1636-1644. http://d.old.wanfangdata.com.cn/Conference/7964304 [6] 孙雪婷, 李磊, 龙光强, 等.三七连作障碍研究进展[J].生态学杂志, 2015, 34(3):885-893. http://d.old.wanfangdata.com.cn/Periodical/sdnyswxb201101016 [7] 刘莉, 刘大会, 金航, 等.三七连作障碍的研究进展[J].山地农业生物学报, 2011, 30(1):70-75. DOI: 10.3969/j.issn.1008-0457.2011.01.016 [8] 张学炜, 黄学森.西瓜连作障碍及其防治方法[J].中国瓜菜, 1993(2):21-23. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK000004993433 [9] 高微微, 张西梅, 田给林, 等.探索中前行——中药材病害研究概况与思考[J].植物保护, 2016, 42(5):15-23. DOI: 10.3969/j.issn.0529-1542.2016.05.003 [10] 吴林坤, 吴红淼, 朱铨, 等.不同改良措施对太子参根际土壤酚酸含量及特异菌群的影响[J].应用生态学报, 2016, 27(1):3623-3630. http://d.old.wanfangdata.com.cn/Periodical/yystxb201611027 [11] 张丽娟, 曲继松, 郭文忠, 等.微生物菌肥对黄河上游地区设施土壤微生物及酶活性的影响[J].中国土壤与肥料, 2014(5):32-36. http://d.old.wanfangdata.com.cn/Periodical/trfl201405007 [12] 王涛, 乔卫花, 李玉奇, 等.轮作和微生物菌肥对黄瓜连作土壤理化性状及生物活性的影响[J].土壤通报, 2011, 42(3):578-583. http://d.old.wanfangdata.com.cn/Periodical/trtb201103013 [13] 逄焕成, 李玉义, 严慧峻, 等.微生物菌剂对盐碱土理化和生物性状影响的研究[J].农业环境科学学报, 2009, 28(5):951-955. DOI: 10.3321/j.issn:1672-2043.2009.05.016 [14] 鲍士旦.土壤农化分析[M].第3版.北京:中国农业出版社, 1999:25-109. [15] 李贺勤, 刘奇志, 张林林, 等.草莓连作土壤酚酸类物质积累对土壤线虫的影响[J].生态学杂志, 2014, 33(1):169-175. http://d.old.wanfangdata.com.cn/Periodical/stxzz201401025 [16] 吴宗伟, 王明道, 刘新育, 等.重茬地黄土壤酚酸的动态积累及其对地黄生长的影响[J].生态学杂志, 2009, 28(4):660-664. http://d.old.wanfangdata.com.cn/Periodical/stxzz200904014 [17] 胡元森, 李翠香, 杜国营, 等.黄瓜根分泌物中化感物质的鉴定及其化感效应[J].生态环境, 2007, 16(3):954-957. DOI: 10.3969/j.issn.1674-5906.2007.03.050 [18] 夏围围, 贾仲君.高通量测序和DGGE分析土壤微生物群落的技术评价[J].微生物学报, 2014, 54(12):1489-1499. http://d.old.wanfangdata.com.cn/Periodical/wswxb201412012 [19] 龙良鲲, 羊宋贞, 姚青, 等. AM真菌DNA的提取与PCR-DGGE分析[J].菌物学报, 2005(4):94-99. http://d.old.wanfangdata.com.cn/Periodical/jwxt200504015 [20] 位小丫, 林煜, 陈婷, 等.田间条件下植物促生细菌缓解太子参连作障碍的效果评价[J].生态学杂志, 2018, 37(2):399-408. http://d.old.wanfangdata.com.cn/Periodical/stxzz201802014 [21] 朱扬玲.采用PCR-DGGE方法研究浙江玫瑰醋酿造过程中的微生物多样性[D].杭州: 浙江工商大学, 2009. [22] 刘鹏飞, 赵丹, 宋刚, 等.变性梯度凝胶电泳技术在微生物多样性研究中的应用[J].微生物学杂志, 2013, 33(6):88-92. DOI: 10.3969/j.issn.1005-7021.2013.06.017 [23] 刘绍雄, 王明月, 王娟, 等.基于PCR-DGGE技术的剑湖湿地湖滨带土壤微生物群落结构多样性分析[J].农业环境科学学报, 2013(7):1405-1412. http://www.cnki.com.cn/Article/CJFDTOTAL-NHBH201307020.htm [24] 梁小刚, 蔡国林, 张中华, 等.利用nested PCR-DGGE技术分析江苏啤酒大麦真菌群落结构[J].食品与发酵工业, 2012, 38(8):1-6. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK201205542302 [25] 牛芳胜, 马志强, 毕秋艳, 等.哈茨木霉菌与5种杀菌剂对番茄灰霉病菌的协同作用[J].农药学学报, 2013, 15(2):165-170. DOI: 10.3969/j.issn.1008-7303.2013.02.07 [26] 王幼珊, 刘润进.球囊菌门丛枝菌根真菌最新分类系统菌种名录[J].菌物学报, 2017, 36(7):820-850. http://d.old.wanfangdata.com.cn/Periodical/jwxt201707004 [27] LI S, MIAO Z, ZHANG Y, et al. Monacrosporium janus sp. nov. a new nematode-trapping hyphomycete parasitizing sclerotia and hyphae of Sclerotinia sclerotiorum[J]. Mycological Research, 2003, 107(7):888-894. DOI: 10.1017/S0953756203008165
[28] NGUYEN H D, NICKERSON N L, SEIFERT K A. Basidioascus and Geminibasidium:a new lineage of heat-resistant and xerotolerant basidiomycetes[J]. Mycologia, 2013, 105(5):1231-1250. DOI: 10.3852/12-351
[29] SMILEY R W, PATTERSON L M. Pathogenic fungi associated with Fusarium foot rot of winter wheat in the semiarid Pacific Northwest[J]. Plant Disease, 1996, 80(8):944-949. DOI: 10.1094/PD-80-0944
[30] 祝英, 彭轶楠, 巩晓芳, 等.不同微生物菌剂对当归苗生长及根际土微生物和养分的影响[J].应用与环境生物学报, 2017(3):511-519. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QKC20172017080200314373 [31] 董占梅, 杜立财, 程林, 等. 4种微生物菌剂对老参地人参产量和质量的影响[J].特产研究, 2016, 38(4):9-13. http://www.cnki.com.cn/Article/CJFDTOTAL-TCYA201604004.htm [32] 罗玉兰, 田龚, 张冬梅, 等.微生物菌剂对连栋大棚土壤养分及硝态氮累积的影响[J].中国农学通报, 2015, 31(13):224-228. DOI: 10.11924/j.issn.1000-6850.casb14120113 [33] 尹淑丽, 张丽萍, 张根伟, 等.复合微生态菌剂对黄瓜根际土壤微生物数量及酶活的影响[J].微生物学杂志, 2012, 32(1):23-27. DOI: 10.3969/j.issn.1005-7021.2012.01.005 [34] 韩忠明, 杨颂, 韩梅, 等.不同菌剂对人参连作土壤酶活性的影响[J].东北农业科学, 2016(1):50-53. http://www.cnki.com.cn/Article/CJFDTOTAL-JLNK201601013.htm [35] 王强, 阮晓, 李兆慧, 等.植物自毒作用及针叶林自毒研究进展[J].林业科学, 2007, 43(6):134-142. http://d.old.wanfangdata.com.cn/Periodical/lykx200706024 [36] 张亚琴, 陈雨, 雷飞益, 等.药用植物化感自毒作用研究进展[J].中草药, 2018, 49(8):1946-1956. http://d.old.wanfangdata.com.cn/Periodical/zcy201808032 [37] KAUR H, KAUR R, KAUR S, et al. Taking ecological function seriously:soil microbial communities can obviate allelopathic effects of released metabolites[J]. Plos One, 2009, 4(3):1-6.
[38] 吴凤芝, 王伟.大棚番茄土壤微生物区系研究[J].北方园艺, 1999(3):1-2. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199900214015 -
期刊类型引用(3)
1. 韩丽红，潘玉君，高敏，黄欢. 贵州耕地生态安全评价及障碍因子分析. 西南林业大学学报(社会科学). 2023(06): 68-73 . 
百度学术
2. 宋林韩，金鹏，吴桐，卫新东. 陕北黄土高原耕地生态安全的时空格局. 江苏农业科学. 2020(12): 249-257 . 百度学术
3. 赵疏航，何刚，朱艳娜，李洁，杜宇. 基于EES-TOPSIS模型的水土生态安全评价与预测. 哈尔滨商业大学学报(自然科学版). 2020(06): 751-759 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 1388
- HTML全文浏览量: 382
- PDF下载量: 16
- 被引次数: 4
