基于土地利用变化的黄河流域景观格局及生态风险评估

黄河流域生态环境脆弱,水土流失严重,诊断其景观格局及生态风险对促进流域高质量发展具有重要意义。该研究采用黄河流域2000、2010和2018年的土地利用数据,计算景观格局指数并构建生态风险评价模型,从网格尺度和县域尺度揭示生态风险的时空演化特征和空间关联性。结果表明:1)草地和耕地面积之和占黄河流域总面积的70%以上,2000—2018年,耕地面积减少14 243 km2,未利用地面积减少9 410.3 km2,建设用地面积增加12 179.8 km2。2)2000—2018年,黄河流域斑块数、斑块密度、最大斑块指数、景观形状指数、香农多样性指数值上升,蔓延度指数值下降,黄河流域景观格局趋于复杂化、破碎化和分散化。3)2000—2018年,黄河流域的生态风险值呈上升趋势,较高风险区和高风险区的面积增加,低风险区和较低风险区的面积大幅度减少,生态风险等级大多由低等级向高等级转移。4)2000—2018年,黄河流域生态风险呈现空间正相关,但空间聚集程度和空间分异性减弱。网格尺度的Moran’s I指数大于县域尺度,表明网格尺度的空间正相关性更强。基于上述结果,该研究将黄河流域划分为生态风险重点管控区、严格管控区和一般管控区,并提出差别化的管控建议,以期为流域管理及风险防控提供定量参考和决策依据。     

黄河流域碳生态安全水平空间格局与动态演进

[目的] 分析中国重要的“能源流域”黄河流域的碳生态安全水平空间差异与动态演进,客观反映黄河流域碳生态安全现状及其演变趋势,为碳生态安全格局优化提供参考。[方法] 在驱动力—压力—状态—影响—响应(DPSIR)分析框架下构建黄河流域碳生态安全评价指标体系,借助逼近于理想解的排序方法(TOPSIS)评估碳生态安全水平并运用核密度估计法分析时空演变特征。[结果] ①2012—2021年黄河流域碳汇量和碳排放量逐年增加,碳排放量增长速度快于碳汇量。从空间范围看,碳汇量具有“西高东低”的特点,而碳排放量呈“西低东高”的特征。②黄河流域碳生态安全水平逐年提升,由2012的0.356增长到2021年的0.639。在空间分布上,上游地区碳生态安全水平高于中、下游地区,呈现出“上游领先,中游追赶,下游超越”的特征。此外,碳生态安全等级均趋于良好。③基于核密度估计结果,黄河流域区域间碳安全水平差异逐渐缩小,且上游区域内碳生态安全水平差异性下降快于中、下游地区。[结论] 基于黄河流域碳生态安全现状,建议上游地区继续增强碳汇能力,中、下游地区发展附加值高的技术密集型和清洁能源产业,推动GDP绿色增长。此外,应建立碳生态补偿机制,协调解决跨区域碳生态安全问题,提升全流域碳生态安全水平。     

基于熵值法和灰色预测模型的耕地生态安全诊断

[目的]弥补已有耕地生态安全诊断中指标体系及评价方法的不足,分析耕地生态安全水平。[方法]在界定耕地生态安全内涵的基础上,构建了基于PSR模型的评价指标体系,采用熵值法和灰色预测模型,对四川省耕地生态安全水平进行了诊断。[结果]1999—2013年四川省耕地生态安全水平总体不断提高,耕地生态安全等级经历了"临界安全—较安全"的演变历程,但2013年"较安全"水平不高;压力指数总体上呈现下降趋势,状态指数呈现波动下降趋势,响应指数呈现上升趋势;单位耕地化肥负荷、人均耕地面积、单位耕地农药负荷、土地垦殖率、水土流失程度等是耕地生态安全等级提升的关键制约因素;2014—2018年四川省耕地生态安全水平呈现稳步上升趋势。[结论]基于PSR模型的评价指标体系能更准确地反映耕地生态安全各要素之间的关系;熵值法和灰色预测模型能有效挖掘耕地生态安全存在的问题,适合用于耕地生态安全评价。     

近14年安徽省耕地资源安全的时空特征

