基于最小数据集的周至县猕猴桃园地土壤质量评价

2019-02-15张福平李肖娟刘洁遥

生态与农村环境学报 2019年1期

张福平,高张,李肖娟,聂硕,李玲,刘洁遥

(陕西师范大学地理科学与旅游学院,陕西西安 710062)

土壤质量是土壤在生态系统边界范围内,维持作物生产能力,保持环境质量及促进动植物健康的能力[1-2]。作为表征土壤条件动态变化最敏感的指标,土壤质量能反映土壤管理水平,对退化土地的恢复、区域土地资源管理和土地持续利用也具有重要意义[3]。对土壤质量进行系统评价,定量分析作物生长过程中的主要限制指标,可为调节作物生产和增产增收提供技术支持。然而,由于土壤质量评价目的和评价对象尺度不同,且评价工作复杂，易受外界因素干扰,国内外对土壤质量评价并没有建立统一标准[4-5]。同时,由于土壤利用方式的多变性、土壤性质的时空变异性和数据获取的高成本等因素,难以获取所有能表征土壤性质指标的数据[6-7]。目前大多数研究是利用一些土壤指标,建立最小数据集(minimum data set,MDS)进行分析[8-9],如GOVAERTS等[10]基于MDS评估墨西哥高原地区小麦和玉米农田土壤质量,RAHMANIPOUR等[11]基于MDS评估伊朗加兹温省土壤质量,VOLCHKO等[12]基于MDS评估生态恢复进程中绿地区域土壤状况。国内学者也基于MDS评价南方地区[13]、黄河三角洲地区[14]和塔里木河流域[15]土壤质量。李桂林等[16]对MDS进行优化,提出将要素在各主成分上的综合载荷作为土壤指标筛选依据。上述研究表明MDS可用于评价土壤质量,但目前对果园土壤鲜有涉及。

陕西省周至县是全球最大的猕猴桃生产基地,生产的猕猴桃畅销全国乃至海外,猕猴桃为当地经济发展和农民增收做出了重要贡献,但土壤质量一直是影响猕猴桃产业发展的重要因素之一。猕猴桃是一种多年生藤本植物,根部在吸收土壤养分的同时会不断分泌化学物质,与土壤之间形成一种耦合关系,使土壤质量处于动态变化中。而猕猴桃对土壤质量要求很高,当土壤条件不适时,往往会出现生理性病害,严重时甚至死苗。近年来,关于猕猴桃园地土壤质量的研究主要集中在土壤养分分布状况、养分空间变异性和污染等理化性质方面[17-19],而对猕猴桃园地土壤质量评价研究较少。因此,笔者以不同年限猕猴桃园地为研究对象,测定18个土壤理化指标,基于MDS对研究区土壤质量进行定量评价,揭示猕猴桃不同生长年限对土壤理化指标和土壤质量的影响,以期为猕猴桃种植研究提供科学支撑,同时为产业调整和种植方式改良提供理论依据。

1 材料与方法

1.1 研究区概况

周至县猕猴桃产业带位于34°06′～34°13′ N,108°00′～108°32′ E,地势西南高、东北低,地形主要为平原、浅山和黄土原区,海拔400～500 m(图1),80%以上区域为砾质潮砂土(潮土),黑垆土、褐土和红土分布区域较少。研究区属温带大陆性季风气候区,年平均气温12～13.6 ℃,作物一年两熟,年降水量600～800 mm,秋季多阴雨,日照少,降温快。周至县猕猴桃种植面积为266.67 km2,占全县耕地总面积的80%以上,为全国提供60%的猕猴桃鲜果和80%的猕猴桃果干,是我国认定的唯一一个猕猴桃标准化管理示范县,素有“猕猴桃之乡”的美称。

图1 周至县猕猴桃产业带及采样点分布

1.2 样品采集与处理

根据周至县行政区划图和土地利用现状图,并结合实地调查,共设33个采样点,于2017年5月分别采用土钻法和环刀法采集土样,环刀尺寸为直径50.46 mm×高50 mm,容积100 cm3。采样深度30 cm,每10 cm 1个土层,每个土层2个土壤样品,每个样品约100 g,共计198个土壤样品。样品采集后封装包好,移入实验室待测定。

