基于生态位理论的富硒土壤资源开发利用适宜性评价及分区方法

实施乡村振兴战略,产业兴旺是关键。立足区域富硒土壤资源禀赋优势,发展富硒农业特色经济,是振兴乡村产业的一条有效途径。以多维超体积生态位理论为指导,以区域富硒土壤资源禀赋和富硒农业发展要求为基础,从土壤自身条件、土地利用和生态保护3个方面明确了富硒土壤资源开发利用的主导因素、区位因素、限制因素;选择土壤硒含量、土壤有机质含量、地形坡度、土地利用类型、土地综合污染指数、生态保护建设等6个关键因子,构建了富硒土壤资源开发利用生态位适宜性评价指标体系;利用生态位适宜性评价相关模型,探讨了生态位适宜性单一指数和综合指数的计算方法;结合土壤利用改良分区、土壤利用功能分区、富硒农业生产分区等,分析了区域富硒土壤资源开发利用适宜性分区问题。研究可为富硒土壤资源开发利用适宜性评价提供科学可行的思路与方法,对于促进区域富硒土壤资源科学合理利用、实现富硒农业产业科学发展具有重要的指导意义和现实意义。     

东北粮食主产区耕地生态安全的时空格局及障碍因子——以吉林省为例

针对东北粮食主产区耕地生态安全状况,以吉林省48个县域为研究对象,构建基于PSR-EES模型的耕地生态安全评价指标体系,综合运用改进的TOPSIS模型、马尔可夫链模型、GIS空间分析方法和障碍度模型对1995-2011年吉林省耕地生态安全时空格局以及障碍因子进行分析.结果表明： 1995-2011年,吉林省耕地生态安全综合指数呈波动上升趋势,耕地生态安全等级从敏感级提高到一般安全级;各县域耕地生态安全水平存在“俱乐部趋同”现象,且空间分布格局差异趋于扩大,“两极分化”趋势日益显著,敏感级和临界安全级由带状趋于分散分布,一般安全级和比较安全级由分散趋于集中成片分布,安全级分布比较零散,耕地生态安全等级水平不稳定趋势越来越明显;农民人均纯收入、经济密度、环境保护投资总额占GDP比重、机耕程度、工业固体废物综合利用率是影响吉林省耕地生态安全水平的主要障碍因子.     

上海市南汇东滩滩涂围垦区农业生态安全评价

针对上海市南汇东滩围垦区农业开发特点, 以压力-状态-响应(PSR)模型框架作为研究方法, 构建了反映和体现滩涂围垦区农业可持续发展对生态环境影响的指标体系。运用层次分析法(AHP)和综合加权评分法对南汇东滩围垦区农业发展进行生态评价, 并对未来农业可持续发展安全进行预测。结果表明, 南汇东滩围垦区围垦初期(2009-2012)农业生态安全综合评价指数范围为4.1-4.8, 其安全状况处于临界安全状态。2015-2050年的生态安全综合评价指数从5.3到7.0, 都高于围垦初期, 其中2020年将达到良好安全状态, 主要原因是社会响应的增加。可见, 在政府引导和政策支持下, 企业自主经营的模式有利于推动围垦区生态农业的建设和发展, 可直接为南汇东滩围垦区农业发展提供科学依据, 其研究方法可为类似围垦区农业生态安全评价借鉴。     

甘南藏族自治州生态足迹与生态承载力分析

运用生态足迹的各种计算模型和评价指数, 对甘南藏族自治州2000-2015 年的生态足迹、生态承载力情况进行动态测算, 并且利用回归分析对2016-2025 年的生态足迹与生态承载力变化趋势进行预测。研究表明:2000-2015 年甘南州人均生态足迹由1.2 hm²·cap^–1 上升到2.2hm²·cap^–1, 万元GDP 生态足迹逐年下降, 人均生态承载力变化较小, 呈微减趋势; 整体上处于生态盈余状态, 但盈余在逐年递减; 生态足迹多样性指数在0.63-0.71 之间平稳波动, 生态经济系统发展能力指数呈快速上升趋势; 生态环境由2000 年的生态很安全强可持续性状态逐渐变为2015 年的较安全弱可持续性状态; 预测显示甘南州自2018 年开始出现生态赤字, 未来发展趋势不容乐观。若不采取有效措施, 其生态可持续发展状况将面临严峻的挑战。     

广东省生态安全状态及趋势定量评价

