黄河三角洲高效生态经济区生态系统服务权衡协同关系

探究黄河三角洲高效生态经济区生态系统服务权衡与协同关系,对区域生态系统服务功能分区和高质量发展具有重要的现实意义。本研究利用InVEST模型、空间自相关和权衡协同度(ESTD)模型,分析黄河三角洲高效生态经济区2000—2020年生境质量、碳储存、土壤保持、水源涵养和水质净化5种生态系统服务在乡镇尺度上的时空变化特征及其权衡协同关系。结果表明：2000—2020年,研究区生境质量、碳储量和氮、磷输出量整体递减,土壤保持量和水源涵养量波动递增;生境质量呈现北高南低的分布格局,碳储存、氮、磷输出、土壤保持和水源涵养大致呈北低南高的格局。研究期间,研究区5种生态系统服务间在时间截面和时段上均以协同关系为主,但仍存在差异,在时间截面上协同关系主要存在于碳储存与其他服务之间,在时段上协同关系主要存在于生境质量与其他服务之间。研究结果可为黄河三角洲高效生态经济区生态系统服务的提升及生态系统功能分区提供理论指导与实践借鉴,并为国土空间格局优化提供基础支撑。     

黄土高原典型退耕区生态系统服务权衡与协同关系研究——以延安市为例

在退耕还林的背景下,从多尺度分析黄土高原生态系统服务的时空变化及权衡协同关系对促进该区经济发展和生态效益多赢具有重要意义。以黄土高原退耕还林典型区延安市为例,运用InVEST模型评估1988-2018年农作物生产、碳储量、生境质量、土壤保持、产水量5种关键生态系统服务物质量,采用偏相关分析法探讨生态系统服务间的权衡与协同关系及其时空变化特征。研究结果表明：（1）受土地利用变化影响,延安市农作物生产水平、土壤保持年际变化趋势波动较大,农作物生产水平呈波动增加的趋势,土壤保持退耕还林前呈减少,退耕还林后呈波动增加的趋势,碳储量和生境质量呈逐渐增加的趋势,产水量呈逐渐减弱的趋势。（2）生态系统服务与土地利用格局联系紧密。碳储量、生境质量的高值区域以及产水量的低值区域随林地分布格局变化而变化。（3）协同关系是延安市生态系统服务间相关关系的主体,主要发生在碳储量、生境质量、土壤保持、农作物生产之间,权衡关系主要存在于产水量和其他生态系统服务之间。（4）市、县域生态系统服务间关系的差异主要发生在生境质量和土壤保持之间。     

南四湖流域关键生态系统服务的时空权衡关系

从多尺度科学地研究流域各生态系统服务间的权衡与协同关系,对流域实现经济发展与生态保护双重目标具有重要研究意义。运用InVEST模型模拟评估产水量、水质净化、土壤保持、碳存储、生物多样性维持和授粉等6类区域关键生态系统服务,并通过热点分析、相关分析和主成分分析,定量分析自改革开放以来（1980-2018年）这些生态系统服务的时空变化及其在不同尺度上的权衡关系。结果表明：（1）南四湖流域产水服务呈上升趋势,而水质净化、土壤保持、碳储存、生境质量、授粉等服务则总体呈下降趋势。（2）研究区东部林地和草地广布,是碳存储、生境质量和授粉服务的高值区,而建设用地分布的地区是产水量和土壤保持服务的高值区。（3）产水量与水质净化服务之间呈权衡关系,水质净化与生境质量、碳存储与授粉之间则存在较强的协同关系。（4）在子流域尺度上主要是存在着供给服务（产水量）与调节服务（水质净化和土壤保持）、供给服务（产水量）与支持服务（生境质量）之间的空间权衡关系;在县域尺度上,供给服务、调节服务与支持服务（产水量、水质净化、碳存储、生境质量和授粉）之间主要是空间协同关系,两个尺度之间具有较大差异。研究结果可为管理者在合适的空间尺度上制定相应的政策提供科学支撑,促进南四湖流域生态环境和社会经济协同发展。     