[目的]探究近14a安徽省耕地资源安全时空特征,为安徽省社会经济发展和耕地保护提供依据。[方法]从耕地资源数量、质量和生态安全3个方面选取27个评价指标,构建耕地资源综合安全评价体系。[结果](1)安徽省耕地数量和质量安全程度分别呈下降和上升趋势,生态安全等级较为稳定,综合安全呈缓慢上升趋势;(2)耕地资源数量安全总体呈皖东、皖北高于皖南的空间格局;质量安全呈现皖中、皖南高于皖北的特征;生态安全表现为皖中部和皖南较低,其余较高的空间格局;综合安全与数量安全空间分布具有相似性,且在时间上,低安全等级地市呈现"集中—分散"的空间特征,高安全程度地市反之。(3)近14a安徽省耕地资源数量、质量和生态安全的主要障碍因子具有明显的时间变化特征。[结论]安徽省耕地资源安全情况不容乐观,各地区需因地制宜,提高安全保障程度。     

1990-2018年河湟谷地耕地的时空演变

[目的] 研究河湟谷地的耕地变化,为研究河湟谷地土地利用现状、农业发展规划、耕地可持续利用和生态环境保护等提供依据,也为全球耕地变化研究提供参考。[方法] 采用遥感影像处理、空间统计学并结合ArcGIS空间分析等方法多角度地对河湟谷地的耕地格局进行分析,得到1990—2018年河湟谷地耕地时空演变规律,分析了耕地空间变化与海拔、坡度等特殊因子的关系。[结果] 1990—2018年河湟谷地耕地格局有明显变化,耕地面积有所增加,1990—2000年耕地增长速度比较平稳;2000年以后,耕地开垦速度加快,向着高坡度地区发展,出现了更多的坡耕地;2010年后耕地继续向着高海拔地区发展,坡耕地的增加速度有减缓的趋势。[结论] 1990—2018年河湟谷地的新增耕地向着高海拔、高坡度地区发展,耕地重心向西北偏移。     

河南省沿黄干流地区社会—生态—生产景观恢复力演化特征及情景模拟

[目的] 定量研究河南省沿黄干流地区2000—2018年社会—生态—生产景观恢复力的演化特征,指导沿黄地市社会生活,生态资源环境与农业生产活动的协调发展。[方法] 通过构建SEPLR （social-ecological-production landscape resilience,社会—生态—生产景观恢复力）评估框架,运用有序加权平均法设置情景偏好,绘制不同政策引导下恢复力预测图。[结果] ①从总体特征看,2000—2018年河南省沿黄干流地区SEPLR值为0.47～0.52之间,呈上升趋势,城市间空间差异逐渐增大,指数相似地区集中分布。②演化特征上,社会系统恢复力和生态系统恢复力缓慢增强,分别在0.14～0.65,0.08～0.30波动,生产系统恢复力变化微弱,在0.12～0.31波动;空间上,SEPLR及各子系统整体呈“中部强两边弱”的分布态势,层级结构具有集聚性。[结论] 河南省沿黄地市在“维持现状型”“生态优先型”“开发优先型”“农业优先型”4种政策情景下,“生态优先型”恢复力最高,未来黄河流域高质量发展应考虑“生态优先型”政策引导。     

宁夏沿黄生态经济带耕地的时空演变及驱动力

[目的] 研究耕地空间分布格局与演变的驱动机制,为严守耕地保护红线及宁夏沿黄生态经济带制定耕地资源保护政策与措施,优化耕地保护格局提供参考。[方法] 利用土地变更调查数据和社会经济数据,采用核密度、空间自相关、土地利用动态度、转移矩阵、二元Logistic回归模型等方法,分析2010—2018年宁夏沿黄生态经济带耕地的空间分布、时空演变及其驱动机制。[结果] ①研究区耕地的空间分布呈现出明显的差异性,整体以中、低密度区分布为主,且由低、高密度区向中密度区演变;②耕地空间分布呈正相性,高—高聚类区域主要集中在银川平原,低—低聚类区域分布在灵武市东部、沙坡头区中西部等生态相对脆弱区;③研究区耕地的双向流转是以未利用地—耕地—建设用地为主的生态转换系统,耕地数量稳中有增;④耕地动态度变化有明显的空间聚类现象,且与距离城市的距离成正比,与城市规模成反比;⑤自然条件因子是耕地增加的主要驱动因素,社会经济因子是耕地减少主要驱动因素。[结论] 宁夏沿黄经济带耕地保护基本实现动态平衡。耕地变化是以一定的自然条件因子为基础,受经济、人口等社会经济因子综合作用下的效益择优过程。     

黄河流域2000-2020年“三生”空间土地利用转型及其生态环境效应