为研究区域生态安全在不同时期的动态变化及发展趋势,在国内外已有的研究基础上,根据压力-状态-响应模型框架,提出了一套完整的区域生态安全水平度量的指标体系和基于熵权法赋权的综合指数评价方法,从时间尺度上（1990～2004年）对广东省生态安全进行定量评估和动态趋势分析.结果表明,广东省生态安全综合指数从1990年的0.125增加至2004年的0.402,生态安全整体水平呈逐年上升趋势,但发展过程中存在的问题是区域生态压力在不断加大及生态安全现状水平仍处于“较不安全”状态.对广东省环境规划的实施效果进行预测评估,得出在实现环境目标的情景下,2010规划年和2020规划年的生态安全综合指数将分别达到0.533（临界安全）和0.691（较安全）.评价方法简单直观,评价结果客观合理.     

吉林省县域环境污染的时空演变特征与影响因素

基于2011年和2015年环境污染物排放和经济社会数据,选取4类污染物的工业源、农业源和城镇生活源排放总量指标,通过构建整体环境污染综合指数和分来源污染指数综合反映吉林省县域尺度单元污染状况,运用Moran I指数、时空跃迁分析和标准差椭圆分析定量研究吉林省环境污染的时空演变特征,并构建空间计量模型对环境污染的影响因素进行分析.结果表明:吉林省环境污染具有空间集聚特征,环境污染主要分布在长春、吉林等中部地区,分来源的环境污染空间分布存在明显的区域差异特征;吉林省县域单元整体环境污染和农业源污染均具有空间相关性,环境污染空间布局的关联结构稳定性高,环境污染具有路径锁定特征;环境污染的空间分布呈现西北-东南分布格局,且环境污染表现出由中心向四周扩散态势,环境污染的重心在43.65-43.66°N,125.83-125.84°E之间变动;地区经济发展水平、城市化水平和农业生产是吉林省县域环境污染和工业源污染的主要驱动因素,产业结构的调整优化改善了整体环境污染与工业源污染状况.化肥的过度使用增加了农业源污染,经济发展水平和城市化水平是城镇生活源污染排放的主要驱动因素.县域整体环境污染和农业源污染均具有空间扩散效应,各自污染状况分别与相邻地区整体环境污染和农业源污染状况密切相关.     

治沟造地背景下延安市农业生态经济系统耦合发展分析

治沟造地是继退耕还林工程后延安市实施的又一项重大生态环境治理工程。研究治沟造地工程背景下延安市农业生态经济系统耦合关系的发展变化,对科学评价治沟造地工程生态经济效益和黄土高原地区农业生态经济协调与可持续发展具有重要意义。本研究通过构建农业生态经济系统耦合协调度模型和耦合度模型,定量分析了2010—2018年延安市农业生态与农业经济系统的综合评价指数变化、耦合协调状态及耦合度演变趋势。结果表明: 2010—2018年间,延安市农业生态系统综合评价指数和农业经济系统综合评价指数均呈现S型增长曲线,两者变化趋势较为一致;农业生态经济系统耦合协调度指数由0.51增长至0.72,耦合协调状态由初级协调水平发展为良好协调水平,说明该区域农业生态与经济关系不断改善,系统耦合关系向相互协调的方向发展;研究期间,延安市农业生态经济系统耦合演变过程较为复杂,系统耦合趋势经历了“衰退耦合、协调耦合、修复耦合、协调耦合”阶段,2018年耦合关系处于经济快速增长的协调发展阶段。治沟造地工程促进了农业生态经济系统的耦合协调发展,但其生态经济效益具有一定的滞后性。     

1990—2015年三江平原生态功能区水禽栖息地适宜性动态