大宁河流域生态系统服务权衡关系梯度效应

深入理解生态系统服务权衡关系是推动生态系统有效管理和实现区域高质量发展的基础。然而,目前涉及生态系统服务权衡关系梯度效应的研究仍相对较少。选取三峡库区大宁河流域为研究区域,对典型平水年份2018年粮食生产、碳固持、土壤保持和产水量4种生态系统服务进行了量化评估和空间制图,探讨了海拔、降雨量、植被覆盖度3种环境因子梯度下生态系统服务的空间分异规律,采用均方根偏差进一步分析了生态系统服务权衡关系梯度效应。研究结果表明：粮食生产在空间上呈现中部、南部高,北部低的趋势,碳固持由西北向东南递减,土壤保持和产水量高值区均分布在西部。海拔、降雨量和植被覆盖度显著影响了生态系统服务,碳固持、土壤保持总体上随着环境梯度的增加而增加,而粮食生产随着3种环境因子梯度的增加呈下降的趋势。从生态系统服务权衡方面来看,碳固持与产水量为弱权衡关系,土壤保持与粮食生产为中度权衡,生态系统服务之间多为强权衡关系。由于多种环境梯度的相互作用,生态系统服务权衡关系随海拔梯度变化无明显规律,植被覆盖度增加显著加强了生态系统服务权衡关系,而降雨量增加显著减弱了生态系统服务权衡关系。因此,今后以流域生态系统服务协调和功能整体提升为目标的生态修复工程中,应更加注重统筹不同环境梯度区域内生态系统服务权衡关系。     

中国快速城市化干燥地区的生态系统服务权衡关系多尺度分析——以呼包鄂榆地区为例

理解生态系统服务间的权衡关系对提高干燥地区人类福祉和实现区域可持续性具有重要的意义。但是,目前人们仍然缺乏对快速城市化干燥地区中生态系统服务权衡关系的多尺度理解。以呼包鄂榆地区这一正在经历快速城市化的干燥地区为例,在城市群、区域和城市3个尺度上来探索生态系统服务间的权衡关系。对呼包鄂榆地区2010年的粮食生产、肉类生产、产水量、土壤保持和碳固持5种关键生态系统服务进行测量。利用相关分析法在3个尺度上对5种服务间的权衡关系进行分析。2010年呼包鄂榆地区产水量和土壤保持以及产水量和碳固持服务表现出显著的权衡关系。其中,产水量和碳固持服务在城市群、农业区和鄂尔多斯市呈现明显的权衡关系,产水量和土壤保持服务在城市群和农业区呈现明显的权衡关系。快速城市化干燥地区中的生态系统服务权衡关系具有明显的尺度效应,同一对生态系统服务在不同尺度上的权衡关系存在较大差异。这种差异主要是由区域人类活动和自然条件的空间异质性所致。研究结果有助于对呼包鄂榆地区生态系统服务权衡关系的认识,为该地区土地系统设计和可持续发展提供了必要的科学依据。     

滇池流域水文生态系统服务权衡与协同时空异质性及其归因分析

城市化扩张影响流域水文生态系统服务供给,导致水文生态环境压力增加。明晰滇池流域水文生态系统服务权衡与协同关系及归因,对于滇池流域社会经济发展与生态可持续管理的双赢具有重要意义。基于InVEST模型估算了滇池流域2000、2010、2018年碳存储、氮输出、土壤保持和产水量4种水文生态系统服务,分析了水文生态系统服务时空格局变化,利用双变量Moran′s I指数分析了滇池流域水文生态系统服务间时空异质性,通过冗余分析和条件效应分析等方法识别水文生态系统服务归因。结果表明：(1)水文生态系统服务的高值分布与土地利用覆盖有密切关系。(2)水文生态系统服务空间变化主要体现在滇池周边人类活动干扰强烈的区域,表现为氮输出、土壤保持和产水量增加,碳存储减小。(3)水文生态系统服务权衡与协同关系变化明显。碳存储和土壤保持为协同;氮输出和产水量先呈协同后呈权衡;产水量、氮输出均与碳存储和土壤保持为权衡,前者空间变化不明显,后者协同范围扩大,权衡范围缩小。(4)水文生态系统服务主要受多个因素共同影响,主要包括：高程、耕地比例、林地比例、建设用地比例。研究结果可为研究区水文生态系统环境治理和土地利用优化配置...     