根据猕猴桃生长习性[20-21]选取18个土壤理化指标进行测定:(1)土壤性质,选择物理性黏粒(SCC)、容重(BD)、毛管孔隙度(CP)、总孔隙度(TPO)、非毛管孔隙度(NCP)、通气度(CA)等指标。(2)水分,选择土壤含水量(SWC)。(3)pH,选择土壤pH值。(4)养分,选择有机质(SOM)、全氮(TN)、全磷(TP)、全钾(TK)、速效氮(AN)、速效磷(AP)和速效钾(AK)含量等养分指标。(5)盐分,选择土壤盐分(SSC)。(6)微量元素,选择Mg2+和SO42-含量。容重、含水量、孔隙度(毛管孔隙度、总孔隙度和非毛管孔隙度)和通气度采用环刀法测定后根据经验公式计算得到;物理性黏粒、土壤盐分、pH值、有机质含量和土壤养分(全氮、全磷、全钾、速效氮、速效磷和速效钾)含量分别采用Mastersizer 2000激光粒度分析仪、离子色谱法、PHSJ-4A型实验室pH计、重铬酸钾氧化-外加热法和HTYF200型微电脑多功能养分速测仪测定。

1.3 研究方法

1.3.1MDS的构建

对选取的指标进行主成分分析,计算各指标所在主成分特征值大于1的载荷,并将在同一主成分载荷大于0.5的土壤指标归为1组;若某个指标同时在2个主成分上的载荷都大于0.5,则将他归入与其他指标相关系数较低的那一组。同时,计算各指标在所在主成分上的综合载荷Norm值,值越大说明该指标解释的土壤质量信息就越大。选择各组Norm值在最高总分值10%范围内的指标[22],进一步进行相关分析,建立MDS。Norm值计算公式为

(1)

式(1)中,Nik为第i个指标在特征值大于1的前k个主成分Norm值;uik为第i个指标在第k个主成分上指标载荷;ek为第k个主成分特征值。

1.3.2土壤质量评价方法

通过土壤质量评价指数进行土壤质量评价,计算公式为

(2)

式(2)中,ISQ为土壤质量评价指数;Wi为第i个指标权重;Ni为第i个指标隶属度。

通过主成分分析法得到各土壤评价指标公因子方差,各指标权重为各指标公因子方差占公因子方差总和的百分数。隶属度由土壤指标所属的隶属度函数确定。隶属度函数可以分为升型和降型两种,计算公式为

(3)

(4)

式(3)～(4)中,f(x)为隶属度;x为土壤指标值;a和b分别为土壤指标阈值,即指标最小值和最大值。

2 结果与分析

2.1 土壤理化性质

18个土壤指标统计特征值见表1。

表1土壤指标统计特征值

Table1Statisticalcharacteristicvalueofsoilindicators

指标最小值最大值极差均值标准差变异系数/%偏度系数峰度系数正态分布 w(SCC)/%18.0362.5244.4939.4310.0025.36-0.26-0.47是 BD/(g·cm-3)0.981.660.681.380.128.91-0.060.40是 w(SWC)/%8.9435.3026.3618.113.3118.281.881.14是 CP/%33.0167.8134.8044.668.8919.900.890.08是 TPO/%35.0674.8839.8249.419.3718.970.870.09是 NCP/%0.1912.9412.744.802.5452.941.021.60是 CA/%4.1650.7346.5724.3510.9344.910.51-0.29是 pH值6.857.460.617.160.111.530.070.43是 w(SOM)/%0.334.454.122.490.8735.12-0.410.17是 w(SSC)/(g·kg-1)0.8726.6225.752.452.53103.070.180.61是 w(Mg2+)/(mg·kg-1)0.5092.0091.5055.000.3531.60-0.751.34是 w(SO42-)/(mg·kg-1)32.00817.50786.00374.002.3531.41-0.303.66是 w(TN)/(g·kg-1)0.012.932.920.680.5275.521.482.64是 w(TP)/(g·kg-1)0.593.342.751.320.5844.091.491.90是 w(TK)/(g·kg-1)24.3031.507.2027.701.806.400.46-0.52是 w(AN)/(mg·kg-1)11.80149.60137.8054.7726.5048.391.201.62是 w(AP)/(mg·kg-1)0.77138.20137.4344.4031.3070.490.880.50是 w(AK)/(mg·kg-1)4.07476.80472.73160.59105.7865.870.870.39是

SCC为物理性黏粒,BD为容重,SWC为含水量,CP为毛管孔隙度,TPO为总孔隙度,NCP为非毛管孔隙度,CA为通气度,SOM为有机质,SSC为盐分,TN为全氮,TP为全磷,TK为全钾,AN为速效氮,AP为速效磷,AK为速效钾。