作为国际重要湿地,三江平原生态功能区是重要的水禽栖息地.随着人类活动干扰、土地利用和全球气候变化,栖息地适宜性逐渐引起生态学家的重视.本文以Landsat TM/ETM+/OLI和HJ-1B为遥感信息源,采用面向对象分类方法提取土地覆被空间信息;采用综合熵值法和层次分析法确定水源状况(湖泊和河流密度)、干扰条件(居民地和道路密度)、遮蔽物(土地覆被类型和坡度)和食物丰富度(NDVI)等因子的权重;根据水禽栖息地适宜性评价系统获得三江平原生态功能区1990、2000、2010和2015年水禽栖息地适宜性结果,并分析其时空分布和变化特征及驱动因素.结果表明: 近25年间,三江平原生态功能区水禽栖息地适宜性最好区域面积减少3.2%,主要由于湿地开垦和退化;适宜性最好区域空间分布特征明显,主要分布于黑龙江、挠力河、乌苏里江、穆棱河以及兴凯湖等水补给充足的沿岸区域.适宜性良好的区域主要分布在饶河县,到2010和2015年,虎林县和抚远县也变为适宜性良好区域的重要分布区,主要由于该区水田面积大幅增加.适宜性一般区域分布较零散,其面积先增加后减少.适宜性差的区域面积在1990—2000年间增加6.7%,2000—2015年间减少3.1%.土地覆被变化是水禽栖息地适宜性等级变化的最重要影响因素;人口和经济增长以及气候的变干、变暖使水禽栖息地适宜性降低;保护区的建立使水禽栖息地得到有效保护.     

基于生态-经济比较优势视角的建设用地空间优化配置研究——以扬州市为例

随着生态空间日益萎缩,生态环境不断恶化,经济发展与生态保护矛盾日趋加剧。为了协调经济发展空间与生态保护空间的矛盾问题,促进生态、经济协调持续发展,综合运用景观生态学方法与经济学原理,从"生态-经济"比较优势视角对建设用地指标进行空间优化配置。运用景观安全格局方法构建区域综合生态安全格局,据此计算空间单元的生态重要性指数,采用最小累积阻力模型对区域建设用地开发经济适宜性进行评价,计算空间单元的建设用地开发经济适宜性指数;其次运用生态重要性指数与经济适宜性指数计算建设用地开发比较优势度,以比较优势约束性满足程度为目标函数构建建设用地空间优化配置模型,并用蒙特卡洛(Monte Carol)法对之求解。以江苏省扬州市为例作了实证研究。结果显示,未来扬州市建设用地向南部经济发达地区集中,北部地区主要承担生态保护功能;建设用地配置结果呈现出经济发展与生态保护的错位格局,有助于扬州市生态保护与经济发展。从"生态-经济"比较优势视角提出了建设用地空间优化配置方法,为土地利用规划及生态保护规划编制提供了科学依据。     

基于PSR 模型的河南省水生态安全综合评价研究

基于河南省水环境现状, 从“压力-状态-响应”(PSR)指标体系框架出发, 建立河南省水环境生态安全评价指标体系, 利用变异系数法确定各指标权重, 运用综合指数评价法对河南省水环境生态安全水平予以评价。结果表明:2007-2013 年河南省水环境生态安全综合指数值整体呈上升趋势, 其生态安全水平已处于预警状态; 压力状态曲线呈逐年下降趋势, 状态曲线与综合指数曲线变化趋势基本一致, 表明河南省水环境综合评价值在很大程度上是状态系统的贡献, 此外响应系统的贡献率逐年增幅显著, 有效促进了河南省水环境生态安全的逐渐改善。     

城市近郊型乡镇农业生态经济系统生态潜力及产业空间布局——以延安市河庄坪镇为例

生态脆弱区抗干扰能力弱，明确其系统的生态潜力，在潜力的基础上开发利用资源对区域的可持续发展具有重要的意义。研究基于生态系统服务价值视角，尝试测度农业生态经济系统生态潜力，同时以三产融合典型代表乡镇为研究对象，对其农业生态经济系统生态潜力和主导产业进行时空分析，最后对生态系统加以功能分区。结果显示：（1）气体调节、气候调节、水源涵养和土壤形成与保护是生态系统提供的主要服务功能。随退耕还林工程实施，林地面积的增加使生态系统服务价值量随之上升，在空间上呈现"中间低两翼高"的格局。（2）农业生态经济系统生态潜力受服务价值与资源开发过程影响，与生态系统服务价值在空间上相对一致，呈现出"两翼高中间低，北高南低"的分布格局，其中，退耕还林工程的推进直接导致林地的生态潜力增加，耕地的生态潜力减少。（3）研究区在退耕过程中形成以大棚种植为主的设施农业、以农林特色销售为主的商品型种植农业，以山地苹果为主的林果产业以及以酒店商贸为主的服务行业。并随城镇一体化的推进，主导产业呈现出由传统粗放型大田种植向集约型设施农业方向发展。在空间上形成了以镇区分布为主的第三产业、依托拐沟坡地的林果产业和沿川道展开的设施农业的产业布局。（4）基于生态系统服务价值，研究区农业生态经济系统功能区划分为城镇生态环境维护区（第Ⅰ类型区）、农业生态保护区（第Ⅱ类型区）、生态安全屏障区（第Ⅲ类区）3个区域。     