[目的] 研究黄河流域生产—生活—生态空间（简称“三生”空间）土地利用转型的结构、空间特征及其生态环境效应,定量权衡土地资源开发与生态环境保护之间的关系,为黄河流域的生态保护和高质量发展提供决策支持。[方法] 以黄河流域为例,利用ArcGIS空间分析与制图技术、生态环境质量指数和土地利用转型对生态环境的贡献指数等模型方法。[结果] ①2000—2020年黄河流域“三生”空间土地利用结构表现为生产用地面积减少,生态用地面积相对稳定,生活用地面积不断增加。②2000—2020年黄河流域生态环境质量指数分别为0.450 2,0.450 5和0.456 7,呈稳定上升趋势。③2000—2020年黄河流域生态环境质量指数在空间分布上有明显的差异性,表现为由东南向北部递增。④2000—2020年黄河流域牧草生态用地对流域生态环境质量贡献率最高。20 a间黄河流域生态环境存在改善和恶化两种趋势,前者大于后者,流域生态环境不断改善。[结论] 城乡生活用地和工矿生产用地对农业生产用地的占用是黄河流域生态环境恶化的主要因素,应合理规划“三生”空间,优化国土资源空间配置,落实黄河流域生态保护和高质量发展国家重大战略。     

松嫩高平原黑土区耕地资源安全及其驱动力分析——以黑龙江省巴彦县为例

耕地资源安全是保护耕地和保障国家粮食安全的基础。选取我国重要的商品粮食基地——松嫩高平原黑土区巴彦县为研究区,从自然、经济、社会和生态共4个方面构建耕地资源安全评价指标体系,运用最优组合赋权法确定指标权重,利用综合指数法定量评价巴彦县2000—2009年间耕地资源安全状态,并采用SPSS软件对耕地资源安全驱动力进行主成分分析,引入Cubic曲线,通过非线性回归分析建立耕地资源安全变化模型,对研究区耕地资源安全值及其驱动力进行结果验证。研究结果表明,2000—2009年,巴彦县耕地资源安全状态存在一定差异。2000年研究区耕地资源安全处于较不安全等级;2001—2003年处于极不安全等级;2004—2005年处于较不安全等级;2006年处于稍不安全等级;2007—2009年处于相对安全等级。研究期间耕地资源安全级别转换呈先负向后正向发展。科技经济社会发展动力、自然干扰力和生态质量影响力是促使巴彦县耕地资源安全变化的驱动力。巴彦县耕地资源安全与3类驱动力表现出非常显著的Cubic曲线关系,说明构建的耕地资源安全变化模型符合研究区实际。     

“双碳”目标下黄河下游耕地绿色利用效率、收敛性与影响因素

[目的] 分析黄河下游耕地绿色利用效率及其收敛性与影响因素,为黄河下游地区耕地生态保护与高质量发展提供参考。 [方法] 将耕地资源利用过程中产生的耕地面源污染和碳排放等非期望产出和总碳汇同时纳入评价指标体系,利用EBM模型、收敛模型、面板Tobit模型系统分析了2007—2020年黄河下游34个地市耕地绿色利用效率、收敛性与影响因素。[结果] ①黄河下游耕地绿色利用效率整体呈上升趋势,但2020年末仍未达到最优状态。黄河下游地区耕地绿色利用效率存在空间非均衡性,豫南及鲁北耕地绿色利用效率值较高,而郑州、威海等地市的效率值一直偏低。 ②从效率分解看,黄河下游地区耕地绿色利用效率提升的主要动力源于规模效率,阻力来自纯技术效率。 ③从收敛检验看,黄河下游耕地绿色利用效率均存在σ收敛和β收敛,河南省耕地绿色利用效率的地市间差距明显低于山东省,但山东省的收敛速度较之更快。 ④从影响因素看,农作物种植结构与政府财政支农对黄河下游地区耕地绿色利用效率起促进作用,而城乡收入差距、经济发展水平和耕地机械投入强度起抑制作用。[结论] 未来黄河下游不但要提高耕地绿色利用的综合管理水平,加强地市间交流合作,更要在推进耕地绿色利用发展的同时,重点关注各地区耕地绿色利用效率的收敛速度与发展差距之间的协调。     

黄河下游城市群生态系统健康时空演变特征