基于生态系统服务权衡的优先保护区选取研究——以南方丘陵山地带为例

南方丘陵山地带作为"两屏三带"国家生态屏障区之一，发挥着保障华南和西南地区生态安全的作用，优先保护区的选取对区域的生态系统服务提升、生态安全和可持续发展具有重要的意义。以南方丘陵山地带为研究对象，采用InVEST模型和森林游憩模型评估产水量、土壤保持、碳储存、生境质量和生态休闲的空间分布，基于有序加权平均算子（OWA）建立了11种情景，通过对比各情景的保护效率筛选出了优先保护区。结果表明：（1）南方丘陵山地带产水量和土壤保持分布不均，平均产水量1108mm，高值区均集中在中部地区，而东部区域和西部区域的生态休闲较高；（2）研究区固碳能力强，平均碳储量为45.58t/hm²，生境质量良好，平均生境质量为0.75，碳储存和生境质量与土地利用/覆被类型关系密切；（3）情景2的保护效率最优，对产水量、土壤保持、碳储存、生境质量和生态休闲的保护效率分别为1.17，1.89，1.32，1.48，1.18。研究结果可为研究区生态系统服务的提升提供决策支持。     

秦岭生态系统服务协同与权衡的时空异质性

秦岭是我国重要的生态功能区,探究秦岭生态系统服务权衡与协同关系对区域的生态保护和高质量发展意义重大。本研究基于遥感、气象、土壤质地等数据,在分析2000、2010和2020年土地利用时空变化的基础上,采用CASA、InVEST和RUSLE模型计算秦岭地区植被固碳量、产水量、生境质量和土壤保持量,并探讨空间异质性特征及生态系统服务关系。结果表明: 2000—2020年,研究区耕地、草地面积持续减少,林地和建设用地面积增长,水域和未利用地变化不明显。研究区固碳量东高西低,呈增长趋势;产水量持续减少,呈现南坡高、北坡低的格局;生境质量稳定增长;土壤保持量先增后减。不同地类的生态系统服务各有差异,同一地类的各个生态系统服务也各不相同。全域尺度下固碳与土壤保持、生境质量与土壤保持、固碳与生境质量皆呈协同关系,固碳与产水、土壤保持与产水、生境质量与产水为权衡关系,但不同条件下服务量间的关系存在明显差异。研究结果有助于降低权衡风险,使总体效益最大化,可为促进生态高质量发展、实现自然资源管理与人类福祉的双赢提供参考依据。     

基于多尺度秦巴山区生态系统服务权衡协同关系的表达