1995-2015年武汉城市湖泊景观生态安全格局演化

以1995、2005、2015年3期遥感影像为数据源,通过景观指数和GIS空间分析方法研究武汉中心城区湖泊系统景观格局的演化特征,进而构建湖泊景观生态安全评价模型,以揭示武汉中心城区湖泊系统景观生态安全格局演变规律,并探讨其驱动因素。研究结果表明：（1）近20年来武汉中心城区湖泊系统总面积呈现不断下降趋势,1995-2005年以年均2.4%的速率缩减了28.95 km²,2005-2015年以年均1.8%的速率缩减了17.47 km²。（2）研究期湖泊系统的斑块密度、分维数呈现不断降低趋势,这表明武汉中心城区湖泊破碎化程度在不断降低,湖泊形状趋于简化;而连接度呈现不断增加趋势,湖泊分布集聚趋势不断加强,空间结构趋于紧凑集中。（3）武汉城市湖泊系统总体景观生态安全呈现不断恶化趋势,而次级湖泊水系景观生态安全呈现不同的演化特征。城市建成区的迅速扩展、湖泊由农业生产对象转化为可利用发展用地、房地产事业成为支柱产业是湖泊系统遭填占最主要的动力,直接危害了其景观生态安全,而市民环境意识提高引发湖泊保护倡议和实践,对湖泊景观生态安全产生着积极影响。对湖泊景观生态安全格局演变特征及驱动因素的研究,可以为武汉制定行之有效的生态管理政策和环境治理措施提供重要支撑。     

中国县域碳汇时空格局及影响因素

以全国1300个县级行政单位作为研究对象，利用全国县域尺度土地利用数据和社会经济数据，核算了1990-2015年间中国县域碳汇总量，并结合标准差椭圆、空间自相关、冷热点分析和地理加权回归分析方法，探析中国县域碳汇的时空分异特征及影响因素，旨在为优化国土空间开发格局，实施差异化减排路径及推动生态文明建设提供参考。研究表明：（1）时空变化上，1990-2015年中国碳汇总量呈波动下降趋势，由13307.79×10⁴t下降至13198.27×10⁴t；林地为主要碳汇类型，其余碳汇类型比例结构基本不变；在空间分布上，中国县域碳汇呈现"西部 > 东北 > 南部 > 中部"的"西高东低"格局。（2）在空间分异和聚集上，碳汇空间分布中心向西南移动，分布范围呈收缩态势，西南地区对整体碳汇空间格局影响作用加强；1990-2015年中国县域碳汇总量的冷热点集聚程度呈现波动稳定特征，空间集聚程度呈现高值与低值聚集，高-低区域零星分布的特征。（3）从影响因素分析来看，2015年经济发展、产业结构、土地利用程度对碳汇产生影响并存在空间异质性。建议通过合理规划县域非建设用地的土地利用方式、差异性制定各区域内策略、控制建设用地规模等方式达到县域的低碳发展。     

基于变权-物元分析模型的老工业基地区域生态安全动态预警研究——以吉林省为例

生态安全预警是生态安全研究的重要内容,对维护区域生态安全具有重要的意义。老工业基地区域生态安全预警研究具有一定的特殊性和典型性,选取吉林省为研究区,基于压力、状态、响应(P-S-R)和生态、环境、经济、社会(E-E-E-S)框架模型构建生态安全预警指标体系,运用变权-物元分析模型对1991—2011年吉林省区域生态安全的预警等级进行了测度,并结合灰色系统GM(1,1)预测模型对2015年、2020年区域生态安全态势进行预测预警。研究结果表明:(1)1991—2011年,吉林省生态安全总体水平呈逐步上升的趋势,区域生态安全预警等级从"巨警"上升为"轻警",指示灯由"红色预警"逐步变为"蓝色预警";(2)2015年,吉林省的生态安全预警等级属于"轻警",指示灯为"蓝色预警",到2020年,吉林省的生态安全预警等级属于"无警",指示灯为"绿色预警",但具有向"轻警"变化的态势;(3)制约吉林省生态安全水平提升的主要因素包括单位面积耕地化肥负荷量、城镇化、人均水资源量、GDP增长率、第一产业占GDP比重等,这些因素应是今后生态环境保护调控的重点,本研究为吉林省经济-生态-社会可持续发展提供参考。     