[目的] 研究城市土地扩张对黄河下游城市群生态系统健康的影响，并构建科学的生态系统健康评价模型以识别生态系统健康状态时空演变特征，为制定有效生态保护策略提供前提条件。 [方法] 以黄河下游城市群为研究区，构建VORS评价模型，以格网为研究单元对2000，2010和2020年黄河下游城市群生态系统健康水平进行动态评价。 [结果] ①2000-2020年，生态系统健康指数平均值从2000年的0.325下降到2020年的0.320，标准差从2000年的0.106上升到2020年的0.112。这表明研究区生态系统健康水平总体上呈现略微降低的趋势，生态系统健康等级高的地区更加健康，而等级低的地区生态系统不断恶化。 ②研究区生态系统健康水平空间异质性显著，各等级空间分布集中，从研究区中部延伸到东北部的大部分区域，生态系统健康程度表现为较差，占比高达约74%，整体生态系统健康水平脆弱。 ③2000，2010和2020年全局莫兰指数分别为0.512，0.522和0.543。这表明研究区生态系统健康分布在整体上具有正向的空间集聚效应，空间依赖性显著。 [结论] 黄河下游城市群生态系统健康水平整体上有所降低，生态系统健康状况不容乐观，生态环境保护与建设亟待加强，未来应着重于提高生态系统服务水平。     

基于景观结构的黄河沿岸土地利用生态风险时空变化分析——以河南省为例

[目的] 探究黄河流域沿岸城市土地利用生态风险时空变化特征,为土地景观资源保护、生态系统优化与生态环境保护、管控措施制定等提供理论依据。[方法] 以河南省黄河沿岸7个城市为研究区,在GIS与Fragstats软件技术支持下,基于2005,2010,2015以及2018年4期土地利用数据,将研究区划分为1 737个6 km×6 km的生态评价单元,根据景观扰动指数与景观脆弱指数构建生态风险评价模型,借助土地转移矩阵与空间自相关分析方法,对研究区内土地利用变化、生态风险时空变化特征及空间关联格局进行评价。[结果] ①2005—2018年,研究区域内建设用地面积增加,且增长趋势较快;其他几类用地面积均有不同程度的减少;耕地是建设用地增加的主要来源,林地、草地和水域主要转化为耕地。②在研究时段内,研究区域土地利用生态风险在空间上呈现出显著的正相关性,存在空间集聚特征,且“高—高”、“低—低”是生态风险主要的空间聚集模式。③2005—2018年,生态风险最小值和最大值大致呈N形变化;4个时期内高风险区分布区域大致相同,呈条状和块状分布,条状高风险区大多分布在黄河沿岸,块状高风险区大多分布在农村居民点零散分布的平原区;低和较低风险区大多分布在研究区域四周以及西南大部分区域。[结论] 研究时段内,低风险区和中等风险区面积在逐年增加,较低风险区面积先减少后增加,高风险区和较高风险区面积在逐年减少。各风险区主要向相邻风险区转化,仅有少量地区出现跨区转移,说明研究区土地利用生态风险变化相对稳定,生态风险急剧变化的区域较小。     

2007-2017年洞庭湖平原耕地生态安全评价及障碍因子诊断

[目的] 研究洞庭湖平原耕地生态安全状况,并诊断其障碍因子,为区域耕地合理利用和耕地污染治理提供参考依据。[方法] 基于DPSIR模型构建洞庭湖平原耕地生态安全评价指标体系,从县域尺度,采用基于熵权的改进TOPSIS模型测度洞庭湖平原2007—2017年耕地生态安全指数,并采用障碍度模型诊断耕地生态安全状态的主要障碍因子。[结果] 2007—2017年洞庭湖平原整体耕地生态安全程度处于临界状态,整体安全指数在研究期内波动上升但增速缓慢;不同县（区、市）耕地生态系统安全变化特征不一;中低产田改良比例、秸秆还田比例、测土配方施肥比例、农药使用强度、地膜使用强度、化肥施用强度、农业自然灾害成灾比例等是制约耕地生态安全的主要障碍因子。[结论] 为改善洞庭湖平原耕地生态安全状况,应大力实施全域土地综合整治,加强农业面源污染治理。     

基于BP神经网络的工业生态安全动态评价及障碍因子诊断