K-S检验结果表明18项理化指标均呈正态分布,说明实验设计科学合理。土壤物理性黏粒、毛管孔隙度、非毛管孔隙度和土壤通气度分别为39.43%、44.66%、4.80%和24.35%;容重为1.38 g·cm-3,且大于1.4 g·cm-3的样品数量占总数40%以上；总孔隙度为49.41%,小于50%,说明土壤紧实度高;含水量为18.11%,表明土壤湿润度适宜;pH值为7.16,土壤呈中性;w(有机质)为2.49%,有机质含量偏高;w(盐分)为2.45 g·kg-1,土壤呈现轻度盐渍化;w(Mg2+)为55.00 mg·kg-1,w(Mg)适宜;w(SO42-)为374.00 mg·kg-1,S含量适宜;w(全氮)为0.68 g·kg-1，土壤全氮含量偏低;w(全磷)为1.32 g·kg-1,土壤全磷含量偏高;w(全钾)为27.7 g·kg-1,土壤全钾含量偏高;w(速效氮)为54.77 mg·kg-1,土壤速效氮含量偏低;w(速效磷)为44.40 mg·kg-1,土壤速效磷含量偏低;w(速效钾)为160.59 mg·kg-1,土壤速效钾含量偏高。容重、pH和全钾3个指标的变异系数小于10%,属于弱变异性;物理性黏粒、含水量、毛管孔隙度、总孔隙度、非毛管孔隙度、土壤通气度、有机质、Mg2+、SO42-、全氮、全磷、速效氮、速效磷和速效钾变异系数介于10%～100%之间,属于中等变异性;盐分变异系数最大,大于100%,为强变异性。变异性强弱反映土壤指标值在空间分布上的离散程度与代表性。

2.2 不同种植年限和土层深度土壤特性的变化

不同种植年限园地土壤指标平均值见表2。随猕猴桃园地种植年限的增加,容重表现出先减小再增加的趋势,最小值出现在种植时间为>10～15 a时。毛管孔隙度、总孔隙度、非毛管孔隙度、通气度和有机质呈逐渐减小的趋势;而含水量、Mg2+、SO42-、全氮、全磷、全钾、速效氮、速效磷和速效钾呈逐渐增加的趋势。单因子方差分析结果表明,除物理性黏粒与园地种植年限相关性不显著外,其他指标均与园地种植年限不同程度呈显著相关,其中pH、盐分、Mg2+、SO42-、全磷、速效磷和速效钾呈极显著相关。

表2不同种植年限和土层深度土壤指标平均值

Table2Theaveragesoilindexesindifferentyearsandwithdifferentdepthsofsoillayer

指标种植年限/a0^5 >5^10>10^15>15^20 >20P显著性 w(SCC)/%36.70±10.0740.12±9.7450.47±6.2544.99±8.2938.25±7.420.469不显著 BD/(g·cm-3)1.40±0.111.36±0.141.34±0.091.41±0.101.43±0.120.099一般显著 w(SWC)/%17.14±2.2817.87±4.3118.49±1.6418.57±2.7819.60±12.880.032显著 CP/%44.87±8.7147.04±10.2743.84±6.3241.23±4.6943.90±9.560.087一般显著 TPO/%50.56±9.2352.19±10.7847.86±5.6545.07±5.1048.10±9.870.090一般显著 NCP/%5.69±2.785.15±2.644.39±2.373.84±1.714.21±2.220.047显著 CA/%26.55±10.0527.96±12.1423.08±7.2019.62±6.7720.00±11.690.042显著 pH值7.16±0.117.15±0.127.19±0.097.12±0.127.19±0.070.003极显著 w(SOM)/%2.68±0.722.50±0.862.42±0.892.37±0.982.41±0.890.096一般显著 w(SSC)/(g·kg-1)1.89±0.583.03±4.202.19±0.362.53±0.362.05±0.640.008极显著 w(Mg2+)/(mg·kg-1)52.00±20.0054.00±17.0056.50±20.0056.50±16.0058.00±12.000.005极显著 w(SO42-)/(mg·kg-1)333.00±135.50374.50±83.50372.00±58.50449.00±116.50337.00±144.500.007极显著 w(TN)/(g·kg-1)0.64±0.480.63±0.460.99±0.760.55±0.410.79±0.410.044显著 w(TP)/(g·kg-1)1.16±0.361.22±0.471.18±0.351.50±0.631.64±0.870.003极显著 w(TK)/(g·kg-1)26.84±1.1828.07±1.8427.21±1.4227.58±1.8228.47±1.850.069一般显著 w(AN)/(mg·kg-1)41.81±16.4155.21±25.0147.57±19.2053.57±27.3879.16±27.680.029显著 w(AP)/(mg·kg-1)45.07±22.1842.53±29.7450.74±9.7352.53±35.1452.96±37.750.003极显著 w(AK)/(mg·kg-1)120.76±73.00147.29±84.37141.35±114.34195.16±93.24219.52±145.890.001极显著