中国农业系统近40年温室气体排放核算

基于排放因子法构建了包含种植业和牲畜养殖业的农业系统温室气体排放核算体系,系统核算了1980-2020年我国全国尺度上的农业系统温室气体排放总量和变化趋势,并在区县级尺度下对1980、2000、2011年的中国农业系统的温室气体排放量进行核算,对比不同阶段农业系统温室气体排放变化的时空异质性规律。研究发现：1980-2020年我国农业系统温室气体排放量呈波动增长趋势,增长了近46%。CH₄是农业系统排放贡献最大的温室气体,占总排放量的47.33%。我国农业系统温室气体排放与不同地区农业生产方式有关,CH₄排放量高的地区主要位于我国主要水稻产区以及旱地作物产区。CO₂排放量高的地区主要位于东北、西北等地区以及华东地区。N₂O排放量较高地区主要位于西北的主要畜牧养殖地区,以及我国农业经济发展水平高的中南部地区。研究有助于揭示我国农业温室气体排放的动态特征,现状规律,以及空间差异性特征,从农业减排角度为实现双碳目标提供科学参考。     

我国农业系统可持续发展协调度分析

农业系统可持续协调发展作为一种发展理念已成为共识,受到我国政府高度重视。采用2004—2015年我国农业生产面板数据,运用熵值法、探索性空间数据分析等方法分析我国农业系统可持续发展协调度情况。研究结果表明:(1)2004—2015年我国各省(市、自治区)农业各子系统可持续能力之间协调度处于不稳定变动态势,且其协调度值较低。(2)从协调度分布来看,协调度处于"比较协调"和"协调"的省(市、自治区)主要分布在我国的东部和东南部地区,其余大多数省(市、自治区)农业各子系统可持续能力之间的协调度处于"不协调"。(3)从空间自相关性来看,各省(市、自治区)农业各子系统可持续能力之间协调度呈正的全局自相关性且相关程度逐渐增大,而大多数省(市、自治区)农业各子系统可持续能力之间协调度的局部自相关性不显著,且其农业各子系统可持续能力之间协调度空间极化不明显,处于较低水平均衡状态。从时空维度、空间相关性方面分析农业各子系统可持续能力之间协调度,为提高我国农业系统可持续发展水平提供政策建议。     

陇中黄土高原农村经济结构变化对退耕还林(草)的响应

黄土高原退耕还林(草)政策的实施有效改善了区域生态环境,然而其对社会经济的影响存在明显的地域差异。以甘肃省定西市安定区近郊凤翔镇和远郊鲁家沟镇为例,采用重复测定方差分析、主响应曲线分析和相关分析等解析2镇1998—2015年农村经济结构及收入变化对退耕还林(草)政策实施的响应特征,确定退耕还林(草)政策实施驱动的土地利用变化与农村收入之间的关系。研究发现:(1)退耕还林(草)政策实施解释了2镇70.4%的农业收入结构差异,因政策的实施推动了小麦等夏粮作物向玉米、马铃薯等秋粮作物的转变,进而影响了农业收入结构的变化;(2)退耕还林(草)政策对农村经济结构的影响弱于对农业收入结构的影响,它解释了2镇54.1%的农村经济结构差异,近郊凤翔镇农村经济收入明显高于远郊鲁家沟镇。退耕还林(草)政策对以农业为主导产业的远郊鲁家沟镇农村经济结构影响相对较小,近郊凤翔镇工业、服务业和商业餐饮业等非农产业占农村经济总收入的比重随着农业收入比重的降低而有所提高;(3)退耕还林(草)政策的实施带动了该区种植业结构的转变,这在一定程度上缩小了凤翔镇和鲁家沟镇农村集体收入之间的差距;(4)退耕还林(草)政策实施过程中,梯田面积、玉米和马铃薯等秋粮作物种植面积的增加提高了农村居民纯收入。研究认为退耕还林(草)在实施过程中需考虑地域差异,依据区域经济发展状况进行相应调整,在推动种植业发展的同时还应兼顾林业、牧业和一些非农产业发展,最终实现生态环境恢复与农村经济的协调发展。     

基于改进TOPSIS-灰色GM(1,1)模型的张家界市旅游生态安全动态预警