[目的] 对安徽省工业生态安全进行评价和预测,为安徽省工业生态安全可持续发展提供科学依据。[方法] 基于压力、状态、响应框架和生态、环境、经济、社会框架构建工业生态安全评价指标体系,采用聚类工具进行安全等级划分,结合熵权法和综合指数法评价2009—2018年安徽省工业生态安全水平,运用BP神经网络模型对2019—2025年工业生态安全水平进行科学预测。[结果] ①从总体上看,2009—2018年安徽省工业生态安全评价指数呈波动上升趋势,安全等级由“临界安全”上升为“较安全”;②从各子系统来看,2009—2018年压力子系统的评价指数呈上升态势,状态和响应子系统的评价指数呈先下降后上升态势;③BP神经网络预测结果显示：2019—2025年安徽省工业生态安全态势总体呈波动上升状态,处于“临界安全”向“较安全”转变态势;④影响安徽省工业生态安全的主要障碍因素包括第二产业从业人员比例、第二产业占GDP的比例、污水集中处理率、森林覆盖率和人口密度,是今后调控的方向。[结论] 研究期间安徽省工业生态安全波动明显,但总体上呈上升趋势,工业生态安全得到明显的改善。     

2000-2018年黄河流域NDVI时空变化及其对气候和人类活动的双重响应

[目的] 分析黄河流域NDVI时空变化特征,探究NDVI变化对气候和人类活动的响应机制,为制定合理的生态工程提供科学依据,为实现黄河流域生态保护和高质量发展提供保障。[方法] 基于MODIS NDVI数据,辅以降水和气温数据,采用一元线性趋势分析、Hurst指数、偏相关分析及残差分析方法,对2000—2018年黄河流域NDVI时空变化特征进行分析,并探讨了NDVI变化对气候和人类活动的双重响应。[结果] ①黄河流域NDVI呈现波动增加趋势,总增速为6.8%/10 a,NDVI在东南部及西部较高,北部及西北部较低,由东南向西北部减少,下游最高,中游次之,上游最低;整体上,以东亚季风生态地理区和西北干旱生态地理区界线（鄂尔多斯市—毛乌素沙地—庆阳市—平凉市—定西市）和青藏高原生态地理区和西北干旱生态地理区及东亚季风生态地理区界线（西宁市—甘南市—定西市）为界呈Ⅴ字型分布。②NDVI整体变化趋势以改善为主,反持续性强于持续性,且表现出较强的弱持续性,其中改善区约占62.32%,主要呈带状集中分布。在中上游地区,在鄂尔多斯市—毛乌素沙地—庆阳市—平凉市—定西市—西宁市—甘南市、太原盆地—临汾盆地—关中平原及运城盆地—三门峡两两界线间的条带区域基本为改善区,而退化区域以中下游为主,呈零散分布。③NDVI与降水和气温呈正相关,且与降水的偏相关强度稍大于气温,流域约76.7%的地区NDVI残差呈增长的趋势,人类活动整体对NDVI增加的影响以促进为主,促进作用大小顺序为：中游>上游>下游。[结论] 黄河流域NDVI变化受气候和人类活动共同影响,具有空间异质性,适度的生态工程及农业生产活动对植被恢复有重要作用。     

基于GIS和模糊数学的川西地区生态安全评价

[目的] 对川西生态安全进行评价研究，以期为该地区生态环境建设、社会经济发展规划提供科学依据。[方法] 基于GIS和遥感技术，从自然和人为影响两个方面构建评价指标体系，通过模糊数学模型对川西地区2000，2005，2010，2015年和2018年的生态安全状况进行评价。[结果] 川西40%的区域生态安全状况处于较危险状态，呈现出“东高西低”的空间分布格局。川西生态安全2000—2018年总体变化趋势为“上升—小幅度下降—上升”。川西地区生态安全Moran’s I指数在0.5左右，西部白玉、巴塘、理塘一带生态安全状况好转，东部九龙、雅江、木里和若尔盖一带生态安全质量有所下降。地形植被是川西生态安全最主要的影响因子。[结论] 川西整体生态安全状况偏差，但呈现出改善趋势，具有较强的空间聚集性，对川西生态安全状况的治理要充分考虑其空间分异特征。     

黄河流域下游城市群生境质量时空演变及其驱动因素

[目的] 对黄河流域下游城市群生境质量及其影响因素开展研究,为该地区生态环境保护和高质量发展提供参考。[方法] 基于2005—2020年土地利用数据并使用InVEST模型分析区域生境质量的演变特征,通过空间分析探索生境质量的集聚效应,采用地理探测器探索生境质量的影响因素。[结果] ①2005—2020年,黄河流域下游城市群生境质量水平整体较低,平均生境质量指数为0.402 2,呈先下降后上升的趋势。②空间自相关分析显示,生境质量在空间上高度相关且呈现出两极集聚特征,生境质量在中原城市群西部、北部及山东半岛城市群中部高—高集聚;在郑州市、青岛市等在城市群中经济发展突出的城市低—低集聚。③地理探测器探测结果显示,土地利用是影响黄河流域下游城市群生境质量的主要因子,不同因子间交互作用强于单个因子对生境质量的影响。[结论] 黄河流域下游城市群生境质量高值区要继续加强对生态环境保护的力度,生境质量低值区要严格控制土地开发强度,杜绝建设用地无序扩张,提高土地的利用效率。     