秦巴山区是中国南北过渡带,具有明显的多维地带性结构和复杂的地貌特征,从大、中、小不同的空间尺度上研究生态系统服务权衡协同关系,为区域间的山水林田湖草生态保护提供基础性参考。从全域、综合分区、典型样区、样点间隔4个尺度分析秦巴山区2005—2015年植被净初级生产力(NPP)、粮食生产、土壤保持、生境质量、水资源供给5种生态系统服务之间权衡协同关系。结果表明:①从全域尺度:水资源供给分别与生境质量、土壤保持、NPP、粮食生产呈现协同关系;NPP分别与生境质量、土壤保持同呈现协同关系;生境质量与土壤保持是协同的关系;粮食生产分别与生境质量、土壤保持呈现权衡的关系。②从综合分区尺度:依据地貌和气候类型划分的综合分区,会出现特定的区域内生态系统服务关系的异同,加上影响生态系统服务的机制不同,促使生态系统服务间的关系在综合分区的体现也有所差别。③典型样区尺度:山地的南北坡在气候类型、植被类型、地貌特征等方面存在明显的差异性,随着时间的推移生态系统服务之间会发生相互转变。④样点间隔尺度:成对服务间的权衡/协同关系与样点间隔尺度存在一定的关联性,不同的间隔尺度权衡/协同关系的发展趋势存在差异性。该研究可以协助决策者在适当的空间尺度制定相应生态管理策略,促进生态-经济协同发展。     

陕北黄土高原景观格局对生态系统服务权衡关系的影响

揭示景观格局对生态系统服务权衡关系的作用机理对于开展多种ES的可持续管理具有重要意义。在评估陕北黄土高原2000—2020年粮食供给、水源涵养、土壤保持、生境质量和景观美学等5种生态系统服务的基础上,利用相关系数和均方根误差测度其权衡强度,结合R指数模型、变异划分法定量表达景观组分和景观配置对生态系统服务权衡关系的影响差异。结果表明：1)2000—2020年陕北黄土高原5种生态系统服务整体呈上升趋势,空间上粮食供给北高南低,水源涵养、土壤保持、生境质量和景观美学则为南高北低的分布格局;2)粮食供给与其他4种服务为权衡关系,而水源涵养、土壤保持、生境质量和景观美学之间呈协同关系;3)斑块丰富度对各组权衡关系具有正向影响,而斑块面积和景观百分比则具有负向影响;4)景观组分对各组权衡关系的单独贡献度高于景观配置,多数情况下,权衡关系更多由组分和配置的联合作用主导。研究从景观管理角度为改善陕北黄土高原生态系统服务间的权衡关系、提升ES的潜在供给能力提供了科学依据。     

子午岭区生态系统服务权衡与协同变化及其影响因素

生态系统服务簇是多种生态系统服务的组合,可指征生态系统服务间权衡或协同关系。辨析服务簇内部的权衡与协同关系,有助于实现基于服务簇的生态分区与管理政策制定。基于InVEST模型,综合运用相关分析、聚类分析、主成分分析和冗余分析等,本研究揭示了子午岭区生态系统服务的空间群聚特征,探究了服务簇内部权衡与协同关系及生态系统服务的影响因素。结果表明: 2000—2017年,子午岭区生态系统服务整体得到改善。除土壤保持高值区在研究区北部较为集中外,其余生态系统服务高值区均在中部和南部较为集中。根据聚类结果将研究区在像元尺度上划分为4个生态系统服务簇,其中,核心保护簇内,产水服务与土壤保持存在较强权衡;粮食生产簇内,生境质量与产水服务、粮食生产存在较强权衡;生态过渡簇和生态脆弱簇内,服务间的协同关系广泛。生态系统服务同时受到自然本底条件和人类活动的影响,植被覆盖、降水和高程对研究区生态系统服务以积极影响为主,而温度、建设用地占比、人口密度等以消极影响为主。     

城市化与生态系统服务的空间交互关系研究——以长株潭城市群为例

探索生态系统服务与城市化的空间交互关系为城市发展和保护提供新的视角。以长株潭城市群为例,以其2015年的土地利用、经济社会等数据为基础,利用格网分析和空间自相关模型解析城市化与生态系统服务的空间交互规律。结果表明:(1)长株潭城市群城市化水平空间分布区域差异性显著,城市化水平由高到低依次为城市化地区、正在城市化地区和生态绿心保护区,表现出从城市化地区向周边地区依次递减的态势。(2)受城市化水平和土地利用变化等影响,生态系统服务空间分布具有明显的差异性,并呈现出"中心-外围"逐步升高的趋势。(3)长株潭城市群城市化与各项生态系统服务之间呈显著负相关,由于驱动机制的不同,城市化地区、正在城市化地区和生态绿心保护区的城市化对生态系统服务的影响存在显著差异。(4)长株潭城市群的各项生态系统服务与城市化之间存在四种空间相关类型。高高集聚区和低低集聚区零散分布在整个研究区域,面积较小;高低集聚区分布在湘潭县西南部、株洲县东南部,面积较大;低高集聚区分布集中在城市化地区,呈现块状分布,面积较大。     

