基于最小数据集的土壤肥力指数评价研究.docxVIP

基于最小数据集的土壤肥力指数评价研究 土壤肥力是指耕地能承受多种养分的能力。科学、合理地评估土壤肥力，是指导农业生产、土地利用规划和土壤肥力的理论依据。 1 材料和方法 1.1 研究区地理位置及土地利用方式 洞庭湖位居我国五大淡水湖第二位，南部平原区为湖南农业主产区，以种植水稻、棉花为主，素有“鱼米之乡”的美称，选定研究区为洞庭湖所在的河湖冲击平原区南部，包括岳阳、益阳、常德、长沙等四个市。研究区的行政面积为57216 km2，其中，耕地面积为21093 km2，占总面积的37%。地理坐标为110°29′23″～114°15′30″E,27°51′13″～30°7′35″N，气候属于中亚热带向北亚热带过渡的湿润气候，年均气温在在16～17℃之间，年均降雨量1398 mm。研究区海拔由中心向四周逐渐升高，地貌由中部的冲积平原向四周过渡到低山红壤丘陵区；土地利用方式以水田为主，由中心向四周过渡到林地（图1）。土壤类型主要以水稻土为主，其次还有潮土、红壤、黄壤、黄棕壤、紫色土（图1）。 1.2 土样的采集与处理 综合考虑研究区内耕地的成土母质、土壤类型、地形部位等基础地理空间属性，采用分层抽样结合网格布点的方法进行样点布设，在水田集中地区加密布点，以满足地统计学插值和肥力评价要求。实际取样时记录每个点位经纬度和剖面性状。通过梅花形采样法采集0～20 cm土样，实际取样3400个点，后续处理使用3381个点位（图2）。 土样经风干、磨碎及过筛后，测定各土样的理化指标。土壤质地采用吸管法；土壤容重采用环刀法；土壤阳离子交换量（CEC）采用乙酸铵交换法；土壤pH采用1∶2.5土水比浸提-电位法；有机质采用重铬酸钾容量法；碱解氮采用碱解扩散法；全氮采用凯氏定氮法；速效钾采用NH 1.3 学习方法 1.3.1 立地条件和土壤理化性质的选取 本研究综合考虑了立地条件、剖面性状、土壤物理性状、土壤养分等对土壤肥力有显著影响的四个方面，共选取了14个指标，立地条件选取了成土母质、地形部位；剖面性状选取了土体构型；土壤物理性状选取了耕层质地和容重；土壤养分选取了pH、有机质、全氮、碱解氮、有效磷、速效钾、CEC、全磷和全钾。 评价指标权重的计算方法如下：通过相关性矩阵对参评指标进行主成分分析，指标的权重为该指标公因子方差占所有参评指标公因子方差和的比值 式中，N 1.3.3 耕地等级标准定量型指标的隶属度函数 数据处理与评价需要对定性型指标进行量化，故结合耕地质量等级标准 定量型指标通过隶属度函数转换 升S型隶属度函数计算如下： 式中，转折点x 抛物线型隶属度函数计算如下： 式中，转折点x 1.3.4 计算全球土壤生产力指数 土壤综合肥力指数（Soil Fertility Index,SFI）的计算如下： 式中，q 1.3.5 处理数据 本研究利用SPSS 16.0软件进行数据统计分析，GS 2 结果分析 2.1 土壤理化性质与结构及养分供应能力分析 对定量型指标进行描述性统计分析，结果如表3所示。土壤容重的均值为1.22 g cm 定量型指标的半方差函数模型及其参数结果如表4所示，土壤容重与有效磷的R2较低，理论模型拟合效果一般，其余指标的拟合效果均比较满意。p H的变程最大，说明土壤母质在研究区内的分布比较均一，容重的变程最小，说明耕作等人为活动对容重的影响比较大。除pH外，其余指标的块基比均小于0.25，说明各指标均属于强空间相关性，主要受结构性因子影响；pH的块基比略大于0.25，属于中等程度的异质性，说明随机因素的影响有所加强。 根据前面得到的各指标的半方差函数模型进行克里格插值，得到了各指标在研究区内的空间分布图，同时对各指标分级并进行面积统计，结果如表5所示。研究区内耕地土壤偏微酸性、酸性和中性的面积占比分别为32.12%、25.40%和22.99%，处于微碱性和碱性的耕地面积占比约20%；研究区内土壤容重主要处于偏紧级别，面积占比74.96%，处于适宜级别的面积占比24.06%，其余较松级别占比小于1%。土壤CEC主要处于1级水平，面积占比71.51%，说明研究区土壤保肥能力强。土壤有机质、全氮和碱解氮主要处于2级水平，面积占比分别为52.46%、51.93%和48.46%，可见研究区内碳氮养分含量较高，土壤肥力供应能力强，供肥潜力大。土壤全磷、有效磷主要处于3级水平，面积占比分别为45.44%和54.86%，其次处于2级水平的耕地面积占比也较大，说明研究区内磷素养分含量较高，磷素供应能力较充足。土壤速效钾和全钾主要处于4级水平，面积占比分别为48.27%和52.12%，其次处于3级水平的耕地面积占比也较大，说明研究区内土壤钾素养分相对缺乏，是作物生长发育的主要限制因子。 2.2 土壤肥力综合评价的最小数据集的构建 2.2.1 载荷