黄河流域生态系统服务与城市化时空交互作用分析——以黄河宁夏段为例.pdf

第54卷第8 期2023年8 月文章编号：10 0 1-417 9(2 0 2 3)0 8-0 0 9 3-0 8引用本文：肖芳，魏文颖，周斌.黄河流域生态系统服务与城市化时空交互作用分析：以黄河宁夏段为例J.人民长江，2 0 2 3,54(8):93-100.人民长江YangtzeRiverVol.54,No.8Aug.，2 0 2 3黄河流域生态系统服务与城市化时空交互作用分析以黄河宁夏段为例肖芳,魏文颖,周周斌？（1.江西科技学院，江西南昌330 0 9 8；2.河海大学水利水电学院，江苏南京2 10 0 9 8）摘要：分析黄河流域生态系统服务与城市化的关系有助于精准化制定流域可持续发展政策。以黄河流域宁夏段为例，量化了198 0 2 0 2 0 年粮食产量、产水量、土壤保持、碳固存和生境质量等关键生态系统服务，分析了城市化变化特征，并结合双变量空间自相关分析模型解析了生态系统服务与城市化之间的时空交互作用。结果表明：198 0 2 0 2 0 年，宁夏粮食产量持续增加，土壤保持、产水量、NPP和生境质量呈增加减少-增加的变化特征。宁夏城市化水平有所提升，呈现南北高、中部低的空间分布格局。1990 2 0 2 0 年，生态系统服务与城市化之间在全局尺度存在着空间正相关向空间负相关转变的特征，并且在局部尺度上以低一低集聚和高一低集聚为主。研究结果可为黄河流域生态保护和城市发展规划提供实践指导。关键词：生态系统服务；城市化；时空变化；交互作用；黄河流域中图法分类号：X171;X826;F291.10引言作为维持人类生存的基础，生态系统服务与人类福密切相关。然而，随着城市化的快速推进，生态系统服务的供给面临着日益加剧的城市化压力，已经威胁到了区域生态系统的可持续发展。在此背景下，理解生态系统服务与城市化间的关系对于有效的生态系统管理至关重要。目前，生态系统服务与城市化之间关系的研究得到了生态学、地理学等领域学者的广泛关注。众多学者在城市化产物对生态系统服务的影响、城市化与生态系统服务的协调关系、生态系统服务对城市化的响应特征及机制等2-6 方面展开了一系列研究,并取得了丰硕的成果。如欧阳晓等7 在洞庭湖的案例中揭示了建设用地扩张对生态系统服务的负向影响;渠玉莲等8 对南京市的研究证实了该地区城市化发展与生态系统服务之间存在着高度协调关系；XIAO等9 在沪杭都市圈采用耦合协调度模型从系统收稿日期：2 0 2 3-0 2-13基金项目：国家自然科学基金项目（52 0 0 90 32）；2 0 19年度南昌市社会科学规划一般项目（YJ201907）作者简介：肖芳，女，讲师，硕士，主要从事公共管理和乡村振兴方面的研究。E-mail：XF2 0 0 2 9 9 0 54 16 3.c o m文献标志码：A1D0I:10.16232/ki.1001-4179.2023.08.013层面和子系统层面揭示了城市化与生态系统服务之间的交互胁迫关系。纵观现有研究发现，在研究尺度上，既包括较小尺度的市域、县域研究，中观尺度的省域研究，也有涉及全国和区域的大尺度研究，如针对各县域10】、各市域】、各省域 等的案例研究明确了行政尺度上两者之间的关系。但由于栅格尺度上城市化相关数据的获取受限，目前对于栅格尺度上生态系统服务与城市化关系的分析较少，这将限制精细化的生态管理实践。此外，在案例区的选择上，研究案例多发生在中国东部发达的城市、城市群及沿海地区，对黄河流域的研究较少,研究成果不足以支撑黄河流域的可持续发展。宁夏回族自治区作为唯一一个全境属于黄河流域的省域,随着“一带一路”战略和黄河流域生态保护和高质量发展战略的实施，城市化水平拥有很大的发展潜力，然而该地区生态系统服务与城市化之间关系的时空动态信息仍然不得而知，这给未来如何协调城市94建设与生态保护带来了挑战。为此，本研究以宁夏为例，在从栅格尺度上解析该地区近40 a来生态系统服务和城市化动态特征的基础上，结合双变量空间自相关分析模型从全局和局部尺度上进一步揭示两者之间的时空交互作用，以期为协调宁夏生态保护和城市化发展的政策制定提供科学依据。1数据与方法1.1石研究区概况宁夏位于黄河流域上游（352 6 3930 N、10 45010740E），面积约6 6 40 0 km，海拔10 8 0 3474m，地势西南高、东北低（见图1）。宁夏属于典型的大陆性季风气候，年平均气温为5.8 8.2，空间上由南向北逐渐升高；年平均降水量为18 0 8 0 0mm,由南向北逐渐减少12 。该地区以草地和耕地为主，随着社会经济的发展和城市化的推进，建设用地呈现明显扩张趋势。与此同时，宁夏面临着包括水资源短缺、水土流失、土地荒漠化、土壤盐渍化、草原退化等一系列问题，对当地可持续发展构成了严重的威胁。因此，有必要分析城市化背景下该地区生态系统服务的变化特征及其与城市化之间的关系,这将增强对“过去当前”状况的了解，并助力“未来”发展政策的制定。黄河宁夏黄河流域0300jkm1980年土地利用类型耕地林地草地建设用地水体未利用地050km1.2石研究数据来源本研究以19 8 0，19 9 0,2 0 0 0,2 0 10 年和2 0 2 0 年为典型研究期，收集整理相关数据，主要包括：由中国人民长江科学院资源环境科学与数据中心（http：/w w w.r e s d a t a.c m a.c n/）提供的气象数据；由中国科学院资源环境科学与数据中心（ht-tp：/w w w.r e s d c.c n/）和国家青藏高原科学数据中心(http：/w e s t d c.w e s t g i s.a c.c n/）提供的NDVI数据；来自国家冰川冻土沙漠科学数据中心（http：/w w w w.r e s d c.c n）的GDP和人口空间分布数据，并采用文献13 提供的方法进行纠正。本研究将所有数据重采样为1km的空间分辨率,进而在1 km1km的栅格尺度上展开相关研究工作。1.3石研究方法1.3.1生生态系统服务量化方法本研究参考黄土高原、西北干旱半干旱区以及研究区的相关研究成果，选择粮食产量、产水量、土壤保持、碳固存和生境质量这5种研究区关键的生态系统服务。相关计算方法如下。（1）粮食产量。结合前人研究14，根据耕地的空间分布，采用NDVI对粮食产量进行空间化，计算公式为G=海拔：m0503474Lkm10802000年土地2020年土地利用类型利用类型耕地耕地林地草地水体建设用地未利用地050km图1研究区概况Fig.1Overview of the study area2023年NDVI,GsumNDVL.um式中：G，为栅格i的粮食产量；Gsum为行政区的粮食产量；NDVI为栅格i的NDVI;NDVIsum为行政区耕地的NDVI总和。（2）产水量。采用InVEST模型量化宁夏的产水林地草地水体建设用地未利用地050km(1)量，模型公式为14WY,=(1-AET,/P,)P;式中：WY,为栅格i的年产水量；AET，为栅格i的年实际蒸散发量；P，为栅格i的年均降水量。（3）土壤保持量。采用RUSLE模型对宁夏的土壤保持服务进行量化，计算公式为15Asc=A,-A.=RKLS(1-C P)(3)式中：A.为土壤保持量;A，为潜在土壤侵蚀量；A。为实际的土壤侵蚀量;R为降雨侵蚀力因子;K为土壤可蚀性因子;LS为坡度坡长因子;C为植被覆盖和管理因子；P为水土保持措施因子。（4）碳固存。以植被净初级生产力（NPP）来表征研究区的碳固存服务。采用CASA模型模拟宁夏的(2)第8 期NPP,模型所涉及的公式为16 NPP,=APAR,8i式中：NPP,为栅格i在t月份的NPP;APAR,为栅格i在t月份所吸收的光合有效辐射；8 i为栅格i在t月份的实际光能利用率。（5）生境质量。采用InVEST模型对宁夏生境质量进行量化，模型涉及公式为17 Q:=Hm (1-式中：Q；为栅格i的生境质量;Hm为生境类型m的生境适宜性;Dm为生境类型m中栅格i的生境退化度;z为模型默认参数；k为半饱和系数。（6）生态系统服务综合指数。由于不同生态系统服务的计量单位不同，参考文献18 ，对5种生态系统服务进行标准化并累加，得到生态系统服务综合指数以表征宁夏生态系统服务的整体水平。具体计算公式为X=maxX,-minX,ESSCI.=X式中：ESSCI为栅格i的生态系统服务综合指数；X为第j项生态系统服务在栅格i上的标准化值;X,为第j项生态系统服务在栅格i上的原始值；maxX,和minX,分别为第j项生态系统服务的最大值和最小值。1.3.2城市化量化城市化一般体现在经济发展、人口增长、建设用地扩张和生活水平提高这4个方面19。由于难以获取空间化的生活水平相关数据，参考他人研究成果3.19-2 0 ，本研究选择经济密度、人口密度和建设用地比例这3个指标来表征城市化。具体而言，对这3个城市化指标采用公式（6）进行标准化，并采用公式（8)进行累加求取平均值得到城市化指数，以反映宁夏的城市化水平。计算公式为Uic+Uip+UicU,=3式中：U,为栅格i的城市化指数；Uic、U i p、U i c 分别为栅格i上标准化后的经济密度、人口密度和建设用地比例。受数据可获性的限制，本研究计算了19 9 0 年以来的城市化指数。1.3.3时空交互作用分析本研究采用双变量空间自相关分析模型2 1 揭示研究期内宁夏生态系统服务与城市化之间交互作用的时空演变特征。具体计算公式为肖芳，等：黄河流域生态系统服务与城市化时空交互作用分析一以黄河宁夏段为例DimDim+k)X,;-minX,95(4)2Zw.(.-)(y.-5)a=I=式中:I为双变量全局MoransI,I0表示存在空间正相关；n为研究区栅格数量；Wa为空间权重矩阵；xa、y 分别为自变量、因变量在栅格、b 上的观测值;s为所有样本的方差。(5)I.=Z.ZWaZ,式中:l。为栅格 的自变量与因变量的局部空间关系；Z。、Z。分别为自变量、因变量在栅格、b 上观测值的方差标准化值。基于I。可形成4种LISA集聚模式2 2 高高集聚，即高生态系统服务综合指数被高城市化指数包围；低低集聚，即低生态系统服务综合指数被低城市化指数包围；低一高集聚，即低生态系统服务综合指数被高城市化指数包围；高一低集聚，即高生态系统服务综合指数被低城市化指(6)数包围。2结果与分析(7)2.1生态系统服务的时空特征图2 显示了宁夏粮食产量、产水量、土壤保持、NPP以及生境质量这5种生态系统服务在19 8 0 2020年的时空变化特征。从时间变化上来看：宁夏的粮食产量自198 0 2 0 2 0 年呈现持续增加趋势，由1980年的6 9.55t/km持续增加至2 0 2 0 年的2 2 0.47t/km,增加了2 倍多，这得益于农业设施发展和现代技术进步所带来的粮食单产能力的提升。宁夏的土壤保持呈现出198 0 1990 年增加、1990 2 0 10 年减少、2 0 10 2 0 2 0 年增加的变化特征，尽管一系列水土保持措施的实施改善了当地的土壤侵蚀状况，但降雨的年际波动也给土壤保持服务的变化带来了很大的不确定性。宁夏的产水量、NPP和生境质量表现出相对一致的时间变化趋势，即呈现出198 0 1990 年增(8)加、1990 2 0 0 0 年减少、2 0 0 0 年以后持续增加的变化特征，这与当地的人类活动密切相关，尤其是2 0 0 0 年以来退耕还林（草)工程的大规模实施，显著改善了当地的生态系统状况。从空间分布上来看，宁夏的粮食产量、产水量、土壤保持、NPP以及生境质量这5种生态系统服务在均存在着明显的空间异质性。其中：粮食产量的高值区主要集中在宁夏北部的引黄灌区,在中南部地区高值呈现零碎分布特征,这与宁夏耕地的空间格局一致。研究期内产水量在空间上大致呈现自南向北递减(9)台台(10)96人民长江2023年1980年粮食产量/tkm330.26050100Ikm1980年产水量/mm391.98050100km1980年土壤保持量/（tkm227164050100km1980年NPP/gCma629.2350.79050100Jkm1980年生境质量1N1990年粮食产量/t:km586.99050100JkmN1990年产水量/mm575.89050100kmN1990年土壤保持量/tkm268628050100JkmN1990年NPP/(gC.m739.1676.53050100LJkmN1990年人生境质量2000年粮食产量/t:km759.28050100Lkm2000年产水量/mm368.920050100JkmN2000年土壤保持量/tkm2194720050100LkmN2000年NPP/(gC m717.830501001kmN2000年生境质量N2010年粮食产量/tkm848.99050100LkmN2010年产水量/mm460.720050100kmN2010年土壤保持量/tkm3561510050100LJkmN2010年NPP/(gCm867.7323.61050100LukmN2010年生境质量1N2020年粮食产量/t:km692.59050100LkmN2020年产水量/mm422.44050100ukmN2020年土壤保持量/tkm4722310050100kmN2020年NPP/(gCm1031.88050100kmN2020年人生境质量1240年平均值220200180160140120100806019801990200020102020年份N50年平均值4540353025201510519801990200020102020年份40年平均值353025201510519801990200020102020年份N3807年平均值34026022018019801990200020102020年份N0.45r年平均值000000.350的空间变异格局,这与当地降雨的空间变化特征高度一致。土壤保持在空间上的分布特征与当地的地形状况密切相关（对比图1），低值区主要位于低海拔地区，这里地势相对平坦，侵蚀潜力小，进而导致较低50100km050100ikmFig.2Spatio-temporal changes of ecosystem services in Ningxia from 1980 to 2020050100Jkm图219802020年宁夏生态系统服务的时空变化050100km的土壤保持；相反，土壤保持的高值区主要位于研究区的高海拔地区，包括宁夏中部丘陵沟壑区和南部山区。NPP呈现“两凸一凹”的空间变化格局，宁夏中部地区的NPP值差于南部和北部地区，这与中部地区干050100km19801990200020102020年份第8 期旱、风大、土地贫等生态基底条件有关。生境质量在空间分布上比较零碎，其空间格局的形成与当地的地形条件和土地利用结构有关，生境质量好的地区主要位于高海拔地区，这里人类活动强度低，主要以林草地为主；而在海拔较低的区域，适宜人类居住和耕种，人类对自然生态系统的干扰频繁,进而导致生境质量较差。2.2城市化的时空特征图3呈现了宁夏城市化在1990 2 0 2 0 年的时空变化特征。从时间变化上来说，宁夏的城市化指数呈现持续增加趋势，这符合中国城市化进程的发展特征。具体来看，宁夏在1990 2 0 0 0 年间城市化指数增长较慢，由19 9 0 年的1.2 5%增长至2 0 0 0 年的1.37%；2000年以后,随着西部大开发战略的实施，引进了大量的人才、资金和技术，宁夏的城市化水平得到了显著提升，城市化指数大幅增加至2 0 2 0 年的3.2 7%，增加了1.39倍，这也与宁夏2 0 0 2 年制定并坚持的“重点突出、梯度推进”城市化发展方针有关。城市化指数/%城市化指数/%56.0563.8200050100km(a)1990年城市化指数/%65.530.05050100km(c)2010年图31990 2 0 2 0 年宁夏城市化的时空变化Fig.3Spatio-temporal changes of urbanization in Ningxiafrom 1990 to 2020从空间变化上来说，宁夏城市化发展水平在空间上总体呈现出南北高、中部低的空间分布格局。其中宁夏北部地区得黄河灌溉之利，是宁夏人类生产生活适宜的区域，也是宁夏社会经济发展水平最高的区域，肖芳，等：黄河流域生态系统服务与城市化时空交互作用分析一以黄河宁夏段为例务的改善也会限制城市化的发展速度。表11990 2 0 2 0 年宁夏生态系统服务与城市化之间的全局MoransI和LISA集聚统计Tab.1 Clobal Morans I and LISA agglomeration statistics betweenecosystem services and urbanization in Ningxia from 1990 to 2020全局年份P高-高高-低低一高低一低050100Lkm(b)2000年城市化指数/%1000.04050100Jkm(d)2020年97得天独厚的发展环境促成了该区域拥有相对较高的城市化水平；相反,在“苦甲天下”的宁夏中部地区，自然条件恶劣，生态本底脆弱，水资源贫乏，社会经济发展水平低，很大程度上限制了城市化进程的推进；而在宁夏南部地区,尽管以山地为主，但降雨丰富,植被茂盛，一定程度上为该区域城市化的发展提供了良好的自然生态支撑，进而导致其城市化水平相对高于中部地区。2.3生态系统服务与城市化的时空交互作用19902 0 2 0 年宁夏生态系统服务与城市化之间时空交互作用的全局MoransI如表1所列。可以发现，在2 0 0 0 年之前，全局MoransI大于0（P=0.01），且由19 9 0 年的0.18 3减小至2 0 0 0 年的0.16 5，这说明在此时期宁夏生态系统服务与城市化之间存在着空间正相关性，即城市化的发展在一定程度上会促进生态系统服务的改善，这可能与当时城市化发展速度较慢有关。随着2 0 0 0 年以后城市化进程的快速推进，两者之间的全局MoransI持续减少至2 0 2 0 年的-0.19 8（P=0.0 1），生态系统服务与城市化之间的关系也由空间正相关逆转为空间负相关，此时城市化的发展是以损害生态系统服务为代价的，反过来说，生态系统服Morans I集聚19900.183N200020100.12720200.1980.01图4的LISA集聚图进一步揭示了宁夏生态系统服务与城市化在局部尺度上的空间集聚特征。由图可知,高一高集聚和低高集聚主要集中在研究区的北部地区，低一低集聚广泛分布在研究区的中部地区，而高一低集聚主要分布在研究区的中南部地区。结合表1分析可以发现，宁夏主要以低低集聚和高低集聚为主,这反映了该地区较低的城市化水平。4种集聚模式随时间变化来看，高高集聚和低低集聚的面积占比在研究期内呈现减少趋势，即高高集聚和低低集聚的面积在缩小，而高低集聚和低高集聚的面积在扩大,表明两者之间的空间正相关性存在着向空间负相关性演变的倾向，这种倾向揭示了宁夏生态系统服务与城市化之间关系的恶化趋势，不利于两者协调可持续发展。面积占比/%集聚集聚0.017.1570.1650.010.01集聚不显著7.1223.4706.5347.6444.1259.8863.35016.18523.85658.3963.48623.47258.8643.85920.51361.6175.14619.55255.76798高-高高-低低-高低-低不显著050100Jkm(a)1990年高-高高-低低-高低-低不显著050100km(c)2010年图41990 2 0 2 0 年宁夏生态系统服务与城市化的LISA集聚图Fig.4 LISA cluster map of ecosystem services and urbanizationin Ningxiafrom1990to20203讨论本研究发现2 1世纪以来宁夏生态系统服务与城市化之间存在着负的空间相关性，即一个地点的城市化发展可能会导致周围地区生态系统服务供给能力的降低，这支持了生态系统服务与城市化之间权衡的观点-2 0.2 1。可能的原因是环境要素（如二氧化碳、废水废气）的变化，很可能会通过大气运动、水文运动等自然过程传递到周围地区，进而导致周围地区生态系统服务供给能力的变化3。从生态系统服务与城市化的交互作用机制来看，一方面，人口增加和经济发展会加大资源的消耗强度，导致大量的废水、废气和废渣等流入生态系统，进而降低了生态系统的承载力2 4；另一方面，城市化进程中大量生态用地向建设用地转变，对生态系统格局、过程和生物地球化学循环产生了强烈的影响，降低了生态系统的供给能力4.7.17 。总体来说，城市化可以通过经济发展、人口增加和建设用地扩张等来影响生态系统的结构和功能4,2 5，进而影响生态系统服务的供给能力。但是，生态系统也可以通过资源赋、环境条件和政策干预来限制城市化发人民长江展 。例如,随着城市化的发展，人类生活水平显著高-高改善，人类逐渐从“求生存”转向“求生态”，这种“诉高-低低-高低-低不显著50100Jkm(b)2000年高-高高-低低-高低-低不显著50100JkmL(d)2020年2023年求”就会很大程度上限制建设用地无序大规模侵占生态用地等行为2 4。那么，是城市建设还是生态保护，对于许多决策者来说至关重要，因为他们面临着如何在不削弱生态系统服务可持续供给能力的情况下开展城市建设活动的难题。尽管生态系统服务与城市化之间整体呈空间负相关，但在局部尺度上并不总是如此，因为研究结果还发现了显著的高一高和低低集聚模式，这也反映了生态系统服务与城市化之间存在着复杂的交互作用机制。也就是说，如果不考虑城市化与生态系统服务之间的交互关系，一些诸如直接将生态系统服务高值区划分为生态保护区的策略可能就会失效,因为城市建设用地的高经济产出能力往往会迫使一些建设活动进人到这些保护区中2 1。鉴于此，有必要将生态系统服务与城市化联系在一起，以探讨其在局部尺度上的交互作用，进而为城市发展和生态保护提出更合理的建议。本研究确定了生态系统服务与城市化在局部尺度上存在着4种集聚模式（即高高、低低、高低、低高集聚），在特定区域采取的措施应该要符合该区域的集聚模式。在生态优先原则下，可以考虑将高一高集聚的区域划分为生态功能区。在这些地区，应严格控制灌区发展规模，可适当鼓励生态友好型的土地利用方式（如林地、草地甚至人工绿地等），以加强对当前生态系统服务高供给能力的维护。高低集聚的区域属于生态优势区，应采取严格的保护措施以避免城市化发展需求下人类活动对生态系统的干扰，城市建设活动（如道路建设、旅游开发等）应被严格地规划以保护区域生态系统的稳定性和可持续性。在低一高集聚的区域应该注意缓解人口压力，集约利用现有建设用地，控制城市扩张规模，合理确定城市发展边界，并可以加快绿地、公园等绿色基础设施建设，以提升生态系统供给能力。在低一低集聚的区域可以将生态建设项目（如土地整治、植树造林等）作为城市化的配套措施，在保证区域生态安全的前提下，适当进行建设活动。此外，研究方法的局限性主要有两点：本研究根据前人相关成果采用了宁夏5种关键的生态系统服务，可能难以充分表征区域的生态系统服务状况；同样，尽管经济密度、人口密度和建设用地比例是被广泛用以反映我国城市化的变量2 0 1,但由于城市发展的多样性和复杂性,可能也存在着难以充分反映区域城市化水平的问题。因此，未来研究中可以考虑增加生态系统服务类型和城市化表征指标。尽管栅格尺度第8 期上的研究结果可以将生态系统服务与城市化之间的关系精确到每个栅格上，对于协调两者关系措施的精准制定具有重要作用，但空间研究中总是存在的空间尺度效应可能会导致当前尺度上的交互作用结果与其他尺度上有所不同16.2 6 1。因此，未来研究中可以开展多尺度分析予以验证。4结论通过对栅格尺度上黄河流域宁夏段生态系统服务和城市化时空动态的分析，进一步揭示了该地区生态系统服务与城市化的时空交互作用,形成的主要结论为：198 0 2 0 2 0 年，宁夏粮食产量持续增加，土壤保持、产水量、NPP和生境质量呈增加减少-增加的变化特征,空间上均存在着明显的空间异质性。19902 0 2 0 年，宁夏城市化指数呈现持续增加趋势，在空间上总体呈现南北高、中部低的空间分布格局。19902 0 2 0 年，宁夏生态系统服务与城市化之间整体存在着从空间正相关向空间负相关转变的特征，并在局部尺度上集聚，表现为高高和低高集聚在宁夏北部地区、低低集聚在宁夏中部地区、高低集聚在宁夏中南部地区。参考文献：1黄敏，杨飞，郑士伟.中国城镇化进程对生态系统服务价值的影响J.水保持研究,2 0 19,2 6（1）：352-359.2杨宜男，李敬，王立，等.长三角地区城市化对典型生态系统服务供需的影响J.自然资源学报，2 0 2 2,37（6）：1555-157 1.3 YANG M,GAO X,SIDDIQUE K H M,et al.Spatiotemporal explora-tion of ecosystem service,urbanization,and their interactive coercingrelationship in the Yellow River Basin over the past 40 yearsJ.Sci-ence of the Total Environment,2022,858:159757.4秦方.河南新郑市“三生”用地转型及其生态服务价值响应J.人民长江,2 0 2 0,51(2)：58-6 3.5 黄卓男，陈颖彪，吴志峰.珠三角城市群地表热环境与生态系统服务价值解耦的时空特征J：应用生态学报，2 0 2 2,33（7）：1993-2000.6 石宇，韩蕊，郭泺.基于地形梯度的南方丘陵山地带生态系统服务价值对城市化的空间响应：以粤北地区为例J.生态学报,2 0 2 1，41(18):7238 7248.7欧阳晓，王坤，魏晓.城乡建设用地关联对生态系统服务的影响：以洞庭湖地区为例J.生态学报,2 0 2 2，42（2 1）：8 7 138 7 2 2.8 渠玉莲，郭宗逵.南京市城市化与土地生态系统服务耦合度分析J.安徽农业大学学报,2 0 16,43（4）：57 6-58 1.9XIAO R,LIN M,FEI X,et al.Exploring the interactive coercing rela-tionship between urbanization and ecosystem service value in the肖芳，等：黄河流域生态系统服务与城市化时空交互作用分析一以黄河宁夏段为例velopment,2021,32(14):4105-4120.22叶璇，康帅直，赵永华，等.陕北黄土高原植被恢复与生态系统服务的时空关系J.应用生态学报,2 0 2 2,33（10）：2 7 6 0-2 7 6 8.23齐静，邓伟，周渝，等.三峡库区生态系统服务价值时空演变及驱动力J.人民长江,2 0 2 0,51(1)：113119.24ZHANG D,HUANG Q,HE C,et al.Impacts of urban expansion onecosystem services in the Beijing-Tianjin-Hebei urban agglomera-tion,China:A scenario analysis based on the Shared SocioeconomicPathwaysJ.Resources Conservation&Recycling,2017,125:115-130.25 1DADASHPOOR H,AZIZI P,MOGHADASI M.Land use change,ur-banization,and change in landscape pattern in a metropolitan areaJ.Science of the Total Environment,2019,655:707-719.26LE PROVOST G,SCHENK N,PENONE C,et al.The supply of mul-tiple ecosystem services requires biodiversity across spatial scalesJ.Nature Ecology&Evolution,2022,7:236-249.（编辑：黄文晋）99Shanghai-Hangzhou Bay Metropolitan RegionJ.Journal of CleanerProduction,2020,253:119803.10刘海龙，丁娅楠，王跃飞，等.山西省城镇化与生态系统服务时空相关性及空间效应分析J.水土保持学报,2 0 2 2,36（1）：12 4-134.【11朱文博，李双成，朱连奇.中国省域生态系统服务足迹流动及其影响因素J.地理研究,2 0 19,38（2）：337-347.12 XU J,WANG S,XIAO Y,et al.Mapping the spatiotemporal heteroge-neity of ecosystem service relationships and bundles in Ningxia,Chi-naJ.Journal of Cleaner Production,2021,294:126216.13 王壮壮，张立伟，李旭谱，等.区域生态系统服务供需风险时空演变特征：以陕西省产水服务为例J.生态学报，2 0 2 0，40（6）：18871900.14 吴树荣，潘换换，姬倩倩，等.基于生态系统服务的山西黄河流域保护优先区识别J.生态学报，2 0 2 2，42（2 0）：8 12 6-8 137.15 杨丽雯，王大勇，王勇智，等.水河流域土壤保持服务供需关系量化评估J.资源科学，2 0 2 0,42（12）：2 451-2 46 2.16郭保祥，聂玉明，张瑞刚，等，山西省生态系统服务多尺度特征及其驱动因素分析J.人民长江，2 0 2 2,53（5）：7 5-8 1.17曾志伟，杨华，宁启蒙，等洞庭湖区土地利用强度演变及其对生态系统服务的影响J.经济地理,2 0 2 2,42(9：17 6-18 5.18张琨，吕一河，傅伯杰，等.黄土高原植被覆盖变化对生态系统服务影响及其阅值J.地理学报，2 0 2 0,7 5（5）：949-96 0.19PENG J,TIAN L,LIU Y,et al.Ecosystem services response to urban-ization in metropolitan areas:Thresholds identification J.Science ofthe Total Environment,2017,607-608:706-714.20陈立恒，龚健，王雄，等.长江中游地区生态系统服务平衡与城镇化的空间关系J.生态学报，2 0 2 1,41（14）：5546-5556.21SHI L,HALIK U,MAMAT Z,et al.Spatiotemporal investigation ofthe interactive coercing relationship between urbanization and ecosys-tem services in arid northwestern ChinaJ.Land Degradation&De-100(1.Jiangxi University of Technology,Nanchang 330098,China;2.College of Water Conservancy and Hydropower,Hohai Univer-sity,Nanjing 210098,China)Abstract:The analysis on the relationship between ecosystem services and urbanization in the Yellow River basin is of greatsignificance for accurate formulation of sustainable development policies in this basin.Taking the Yellow River reach in Ningxia asan example,we quantified the grain yield,water yield,soil conservation,carbon sequestration and habitat quality from 1980 to2020,analyzed the change characteristics of urbanization,and revealed the spatio-temporal interaction between ecosystem serv-ices and urbanization with the bivariate spatial auto-correlation analysis model.The results showed that the grain yield in Ningxiacontinued to increase,and soil conservation,water yield,NPP and habitat quality showed the characteristics of increase-decrease-increase during the study period.The urbanization level of Ningxia had improved,showing a spatial distribution pattern of highin the north and south and low in the middle.From 1990 to 2020,there was a shift from positive spatial correlation to negative spa-tial correlation between ecosystem services and urbanization at the global scale,and low-low cluster and high-low cluster weredominant at the local scale.The research results can provide practical guidance for ecological protection and urban project devel-opment of Yellow River basin.Key words:ecosystem services;urbanization;spatio-temporal change;interaction;Yellow River basin人民长江Spatio-temporal interaction between ecosystem services and urbanizationalong the Yellow River reach in NingxiaXIAO Fang,WEI Wenying,ZHOU Bin?2023年(上接第8 4页）Research progress and hot spot analysis of aeration technology in river ecological r