中巴经济走廊巴基斯坦段产水服务供需风险评价及分区.pdf

第4 3卷第3期2 0 2 3年6月水土保持通报B u l l e t i no fS o i l a n dW a t e rC o n s e r v a t i o nV o l.4 3,N o.3J u n.,2 0 2 3 收稿日期:2 0 2 2-0 7-0 6 修回日期:2 0 2 2-1 0-0 2 资助项目:教育部规划基金项目“基于充分均衡发展视角下的彝汉交错深度贫困区生态稳定脱贫机制研究”(1 8 X J A 6 3 0 0 0 5)第一作者:钟孟君(1 9 9 7),女(汉族),四川省通江县人,硕士研究生,研究方向为土地生态评价。E m a i l:1 3 7 6 9 3 4 5 4 5q q.c o m。通信作者:蒋贵国(1 9 7 1),男(汉族),江苏省靖江县人,博士,教授,主要从事土地评价等方面研究。E m a i l:d l x y j g g s i c n u.e d u.c n。中巴经济走廊巴基斯坦段产水服务供需风险评价及分区钟孟君,蒋贵国,陆忆文,程珍旎,董威威(四川师范大学 地理与资源科学学院,西南土地资源评价与监测教育部重点实验室,四川 成都6 1 0 0 0 0)摘 要:目的研究中巴经济走廊巴基斯坦段产水服务供需状况,科学识别产水服务供需风险区并制定合理的风险防范措施,为中巴经济走廊实现绿色发展提供重要依据。方法以中巴经济走廊巴基斯坦段为研究区,基于2 0 0 0,2 0 0 5,2 0 1 0,2 0 1 5,2 0 2 0年遥感影像数据和社会经济数据,采用I n V E S T模型和以人定需法估算产水服务供给量与需求量,结合 供需 指数 及 趋势 评估 产 水服 务供 需风 险 并进 行风 险 划区。结果2 0 0 02 0 2 0年总产水量呈先增后减的趋势,产水服务供给能力中部高南北低;2 0a来产水需求总量增幅达5 8.4 4%。2 0 0 02 0 2 0产水服务供需指数变化趋势较弱,整体稳定。严重赤字区主要分布于印度河平原与俾路支省西部,高度盈余区集中分布于吉尔吉特巴尔蒂斯坦与俾路支省东部。产水服务供需高风险区面积占2 2.6 9%,其供需匹配类型以低供给低需求为主(4 7.5 7%)。结论动态评估中巴经济走廊巴基斯坦段产水量,制定区域水资源保护措施,构建生态补偿机制,使产水服务供需达到平衡,有助于促进中巴经济走廊的绿色发展。关键词:产水服务;供给与需求;供需风险;I n V E S T模型;中巴经济走廊文献标识码:A 文章编号:1 0 0 0-2 8 8 X(2 0 2 3)0 3-0 2 2 5-0 9 中图分类号:X 1 7 1.1,K 9 0 3文献参数:钟孟君,蒋贵国,陆忆文,等.中巴经济走廊巴基斯坦段产水服务供需风险评价及分区J.水土保持通报,2 0 2 3,4 3(3):2 2 5-2 3 3.D O I:1 0.1 3 9 6 1/j.c n k i.s t b c t b.2 0 2 3.0 3.0 2 7;Z h o n gM e n g j u n,J i a n gG u i g u o,L uY i w e n,e ta l.R i s ka s s e s s m e n ta n dz o n i n go ns u p p l ya n dd e m a n do fw a t e ry i e l ds e r v i c e si nP a k i s t a nS e c t i o no fC h i n a-P a k i s t a ne c o n o m i cc o r r i d o rJ.B u l l e t i no fS o i la n d W a t e rC o n s e r v a t i o n,2 0 2 3,4 3(3):2 2 5-2 3 3.R i s kA s s e s s m e n t a n dZ o n i n go nS u p p l ya n dD e m a n do fW a t e rY i e l dS e r v i c e s i nP a k i s t a nS e c t i o no fC h i n a-P a k i s t a nE c o n o m i cC o r r i d o rZ h o n gM e n g j u n,J i a n gG u i g u o,L uY i w e n,C h e n gZ h e n n i,D o n gW e i w e i(K e yL a b o r a t o r yo fL a n dR e s o u r c e sE v a l u a t i o na n dM o n i t o r i n gi nS o u t h w e s t,M i n i s t r yo fE d u c a t i o n,S c h o o l o fG e o g r a p h ya n dR e s o u r c e sS c i e n c e,S i c h u a nN o r m a lU n i v e r s i t y,C h e n g d u,S i c h u a n6 1 0 0 0 0,C h i n a)A b s t r a c t:O b j e c t i v eT h es u p p l ya n dd e m a n do fw a t e ry i e l ds e r v i c e si nt h eP a k i s t a ns e c t i o no fC h i n a-P a k i s t a ne c o n o m i cc o r r i d o rw e r ea n a l y z e d,t h er i s ka r e a s f o rw a t e ry i e l ds e r v i c e sw e r e i d e n t i f i e ds c i e n t i f i c a l l y,a n dt h er i s kp r e v e n t i o nm e a s u r e sw e r e f o r m u l a t e dr e a s o n a b l yi no r d e rt op r o v i d ea ni m p o r t a n tb a s i sf o rt h eg r e e nd e v e l o p m e n to f t h eC h i n a-P a k i s t a ne c o n o m i cc o r r i d o r.M e t h o d sB a s e do nr e m o t es e n s i n g i m a g ed a t aa n ds o c i o-e c o n o m i cd a t a f o r2 0 0 0,2 0 0 5,2 0 1 0,2 0 1 5,a n d2 0 2 0f r o mt h er e s e a r c ha r e a(t h eP a k i s t a ns e c t i o no f t h eC h i n a-P a k i s t a ne c o n o m i cc o r r i d o r),t h eI n V E S T m o d e l a n dt h eh u m a n-d e m a n dm e t h o dw e r eu s e dt oe s t i m a t e t h es u p p l ya n dd e m a n do fw a t e ry i e l ds e r v i c e s.T h es u p p l ya n dd e m a n d i n d e xa n dt r e n d sw e r eu s e dt oa s s e s s t h es u p p l ya n dd e m a n dr i s k so fw a t e ry i e l ds e r v i c e sa n dt od e l i n e a t er i s ka r e a s.R e s u l t s T o t a lw a t e ry i e l df r o m2 0 0 0t o2 0 2 0 i n i t i a l l y i n c r e a s e da n dt h e nd e c r e a s e d.T h ew a t e ry i e l ds e r v i c es u p p l yc a p a c i t yw a sh i g h i nt h ec e n t e ro f t h es t u d ya r e aa n dl o wi nt h en o r t ha n ds o u t h.T h et o t a ld e m a n df o rw a t e ry i e l di n c r e a s e db y5 8.4 4%o v e r t h ep a s t2 0y e a r s.F r o m2 0 0 0t o2 0 2 0,t h es u p p l ya n dd e m a n di n d e xo fw a t e ry i e l ds e r v i c e sh a daw e a k t r e n da n dw a sg e n e r a l l ys t a b l e.S e v e r ed e f i c i t a r e a sw e r em a i n l y l o c a t e d i n t h e I n d u sP l a i na n d W e s t e r nB a l o c h i s t a n,a n dh i g hs u r p l u sa r e a sw e r ec o n c e n t r a t e di nG i l g i t-B a l t i s t a na n dE a s t e r nB a l o c h i s t a n.T h eh i g h-r i s ka r e a sa c c o u n t e df o r2 2.6 9%o ft h en a t i o n a la r e a,a n dt h et y p eo fs u p p l ya n dd e m a n dm a t c hw a sm a i n l y l o ws u p p l ya n d l o wd e m a n d(4 7.5 7%).C o n c l u s i o nD y n a m i c a l l ya s s e s s i n gw a t e ry i e l di nt h e C h i n a-P a k i s t a ne c o n o m i cc o r r i d o r(P a k i s t a ns e c t i o n),f o r m u l a t i n gr e g i o n a lw a t e rr e s o u r c ep r o t e c t i o nm e a s u r e s,a n dc o n s t r u c t i n ga ne c o l o g i c a lc o m p e n s a t i o n m e c h a n i s mt ob a l a n c et h es u p p l ya n dd e m a n do fw a t e ry i e l ds e r v i c e sa r ec o n d u c i v et op r o m o t i n gt h eg r e e nd e v e l o p m e n to ft h eC h i n a-P a k i s t a ne c o n o m i cc o r r i d o r.K e y w o r d s:w a t e ry i e l ds e r v i c e;s u p p l ya n dd e m a n d;s u p p l y-d e m a n dr i s k;I n V E S T m o d e l;C h i n a-P a k i s t a ne c o n o m i cc o r r i d o r 产水服务作为生态系统服务的一项重要分支,对维持生态系统的稳定与和谐,促进区域经济的可持续发展具有极其重要的作用1。近年来随着全球生态环境的变化与人类活动的加剧,产水服务供需矛盾日益凸显。I P B E S发布的 全球生物多样性和生态系统服务评估报告 指出,包括产水服务在内的绝大多数生态系统服务功能迅速下降2。在此背景下,研究区域生态系统产水服务供需状况、识别产水服务供需风险区并制定合理的风险防范措施对促进区域可持续发展显得尤为重要。目前国内外学者大多采用模型模拟法研究区域产水量,评估产水服务供需风险。模型模拟法包括A R I E S模型和I n V E S T模型等,其中I n V E S T模型操作简便在国内外应用广泛。包玉斌等3以黄土高原为研究区,采用I n V E S T模型定量评估水源涵养功能,开展水源涵养的空间分区研究。吕乐婷等4利用I n V E S T模型的产水量模块对大连市产水量的时空分布特征及其成因进行研究。M a n s o o rD.K.L e h等5利用I n V E S T模型对西非加纳和科特迪瓦河流域产水量进行评估,取得良好效果。产水量评估方法的不断发展为产水服务生态风险的评估奠定了良好的基础6。现阶段众多学者对产水服务生态风险评估的研 究 更 多 考 虑 到 了 供 给 与 需 求 间 的 耦 合 关系7-8。B a r oF等9,B o i t h i a sL等1 0以欧洲的5个城市、西班牙E b r o盆地为研究区,将流域产水服务供需结合共同评估在空间上的匹配度来确定区域产水服务生态风险。王壮壮等1 1,L iD e l o n g等1 2分别对陕西省、京津冀地区产水服务进行研究,揭示产水服务供需风险时空变化格局、量化淡水生态系统服务流量对下游地区的影响。2 0 1 3年5月,中国国家领导人访问巴基斯坦时提出中巴经济走廊远景规划,得到国际社会的广泛关注1 3。由于中巴经济走廊地缘价值的重要性,该走廊也成 为“一 带 一 路”倡 议 的 样 板 工 程 和 旗 舰 项目1 4。然而一系列的生态环境问题制约着中巴经济走廊项目的发展,而且项目的建设也会对巴基斯坦脆弱的生态环境产生一定的负面影响1 5。A z i z等1 6研究发现巴基斯坦产水服务价值由于开放水域的减少等原因在2 0 1 5年和未来土地使用方案中急剧下降。为秉持“绿色发展”理念,在建设中巴经济走廊过程中就需充分考虑生态风险的问题。鉴于此,本研究选取2 0 0 0,2 0 0 5,2 0 1 0,2 0 1 5,2 0 2 0年5个时间节点,采用I n V E S T模型和以人定需法研究巴基斯坦产水服务供需状况,识别产水服务供需风险区并提出风险防范措施为优化区域生态系统管理提供科学依据。1 研究区概况中巴经济走廊覆盖中国新疆维吾尔自治区和巴基斯坦全境,北至新疆的喀什,南至巴基斯坦的西南港口城市瓜达尔港1 7。结合中巴经济走廊项目的实际分布以及已有研究,选取巴基斯坦全境(包括巴控自由克什米尔地区)作为研究区(图1)。注:本图地理底图源于全球行政区划数据库,并基于自然资源部标准地图服务系统G S(2 0 2 0)4 3 9 3号标准地图校准制作。下同。图1 中巴经济走廊巴基斯坦段地形图F i g.1 T o p o g r a p h i cm a po fC h i n a-P a k i s t a ne c o n o m i cc o r r i d o r(P a k i s t a ns e c t i o n)622 水土保持通报 第4 3卷巴基斯坦地处南亚次大陆西北部(2 3 3 7 N,6 0 8 0 E),总面积约8.8 01 05k m2。巴基斯坦位于亚热带,气候干旱、降水少,植被稀疏,森林面积仅占巴基斯坦总面积的4.8 0%。巴基斯坦设旁遮普省、开伯尔普什图省(简称开普省)、俾路支省、信德省、伊斯兰堡首都特区和联邦直辖部落区等6个省级行政区,此外还包括吉尔吉特巴尔蒂斯坦、自由克什米尔地区,各省区的自然环境和社会经济发展具有较大的差异。缺水、土壤侵蚀、森林的过度砍伐等是巴基斯坦目前面临的主要生态问题。2 数据来源与方法2.1 数据来源及预处理本研究所需数据来源如表1所示。表1 数据来源及预处理T a b l e1 D a t a s o u r c ea n dp r e p r o c e s s i n g数据名称 来 源 预处理 降水数据NA S A的T RMM(TMP A/3 B 4 3)数据集结合N D V I数据降尺度为1k m分辨率N D V I数据T e r r a卫星MO D I S传感器的植被指数数据MO D 1 3 A 3数据集通过MR T软件处理成1k m分辨率参考作物潜在蒸发量数据NA S A的MO D 1 6 A 3G F数据集通过MR T软件和A r c G I S软件处理成1k m分辨率植物有效水分含量世界土壤数据库(H a r m o n i z e dW o r l dS o i lD a t a b a s e,HWS D)由S P AW的S o i lW a t e rC h a r a c t e r i s t i c s模块计算土壤数据世界土壤数据库(H a r m o n i z e dW o r l dS o i lD a t a b a s e,HWS D)土地利用数据NA S A的数据集MC D 1 2 Q 1人口密度数据世界人口网(www.w o r l d p o p.o r g/)通过A r c G I S软件处理成1k m分辨率人均需水量数据联合国粮食及农业组织(h t t p:www.f a o.o r g/h o m e/e n)D EM数据S R TM数据集2.2 研究方法2.2.1 产水服务供给量化方法 I n V E S T模型产水模块以水量平衡原理为基础,产水量为降水量减去蒸散发量1 8,公式如下:Y(x)=1-A E T(x)P(x)P(x)(1)A E T(x)P(x)=1+P E T(x)P(x)-1+P E T(x)P(x)1(2)(x)=ZAWC(x)P(x)+1.2 5(3)P E T(x)=k c(l x)E T0(x)(4)式中:A E T(x)为栅格的年度实际蒸散发量;P(x)为栅格的年降水量;P E T(x)为潜在蒸散发量;(x)是表征自然气候土壤特性的非物理参数。AWC(x)为植物可用水分的体积(毫米)。E T0(x)为该位置生长的参考植被(例如草或苜蓿)的蒸散量。k c(l x)取决于该像素栅格上的土地覆被营养特征。Z是一个经验常数,与N(每年降水事件的数量)呈正相关。2.2.2 产水服务需求量化方法 每个像元的需水量是由人口密度与人均水需求量的乘积求得,人均水需求量包括工业用水量、农业用水量和生活用水量。人均水需求量一般由人均用水量表示,但由于2 0 0 0年后巴基斯坦面临着非常严重的缺水问题,人均用水量不可替代为人均需水量,因此,本研究将2 0 0 5,2 0 1 0,2 0 1 5,2 0 2 0年人均需水量统一为2 0 0 0年的人均需水量(12 2 3.1 0m3/人)。E Sd e m a n d,i=iDi(5)式中:E Sd e m a n d,i为像元栅格上的水资源需求量;i为像元栅格的人口密度(人/k m2);Di为人均水资源需求量。2.2.3 产水服务供需指数及趋势 本研究采用供需指数s u p p l y-d e m a n di n d e x(S D I)1 9-2 0来刻画生态系统服务供给与需求在空间上的矛盾和差异。S D I的取值范围为-1至1,当S D I 0时,供给量大于需求量,处于盈余状态;当S D I=0时,供给量等于需求量,处于供需平衡状态;当S D I 0,供需指数呈上升趋势,供需矛盾缓解;当k0,p0,p0.0 5)、不显著下降区(k0.0 5)、显著下降区(k0,p0.0 5)。拟合五期供需指数得出巴基斯坦产水服务供需指数趋势(图6),结果表明巴基斯坦产水服务供需指数变化趋势以不显著上升和下降为主(9 4.3 0%),显著上升区占2.2 6%,显著下降区占3.4 4%。显著上升区集中分布于吉尔吉特巴尔蒂斯坦和俾路支省北部。显著下降区集中分布于开普省、旁遮普省北部。不显著上升区广泛分布于俾路支省、旁遮普省、信德省等地。不显著下降区集中于巴基斯坦北部和西部,其他区域较分散。总体来看,巴基斯坦大部分地区变化趋势较弱,产水服务相对稳定。图6 2 0 0 02 0 2 0年中巴经济走廊巴基斯坦段产水供需指数空间变化趋势F i g.6 S p a t i a lc h a n g et r e n do fw a t e ry i e l ds e r v i c e ss u p p l ya n dd e m a n di n d e xo fC h i n a-P a k i s t a ne c o n n o m i cc o r r i d o r(P a k i s t a ns e c t i o n)f r o m2 0 0 0t o2 0 2 03.2 产水服务供需风险分区结合产水服务供需指数和趋势划分出产水服务供需风险(图7)。由表6可知,2 0 2 0年巴基斯坦产水服务高风险区面积比例为2 2.6 9%,主要分布于俾路支省西部、开普省中部、吉尔吉特巴尔蒂斯坦西部、东部省区的印度河平原地区。低风险区面积占比为3 6.1 3%,集中分布于各省份的中部和边界地区。无风险区面积占比为4 1.1 9%,分布于吉尔吉特巴尔蒂斯坦大部分地区、联邦直辖部落区南部、俾路支省东部、旁遮普省东南部和信德省东部。高风险区土地利用类型主要为农耕区、城市建成区、荒漠区等,农耕区和城市建成区需大量水源,产水服务供不应求,风险等级 高。荒 漠 区 降 水 稀 少,产 水 服 务 供 需 矛 盾突出。图7 中巴经济走廊巴基斯坦段产水服务供需风险F i g.7 S u p p l ya n dd e m a n dr i s k so fw a t e ry i e l ds e r v i c e so fC h i n a-P a k i s t a ne c o n n o m i cc o r r i d o r(P a k i s t a ns e c t i o n)将高风险区基于Z-s c o r e标准化后得出2 0 2 0年产水服务高风险区供需匹配空间类型(图8)与分布(表7)。2 0 2 0年巴基斯坦产水服务高风险区供需匹配类型以低供给低需求为主,其比例为4 7.5 6%,主要分布于俾路支省西部和吉尔吉特巴尔蒂斯坦,这些地区虽人口密度小,需求量也低,但产水供给满足不了需求导致供需矛盾恶化。高供给高需求区(占2 5.1 9%)分布于开普省南部和旁遮普省北部,这些地区产水供给虽高,但人口密集,水需求量远远超过供给量。低供给高需求区(占2 6.7 4%)集中于印度河中下游平原,这些地区降水稀少且人口密度大,水需求量大,供需矛盾严重。132第3期 钟孟君等:中巴经济走廊巴基斯坦段产水服务供需风险评价及分区表6 2 0 2 0年中巴经济走廊巴基斯坦段产水服务供需风险等级比例T a b l e6 P r o p o r t i o no f r i s kl e v e l so f s u p p l ya n dd e m a n df o rw a t e ry i e l ds e r v i c e so fC h i n a-P a k i s t a ne c o n n o m i cc o r r i d o r(P a k i s t a ns e c t i o n)i n2 0 2 0风险划区风险等级比例/%主要分布区域无风险区富 足2 3.9 8吉尔吉特省巴尔蒂斯坦大部分地区、联邦直辖部落区中部南部、俾路支省北部南部、旁遮普省东南部、信德省东部安 全1 7.2 1低风险区脆 弱3 1.4 9各省份中部和省边界周围地区供给不足4.6 4危 险4.3 7高风险区中度危险5.4 7俾路支省西部、开普省中部、吉尔吉特巴尔蒂斯坦中部、旁遮普省东部和信德省中部的印度河平原地区高度危险1 2.8 5表7 2 0 2 0年中巴经济走廊巴基斯坦段产水服务高风险区供需匹配空间类型分布T a b l e7 S p a t i a l d i s t r i b u t i o no f s u p p l ya n dd e m a n dm a t c h i n g i nh i g h-r i s ka r e a so fw a t e ry i e l ds e r v i c e so fC h i n a-P a k i s t a ne c o n n o m i cc o r r i d o r(P a k i s t a ns e c t i o n)i n2 0 2 0供需匹配空间类型面积比例/%分布地区 低供给高需求2 6.7 4印度河中下游平原高供给高需求2 5.1 9开普省南部和旁遮普北部低供给低需求4 7.5 6俾路支省西部和吉尔吉特巴尔蒂斯坦以及印度河平原两侧高供给低需求0.5 0 图8 2 0 2 0年中巴经济走廊巴基斯坦段产水服务高风险区的供需匹配空间类型F i g.8 T y p e so f s u p p l ya n dd e m a n dm a t c h i n gs p a c e s i nh i g h-r i s ka r e a s o fw a t e ry i e l d s e r v i c e s o fC h i n a-P a k i s t a ne c o n n o m i cc o r r i d o r(P a k i s t a ns e c t i o n)i n2 0 2 04 讨 论产水服务作为生态系统服务的一项重要分支,对维持生态系统的稳定与和谐、促进区域经济的可持续发展具有极其重要的作用。产水量的变化受降水量、植被分布状况和土地利用/覆被等因素的影响。近年来多位学者对降水量、土地利用/覆被与区域生态系统产水量之间的关系进行了大量研究。杨洁等2 5在对黄河流域产水量的时空格局研究中发现降水对产水量的影响比较显著。李怡颖等2 6通过研究发现影响张家口市产水量的主要因素是降水,其次为土地利用方式。谢余初等2 7在白龙江流域的研究也得到相似结论。由于不同土地利用类型的土壤含水量、蒸散发能力、凋落物持水能力及冠层截留量均存在差异,因此不同土地利用类型产水能力不同2 5。本研究的5个时间节点中2 0 1 0年降水量最大,产水量也最高;同时本研究区内产水量较高的地区地类以耕地、林地为主,植被覆盖率高。另外产水服务需求量的变化与人口密度高度相关,随着人口的增长,水需求量也逐年增加;在空间格局上,研究区内水需求量高值区通常为人口密集区。本研究结合供需指数与最小二乘法拟合的供需指数趋势动态评估了巴基斯坦产水服务供需风险,高风险区生态脆弱、生态敏感性高,在未来发展中要密切关注风险区内项目的开发建设,监测产水服务动态变化,防止其进一步退化;加强水源地生物多样性保护,依托水土保持等工程措施进行生态保护;推行节水农业,采取跨流域调水等方式弥补用水需求;借鉴中国的生态补偿机制,结合供给与需求确定补偿主体、补偿受体、补偿金额,使产水服务供需相对平衡。由于数据的可获得性以及研究水平的限制等原因,本研究还存在一些局限和不足。例如未考虑未来土地利用、气候变化以及区域间产水服务流动等,因此在未来研究中可从这些方面进行完善,将多情景模拟纳入产水服务供需风险研究,考虑生态系统服务的流动性,将其量化以便动态识别产水服务在空间上的供需匹配。232 水土保持通报 第4 3卷5 结 论(1)在 产 水 服 务 供 需 时 空 变 化 方 面,2 0 0 02 0 2 0年总产水量呈先增后减的趋势,产水服务供给能力 中 部 高 南 北 低;2 0a来 水 需 求 量 增 幅 高 达5 8.4 4%,水需求量高值区集中分布于人口密集的城市建成区和耕地区。(2)在产水服务供需指数方面,供需指数趋势以不显著上升和下降为主(9 4.3 0%),总体来看巴基斯坦大部分地区变化趋势较弱,产水服务相对稳定;严重赤字区主要分布于印度河平原与俾路支省西部,高度盈余区集中分布于吉尔吉特巴尔蒂斯坦与俾路支省东部。(3)在产水服务供需风险方面,2 0 2 0年产水服务供需高风险区面积比例为2 2.6 9%,低风险区面积比例为3 6.1 3%,无风险区面积比例为4 1.1 9%,其中高风险区供需匹配类型以低供给低需求为主(4 7.5 6%),其次为 低 供 给 高 需 求(2 6.7 4%)、高 供 给 高 需 求(2 5.1 9%),高供给低需求(0.5 0%)。参考文献1 戴尔阜,王亚慧.横断山区产水服务空间异质性及归因分析J.地理学报,2 0 2 0,7 5(3):6 0 7-6 1 9.2 吴会慧,赵艳玲,李素萃.基于G E E的淮南市土地利用变化对生态系统服务价值的影响J.生态经济,2 0 2 1,3 7(2):1 4 6-1 5 2.3 包玉斌,李婷,柳辉,等.基于I n V E S T模型的陕北黄土高原水源涵养功能时空变化J.地理研究,2 0 1 6,3 5(4):6 6 4-6 7 6.4 吕乐婷,任甜甜,李赛赛,等.基于I n V E S T模型的大连市产水量时空变化分析J.水土保持通报,2 0 1 9,3 9(4):1 4 4-1 5 0.5 L e h M D K,M a t l o c k M D,C u mm i n g sE C,e ta l.Q u a n t i f y i n ga n d m a p p i n g m u l t i p l ee c o s y s t e m s e r v i c e sc h a n g e i n W e s tA f r i c aJ.A g r i c u l t u r e,E c o s y s t e m s&E n v i r o n m e n t,2 0 1 3,1 6 5:6-1 8.6 W e iH e j i e,F a nW e i g u o,W a n gX u e c h a o,e ta l.I n t e g r a t i n gs u p p l ya n ds o c i a l d e m a n d i ne c o s y s t e ms e r v i c e s a s s e s s m e n t:Ar e v i e wJ.E c o s y s t e mS e r v i c e s,2 0 1 7,2 5:1 5-2 7.7 白杨,王敏,李晖,等.生态系统服务供给与需求的理论与管理方法J.生态学报,2 0 1 7,3 7(1 7):5 8 4 6-5 8 5 2.8 赵文武,刘月,冯强,等.人地系统耦合框架下的生态系统服务J.地理科学进展,2 0 1 8,3 7(1):1 3 9-1 5 1.9 B a r F,H a a s eD,G m e z-B a g g e t h u nE,e ta l.M i s m a t-c h e sb e t w e e ne c o s y s t e ms e r v i c e ss u p p l ya n dd e m a n di nu r b a na r e a s:Aq u a n t i t a t i v ea s s e s s m e n t i nf i v eE u r o p e a nC i t i e sJ.E c o l o g i c a l I n d i c a t o r s,2 0 1 5,5 5:1 4 6-1 5 8.1 0 B o i t h i a sL,A c u aV,V e r g o sL,e ta l.A s s e s s m e n to f t h ew a t e r s u p p l y:D e m a n dr a t i o s i naM e d i t e r r a n e a nB a s i nu n d e rd i f f e r e n tg l o b a l c h a n g es c e n a r i o sa n dm i t i-g a t i o na l t e r n a t i v e sJ.S c i e n c eo ft h eT o t a lE n v i r o n-m e n t,2 0 1 4,4 7 0/4 7 1:5 6 7-5 7 7.1 1 王壮壮,张立伟,李旭谱,等.区域生态系统服务供需风险时空演变特征:以陕西省产水服务为例J.生态学报,2 0 2 0,4 0(6):1 8 8 7-1 9 0 0.1 2 L iD e l o n g,W uS h u y a o,L i uL a i b a o,e ta l.E v a l u a t i n gr e g i o n a lw a t e r s e c u r i t yt h r o u g ha f r e s h w a t e r e c o s y s t e ms e r v i c ef l o w m o d e l:Ac a s es t u d yi nB e i j i n g-T i a n j i a n-H e b e i r e g i o n,C h i n aJ.E c o l o g i c a l I n d i c a t o r s,2 0 1 7,8 1:1 5 9-1 7 0.1 3 玛依努尔,赵军.中巴农产品贸易增长特质及影响因素研究:基于贸易引力模型的实证分析J.价格月刊,2 0 1 9(1 1):3 2-3 7.1 4 楚新元.中巴经济走廊研究的热点、趋势与动态演进:基于C i t e S p a c e的文 献 计量 分析J.金 融发 展 评 论,2 0 2 0(9):2 6-4 1.1 5 M u h a mm a dS a l i m.“中巴经济走廊”开启中巴合作新模式J.中国投资,2 0 1 5(7):5 4-5 7.1 6 A z i zT.C h a n g e si nl a n du s ea n de c o s y s t e m s e r v i c e sv a l u e si n P a k i s t a n,1 9 5 02 0 5 0J.E n v i r o n m e n t a lD e v e l o p m e n t,2 0 2 1,3 7:1 0 0 5 7 6.1 7 田丰.中巴经济走廊新愿景J.中国外资,2 0 1 8(3):2 8-3 0.1 8 Z h a n gL i w e i,F uB o j i e,L Y i h e,e t a l.B a l a n c i n gm u l-t i p l ee c o s y s t e ms e r v i c e s i nc o