2005-2018年澜沧江下游流域景观生态安全时空变化及其驱动因素

[目的] 分析澜沧江下游流域景观生态安全时空变化规律和驱动机制，为该流域及流经国家的人与自然协调可持续发展提供科学依据。[方法] 以2005，2010，2015和2018年4期遥感影像为数据源，通过景观指数和GIS空间分析方法评价流域景观生态安全，采用双变量空间自相关方法明确城市、道路、河流、高程、坡度、气温和降水等因素对景观生态安全的驱动作用。[结果] ①2005—2018年，澜沧江下游流域的整体景观生态安全度呈先升后降的变化特征。其中，北部和中西部景观生态安全较差，东部、南部和东南部较好，但也有恶化趋势；②各景观类型中，景观生态安全度从大到小为：水域 > 林地 > 草地 > 耕地 > 建设用地 > 裸地。其中，林地和草地的景观生态安全度存在波动，耕地和水域的景观生态安全度一直上升，而建设用地和裸地的景观生态安全度持续下降；③距城市距离、距道路距离、气温和降水等社会和气象因素对景观生态安全时空变化的驱动作用最强，距河流距离也对景观生态安全的时空变化有一定解释作用，而高程和坡度等地形因素并非澜沧江下游流域景观生态安全时空变化的主要驱动因素。[结论] 应尽量降低由人类活动引起的生态环境负荷，积极发挥气象因素对生态安全的正向效应，促进澜沧江下游流域的人与自然协调可持续发展。     

基于“贡献力—活力—组织力—恢复力”的武汉市自然资源生态安全评价

[目的] 构建生态系统健康与生态系统服务关联的自然资源生态安全评估框架,分析城市化过程中自然资源生态安全的演变规律,为区域国土空间优化和生态文明建设提供理论依据。[方法] 基于生态系统健康与生态系统服务关联的视角,构建“贡献力—活力—组织力—恢复力”的自然资源生态安全评估框架,综合运用InVEST模型、景观指数和德尔菲法评估武汉市2000和2021年的自然资源生态安全,并分析其时空演变规律。[结果] ①2000和2021年武汉市自然资源生态安全指数分别为0.394 4和0.374 5,生态安全等级在空间上呈现中部低、四周高的分布格局,且黄陂区生态安全水平最高; ②2000—2021年武汉市自然资源生态安全整体呈下降趋势,下降的地区主要沿中心城区边缘呈环状分布,“低安全”范围有从中部核心区逐步向外扩张的趋势,且“中等安全”降为“较低安全”,“较高安全”降为“中等安全”的分布范围较广。[结论] 城市扩张过程中面对自然资源生态安全不断下降的压力,应进一步加强国土空间规划的引领作用,严惩非法侵占生态用地,加强生态空间的保护。     

粮食主产区耕地生态安全障碍诊断及调控策略

生态安全作为耕地“三位一体化”的核心保护任务,越来越得到政府、公众及学术界广泛关注,耕地生态安全评价视为关键技术及重要依据。基于PSR模型构建耕地生态安全评价指标体系,以国家粮食主产区安徽省为例,综合运用改进熵权法、TOPSIS模型、ESDA模型及障碍诊断模型,评价安徽省2000年以来,以及分市级2000年、2010年、2018年3个年份的耕地生态安全度,剖析耕地生态安全时空分异特征,进一步诊断耕地生态障碍因子及对应调控策略。结果表明：(1)安徽省耕地生态安全度具有一定波动性,总体呈现上升趋势,由2000年的0.234 1增至2018年的0.724 7,大致分为低位震荡阶段、快速增长阶段及高位爬坡阶段;(2)近20年来耕地生态安全状态进一步呈现不平衡、不均衡演变趋势,导致全省耕地生态安全空间格局差异显著,尤其是以高风险、安全两类区域分布空间集中区域发生很大改变,基本呈现皖南地区安全程度高并由皖中地区过渡至皖北地区风险程度高的空间梯度格局,耕地生态安全空间集聚情况并未显著;(3)从全省维度分析,障碍因子以建成区绿化覆盖率、污水处理率等响应准则层为主,逐渐转变至复种指数、单位耕地粮食产量、...