1、第 卷 第 期 年 月中 国 水 土 保 持 科 学 近 年东北沙地生态系统服务动态张 丽，刘佳雨（吉林师范大学旅游与地理科学学院，吉林四平；吉林师范大学附属小学，吉林四平）摘要：东北沙地位于北方农牧交错带，揭示该区长时间序列生态系统服务功能时空演化特征，可为沙地生态恢复与水土保持提供科学依据。基于 模型，对 年该区土壤保持、水源涵养、碳储存量和生境质量功能进行评估。结果表明：）东北沙地土壤保持总量呈现波动增加的变化趋势，增幅.亿 ，高值区集中于科尔沁沙地的西南地区、松嫩沙地的东部和西部以及呼伦贝尔沙地的东部；）水源涵养量总体呈现增加趋势，高值区集中于科尔沁沙地东南部、松嫩沙地东北部和呼伦贝尔
2、沙地的东部；）碳储存量呈增加趋势，增幅.亿 ，高值区主要分布于科尔沁沙地东北部、松嫩沙地的东部、西部和呼伦贝尔沙地的东部；）平均生境质量指数总体变化不大，呈现先上升后下降趋势，高值区域集中于松嫩沙地东部和西部、呼伦贝尔沙地的东部。东北沙地土壤保持、水源涵养、碳储存量和生境质量 大生态系统服务功能总体呈现增加趋势，增加的区域在空间上表现一定聚集性。关键词：生态系统服务功能；模型；风沙地貌；东北沙地中图分类号：.文献标志码：文章编号：（）：收稿日期：修回日期：项目名称：吉林省社会科学基金“乡村振兴背景下的智慧农业创新模式与政策协同机制研究”（）第一作者及通信作者简介：张丽（），女，博士，教授。主要
3、研究方向：资源与环境信息系统。：，（，；，）：（，），）.，）.，中国水土保持科学 年 ）.，），：；生态系统是在一定的区域范围内，生物与自然地理环境之间相互作用而形成的动态综合体，是人类生存发展的基础，是无法替代的自然资源和自然资产。生态系统及其生态过程不断地提供生态系统物品和服务，从而形成与维持着人类赖以生存的环境条件和物质基础。生态系统服务是在生态系统与生态过程中所能维持和满足的人类赖以生存的各种环境条件与效用，是生态系统在形成过程中产生对人类所需的功能和价值，包括气候调节、净化环境等。生态系统服务对区域生态系统服务功能进行长时间序列分析，可为水土保持与区域可持续发展提供科学依据。沙地是
4、一种典型的荒漠生态系统。东北沙地（科尔沁沙地、松嫩沙地和呼伦贝尔沙地）属于半干旱、半湿润区向半干旱区的过渡地带，为气候变化敏感带、生态环境脆弱带和农牧交错带。该区沙地土地沙化严重，是国家防风治沙的重点区域，是开展相关研究的理想靶区。国内外相关学者 基于生态系统服务能力综合估价和权衡得失评估模型（，），针对生态系统服务功能进行了大量的研究和应用。如 等基于 和遥感数据，利用 生境质量模型评估 年 流域不同土地利用类型下生境质量的时空演替特征并探讨了影响因子与生境质量的相关性。等利用马尔可夫链和模型对、和 年不同情景下基于总碳变化的碳固存进行评估和预测。笔者对 年东北沙地 种主要生态系统服务功能进
5、行分析，以期为东北沙地生态系统保护提供科学参考。研究区概况东北沙地分布在内蒙古自治区、辽宁省、吉林省和黑龙江省等省（区）的部分地区（，），地处内蒙古高原和东北平原的过渡地带并且处于亚洲东部温带季风气候区。该区风沙地貌发育受风水交互侵蚀作用。以水系流域分割方法，确定东北沙地的范围，面积.万（图）。其中，科尔沁沙.万，年均气温.，年降水量 ，年蒸发量 ；松嫩沙地位于松嫩平原的中西部，面积.万，年均温.，年降水量 ，年蒸发量 ；呼伦贝尔沙地位于大兴安岭西麓的呼伦贝尔高平原，以固定、半固定沙地为主，面积约.万，年平均气温 ，年平均降水量 ，年蒸发量 。方法与数据.生态系统服务模型）土壤保持模型。土壤保
6、持模型是基于修正的通用土壤流失方程，通过潜在土壤流失量减去实际土壤流失量来评估不同土地利用 覆被类型的土壤保持服务，公式见参考文献。）土壤可蚀性因子 的计算。因子代表土壤抵抗侵蚀的能力，反映土壤对侵蚀的敏感性，是土壤渗透率、渗透量、容水量以及土壤颗粒抵抗分散和输送能力的综合反映。值的计算方法很多，结 第 期张丽等：近 年东北沙地生态系统服务动态图 东北沙地所在位置 合研究区特点采用 模型中提出的 值计算方法。）水源涵养模型。利用 模型中的水源涵养模块来量化产水量。该模块基于水平衡原理，根据年平均降水量与年实际蒸散发之差来评估各个土地利用类型 上每个像素的年水源供给量，计算式见参考文献。）碳储存
7、量模型。碳储存量模型把生态系统碳储存量分为 个子碳库，即地上生物碳、地下生物碳、土壤碳和死亡有机碳，计算式见参考文献。）生境质量模型。生境质量结果描述值范围为，代表生境适宜度最高。模型生境质量（）模块假设生境质量与生物多样性正相关，用生境退化程度来计算生境质量。.数据来源与预处理遥感影像选取 年 个时期 影像，数据来源于地理空间数据云（：）。基于 软件对 系列影像进行多波段合成、几何精校正（误差控制在一个像元内）、图像拼接和按研究区范围裁剪等预处理。根据东北沙地立地条件，划分 类，即耕地、林地、草地、水域、建设用地、沙地和未利用地等。为提高解译精度，基于 .软件平台，采用面向对象分割方法进行遥
8、感影像的解译，对不明确的点位，采用实地调研，获取野外地面验证点。气象数据下载于国家气象信息中心（：），土壤数据下载于寒区旱区科学数据中心（：）。.数据分析.土壤保持主要包括高程数据，降水量，黏粒、粉粒、砂粒、有机质含量，植被覆盖管理因子，水土保持措施因子，并利用这些数据进行如下计算：）高程数据来源于地理空间数据云（：），根据研究区范围进行镶嵌和裁剪，可得到研究区域高程图。）降雨侵蚀因子。该因子受研究区降雨强度和历时的影响，单次降雨强度越大和历时越长，降雨侵蚀力（）越大。降雨侵蚀因子（）参照 等的研究计算。）土壤可侵蚀性因子。值越大，说明土壤较易受到侵蚀。土壤可侵蚀性因子 参照王蓓等的研究计算。
9、）生物物理参数表。通过参考 模型手册和相关文献 获得。.水源涵养 模型水源涵养模块（）运行所需数据包括降水量、潜在蒸发量、土壤质地和有机质数据、土地覆被、土壤最大根系深度、植被蒸散系数和高程数据，并利用数据进行如下计算：）年均潜在蒸散量。计算采用 法。）植物可利用水含量（，）受土壤结构、有机质含量和土壤密度影响。选择 等的“”公式进行计算。）生物物理参数表。包括土地覆盖、土壤最大根系深度和植被蒸散系数（表）。）系数。通过多次模拟与校正，值定为.时，模型评估结果最优。）子流域。模型中产水量的计算是在流域尺度上进行的，根据研究区 数据进行流域提取及次一级小流域的划分。.碳储存量本研究碳密度表数据采
10、用前人研究中的数据，通过单位换算，得到碳密度表（表）。.生境质量 模型生境质量模块运行所需数据主要包括土地利用类型、胁迫因子、最大距离、权重、衰减率，并利用数据进行如下计算：）威胁源。选取耕地、建设用地、沙地和未利用地作为威胁源。）生境敏感型参数见表。）半饱和常数。半饱和常数参考 模型手册设置为.。中国水土保持科学 年表 生物物理参数表 土地利用类型 土地覆被 土壤最大根系深度 植被蒸散系数 耕地 .林地 .草地 .水域.建设用地 .沙地 .其他.表 土地利用覆盖碳密度 土地利用类型 地上生物碳密度 地下生物碳密度 土壤碳密度 死亡有机物碳密度 耕地 .林地.草地 .水域.建设用地 .沙地 .
11、其他.表 生境敏感性参数 土地利用类型 生境威胁源耕地 威胁源建设用地 威胁源沙地 威胁源未利用地 耕地 .林地.草地 .水域.建设用地 沙地 其他 生态系统服务功能.土地利用变化根据土地利用分类结果，基于 得到各地类面积大小依次为耕地、草地、林地、其他未利用土地、水域、建设用地和沙地。年间，耕地面积增幅最大，共增加 万 .；林地面积处于减少的趋势，共减少 .；草地面积减幅最大，共减少 万 .，年减少最快，减少 万.；水域面积也在不断减少，共减少 .；建设用地面积呈现稳定的增长趋势，共增加 .；沙地面积呈现先下降再增长的趋势，年减少.，后又增加到 万.；其他未利用土地呈现增长趋势，共增加 万.
12、。.土壤保持.年际变化东北沙地土壤保持量总体呈现先减少后增加的变化趋势（图）。年达到最大 第 期张丽等：近 年东北沙地生态系统服务动态图 年东北沙地土壤保持总量 值，为.亿。年土壤保持总量共增加.亿，年均增加量为.亿。年土壤保持总量呈减少趋势，减少.亿，年均减少量为.亿。年土壤保持总量呈现不断增长的趋势，总共增加.亿，年均增加量为.亿。.空间分布变化 年间，东北沙地的土壤保持空间分布特征非常明显（图），土壤保持量高值区域主要位于科尔沁沙地的西南地区、松嫩沙地的东部和西部以及呼伦贝尔沙地的东部，而其余地区土壤保持量均较低。图 年东北沙地土壤保持服务空间分布 .水源涵养.年际变化东北沙地在水源涵养
13、量呈波动增加趋势（图），净增加.亿，年均增加量为.亿。年水源涵养量呈减少趋势，减少.亿，年均减少量为.亿。从 年起水源涵养量呈现不断增加的趋势，至 年总共增加.亿，年均增加量为.亿。.空间分布变化 年间，东北沙地水源涵养量呈现增加的趋势，且出现地域集聚性特点（图）。其中，年水源涵养量高值区域主要集中于科尔沁沙地东南部、松嫩沙地东北部和西部部分地区以及呼伦贝尔沙地的东部，低值区域集中于松嫩沙地中部以及呼伦贝尔沙地西部地区；年水源涵养量高值区域主要集中于松嫩沙地东部和呼伦贝尔沙地的东部；年水源涵养量高值区域集中于科尔沁沙地东南部和松嫩沙地的东北部，其余地区水源涵养量均较低；年除呼伦贝尔沙地西部以及
14、科尔沁沙地西部以外，其余地区水源涵养量均处于中高水平。.碳储存量.年际变化东北沙地碳储存量呈现波动增中国水土保持科学 年图 年东北沙地产水量 加趋势（图），其值在.亿.亿 之间波动，研究期净增加.亿，年均增加量为.亿。东北沙地碳储存量在 年达到最大值，为.亿。.空间分布变化 年间，东北沙地主要以中高水平碳储量分布为主（图），碳储量高值区主要分布于科尔沁沙地东北部、松嫩沙地的东部、西部和呼伦贝尔沙地的东部。中等水平的碳储量主要集中于松嫩沙地中部和呼伦贝尔沙地西部；碳储量低值区主要集中于科尔沁沙地西南部，松嫩沙地和呼伦贝尔沙地少量区域。图 年东北沙地产水量空间分布 图 年东北沙地碳储存量 .生境质
15、量.年际变化 东北沙地在 年的平均生境质量指数变化不大，呈现先上升后下降的变化趋势。林地的生境质量最高，其次是草地、水域和耕地。林地的生境质量最高，其在 个年份的值分别为.、.、.和.，其次是草地，其在 个年份的数值分别为.、.、.和.。.空间分布变化 年间，东北沙地生境质量高值区域主要分布于松嫩沙地东部和西部、呼伦贝尔沙地东部（图）。中值区域主要集中 第 期张丽等：近 年东北沙地生态系统服务动态图 年东北沙地碳储存量服务空间分布 于科尔沁沙地、呼伦贝尔沙地西部以及松嫩沙地中部的部分区域。低值区域主要集中于松嫩沙地中部。讨论从东北沙地近 年生态系统服务功能的结果来看，区域内各项生态系统服务功能
16、量变化趋势大致相同。土壤保育与水源涵养变化趋势一致，在 年达到最大；碳储存量在 年达到最大；生境质量在 年最好。年水源涵养量特别小，年碳储存量最大。这与陈童尧等研究一致。这可能与土地利用方式与国家恢复政策有关，此外，地理环境因子对区域内各项功能也有显著影响。对比结果分析可知，区域内的土壤保持量、水源涵养量、碳储量和生境质量指数的变化均以 年为节点，大体上呈现先下降后上升的变化特征。其中土壤保持总量增幅为.亿 ，水源涵养量增幅为.亿 ，碳储存量增幅为.亿，生态系统服务功能整体变化趋势与张华等在科尔沁沙地的研究结果具有相似性。从空间变化来看，各项生态系统服务功能的高值区空间分布大体上保持一致，并呈现空间的集聚性，表明生态系统服务功能变化受控于土地利用变化。这与宋洁等和孙兴辉等的研究结果较为一致。东北沙地地区处于干旱半干旱区农牧交错带，生态环境比较脆弱，不同生态系统服务会受到相同影响因子的作用而发生改变，其时间变化的关联性在不同影响因子变化过程中甚至会发生显著变化。在本次研究阶段内各项生态系统服务功能量的变化均以 年为节点，这与国家的生态保护修复措施有关，如“三北”防护林工程、封沙育林育草、