快速城市化生态系统服务格局分析与空间管控——以厦门市为例

快速城市化对生态系统服务产生消极影响，科学识别二者的空间规律并进行管控对国土空间资源可持续开发意义重大。基于土地利用、遥感影像、POI等多源数据，借助ArcGIS、InVEST平台分析厦门市六类生态系统服务的空间分布及权衡协同格局，并引入人工神经网络算法与土地利用程度指数模型定量识别城市化空间格局，通过决策矩阵提出空间管控策略。结果表明：①厦门生态系统服务分布在"陆-海"方向上分异显著，岛外地区综合生态系统服务指数总量占全域的91.78%，其中同安区的总量和均值都处于首位；②生态系统服务之间以协同作用为主，在全局关系中生境质量与雨涝调蓄以及与温度调节的协同水平最高，在局部空间关系中水源涵养与雨涝调蓄的协同水平最高，休闲游憩服务与其他服务的权衡效应最强；③翔安区开发潜力最大，同安区过度城市化问题突出，市域未开发单元中生境质量、水源涵养和雨涝调蓄服务较强，建成单元仅在休闲游憩服务具有优势；④叠加空间细分单元的城市化与生态系统服务特征，得到厦门市生态管控分区和决策矩阵，提出四类单元八种情景的具体管控策略。     

基于生态系统服务的宁夏回族自治区自然-社会经济协调性分析

近年来,受全球气候变化和人类活动的双重影响,我国干旱半干旱区生态环境保护和社会经济发展面临着严峻挑战。生态系统服务是链接自然生态系统和人类福祉的桥梁,其持续供给对区域可持续发展具有重要意义。在具有生态重要性、脆弱性的地区,如何协调生态保护与社会经济发展的关系以达到可持续发展是重要的科学问题。以宁夏回族自治区为例,基于InVEST模型及相关地理图件,分析研究区生态系统服务变化及其权衡关系,通过耦合协调度模型分析生态系统服务与社会各项经济指标的耦合关系及空间差异,辨析区域耦合协调度的主要影响因素。结果表明：（1）研究区生态系统服务空间分异明显,土壤保持、水质净化、水资源供给和碳储存均呈现南高北低的特征。（2）县域尺度上各项生态系统服务之间显著正相关,呈现为协同关系;其中水质净化与水资源供给的相关性最高（相关系数是0.73）,其次是碳储量和水资源供给（相关系数是0.60）;而土壤保持与水质净化、水资源供给的相关性较低（相关系数均在0.50以下）。（3）宁夏回族自治区耦合协调度与城镇人口比重、城镇居民人均可支配收入负相关,与农村常住居民人均可支配收入正相关;产业分配模式与耦合协调度显著相关。（4）宁夏回族自治区约50%的县（区）处于勉强协调发展状态,其中73%的县（区）更注重经济发展;约41%的县（区）属于协调发展类,其中2/3的县区环境保护相对超前;中部地区的发展优势较强。研究结果对于准确把握经济社会发展过程中生态系统服务特征具有重要理论和实践意义,为促进宁夏生态保护与经济发展的统一协调提供科学辅助依据。     

区域一体化视角下长三角地区生态用地时空演变影响因素