旅游生态安全预警对旅游地协调生态环境保护与旅游产业快速发展的矛盾,推进区域社会经济持续健康发展具有重要意义.本文基于DPSIR模型,从驱动力、压力、状态、影响、响应5方面构建了张家界市旅游生态安全预警指标体系,在此基础上,运用改进TOPSIS法对2001—2014年张家界市的旅游生态安全警情格局进行评估,并运用灰色GM(1,1)模型对其2015—2020年的警情演变趋势进行预测.结果表明: 总体上,2001—2014年,张家界市的旅游生态安全贴近度呈略微上升态势,警度处于“中警”.就各子系统而言,2001—2014年,张家界市旅游生态安全驱动力系统和压力系统的警度呈上升趋势,由“轻警”演变成“重警”;状态系统和影响系统的警度变化不大,一直处于“中警”;响应系统的警度呈下降趋势,由“巨警”转变成“无警”.按现状发展态势,2015—2020年,张家界市的旅游生态安全贴近度将进一步上升,警度将由“中警”向“轻警”转变,但协调旅游发展与生态建设和环境保护关系的任务依然较艰巨.     

供给侧改革视角下中国省域生态产品有效供给及影响因素

生态产品是事关人类生命健康、经济社会建设的世界性稀缺产品.为了更好地分析生态产品有效供给问题,找出生态产品有效供给的规模效率的变化规律,本研究利用2004—2015年的面板数据,以劳动力、资本、能源和技术为投入变量,以GDP、生态产品和非期望产出为产出变量,运用SBM-Undesirable模型和泰尔指数对我国生态产品有效供给的规模效率及其差异性进行测算,并运用Tobit模型进一步对其影响因素进行实证分析.结果表明: 我国生态产品供给规模效率在2011年以前呈下降趋势,以2012年为转折点呈现上升态势,但地区间仍有较大差距.全国及各地区生态产品有效供给效率的影响因素不尽相同,从全国总体情况看,经济发展水平、产业结构(第二产业增加值占GDP比重)、科技投入、金融资本投入(环境污染治理占GDP比重)、能源消费结构(煤炭消费量占能源消费总量比重)和农业面源污染的影响为负,而生态空间及外贸结构的影响为正;在东部地区,产业结构、科技投入、金融资本、能源消费结构和农业面源污染的影响为负,经济发展水平、生态空间及外贸结构的影响为正;在东北地区,经济发展水平、农业面源污染和外贸结构的影响为负,产业结构、科技投入、能源消费结构和生态空间的影响为正;在中部地区,经济发展水平、科技投入、金融资本、能源消费结构、农业面源污染和生态空间的影响为负,产业结构和外贸结构的影响为正;在西部地区,产业结构、能源消费结构和农业面源污染的影响为负,经济发展水平、科技投入、金融资本、生态空间及外贸结构的影响为正.在进一步对以上实证测算进行理论分析的基础上,从供给侧改革的角度提出有效提高生态产品供给规模效率的政策启示,包括注重转变经济发展方式、优化生态系统结构和遵循生态法则进行整体治理等.     

种养系统可持续发展指数的空间格局及其演变趋势——以山东省为例

农业生产专业化、集约化背景下,种养分离越发严重,畜禽粪便排放与环境之间矛盾日益突出,农业可持续发展面临重大挑战.本研究以山东省为例,基于1999-2015年种养系统投入产出资料,利用能值的方法,定量分析了山东省域及市域种养系统可持续发展指数的空间格局及其演变趋势.结果表明:1999-2015年,山东省种养系统的可持续性减弱.可持续发展指数呈极显著下降趋势,2015年较1999年下降了22.0%;净能值产出率和环境负载率均呈极显著增加趋势,单位经济投入获得的利益显著增加,同时种养生产对环境的压力明显加大,这主要与电力、复合肥、农业机械等工业资源投入量的增加密切相关.山东省各地市种养系统可持续发展水平具有差异.大部分地市种养系统的可持续发展指数较高,而沿海地区(威海、烟台)和工业城市(淄博)则较低.同时,各地市种养系统的可持续发展态势也具有差异.鲁中和鲁南地区可持续性逐年减弱,鲁北地区可持续性逐年增强.以2015年为本底数据进行情景分析,发现种养高度结合(100%粪便有机肥替代化肥)情景的可持续发展指数可达8.4,是种养结合现状(30%粪便有机肥替代化肥)的2.6倍.