SCC为物理性黏粒,BD为容重,SWC为含水量,CP为毛管孔隙度,TPO为总孔隙度,NCP为非毛管孔隙度,CA为通气度,SOM为有机质,SSC为盐分,TN为全氮,TP为全磷,TK为全钾,AN为速效氮,AP为速效磷,AK为速效钾。P≥0.1表示呈不显著相关;0.05≤P<0.1表示呈一般显著相关;0.01≤P<0.05表示呈显著相关;P<0.01表示呈极显著相关。

2.3 MDS的指标选择

前7个主成分特征值大于1,累积贡献率为72.27%,可将18个土壤指标分为7组(表3)。研究区土壤质量评价的MDS为总孔隙度、全磷、速效氮、Mg2+、全钾、有机质、含水量、速效磷和物理性黏粒9个指标,且指标间相关系数均小于0.5(表4)。这些指标中,总孔隙度包含大孔隙和小空隙,大孔隙通水透气,小空隙持水保水;全磷、速效氮、全钾、有机质和速效磷为猕猴桃生长必需的养分;Mg2+为猕猴桃生长所需的微量元素,同时能调节根部渗透作用;含水量直接影响土壤养分有效性、无机阴阳离子吸收和各种物理化学反应的顺利进行;物理性黏粒能够防止土壤沙化,保持土壤活性。

表3主成分载荷矩阵与MDS的确定

Table3Principalcomponentloadmatrixandthedeterminationoftheminimumdataset

指标主成分1234567Norm值公因子方差分组MDS SCC-0.14-0.110.160.430.110.090.771.030.867是 BD-0.660.03-0.170.030.16-0.180.041.310.531 SWC0.170.350.290.36-0.59-0.34-0.271.120.905是 CP0.83-0.350.190.120.20-0.14-0.121.710.931 TPO0.90-0.300.150.060.08-0.08-0.041.780.941是 NCP0.560.13-0.09-0.22-0.420.230.271.110.691 CA0.85-0.420.05-0.110.240.110.071.760.991 pH0.22-0.090.050.58-0.300.310.120.790.594 SOM0.240.410.31-0.56-0.30-0.060.121.130.754是 SSC-0.53-0.150.10-0.050.260.100.070.750.401 Mg2+0.290.21-0.72-0.030.210.15-0.211.310.763是 SO42--0.33-0.260.55-0.220.270.11-0.191.200.653 TN0.160.570.180.430.250.27-0.121.170.712 TP0.270.69-0.18-0.240.18-0.120.251.300.752是 TK0.240.30-0.800.050.100.100.001.380.813是 AN0.120.640.410.380.240.12-0.131.320.832是 AP0.210.180.010.040.25-0.770.211.020.786是 AK0.130.480.51-0.440.200.230.061.260.803 特征值3.562.422.281.391.191.151.02 贡献率/%19.7713.4412.657.726.616.415.68 累计贡献率/%19.7733.2145.8653.5860.1966.6072.28

表4主成分不同分组内土壤指标之间Pearson相关系数

Table4PrincipalcomponentsPearsoncorrelationcoefficientbetweensoilindicatorsindifferentgroups

分组指标Pearson相关系数BDCPTPONCPCA1CP-0.14TPO-0.18 0.32∗∗NCP0.00-0.07 0.28∗∗CA-0.140.130.31∗∗0.13SSC-0.12-0.21∗0.030.160.08分组指标Pearson相关系数TNTP2TP 0.73∗∗AN-0.26∗∗-0.11分组指标Pearson相关系数Mg2+SO42-TK3SO42- 0.11TK-0.06-0.15AK-0.14-0.08-0.29∗∗分组指标Pearson相关系数pH4SOM 0.45∗∗

BD为容重,CP为毛管孔隙度,TPO为总孔隙度,NCP为非毛管孔隙度,CA为通气度,SOM为有机质,SSC为盐分,TN为全氮,TP为全磷,TK为全钾,AN为速效氮,AK为速效钾。*表示P<0.05,**表示P<0.01。

2.4 土壤质量评价

MDS指标权重见表5。有机质、全氮、全钾、速效氮和速效磷都反映土壤养分,在土壤质量评价中权重总占比达53%,说明养分是研究区土壤质量的重要限制指标。

表5MDS指标权重

Table5Minimumfactorweightsofdatasets