区域一体化在推动要素流动、整合地区资源、提高区域竞争力的同时,也为生态用地带来一定风险与压力,造成生态用地规模缩减、景观格局破碎、生态功能降低等问题。探讨区域一体化背景下生态用地的时空演变特征与影响因素对于平衡区域发展与生态保护的关系、制定合理生态用地保护策略、实现区域协调可持续发展等具有重要意义。以长三角地区为例,在分析生态用地时空演变的基础上分别构建两个地理加权回归（GWR）模型,一个模型为考虑区域一体化因素的实验组,一个模型则是不考虑区域一体化因素的对照组,通过结果对比分析区域一体化对长三角地区生态用地演变的影响,结果表明：（1）2000-2020年间,长三角地区生态用地规模缩减显著,减少166.56×10⁴ hm²,生态用地急剧减少时间集中在2015-2020年间;（2）长三角地区生态用地减少面积的时空演变具有以上海、江苏南部和浙江北部为中心向周边地区扩散的趋势,扩散方向覆盖长三角全部地区,生态用地集聚性快速减少区域则以安徽中北部、江苏北部地区为主;（3） GWR回归对照试验分析结果表明区域一体化对生态用地演变具有显著促进作用,其作用强度和作用范围伴随时间增加;（4）根据生态用地演变影响机制,长三角生态用地演变类型可分为内生驱动型、外生驱动型以及发展监测型,可为平衡国土空间开发与生态用地保护提供决策支持。     

基于DPSIR模型的长三角城市群生态安全评价研究

长三角城市群作为中国经济发展最具活力的区域之一,城镇化水平高,同时面临巨大的环境压力。以长三角城市群41个城市为研究对象建立DPSIR生态安全指标评价模型,利用熵权法与均方差决策法确定DPSIR权重,采用主成分分析法和综合指数法得到长三角城市群生态安全综合指数YDESI与生态安全等级。结果表明,长三角城市群41个城市依据生态安全水平可划分为4类。总体水平属于临界安全偏较不安全水平,上海市与江苏省生态安全情况较好,浙江省生态安全水平一般,安徽省生态安全情况稍差。DPSIR模型显示,驱动力与响应指标对长三角城市群生态安全建设有很大贡献。长三角城市群生态安全主要影响因素为农业发展、环境空气质量、环保投资、城市绿化建设。最终提出相应的策略:即提高区域人口素质、强调新型城镇化、加强城市群基础设施建建设,为长三角城市群生态安全建设提供思路与依据。     

柑橘园扩张对生态系统服务和景观格局的影响——以江西省信丰县为例

柑橘园扩张是我国柑橘主产区较典型的土地利用/覆盖变化过程之一。但柑橘园扩张对区域生态系统服务和景观格局的影响尚不清楚,探明柑橘园种植的生态效应对柑橘主产区实现生态-经济协同发展具有重要的理论和现实意义。以江西省信丰县为研究区,基于第二次、第三次全国国土调查数据,借助InVEST模型、景观格局指数、图谱分析等方法,结合地形梯度和相关性分析,对柑橘园扩张所致的生态系统服务和景观格局变化及其权衡与关联进行分析。结果表明：①2010-2019年,柑橘园扩张方向以占用林地和耕地为主,柑橘园面积从9961.29 hm²增至12633.84 hm²,其中新增柑橘园7504.72 hm²,退出柑橘园4832.16 hm²,区域内柑橘园转换频繁。②柑橘园扩张导致信丰县景观格局总体呈现破碎化趋势,景观连通性降低,并对土壤保持服务和碳储量服务产生负向影响,对柑橘供给服务产生正向影响,对产水服务几乎无影响。③不同类型柑橘园在不同地形梯度上的扩张对生态系统服务的影响方向和程度存在差异。④在权衡与协同分析中,供给服务与支持服务呈权衡关系,碳储量-产水、碳储量-土壤保持等支持服务之间呈协同关系;景观破碎化与生态系统服务显著关联。研究结果可为经济作物大规模种植区的园地合理开发利用和生境保护提供科学依据。     

