人类活动影响下蒙辽农牧交错区的生态系统服务空间格局_徐梦冉.pdf

第 25 卷第 1 期2023 年 1 月大 连 民 族 大 学 学 报Journal of Dalian Minzu UniversityVol25，No1January 2023收稿日期:20220804;最后修回日期:20221128基金项目:辽宁省联合基金项目(2020MZLH11)。作者简介:徐梦冉(1997)，女，内蒙古锡林郭勒盟人，大连民族大学环境与资源学院硕士研究生，主要从事生态系统评估研究。通讯作者:张靖(1981)，男，内蒙古赤峰人，讲师，博士，主要从事生态系统评估等方面研究，Email:zhangjing dl-nueducn。文章编号:20961383(2023)01002107人类活动影响下蒙辽农牧交错区的生态系统服务空间格局徐梦冉，张靖，鲍雅静，李政海，刘睿，张贻龙(大连民族大学 环境与资源学院，辽宁 大连 116605)摘要:以蒙辽农牧交错区为研究区，通过 InVEST 模型，计算其 2000 年2019 年的产水量、土壤保持、碳固持和生物多样性保护等 4 种生态系统服务指标，分析各生态系统服务指标的空间分布格局及空间自相关性，探讨人类活动对蒙辽农牧交错带生态系统服务空间格局的影响。结果表明，2000 年2019 年研究区产水量空间差异明显，呈现出东南高西北低的分布格局;20 年间土壤保持呈现增加趋势;碳固持量整体偏低，生态系统多样性保护服务量较高，但也存在空间异质性。4 种服务均呈现不同程度的空间自相关性，具有空间聚集性;2000 年和 2010 年产水量的 Morans I 分别为 0676 5、0682 3，空间聚集性最强;生物多样性保护的 Morans I 提高最明显，从 2010 年增长至大于 07，具有较强的空间聚集性。关键词:蒙辽交错区;生态系统服务;人类活动;空间格局中图分类号:Q148文献标志码:ASpatial Pattern of Ecosystem Services under the Influence of HumanActivities in Ecotone of Inner Mongolia and LiaoningXU Mengran，ZHANG Jing，BAO Yajing，LI Zhenghai，LIU ui，ZHANG Yilong(School of Environment and Bioresources，Dalian Minzu University，Dalian Liaoning 116605，China)Abstract:In this study，four ecosystem service indicators(water yield，soil conservation，car-bon sequestration and biodiversity conservation)were calculated by InVEST model in the Mon-goliaLiaoning egion from 2000 to 2019，and the spatial distribution pattern and spatial auto-correlation of each ecosystem service index were analyzed，to explore the effects of human activi-ties on the spatial pattern of ecosystem services in the Mengliao agropastoral ecotone The re-sults showed that the water yield in the study area was significantly different from 2000 to 2019，showing a distribution pattern of high water yield in southeast and low water yield in northwestSoil conservation showed an increasing trend in the past 20 years The overall carbon sequestra-tion was low，and the ecosystem diversity protection services were high，but there was also spa-tial heterogeneity The four services showed different degrees of spatial autocorrelation and spatialaggregation Morans I of annual water volume in 2000 and 2010 were 0676 5 and 0682 3，re-spectively，showing the strongest spatial aggregation Morans I of biodiversity conservation in-creased significantly from 2010 to more than 07，with strong spatial aggregationKey words:ecotone of Inner Mongolia and Liaoning;ecosystem services;human activity;spa-tial patternDOI:10.13744/21-1431/g4.2023.01.009生态系统服务是指生态系统形成和所维持的人类赖以生存和发展的环境条件与效用1，它支撑与维持了地球的生命过程，为人类提供生活、生产所需的原料2。在过去几十年，由于人类对生态系统服务及其重要性缺乏充分认识及科学而有效的管理3，如对供给服务(如粮食)的过分追求4，导致除供给服务以外的大多数生态系统服务呈不同程度的降低5，6。生态系统服务的退化和丧失必将对人类的福祉产生不利的影响，直接威胁区域和国家，甚至全球的生态安全2。在这种情况下，生态系统服务作为生态系统管理的有效工具，受到学术界越来越高的关注。而生态系统服务管理的关键在于明确其空间分布特征及聚集程度，以把握服务供给、明晰空间异质性。因此，研究人类活动影响下的生态系统服务空间格局，对人类生存空间与生态环境可持续发展具有重要意义。国内外学者评估生态系统服务的方法主要分为物质量评估法7 和价值量评估法8，传统的价值量评估虽然简单实用，但科学性较低，例如，仅估算生态系统服务的价值，则会忽略其空间异质性9;而物质量评估法可以根据生态系统服务的形成机制进行评价，结果可靠，解释力高，因而应用广泛10。随着对生态系统服务研究的深入和3S 技术的发展，众多模型应运而生。InVEST 模型(Integrated valuation of ecosystem services andtradeoffs)是由美国资本项目组研发的生态系统服务评估模型，它广泛应用于各种尺度的评估，其结果空间可视化、评估精度方面更具优势，可用于量化多种生态系统服务功能(如产水量、土壤保持、碳固持和生物多样性)1114。例如，Leh 等运用InVEST模型中的水源涵养、碳固持、土壤保持和生境质量模块，分析了科特迪瓦和加纳地区的生态系统服务功能空间分布和时间变化情况15。Himlal Baral 等人利用 InVEST 模型评价了澳大利亚维多利亚州不同土地覆盖情景下生物多样性的变化情况16。国内学者对 InVEST 模型的应用多位于长江上游、北京等地，基于农牧交错区的研究较为罕见。例如，方露露等人对比长江和黄河流域 2000 年2015 年植被净初级生产力(Net PrimaryProductivity，NPP)、产水服务和土壤保持的时空变化规律17。吴丹等人基于遥感影像和 InVEST模型，探究了长江经济带的产水量、土壤保持、碳固持与生物多样性的生态系统服务功能18。如今，基于 InVEST 模型关于生态系统服务功能的研究非常成熟。农牧交错带是农业和牧业两种生产方式间的交错区。近年来，在不合理人类活动和全球变化的双重影响下，蒙辽农牧交错带的草地退化明显、土地荒漠化现象加剧，给当地人民生产、生活带来了极大危害，并对我国东部地区的生态环境和经济发展带来了不良影响19。为此，分析和评估人类活动对生态系统服务的影响是十分必要的。因此，本文以蒙辽农牧交错区为研究对象，采用InVEST模型定量评估该地区产水量、土壤保持、碳固持及生物多样性保护等 4 项生态系统服务指标，并利用 Morans I 指数描绘出各生态系统服务的空间聚集特征，探究人类活动对生态系统服务空间格局的影响，以期为蒙辽农牧交错区及其他区域的可持续发展提供科学的数据。1研究区概况本文依据史文娇等20 划定的北方农牧交错带界线确定蒙辽农牧交错区位于内蒙古自治区和辽宁省的交界处，坐标范围:E11905E12231，N4137N4436，具体包含辽西北的 10 个半农半牧县(康平、彰武、阜蒙、建平、北票、喀左、义县、凌源、朝阳、建昌)，内蒙古通辽的 7 个旗县(开鲁、科左中旗、科左后旗、库伦、通辽县、扎鲁特、奈曼)和赤峰的 9 个旗县(翁牛特、巴林左、巴林右、林西、敖汉、阿鲁科尔沁、喀喇沁、克什克腾、宁城)。研究区 总 面 积 13 9 104km2，其 中 草 地 面 积57104km2。2数据与方法21数据来源与处理以蒙辽农牧交错区为研究对象，通过 InVEST模型计算 2000 年2019 年逐年蒙辽农牧交错区的产水量、土壤保持、碳固持和生物多样性保护等4 种生态系统服务。数据来源如下:归一化植被指数(Normalized Difference Vegetation Index，ND-VI)，采用 MODIS MOD13Q1 NDVI(2000 2019年，时间分辨率16 d、空间分辨率250 m，数据获取地址:https:/ladsweb modaps eosdis nasa gov/search/);年降水量、10 积温、降水侵蚀力()、潜在蒸散发(Potential Evapotranspiration，PET)和陆表太阳辐射(s)等气候参数来源于中国地面气候资料日值数据集(V30)，下载地址:22大连民族大学学报第 25 卷http:/datacmacn/，经在 MySQL80 和 ArcGIS 中处理获得;土壤参数，如沙粒、粉粒、黏粒、有机碳、碳酸盐的比率或含量来源于 Harmonized WorldSoil Database V12，下载地址 http:/westdcwest-gis accn/;土地利用图，基于中国科学院资源环境科学与数据中心每 5 年更新的土地利用数据，分辨率 30 m。22生态系统服务计算与空间格局分析221生态系统服务计算(1)产水量，由 InVEST 产水模型计算获得，公式如下:Yieldxj=(1 AETxj/Px)Px;(1)AETxj/Px=1+PETxj/Px1+(PETxj/Px)x1/x。(2)式中:Yieldxj为某一栅格中土地覆盖类型年单位区域产生的径流量，mm;Px 为栅格的一年中平均降水量(mm)，AETxj为栅格中土地覆盖类型包括植被蒸腾、土壤水分蒸发等真正失去水分量，mm;x为自然气候等非物理参数;PETxj为栅格中土地覆盖类型 j 的年有可能会蒸散量，mm21。(2)土壤保持服务，采用土壤保持服务量来表示，计算公式如下:SC=K L S 1 P C()。(3)式中:SC为土壤保持服务量，t(hm2a)1，是降雨侵蚀力因子，MJ(mmhm2ha)1，K 是土壤可蚀性因子，MJ(mmhm2ha)1，L 为一定坡度下单位面积土壤流失量与标准坡度流失量的比值;S 为定坡度下单位面积土壤流失量与标准坡度流失量的矢量之比(无量纲)，C 为作物覆盖的比率和人为管理的因素(无量纲);P 为水土保持所做的措施的数量因子(无量纲)。降水侵蚀力()指的是一个时间跨度内降水侵蚀力，通过月尺度的降雨量可以估算出一个降雨侵蚀力的模型22。土壤可蚀性(K)反映土壤遭受侵蚀的敏感程度，本文利用土壤可蚀性因子 K 值的估算方法来测算23。(3)碳固持，与 NPP 有关，计算采用 CASA 模型，其计算公式24 如下:NPP x，t()=APA x，t()x，t();(4)APA x，t()=SOL x，t()FPA x，t()05;(5)(x，t)=T1(x，t)T2(x，t)W(x，t)max。(6)式中:APA(x，t)为在 t 时间内 x 中能被植被光合作用吸收的有效光谱成分，g C(m2month)1，(x，t)为植被对光能资源利用效率的指标数值，gCMJ1;SOL(x，t)是 t 时间内 x 平面上太阳照射量的累计值，MJ(m2month)1;05 是指植被吸收有效的光谱成分与总太阳照射值的比率;FPA(x，t)为植被利用太阳辐射的利用比例;max为植被可以利用的太阳辐射的最大值的比率，gCMJ1;T1(x，t)和 T2(x，t)表示低高两种极端温度情况下光能利用的应力影响;W(x，t)为水分威胁程度常数。(4)生物多样性保护，采用 NPP 定量指标的方法来评估该地区生物多样性保护服务，计算公式25 如下:Sbio=NPPmean Fpre Ftem1 Falt()。(7)式中:Sbio为评价生物多样性维护服务的重要指标;NPPmean为植被通过光合作用所吸收的碳除去自身呼吸消耗的量的平均值;Fpre为所研究年份降水均值;Ftem为所研究年份当地的一年气温的均值;Falt为当地的海拔高度对生态环境影响的值，采用极限大小值的方法将这三个数归一到 0 1之间。222生态系统服务变化趋势分析采用 TheilSen median 趋势分析，分析 20 年间蒙辽农牧交错区生态系统服务变化趋势，这种方法是一种非常合理的用来统计趋势计算方法26，计算公式如下:=medianNDVIi NDVIji j，j i()。(8)式中:1jin。是计算 n(n1)/2 个数据排列组合之后绘制出斜率的中位数，用来把一个趋势用数字表现出来。当 0 时，说明在固定时间以内的植被数据或生态系统服务的各项指标呈增强趋势;反之，则说明该区域生态系统往相反的方向发展。223生态系统服务空间关系采取使用全局 Morans I 测算出各生态系统服务的空间自相关性以及各生态系统服务之间彼此功能的影响，公式如下:I=ni=1nj=1wij(xi xm)(xj xm)S2ni=1nj=1wij。(9)式中:wij为斑块 i 与斑块 j 之间不同研究区间的权32第 1 期徐梦冉，等:人类活动影响下蒙辽农牧交错区的生态系统服务空间格局重，xi、xj为一个区间内的属性值，n 为研究区内划分出的区间数量，S2为方差。莫兰指数值的范围是1010 之间。Morans I0 意味着研究区内各空间内成正相关，莫兰指数越大，空间之间的相关性就越明显，Morans I 0 表示研究区内各空间内成负相关，Morans I=0，意味着空间无规律性。224人类活动与生态系统服务空间格局分析在评价人类活动强度时，通常考虑人类足迹指数。探索人类足迹指数与生态系统服务之间的相互关系，计算人类活动强度同各生态系统服务间相关性。人类足迹指数可以将指数将人类影响与影响区域的相互作用考虑在内，计算公式如下:HFi=1 nj=1(1 hij)。(12)式中:HFi为人类足迹指数;hij为人类活动项得分;n 为人类活动类型数。根据连喜红等人27 的赋值标准对多期土地利用类型进行赋值见表 1。将人类活动足迹指数按照 0、00.4、0.40.6、0.60.85、085 划分为无、低、中、高、极高五个等级，应用 ArcGIS 10.3 随机采样工具获得采样点，在 SPSS 22.0 中采用皮尔逊相关分析，统计出人类活动强度与各个生态系统服务之间的相关程度。表1各人类活动类型得分标准地类得分地类得分城镇用地100土路070矿山090林地045硬化路085耕地080旅游080草地(500m 缓冲区外)030农村居民地070草地(500m 缓冲区内)0453结果与分析31生态系统服务空间格局由图1、图2 可知，2000 年2019 年蒙辽农牧交错区产水量空间分布差异明显，东南部地区较高、西北部地区较低;全域单位面积的产水量服务年均值为 36939mmm2，平均生态系统服务量为 5.081016t。研究区南部产水服务能力较高，西部地区产水服务能力较低。其中，产水服务分布较高的区域以林地为主，而产水服务较低的区域以草地为主。研究区产水量总体上呈现增加趋势，生态系统服务量由 2000 年的 4.481016t 增加至 2010 年的5.681016t。其中蒙辽农牧交错区的东南部产水量出现降低的情况，而东北部产水量增加量最高，其它区域产水量自东向西增加量减少。图 1蒙辽农牧交错区不同年份生态系统服务空间分布图图 22000 年2019 年蒙辽农牧交错区生态系统服务变化趋势20a 间该区单位面积的土壤保持服务均值为42大连民族大学学报第 25 卷7235 t(hm2a)1，平均生态系统服务量为 0991015t，空间分布差异明显。土壤保持服务较高的区域集中于研究区的南部、西南部及北部，地形多为山地、丘陵，虽然坡度较大土壤易被侵蚀，但地表植被茂密，林草地居多，对土壤有更好的持留作用。与土地利用图对比分析发现林地的土壤保持量最高，而耕地、草地和未利用土地的土壤保持较低，表明农田开垦、过度放牧等人类活动会引起土地退化和沙化，导致土壤保持服务减少。研究区土壤保持总体上呈现增加趋势，2000 年土壤保持量为 5782 t(hm2a)1，2010 年的土壤保持量为 86.89 t(hm2a)1。同时生态系统服务量也是增长的趋势，由 7921014t 增长至 1191015t。土壤保持增长的区域主要为蒙辽农牧交错区的南部、西北部及西南部。其余区域均无明显变化，面积占比为 5467%。研究区碳固持量整体偏低，年均值为 19059 gCm2，平均生态系统服务量为 2611013t C。碳固持的最小值集中于研究区西部、北部和中部，东部地区碳固持值相对较高。研究区的碳固持整体呈减少趋势，由 2000 年的 20436 g C/m2减少为2010 年的17681 g Cm2，生态系统服务的量也出现减少，由 2821013t C 减少至 2421013t C。其中，研究区的东北部碳固持明显减少，主要原因为该区域为草地生态系统，草地的退化、沙化可能导致了碳固持服务的减少。其西北部和南部的部分区域也有减少的情况。而蒙辽交界处的多个城市如建平县、敖汉旗等地碳固持量则为增加趋势。2000 年2019 年蒙辽农牧交错区生态系统多样性保护服务由东南向西北部递减，年均值为7541(无量纲)，平均生态系统服务为 1031019。其空间分布与产水量相似，其高值集中于研究区的南部及东部，土地利用以林地为主，加上研究区处于暖温带半干旱半湿润气候区，夏季炎热多雨，东南部地区雨量充沛，水热条件好，而水热条件又是限制植被生长的重要因素，所以该地区草地生态系统完整，生物种类丰富，生物多样性功能较高。生物多样性保护较低的区域大多为林草交界处。该服务上呈现降低的趋势，由 8534 降低为 65.49，服务量由 1 17 1019降低至 0 90 1019。研究区东部生态系统多样性保护服务减少最为显著，这可能是由于该地区草地的退化，使生态环境恶化，植被减少，其南部的部分地区也出现了生态系统多样性保护减少的情况，蒙辽交界处的多个城市如建平县、敖汉旗等地生态系统多样性保护服务增加并不明显，研究区的西北部无明显变化。32生态系统服务空间关系2000 年2019 年蒙辽农牧交错区 4 项生态系统服务呈现不同程度的空间正向自相关性(图3)。2000 年和 2010 年产水服务具有较强的空间聚集性(全局 Morans I 分别为0697 9 和0.680 7)，但在时间上该服务聚集性有减少的趋势;土壤保持服务聚集性 2010 年比 2000 年有所增加，分别为0676 5 和 0682 3;碳固持的全局 Morans I 则由2000 年的 0557 9 增加至 2010 年的 0649 9;生物多样性保护的全局 Morans I 增加明显，由2000 年的 059 增加至 2010 年的 0731。图 3不同年份蒙辽农牧交错区生态系统服务局部 Morans I 指数通过局部 Morans I 分析，产水量在研究区北部出现低低聚集，在南部喀喇沁左翼蒙古族自治52第 1 期徐梦冉，等:人类活动影响下蒙辽农牧交错区的生态系统服务空间格局县出现高高聚集，在 2010 年新增了彰武县为高高聚集，相比 2000 年新增了赤峰市市辖区为低低聚集区，而克什克腾旗不再有低低聚集的情况。土壤保持 2000 年与 2010 年相比并无变化，在该区南部出现高高聚集，东北部出现低低聚集。碳固持聚集情况变化较大，高高聚集除了南部的三个城市之外还新增了蒙辽交界处的喀喇沁旗、宁城县和建平县，另外，科尔沁旗左翼后旗由原来的高高聚集变为低低聚集，低低聚集由原来的林西县、巴林左旗、巴林右旗和翁牛特旗变为东北部的扎鲁特旗、科尔沁左翼中旗、开鲁县、通辽市，而阿鲁科尔沁旗还是为低低聚集没有变化。生物多样性保护聚集区域有增加，相较 2000 年而言 2010 年低低聚集增加了扎鲁特旗和开鲁县，减少了赤峰市，高高聚集增加了凌源市和朝阳县，而蒙辽交界处的建平县出现了特殊的低高异常区。说明辽农牧交错区南部生态在恢复，北部在退化。33人类活动与生态系统服务空间格局分析人类活动强度同各生态系统服务间相关性见表 2。表 2不同年份蒙辽农牧交错区生态系统服务与人类足迹指数相关系数年份统计指标产水量土壤保持碳固持生物多样性保护2000 年相关系数 r002300140077*0052*P 值03780601000300442010 年相关系数 r0035001300090013P 值0176062007230620注:表示 P001，极显著相关;表示 P005，显著相关。结果表明，人类足迹指数与产水量和土壤保持在 2000 年2019 年没有显著相关性。碳固持与人类足迹指数相关性在 2000 年有极显著相关(相关系数为0077)，在 2010 年人类足迹指数同碳固持呈弱负相关，即碳固持量随人类活动强度的增加呈递减趋势，人类活动较强的区域存在大量耕地，且人类活动最强的区域存在城镇和工业用地，植被覆盖较少，碳固持低。生物多样性保护与人类足迹指数相关性在 2000 年呈现显著负相关，人类活动越强的区域对生物多样性保护越弱。人类活动较强的区域会进行城镇建设、耕地种植等造成生境破坏，一定程度上会导致生物多样性降低。4结论本文计算 2000 年2019 年逐年蒙辽农牧交错区的产水量、土壤保持、碳固持和生物多样性保护等 4 种生态系统服务，分析了各种生态系统服务的相关程度、以及人类活动强弱对生态系统服务空间格局的影响，得出以下结论:(1)2000 年2019 年蒙辽农牧交错区各项生态服务的空间分布相似，产水量和生物多样性保护服务整体呈现出“东南高、西北低”的分布格局，全区大部分碳固持量的空间分布较为平均，高值区主要分布在林地，研究结果对农牧交错区的减排增汇政策的制定以及生态系统碳库管理提供了一定的参考。(2)蒙辽农牧交错区的产水量、土壤保持、碳固持、生物多样性保护服务均呈现出不同程度的空间正向自相关性。2000 年的碳固持、生物多样性保护的生态系统服务量随人类活动的增强呈现显著递减趋势，因此协调人类活动与生态环境之间的关系，对蒙辽农牧交错区的生态保护具有重要意义。该研究是分析人类活动需求与生态系统服务供给之间反馈机制的关键，但缺乏生态系统服务空间异质性的归因分析，未来研究将着力于分析蒙辽农牧交错区各生态系统服务空间异质性的主要影响因素，为该区域制定生态保护政策提供依据。参考文献:1DAILY G C Natures services:societal dependence onnatural ecosystems M Washington DC:Island Press，1997:1373 2傅伯杰，吕一河，高光耀 中国主要陆地生态系统服务与生态安全研究的重要进展 J 自然杂志，2012，34(5):261272 3DAILY G C，POLASKY S，GOLDSTEIN J，et al Eco-system services in decision making:time to deliver J Frontiers in Ecology and the Environment，2009，7(1):2128 4LEMAITE D C，MILTON S J，JAMAIN C，et alLinking ecosystem services and water resources:land-scapescale hydrology of the Little KarooJ Frontiersin Ecology and the Environment，2007，5(5):261270 5FU B，FOSIUS M，LIU J Ecosystem services:climatechange and policy impacts J Current Opinion in Envi-ronmental Sustainability，2013，5(1):13 6KAEIVA P，WATTS S，MCDONALD，et al Domes-ticated nature:shaping landscapes and ecosystems for62大连民族大学学报第 25 卷human welfare J Science，2007，316:18661869 7陈红，江旭聪，任磊，等基于 USLE 模型的淮河流域土壤侵蚀定量评价J 土壤通报，2021，52(1):165 176 8谢高地，张彩霞，张雷明，等 基于单位面积价值当量因子的生态系统服务价值化方法改进 J 自然资源学报，2015，30(8):12431254 9张舟，吴次芳，谭荣 生态系统服务价值在土地利用变化研究中的应用:瓶颈和展望 J 应用生态学报，2013，24(2):556562 10王文萱，李明孝基于 DPSI 的湖南省土地生态安全时空变化J 生态学杂志，2020，39(8):27242736 11常虹，杨武，石磊，等 基于土地利用的北方农牧交错区碳固定和土壤保持时空变化 J 农业资源与环境学报，2021，38(3):484493 12肖建设，乔斌，陈国茜，等 黄河源区玛多县土地利用与生态系统服务价值演变研究J 生态学报，2020，40(2):510521 13梁潇洒，张立奇，闫庆忠，等 围栏封育恢复和提升辽西退化草地的碳固持功能J 中国草地学报，2019，41(2):6570 14郑华，欧阳志云，赵同谦，等 人类活动对生态系统服务功能的影响J 自然资源学报，2003，18(1):118126 15LEH M D K，MATLOCK M D，CUMMINGS E C，etal Quantifying and mapping multiple ecosystem serviceschange in West AfricaJ Agriculture Ecosystems Environment，2013，165(1751):618 16BAAL H，KEENAN J，SHAMA S K，et al Spa-tial assessment and mapping of biodiversity and conser-vation priorities in a heavily modified and fragmentedproduction landscape in northcentral Victoria，Australia J Ecological Indicators，2014，36(1):552562 17方露露，许德华，王伦澈，等 长江、黄河流域生态系统服务变化及权衡协同关系研究J 地理研究，2021，40(3):821838 18吴丹，邹长新，林乃峰，等 长江经济带生态系统服务权衡与协同关系研究J 环境生态学，2021，3(9):17 19赵哈林，赵学勇，张铜会，等 北方农牧交错带的地理界定及其生态问题J 地球科学进展，2002，17(5):739747 20史文娇，刘奕婷，石晓丽 四个历史时期中国北方农牧交错带界线数据集(19702010)J 全球变化数据学报(中英文)，2018，2(4):419425+537543 21DOGES P，ALLEN Estimating reference evapo-transpiration under inaccurate data conditions J Irri-gation and Drainage Systems，2002，16(1):3345 22WISCHMEIE W H，SMITH D D Predicting rainfallerosionlossesfromcroplandeastoftheockyMountains:guide for selection of practices for soil andwater conservation M Washington，D C:United StatesDepartment of Agriculture，1965:2644 23WILLIAMS J Epicerosion/productivity impact cal-culator:1 model documentation M Washington，DC:United States Department of Agriculture，1990，4(4):206207 24ALLEN G，PEEIA L S，AES D，et al Crop evap-otranspiration guidelines for computing crop water re-quirements FAO Irrigation and drainage paper 56 M ome:FAO Food and Agriculture Organization of theUnited Nations，1998:56 25杨文仙，李石华，彭双云，等 顾及地形起伏的 In-VEST 模型的生物多样性重要区识别 以云南省为例 J 应用生态学报，2021，32(12):43394348 26刘洋，李诚志，刘志辉，等 1982 2013 年基于GIMMSNDVI 的新疆植被覆盖时空变化 J 生态学报，2016，36(19):61986208 27连喜红，祁元，王宏伟，等 人类活动影响下的青海湖流域生态系统服务空间格局 J 冰川冻土，2019，41(5):12541263(责任编辑冉春秋)72第 1 期徐梦冉，等:人类活动影响下蒙辽农牧交错区的生态系统服务空间格局