MDS指标公因子方差权重 SCC0.860.12 SWC0.900.12 TPO0.940.13 SOM0.750.10 Mg2+0.760.10 TP0.750.10 TK0.810.11 AN0.830.11 AP0.780.11

SCC为物理性黏粒,SWC为含水量,TPO为总孔隙度,SOM为有机质,TP为全磷,TK为全钾,AN为速效氮,AP为速效磷。

速效钾与全钾呈极显著相关(表4),说明研究区范围内全钾能代表土壤钾元素信息。周至县猕猴桃果园土壤物理性黏粒、含水量、总孔隙度、微量元素和养分含量都略偏低,因此MDS指标对土壤质量均表现为正效应,属于升型隶属度函数。不同种植年限猕猴桃园地土壤质量评价指数随种植年限增加而增加,大小依次为0.68(>20 a)、0.64(>15～20 a)、0.64(>10～15 a)、0.64(>5～10 a)和0.60(0～5 a),表明种植年限越高,园地土壤质量越好。种植年限<10 a的园地土壤指标值均有所提升,土壤质量提升较快;种植年限>10～15 a的园地土壤孔隙度、通气度和有机质含量开始降低,而含水量、微量元素和养分持续增加,土壤质量相对稳定;种植年限>15～20 a园地土壤孔隙度、通气度和有机质保持稳定,其他指标增加明显,土壤质量再次提升。随猕猴桃园地种植年限的增加,土壤质量呈现先提升后稳定再提升的趋势。

2.5 MDS的检验

如表6所示,MDS指标与其他指标均呈极显著相关(P<0.01)或显著相关(P<0.05),说明MDS指标一定程度上可反映其他指标代表的土壤质量信息,采用MDS指标评价土壤质量具有可行性。

表6MDS指标与其他土壤指标间的Pearson相关系数

Table6MinimumdatasetfactorandotherPearsoncorrelationcoefficientbetweensoilindexes

MDS指标BDCPNCPCApHSSCSO42-TNAK SCC0.11-0.15-0.03-0.060.04-0.190.130.120.14 SWC0.100.20-0.02-0.15-0.06-0.19-0.20-0.29∗∗-0.09 TPO-0.180.32∗∗0.28∗∗0.31∗∗0.150.03-0.24∗-0.35∗∗0.39∗∗ SOM-0.060.22∗0.050.45∗∗0.45∗∗0.45∗∗-0.010.18-0.17 Mg2+-0.060.010.11-0.060.080.42∗∗0.11-0.06-0.14 TP-0.04-0.27∗∗-0.14-0.21-0.08-0.060.66∗∗0.73∗∗-0.12 TK0.29∗∗-0.16-0.11-0.31∗∗-0.25∗-0.16-0.150.18-0.29∗∗ AN0.03-0.05-0.110.00-0.20-0.05-0.25∗-0.26∗∗0.10 AP-0.160.180.00-0.03-0.21∗-0.21-0.22∗-0.35∗∗0.09

SCC为物理性黏粒,BD为容重,SWC为含水量,CP为毛管孔隙度,TPO为总孔隙度,NCP为非毛管孔隙度,CA为通气度,SOM为有机质,SSC为盐分,TN为全氮,TP为全磷,TK为全钾,AN为速效氮,AP为速效磷,AK为速效钾。*表示P<0.05,**表示P<0.01。

3 讨论

采用主成分分析和相关分析等统计学方法定量分析不同种植年限猕猴桃园地对土壤因子的影响,并建立土壤质量评价最小数据集。主成分分析法被广泛用于土壤质量评价,但是仅靠变量在单个主成分中的载荷值作为选入MDS的依据,会剔除过多指标,而造成具有重要价值的土壤质量信息缺失,笔者引入Norm值则可避免这一问题[1,13,23],提升评价结果准确性。然而,土壤质量还受到其他因素影响,该研究仅依据测定的18项土壤理化性质指标进行分析,而未考虑土地利用方式、生态系统类型、地理区位等其他因素对土壤质量的影响,这将是下一步研究工作的重点。

从不同种植年限猕猴桃园地土壤指标的变化趋势来看,随种植年限增加,物理性黏粒先增加再减小,而容重先减小再增加,毛管孔隙度、总孔隙度、非毛管孔隙度、土壤通气度和有机质则逐渐减小。这说明随种植年限增加,土壤逐渐形成板结。另外,制约土壤质量的决定指标是养分,养分在土壤质量评价中的总权重占比达53%,而研究区土壤养分含量总体偏低,因此提高猕猴桃产量与质量,首先需要提高土壤养分质量,改变耕作方式,改善土壤结构。