长三角地区休闲游憩服务供需关系研究

依托生态空间的休闲游憩服务是提升人类福祉至关重要的贡献者。但由于经济发展、城市扩张、人口增长等因素,区域休闲游憩服务供需矛盾突出,且供给与需求在空间上错位、不匹配等问题严重。研究休闲游憩服务供需关系对合理规划、管理区域生态用地,促进生态文明建设等方面具有重要的理论和现实价值。以长三角为例,基于土地利用与人口数据,提出不同可达半径下区域休闲游憩服务供需关系的评价方法,研究41个地级市不同可达半径下供需关系的空间分异规律。研究结果表明：（1）2015年长三角地区休闲游憩服务的供需比为1028%,供给远大于需求,但供需关系空间异质性明显,高供给区与高需求区明显错位,空间不匹配问题严重。（2）随着可达半径由0 km扩大至50 km,区域赤字区面积占比由43.14%减少至7.83%,休闲游憩服务需求未得到满足的人口占比由65.47%下降至10.64%。（3）根据41个地级市不同可达半径下休闲游憩服务供需关系的演变规律,将其划分为5种类型区,提出相应的发展建议。研究表明,快速城市化和城市扩张过程中,应强化城区内绿色空间的分散化配置;同时,在长三角经济一体化发展过程中,休闲游憩服务供需的盈余区、赤字区可以建立生态补偿和经济合作关系,促进长三角绿色发展的一体化。     

Optimal conservation planning of multiple hydrological ecosystem services under land use and climate changes in Teshio river watershed,northernmost of Japan

Most anthropogenic activities impacted on water quality and quantity, and further impacted on ecosystem services (ESs) in watershed are related to land use and climate changes those may cause losses of ecosystem functions. Effective information regarding ESs and their optimal priority conservation planning responded to land use and climate changes provide useful support for diverse stakeholders in ESs planning, management and policies. This study integrated the approach of spatially explicit ESs (water yield, inorganic nutrient, organic nutrient and sediment retentions) by using hydrology and material flow model (Soil and Water Assessment Tools, SWAT model) into systematic conservation of hydrological ESs according to land use and climate changes in Teshio watershed located in the north of Hokkaido, Japan. We investigated the spatial patterns and the hotspots of ESs changes to determine the spatial pattern of changes in systematic conservation optimal area of ES protection in terms of ESs protection targets. Under the land use and climate change scenarios, the forest land use significantly affected on the water yield, sediment, organic-Nitrogen (N) and organic-Phosphorous (P) retentions. The agricultural land (paddy and farmland fields) impacted on the inorganic-N and inorganic-P retentions. We applied the systematic conservation model (MARXAN model) to optimize the area for management of hydrological ESs satisfied the protection targets (30% and 50% of potential maximum ESs values among all scenarios) in all and individual ecosystem services, respectively. The simulated results indicated that the areas of spatial optimal ESs protection for all hydrological ESs were totally different from those for individual ESs. For bundles of ESs, the optimal priority conservation areas concentrated in southwest, north, and southeast of this watershed, which are related to land use, topography and climate driving factors. These places could guarantee ESs sustainability from both environmental protection and agricultural development standpoints. The priority conservation area turned more compact under climate change because the increased precipitation and temperature increased ESs amount. For individual ESs, the optimal priority conservation areas of water yield, sediment retention and organic nutrient retention were traded off against those of inorganic nutrient retention (lower Jaccard's indexes and negative correlations of selection times). Especially, the negative correlation of selection times increased as the conservation target increased from 30% to 50%. The proposed approach provided useful information for assessing the responses of ESs and systematic conservation optimal planning to the land use and climate changes. The systematic conservation optimal areas of hydrological ESs provided an effective trade-off tool between environmental protection (sediment and organic nutrient retentions) and economic development (water yield and inorganic nutrient retention).