干旱区耕地质量等级时空变化及其评价-以西宁市为例.docx

    干旱区耕地质量等级时空变化及其评价以西宁市为例    潘 雪，关宇淇，潘占东，刘 杰，蔡立群，董 博，杜 健（1.黑龙江省黑土保护利用研究院，黑龙江 哈尔滨 150000；2.甘肃农业大学资源与环境学院，甘肃 兰州 730070；3.甘肃省农业科学院旱地农业研究所，甘肃 兰州 730070）耕地是保障农业生产和农产品供应的基础，随经济和工业化的快速发展，中国耕地出现不同程度的退化，影响作物安全生产和国家农业生产结构的总体战略布局［1-4］现阶段，消费者对农产品的需求由数量转变为安全、绿色、优质［5］；因此，保护耕地、提升耕地质量刻不容缓摸清现有耕地资源状况对于中国生态、经济、民生发展均起到至关重要的作用，而耕地质量评价是调查和摸清现有耕地资源状况的重要手段［6］为保障粮食供应安全，国家划定了1.2×108hm2耕地红线，出台有关加强耕地保护和征地-补贴平衡政策［7］2016年，农业农村部下发《耕地质量等级》（GB/T33469-2016）（简称“标准”），规范和统一了现阶段耕地质量评价的方法和步骤，并要求全国开展新时期耕地质量评价工作，当前耕地质量评价主要从耕地质量等级、时空变化和主要养分特征3个方面开展工作。

姚东恒等［8］对东北典型黑土区2008 年和2018 年耕地质量进行评价，探讨该地区耕地质量等级变化和空间分布特征，10 a间，耕地质量等级提升了1.22个等级；康钦俊［9］对全国九大农区之一的黄土高原区2018 年耕地质量等级及土壤养分空间变化进行评价，初步掌握该区域的土壤养分和耕地质量状况，分析影响土壤养分空间差异的主要因素并给出了该区域未来耕地利用的建议；戴文举等［10］以吴川县为例，对县域耕地质量和产能展开了评价近年来，耕地质量评价在全国逐步开展，但依然存在一些问题，如：目前已有的评价成果，缺乏统一标准，并不能很好地反映实际耕地质量现状，需要及时更新和补充相应的数据资料另外，中国地形、气候复杂多样不利于采取大面积、相同方式的耕地调查就该领域目前研究现状而言，东部沿海等经济发达地区耕地质量提升、土壤培肥等研究较多，大多数研究体系以大区域或者省级为基础，关于西部地区、小范围市级、县级耕地质量时空变化的研究较少而青海省西宁市是青藏高原地区人口超过百万的省会城市，也是国务院批复的西北地区重要的中心城市［11］，自国家“标准”颁布和辖区更名以来，西宁市全域范围内耕地质量评价工作尚未见报道。

同时，青海省耕地资源数量少，主要集中在西宁市和海东地区，整体质量偏低，且水土流失严重，因此准确掌握研究区域耕地质量时空变化及养分变化特征对青海省农业发展尤为重要［12-14］本文以青海省西宁市作为研究对象，以全市的耕地质量监测点为坐标，依据“标准”，结合GIS技术，构建耕地质量评价指标体系，对其耕地质量等级进行科学、合理评价，为该区域进一步调整农业结构、提升耕地质量等提供技术支撑［15］1 研究区概况青海省西宁市（101°52＇～101°54＇E，36°13＇～37°28＇N）包括城中区、城东区、城西区、城北区、湟中区、湟源县、大通回族土族自治县（简称大通县）5区2 县，地处青海省东部、青藏高原东北部，地势西南高、东北低，东西向呈条带状分布，是西部大发展的重要城市之一该区域属大陆性高原半干旱、高原高山寒温性气候全年平均日照时数2510.1 h；年平均气温5.5 ℃，最高气温34.6 ℃，最低气温-18.9 ℃；年平均降水量493.4 mm，蒸发量1363.6 mm2018年，西宁市总面积76.60×104hm2，耕地面积14.47×104hm2，其中水浇地2.72×104hm2，旱地11.75×104hm2；常住人口238.71×104人，农村人口64.81×104人，占总人口的27.15%［16］。

2 数据与方法2.1 数据来源与处理本研究共布设调查样点310 个（图1）资料主要包括西宁市土地利用现状栅格数据图、土壤类型栅格数据图、行政区划矢量数据图，以及2015 年和2018 年土壤样品采集分析数据土壤理化指标参考《土壤农化分析》、“标准”规定方法检测，并建立耕地质量等级评价数据库［17］将土壤类型栅格数据图、行政区划矢量数据图和土地利用现状栅格数据图叠加形成评价单元图，对评价单元属性赋值，建立耕地资源管理信息系统［18-20］图1 采样点分布Fig.1 Distribution of sampling points2.2 评价单元确定评价单元作为评价耕地质量的最小单元，评价单元的合理性直接关系到评价结果的准确性和合理性；因此评价单元的基本条件、单元属性应基本一致，又要存在差异本研究采用叠置法和地块法，将土壤类型栅格数据图、行政区划矢量数据图和土地利用现状栅格数据图三者叠加形成图斑底图，将底图按照相关技术标准进行处理，共生成36762个评价单元，经过处理后的图斑属性唯一［21］2.3 评价指标体系及权重确定评价指标的选取直接影响评价结果的准确性和实际性，在选择的评价指标中，要既能够反映本区域特点，在区域内存在较大变异，又在时间序列上能够保持相对稳定，并且容易获得［22］。

