南四湖位于山东省西南部，北纬34˚15'~34˚53'和东经114˚50'~117˚48'之间，全流域跨山东省、河南省、江苏省和安徽省四个省份。作为华北地区最大的淡水湖，南四湖在地区生态系统中扮演着重要的角色，同时也是山东省的关键水源地和生态宝库。南四湖流域作为典型的农业流域，耕地面积占比高达69%，是流域内主要的土地利用类型。地形上，南四湖流域呈现出多样的特征，地势东高西低(图1)。气候方面，该流域属于暖温带半湿润地区，表现出明显的暖温带大陆性季风气候，四季分明，雨热同期。南四湖流域是我国重要的湿地保护区，其丰富湿地环境为各类生物提供了良好的栖息地。流域地下矿产资源丰富，在能源需求方面发挥着关键作用。此外，南四湖流域还以其广泛的耕地面积和适宜的气候条件成为粮棉生产基地，农业活动的发展不仅为当地经济注入活力，也为国内粮食和棉花产量做出了贡献。总之，南四湖流域多元化的资源特征使其在湿地保护、农业生产、渔业和能源方面发挥着重要的战略作用。

图1. 南四湖流域研究区概况

本研究所采用的数据涵盖多个方面，包括土地利用地图、气象数据、土壤分布数据、数字高程模型(DEM)等社会数据(表1)。研究采用的流域2018年土地利用数据(1:10万)，分辨率为30 m，由中国科学院资源环境科学与数据中心(http://www.resdc.cn/)提供。本研究参考刘纪远提出的中国土地资源分类系统，将土地利用划为六个一级地类。数字高程模型(DEM)数据的来源是美国马里兰大学全球土地覆被数据库(https://glad.umd.edu/)，分辨率为30 m。气象数据来自中国气象数据网(http://data.cma.cn/)。土壤数据包括了最大根系深度、土壤沙粒、粉粒、黏粒和有机质含量等参数，来自世界土壤数据库(HWSD)中的中国土壤数据集(1:100万)。行政区划矢量数据来自全国基础地理数据(1:400万)。

表1. 本研究数据来源

在InVEST模型中测算南四湖流域的六种生态系统服务，评价方法见表2。参数的设置和校验参照了InVEST模型使用手册和课题组前期成果 [ 6 ] 。通过合理的参数设置，确保了模型对于南四湖流域的适用性和准确性。

表2. 本研究选取的生态系统服务及其评价指标

为了深入了解生态系统服务对区域可持续发展的贡献，本研究构建了一套流域生态系统服务效益评估方法。即对数据进行标准化，这一步骤有助于消除不同生态系统服务数据之间的单位差异，使得评估更具科学性，为决策者提供更为可靠的信息。

本研究考虑到城市化进程和现有土地政策的影响，故而采用CLUE-S模型进行土地利用的模拟。基于研究区1975年和2018年的土地利用数据，预测了南四湖流域在城市自然扩张背景下2050年的土地利用情况。

本研究采用多种土地利用情景模拟，设定2050年社会经济发展现状和生态保护目标，通过模拟相应的土地利用情景作为背景变量，进行五种情景的设定，旨在更为全面地考虑不同土地利用决策对生态系统服务的潜在影响，进而深入了解各种土地利用模拟情景下生态系统服务之间的相互作用关系，情景设定见表3。

表3. 五种情景的设定及模拟方法

本研究模拟了五种土地利用情景，其土地利用分布如图2所示。表4介绍了五种情景下的土地利用结构。

表4. 五种情景下的土地利用结构(km²)

图2. 五种情景下土地利用分布图

基于CLUE-S模型，我们模拟了2050年的土地利用模式(情景1)。在这一情景中，耕地作为流域内主要的用地类型，在情景1、情景3和情景4下被大量用作补充其他用地，因而相比2018年呈现大幅下降趋势。林地占地面积在情景3和情景5下呈现明显的上升趋势，分别增加了181.75 km²和40.13 km²。林地在情景1和情景2下有所上升，但增加量较少。草地和湿地在情景1下主要用于弥补其他用地需求，导致这两种地类的面积分别减少了10.43 km²、61.59 km²。此外，在情景1下，由于城市化的推进，越来越多的土地被开发为建设用地，建设用地面积相较2018年增加了543.83 km²。在情景4下，专门针对湿地的生态保护政策导致湿地面积大幅度增加，增加了181.75 km²。

使用InVEST模型预测上述五种情景下六种生态系统服务的变化，结果见图3和图4。在情景1下，产水量增加，而水质净化服务下降，说明水体中存在较多农业化学物质、工业废水或城市污水等，降低了水体的自净能力；碳存储和授粉服务的下降与林地转化为其他用途、土地退化、生境破坏有关。情景2与情景3下的变化相似，其中授粉服务和生境质量变化最为显著。具体而言，在情景2和3下，生境质量的值分别增加了11.50%和12.50% (图3(f))，由此可见，退耕还林、保护自然植被等措施有助于区域生境质量的提高；与此相反，授粉服务则呈现出下降趋势。在情景4下，产水功能、碳存储和授粉呈现下降趋势，这与湿地恢复过程中植被减少和花粉传播者的栖息地减少有关；而水质净化、土壤保持和生境质量出现了不同程度的提高，退耕还湿导致湿地面积增加，有助于改善生境的多样性；生境质量的变化最为显著，增加了13.50%，是五种情景中增加最多的(图3(f))，而授粉服务减少了23.76% (图3(g))。在情景5下，产水量、土壤保持和授粉服务下降，这与建设用地特征、土地利用结构变化等因素有关；而水质净化、碳存储和生境质量有所改善，说明植被的增加和多样性的提高可以提高碳的储存能力，并为各种生物提供更适宜的生存条件。

图3. 五种情景下的各类生态系统服务价值及其变化情况

总之，在这五种情景下，与其他生态系统服务相比，生境质量和授粉服务的变化更加显著。此外，生境质量在各个情景下都有明显的改善(图3(f))，而授粉服务则显著下降(图3(g))。此外，在情景1下，产水量和土壤保持服务变化显著，然而，除了授粉服务和生境质量之外，其他生态系统服务的变化方向并不一致。

图4. 五种情景下生态系统服务分布

通过相关分析，获得五种模拟情景下研究区生态系统服务相关性热图(图5和图6)。两幅图展示了在不同土地利用模拟情景下，六种生态系统服务之间的相互关系及其强度，并与2018年的情况进行了比较。

图5. 2018年及五种情景下各类生态系统服务相关性

图6. 五种情景下各类生态系统服务相关性变化

根据2018年以及五种不同情景下的数据，我们观察到水质净化(N输出)与水质净化(P输出) (WppWpn，图6(b))、水质净化(N输出)与生境质量(WpnH，图6(b))、水质净化(P输出)与生境质量(WppH，图6(c))、碳存储与授粉(CP，图6(d))之间存在着明显的协同作用，且这些生态系统服务之间的相关系数都高于0.5 (图5)。另外，在2018年以及五种不同情景下，产水量与土壤保持(WyS，图6(a))、水质净化(N)与碳存储(WpnC，图6(b))、水质净化(P输出)与授粉(WpnP，图6(b))、水质净化(P输出)与碳存储(WppC，图6(b))、碳存储与生境质量(CH，图6(d))之间协同关系较弱，相关系数约为0.1 (图5)。然而，在2018年以及五种不同情景下，产水量与水质净化以及产水量与生境质量之间权衡关系较强(图6(a))，相关系数均处于0.4~0.5之间。值得注意的是，其他生态系统服务之间的相关性很弱，且不显著。

综上所述，城市自然增长情景与四种生态保护情景相比，各服务间的相互关系情况存在显著差异，但是在四种生态保护情景下差别相对较小。此外，在城市自然增长情景下，大多数生态系统服务之间的权衡作用更强，而协同作用相对更弱。

将五种情景下的生态系统服务数据进行标准化，处理结果见表5。通过比较总体效益，得出以下排名：情景3 > 情景5 > 情景2 > 情景4 > 2018年> 情景1。此外，相对于情景1而言，各种生态保护情景均呈现出较高的总体效益，生态保护情景的高效益可能与其有助于改善水质、提高生境质量、增加碳存储等生态系统服务相关。相反，城市自然增长情景可能受到城市扩张、资源开发等负面因素的影响，导致生态系统服务效益的降低。

表5. 不同情景下的生态系统服务效益

近年来，通过制定生态环境保护政策和调整土地利用规划，人们希望能最大化地提高生态系统服务对人类社会的效益在管理生态系统时，必须要考虑多种生态系统服务间的相互作用，以避免偏好某一服务导致其他服务的损害和丧失 [ 20 ] 。然而，许多研究过于强调生态系统服务之间的相互作用或土地利用对生态系统服务的影响，而忽略了在不同生态政策下预测多个生态系统服务之间相互作用以优化某一地区土地利用的研究 [ 27 ] [ 28 ] 。本研究的独特之处在于不仅考虑了生态系统服务之间的相互作用，还关注了在不同土地利用情景下，生态政策对多个生态系统服务的影响。这一综合性的研究设计旨在更全面地评估土地利用的优化可能性，为制定有效的生态政策提供科学依据。

根据土地利用的潜在变化，本文提出了5种土地利用情景，以便为南四湖流域土地利用的未来发展提供一个框架，而不是侧重于单一的预测结果。研究结果表明，在各生态保护情景下，生态系统服务之间的权衡协同关系呈现出健康的发展趋势。因此，决策者进一步落实生态保护政策，更有利于城市生态系统的健康发展，有利于生态安全问题的解决。

本研究通过深入研究生态系统服务、土地利用情景和相互关系的变化，旨在为制定更有效的生态保护政策和可持续发展战略提供指导。然而，由于本人专业知识以及研究水平等方面有限，论文存在一定局限性。例如，本研究仅构建了五种场景来研究生态系统服务，最优场景仅限于这五种场景之中。随着城市化的快速发展，有必要构建大量的方案，以分析未来不同土地利用结构下生态系统服务之间的权衡和协同效应，进一步提高研究的深度和广度，以便根据地区发展的需要，选择最为适宜的方案，为未来的规划提供坚实的基础。

(1) 在五种模拟的情景中，生境质量和授粉服务的变化最为显著。与其他四种生态保护情景相比，在城市自然增长情景下，除了碳存储和生境质量外，其他四种生态系统服务的变化程度不如前者大。

(2) 在城市自然增长情景下，生态系统服务之间的相互关系表现出比其他四种情景更为错综复杂的变化。在该情景下，大多数生态系统服务之间的权衡作用呈上升趋势，而协同作用则显示减弱的趋势，这意味着城市扩张可能对流域内的生态平衡和生态系统服务产生负面影响。

(3) 四种生态保护情景的整体效益均高于城市自然增长情景。此外，在五种情景中，滨岸带退耕还林生态保护情景下总体效益最佳，这表明这一情景对提高生态系统服务效益具有显著的正面影响。因此，为了可以更有效地平衡经济发展和生态环境保护的关系，未来规划应偏向滨岸带退耕还林生态保护情景，优化土地利用结构，切实保护好研究区的生态用地资源。