耕地质量受很多因素的综合影响，其中自然因素是影响耕地质量最重要的指标［23］如表1 所示，本研究根据“标准”中规定的“N+X”指标体系结构，结合本研究区域的自然条件和实地状况，建立西宁市耕地质量等级评价指标体系［24-26］其中，包括目标层、准则层和指标层3级结构，采用特尔斐法构造各指标的判断矩阵，后经计算得出各指标权重，将各指标权重按照准则层的分类标准分别相加后得出组合权重［27-28］表1 西宁市耕地质量评价指标及权重Tab.1 Cultivated land quality evaluation indices and weights in Xining City2.4 确定评价指标的隶属度根据模糊数学的基本原理，并结合西宁市耕地质量评价过程中部分评价指标的数值特性，拟合出峰型、戒上型、戒下型3 种隶属函数，其他概念型函数直接采用特尔斐法给出相应隶属度［29-30］其结果如表2、表3所示表2 西宁市耕地质量评价指标隶属函数Tab.2 Membership functions of cultivated land quality evaluation indices in Xining City表3 西宁市耕地质量评价指标属性及其隶属度Tab.3 Attribute and membership degree of cultivated land quality evaluation indices in Xining City2.5 耕地质量等级划分及综合指数计算根据“标准”，采用累加法计算各评价单元的耕地质量综合指数，计算公式［31-32］如下：式中：IFI为耕地质量综合指数；Ci为第i个评价指标的组合权重；Fi为第i个评价指标的隶属度。

按照“标准”的规定，等间距将耕地质量等级划分为10个等级IFI越大，耕地质量越高，西宁市耕地质量综合指数范围在二至十等之间为确保评价结果的与实际现状相符，使其能够对当地农业生产结构调整起到指导作用，对评价结果进行了实地验证和专家论证；最终结果如表4 所示研究区域平均耕地质量等级采用面积加权平均法计算，计算公式［27］如下：表4 西宁市耕地质量综合指数划分标准Tab.4 Classification standard for comprehensive index of cultivated land quality in Xining City式中：AI为耕地质量平均等级；I为耕地质量等级；A为各等级耕地所占面积（hm2）；SA 为研究区总耕地面积（hm2）3 结果与分析3.1 耕地质量等级时空变化3.1.1 耕地质量等级面积变化由表5可知，西宁市2015年耕地总面积为144707.79 hm2，高、中、低等地面积由低到高依次为6237.12 hm2、51460.56 hm2、87010.11 hm2，占比分别为4.31%、35.56%、60.13%；2018 年高、中、低等地面积由低到高依次为10418.15 hm2、49428.85 hm2、84860.79 hm2，占比分别为7.20%、34.16%、58.64%；2018 年高等地面积较2015 年增加2.89%，中、低等地面积分别减小1.40%、1.49%。

高等地中，三等地2015、2018年占比分别为97.84%、99.48%，2018年面积较2015年增加2.94%中等地中，四等地2015、2018年占比分别为28.66%、35.72%，五等地2015、2018 年占比分别为43.33%、52.03%，六等地2015、2018 年占比分别为28.02%、12.25%；2018 年四、五等地面积较2015 年分别增加2.01%、2.36%，六等地面积下降5.78%低等地中，七等地2015、2018 年占比分别为20.72%～44.70%，八等地2015、2018 年占比分别为57.34%～45.22%，九等地2015、2018 年占比分别为18.90%～10.08%；2018 年八、九等地面积较2015 年分别减小7.96%、5.46%，七等地面积增加13.76%西宁市无一等地，二等地和十等地占比低于2%，影响极低西宁市平均耕地质量等级进行计算可知，2015年和2018 年平均耕地质量等级分别为6.74 和6.33，2018年较2015年质量等级提升0.41表5 西宁市耕地质量等级面积变化Tab.5 Area change of cultivated land quality grades in Xining City3.1.2 主要区县耕地质量等级变化由图2可知，西宁市2015 年的高等地主要集中在大通县东南部的长宁镇和城北区大堡子镇，湟中区中部区较少；2018年主要分布在湟中区汉东乡及城北区，4 a间高等地的面积减小。

中等地2015年和2018年在研究区均有分布，4 a 间区域变化不大，主要集中在中部的湟中区和北部大通县的中南部低等地分布较广，各县（区）均有分布，主要在高、中等地四周及边缘靠山区，其中2015 年以湟中区和大通县居多，2018年湟中区和湟源县分布较多图2 2015年和2018年西宁市耕地质量等级空间分布Fig.2 Spatial distributions of cultivated land quality grade in Xining City in 2015 and 2018由图3和表6可知，2018年与2015年相比，西宁市耕地质量提升、下降、不变的面积分别占总面积的43.585%、7.231%、49.180%耕地质量等级最高可提升3个等级，仅出现在城中区；等级提升1的耕地面积最大，为56560.41 hm2，占总耕地面积的39.086%，主要集中在中部的湟中区、北部的大通县及西部的湟源县耕地质量等级最多下降2 个等级，出现在大通县和城北区；等级下降1的耕地面积达10427.24 hm2，占总耕地面积的7.206%等级下降2 和提升2、3 的耕地面积较小，仅占评价区总面积的4.524%。

等级不变的耕地面积在7 个县（区）均有分布表6 西宁市耕地质量等级变化面积分布Tab.6 Area distribution of cultivated land quality grade change in Xining City图3 西宁市耕地质量等级变化分布Fig.3 Distribution of cultivated land quality grade change in Xining City3.2 耕地质量等级土壤主要养分变化由表7 。