退耕还林还草工程对干旱区内陆河流域生态系统服务的影响

探讨人类活动尤其是大型工程、政策对生态系统服务的影响是国际上地理学和生态学研究的前沿与热点领域。利用土地利用数据、气象数据、遥感数据等, 结合InVEST产水量模型、CASA模型和RUSLE模型, 对石羊河流域退耕还林还草工程实施期间的生态系统产水、净初级产品生产和土壤保持服务进行精准测度, 分析其演变特征, 并进一步细化对比了工程实施区与非工程实施区生态系统服务变化的差异。结果表明: (1)石羊河流域退耕还林还草工程实施期间(2000—2015年), 林、草地面积有所增加, 耕地和未利用地面积有所减少; 面积变化比例最大的土地转化类型为“未利用地—草地”、“未利用地—林地”、“耕地—草地”和“耕地—林地”。(2)石羊河流域在退耕还林还草工程实施期间整体生态系统产水服务和土壤保持服务有不同程度衰退, 但是NPP有较为明显增加; 各服务在空间分布上具有一致性, 均呈现出南高北低, 由西南向东北阶梯递减的趋势。(3)石羊河流域“退耕区”产水、NPP和土壤保持总量和平均量均呈现出明显的增加趋势, 而“非退耕区”除了NPP总量和平均量有所增加以外, 其余两类服务均有不同程度的减少; 各退耕模式中, “荒地造林”区对各类服务的提升最为明显, “退耕还林”区产水总量和平均量均有所减少。     

滦河流域景观格局变化对水沙过程的影响

以滦河流域为研究区域,基于SWAT模型模拟1976—2012年滦河流域的水沙过程,分析2000年京津风沙源治理项目实施前后流域产水产沙时空格局变化;研究1980—2010年流域景观格局变化特征,揭示景观格局变化的水沙响应;应用Spearman相关分析法分析流域景观格局变化对水沙过程的影响。结果表明:与20世纪80年代相比,2010年流域林地和建设用地增加,其他用地类型减小;流域景观集中程度提高、连通性变优、优势斑块显著、形状趋于规则、多样性减少、破碎化程度降低、景观类型向非均衡方向发展;流域年均地表径流减少9mm,产水量增加5.44mm,产沙减小1.59t/hm~2;地表径流减少区域占全流域89.32%,产水量增加区域占76.71%,产沙量减少区域占93.89%;地表径流、产水、产沙与林地面积呈负相关,产水与草地面积呈正相关,地表径流、产沙与农业用地面积呈正相关;地表径流、产水、产沙与景观形状、Shannon′s均匀度、景观分离度呈正相关,与蔓延度、最大斑块指数呈负相关;产水、产沙与斑块密度和Shannon′s多样性指数呈正相关;工程治理后,流域年均径流量与产沙量显著下降,产水产沙高值区显著缩小,产沙关键区域仍需治理。     

清江流域生态系统服务权衡时空效应及驱动因素

识别流域生态系统服务权衡协同关系的时空分异特征及其影响因素,是生态系统管理和调控的重要基础,对高效配置环境资源和合理制定生态环境政策具有重要意义。本研究运用相关性分析、均方根偏差等方法分析清江流域2000—2020年粮食供给、净初级生产力(NPP)、土壤保持和产水服务间的权衡协同关系,并采用地理探测器分析影响生态系统服务权衡的决定因素。结果表明：2000—2020年,清江流域粮食供给服务呈下降趋势,NPP、土壤保持和产水服务呈递增态势;粮食供给与土壤保持服务、NPP与产水服务间权衡强度呈下降趋势,其他服务间权衡强度均呈增加趋势。粮食供给与NPP、土壤保持、产水服务间呈现东北权衡、西南协同的格局;NPP与土壤保持、产水服务在中部为协同关系、四周为权衡关系;土壤保持与产水服务表现为高度协同。土地利用和归一化植被指数为粮食供给与其他生态系统服务间权衡强度的主导因素,降水、气温和高程为产水服务与其他生态系统服务间权衡强度的主导因素,生态系统服务间权衡强度并非受单一因素的影响,两种服务间的相互作用或共同影响因素的驱动是其决定因素。研究结果可为国土空间生态修复规划策略的制定提供参考。     

基于供需关系的西南喀斯特区生态系统服务空间流动研究

生态系统服务的持续供给是社会和自然可持续发展的基础,人类通过对生态系统服务的消费来满足需求和提高福祉。科学理解生态系统服务从产生、传递到使用的全过程,明确区域生态系统服务供给与需求的平衡状况,对于实现区域可持续发展与提高人类福祉具有重要意义。基于RUSLE、InVEST、CASA等多种模型和方法,量化西南喀斯特区2000-2015年土壤保持服务、产水和碳固定服务的供给量和需求量,探讨供需空间盈余变化特征,确定服务流传输路径和流量。结果表明：（1）研究区土壤保持、产水和碳固定服务供给量和需求量均呈增加趋势,供给量分别增长679.52 t/hm²、2.2×10⁵ m³/km²、72.91 t/km²,需求量分别增加298.6 t/hm²、0.04×10⁵ m³/km²、168.36 t/km²。（2）整个研究区上三种生态系统服务供给尚能满足需求,土壤保持和碳固定服务供需比均呈减少趋势,而产水服务的供需比呈增加趋势。（3）土壤保持高服务流主要分布在云南的黑水河流域、四川的雅砻江流域、广东的北江和东江流域;广西的红水河流域在2000年为产水服务需求区,而到2015年转变为高服务流;碳高服务流主要分布在云南的西南部,方向由西南至东北。     

大宁河流域生态系统服务权衡关系梯度效应

深入理解生态系统服务权衡关系是推动生态系统有效管理和实现区域高质量发展的基础。然而,目前涉及生态系统服务权衡关系梯度效应的研究仍相对较少。选取三峡库区大宁河流域为研究区域,对典型平水年份2018年粮食生产、碳固持、土壤保持和产水量4种生态系统服务进行了量化评估和空间制图,探讨了海拔、降雨量、植被覆盖度3种环境因子梯度下生态系统服务的空间分异规律,采用均方根偏差进一步分析了生态系统服务权衡关系梯度效应。研究结果表明：粮食生产在空间上呈现中部、南部高,北部低的趋势,碳固持由西北向东南递减,土壤保持和产水量高值区均分布在西部。海拔、降雨量和植被覆盖度显著影响了生态系统服务,碳固持、土壤保持总体上随着环境梯度的增加而增加,而粮食生产随着3种环境因子梯度的增加呈下降的趋势。从生态系统服务权衡方面来看,碳固持与产水量为弱权衡关系,土壤保持与粮食生产为中度权衡,生态系统服务之间多为强权衡关系。由于多种环境梯度的相互作用,生态系统服务权衡关系随海拔梯度变化无明显规律,植被覆盖度增加显著加强了生态系统服务权衡关系,而降雨量增加显著减弱了生态系统服务权衡关系。因此,今后以流域生态系统服务协调和功能整体提升为目标的生态修复工程中,应更加注重统筹不同环境梯度区域内生态系统服务权衡关系。     

基于地理探测器的喀斯特山区生态系统服务关系分异特征及驱动力解析——以贵州省为例

生态系统服务关系的空间异质性及其驱动机制分析一直都是喀斯特生态系统服务研究的关键问题。以典型喀斯特山区--贵州省为例,选取植被净初级生产力(NPP),土壤保持(SC)和产水(WY)3类生态系统服务,在对其进行准确估算的基础上,采用相关分析法、均方根误差法从多角度剖析生态系统服务关系变化及其空间分异规律,借助地理探测器对生态系统服务关系变化的驱动机理进行解析。结果表明：①2000-2015年贵州省NPP、土壤保持、产水服务均值均呈增长趋势,其中,NPP和土壤保持服务在研究区东北部及西北角增长最显著,产水服务在研究区东南部增加最明显。②研究时段内,区域生态系统服务关系以协同为主,只是协同关系在减弱并有向权衡发展的趋势。③不同背景条件下区域生态系统服务关系具有明显的分异特征。相较于白云岩,NPP与土壤保持的权衡度在石灰岩上更大;土壤保持与NPP的权衡度在地势平缓的岩溶盆地上较大,产水服务与NPP、土壤保持间的权衡度在陡峭的岩溶槽谷(峡谷)上更大;产水服务与土壤保持间的权衡度随着植被覆盖度的增加而逐渐增大。④地理探测器结果表明：造林总面积是影响NPP与土壤保持间权衡关系变化的主导因子,降水总量已成为制约NPP与产水服务间协同的最主要因子,高程和人均国内生产总值(GDP)是影响土壤保持与产水服务间关系的最主要因子;因子间的交互影响均高于单因子的影响。研究深化了喀斯特地区生态系统服务权衡/协同关系变化的机理解析研究,为喀斯特生态系统有效保护和管理提供了科学参考。     

新疆生态系统服务时空变化及权衡协同关系

基于CASA模型、RUSLE模型和InVEST模型定量测算了新疆1982-2015年的净初级生产力(NPP)、土壤保持和产水三种生态系统服务物质量。使用空间分析和热点分析方法对生态系统服务的分区差异进行了刻画,并使用Sen趋势度方法分析了生态系统服务的时空变化趋势,进一步基于逐像元偏相关分析,对三种生态系统服务之间的权衡与协同关系及时空变化进行了分析。结果表明:1982-2015年,新疆NPP年际变化较为稳定,呈现逐渐增加的趋势,土壤保持和产水量的年际变化波动较大;NPP与土壤保持、产水之间为协同关系;土壤保持与产水之间呈现明显的协同关系,且主要是强协同关系;NPP与土壤保持的协同关系在全省范围内有弱化趋势,NPP与产水之间由权衡关系过度到协同关系,土壤保持和产水协同关系先强化后弱化。     

新疆生态系统服务时空变化及权衡协同关系

基于CASA模型、RUSLE模型和InVEST模型定量测算了新疆1982—2015年的净初级生产力(NPP)、土壤保持和产水三种生态系统服务物质量。使用空间分析和热点分析方法对生态系统服务的分区差异进行了刻画,并使用Sen趋势度方法分析了生态系统服务的时空变化趋势,进一步基于逐像元偏相关分析,对三种生态系统服务之间的权衡与协同关系及时空变化进行了分析。结果表明:1982—2015年,新疆NPP年际变化较为稳定,呈现逐渐增加的趋势,土壤保持和产水量的年际变化波动较大; NPP与土壤保持、产水之间为协同关系;土壤保持与产水之间呈现明显的协同关系,且主要是强协同关系; NPP与土壤保持的协同关系在全省范围内有弱化趋势,NPP与产水之间由权衡关系过度到协同关系,土壤保持和产水协同关系先强化后弱化。     

土地利用变化对生态系统服务的影响——以张家口-承德地区为例

生态系统服务与人类福祉密切相关,人类活动引起的土地利用变化是生态系统服务变化的主要驱动力之一。运用InVEST模型估算张家口-承德地区产水量和土壤保持量,并分析不同海拔和坡度条件下,2011—2015年土地利用以及产水和土壤保持服务的变化,研究结果表明:建设用地和乔木林地分别为面积增加量和减少量最多的用地类型,陡坡地段建设用地增加的同时水体减少;灌木林地是产水和土壤保持能力均最高的生态用地类型;5年间产水量和土壤保持量分别增加了32.16×10~6 m~3和17.69×10~6 t;产水量在海拔900 m地段,海拔越高增加量越少,土壤保持量在较低海拔有所增加;坡度越陡产水增加量越少,土壤保持量在各坡度均有增加。基于此研究结果,有助于选择优先建设的用地类型,使当地在保护源头水系和维持生态功能的同时发展冰雪运动。     

森林固碳释氧服务价值与异养呼吸损失量评估

固碳释氧是森林最重要的生态系统服务之一,将森林碳收支与固碳释氧服务价值评估相结合对于准确评估生态系统服务价值具有重要意义。应用森林碳收支模型(CBM-CFS3),分别基于净初级生产力(NPP)和净生态系统生产力(NEP)评估了2009—2030年湖北省兴山县森林生态系统总、净固碳释氧服务价值的时空动态,量化了异养呼吸造成的固碳释氧服务价值损失。模拟期间兴山县森林生态系统NPP逐渐增加(0.46—0.70 Tg/a),NEP由0.12 Tg/a先增加至0.21 Tg/a,然后逐渐下降至0.18 Tg/a;所对应的森林总、净固碳释氧服务价值范围分别为7.59—11.53亿元/a和2.21—3.70亿元/a。异养呼吸逐年增加,导致固碳释氧价值每年损失平均值为7.29亿元/a或4509元hm~(-2) a~(-1),约占总价值的68.6%。兴山县东南部异养呼吸造成的森林固碳释氧服务价值损失较高,而中部及西南部森林净固碳释氧价值较高。模拟期间兴山县森林为碳汇,稳定地提供固碳释氧服务。与NPP相比,使用NEP评估固碳释氧服务价值更为合理。忽视异养呼吸将严重高估森林生态系统固碳释氧服务价值;因而必须将物质循环过程与生态系统服务评估相结合,以降低评估结果的不确定性、提高生态系统服务的评估能力。     

Spatial interactions among ecosystem services and the identification of win-win areas at the regional scale

Although the quantification and evaluation of ecosystem services has been discussed on the Loess Plateau of China, few studies have specifically focused on quantifying the ecosystem services of trade-off and finding win-win areas. Based on previous studies about the evaluation of ecosystem services, we proposed a comprehensive spatial analysis method based on the InVEST and CASA models to analyse the spatial interaction between the four ecosystem services of water yield, sediment retention, NPP and habitat quality in the study area. The results showed that soil retention, NPP, and habitat quality gradually increased from 2000 to 2010 as water yield increased. Relationships among sediment retention, NPP and habitat quality were synergistic, but trade-off occurred between those three ecosystem services and water yield. Hotspots accounted for 3.27% of the basin, and these areas were mainly located in the downstream area of the basin and coincided with the forest and grassland areas. Coldspots accounted for 3.64% of the basin and were mainly in line with the upstream area and the urban road area in the basin. Win-win areas only accounted for 1.6% of the basin. The existence of win-win zones was positively correlated with soil silt content and was negatively correlated with distance to grassland, distance to forest, soil erodibility and slope. This research may be helpful for managing water resources and the sustainable development of the basin on the Loess Plateau of China.     

藏北草地退化与生态服务功能价值损失评估——以申扎县为例

藏北草地是我国重要的畜牧业生产基地和生态安全屏障。基于遥感和GIS技术,利用1990、2000、2010年3期不同时相的TM、ETM+和CEBERS遥感影像,对申扎县草地资源退化状况进行了遥感监测,并采用生态经济学评估模型对草地生态系统8个方面的服务功能价值损失进行了评估测算。结果表明:1990—2010年,申扎县草地退化面积增加了47.40×10~4hm~2,生态系统服务功能价值损失高达5.20×10~8元;其中1990—2000年,草地退化较严重,该时段也是生态系统服务功能价值损失较多的时期;2000—2010年,草地退化趋势变缓。藏北草地提供生物量价值仅约占生态系统服务功能总价值的7.0%,草地生态服务功能远大于其提供的生物量价值,因此必须从生态服务功能的的理念出发去经营草地,从而实现草地的可持续发展。     

京津冀地区生态系统服务时空变化及驱动因素

京津冀地区地处蒙古高原向华北平原过渡地带,是华北平原关键的生态屏障,在区域生态安全格局上具有重要的地位。结合地面观测、遥感监测、模型模拟和GIS空间分析技术,构建了一套京津冀地区生态系统关键服务的生态本底图谱,分析和评估了京津冀地区生态系统结构与防风固沙、水源涵养、土壤保持等关键服务的时空变化特征,并对气候变化和生态工程等驱动因素进行了讨论。结果表明:(1)京津冀地区生态系统类型以农田为主,占全区面积的49.7%。2000—2015年,聚落面积明显增加,农田面积大量减少,森林、草地、水体与湿地面积微弱减小。(2)多年平均生态系统防风固沙量、水源涵养量和土壤保持量分别为5.61×10~8 t、74.58×10~8 m~3、7.98×10~8 t,燕山-太行山森林生态区的生态系统关键服务量最高,西北草原生态区其次,京津唐和华北农业生态区较低。(3)2000—2015年,防风固沙量、水源涵养量均呈不显著上升趋势,增幅分别0.11 t hm~(-2) a~(-1)、0.03×10~4 m~3 km~(-2) a~(-1);土壤保持量呈显著上升趋势,增幅为1.08 t hm~(-2) a~(-1)。森林生态区各生态系统服务量均呈上升趋势(部分地区P0.05);草原生态区的防风固沙量呈不显著下降趋势,其他服务量整体变化幅度较小;京津唐农业生态区和华北农业生态区东部的防风固沙量呈极显著下降趋势,土壤保持量整呈不显著上升趋势。(4)研究区生长季气候趋于暖湿,叠加实施了多项重点生态修复工程,显著提高了区域植被覆盖度,进而提升了生态系统服务量。而在城市化发展过程中,建设用地不断扩张,导致部分农业生态区植被覆盖度下降,生态系统服务量有所降低。     

用RS和GIS技术评价福建省长汀县土壤保持功能对生态系统变化的响应

人类活动直接影响生态系统的变化,而生态系统变化又改变着其土壤保持功能。通过对生态系统及其土壤保持量变化的对比分析可以发现两者间响应机制,分析水土保持工作的得失,并进而为下一步治理提供指导。基于e Cognition平台利用landsat影像解释出长汀县2000、2010和2016年的生态系统类型,并利用RULSE模型定量评估土壤保持功能,探讨土壤保持量相对于生态系统变化的响应规律。结果表明:2000—2016年间,长汀县生态系统变化主要表现为两大特征:其一是低覆盖类型向高覆盖类型转变,表现为稀疏林大量地转为常绿阔叶林和常绿针叶林,研究期内稀疏林减少14729.38hm~2(减幅56.75%),而常绿阔叶林和常绿针叶林分别增加6791.50hm~2(增幅23.21%)、10228.07hm~2(增幅5.80%);其二是受经济利益驱动,大量耕地和常绿针叶林转为城镇和茶果园,研究期内城镇、茶果园分别增加4469.14hm~2(增幅85.43%)、1949.32hm~2(增幅69.38%)。研究期内全县土壤保持量显著提高,单位面积土壤保持量从2000年217.32t hm~(-2)a~(-1)提高2016年246.02t hm~(-2)a~(-1)。各生态系统中,常绿阔叶灌丛和常绿阔叶林土壤功能最强,单位面积土壤保持量大于277t hm~(-2)a~(-1),而城镇土壤保护功能最弱,单位面积土壤保持量低于65t hm~(-2)a~(-1);研究期内大部分生态系统的土壤保持能力有不同程度提高,只有城镇和茶果园下降,单位面积土壤保持量分别减小21.03t hm~(-2)a~(-1)和16.74t hm~(-2)a~(-1),这说明城镇和茶果园迅速扩张可能引发新的水土流失。从空间上看,土壤保持量提高的区域占全县96.5%,其中东北部森林区和中南部增幅最显著,说明了十几年的水土流失综合治理工作极具成效;另一方面土壤保持量下降的区域散布于人类活动较频繁的各个山间盆谷,其中汀州、大同最为集中,迅速的城镇化与茶果园开发是造成保土功能下降的主要原因。     

生态系统服务权衡与协同关系及驱动力——以江苏省太湖流域为例

生态系统服务权衡与协同关系研究是生态系统综合管理的前提。以太湖流域江苏省为例,通过空间显示的生物物理模型计算氮、磷净化、水量供给及土壤保持4种服务指标,借助GIS空间分析表征氮、磷指标与其他指标之间关系,并利用多元Logistic模型定量识别主导驱动力。结果表明:2000到2010年间,各指标的单位面积年均值呈现不同程度的变化,氮输出指标先增加后略下降,磷输出指标逐渐增加,水量供给先下降后上升,土壤保持逐渐下降,并且服务指标的空间格局显著差异。氮、磷净化与水量供给关系在空间上表现为广泛分布的权衡及协同变化区,但与土壤保持的关系不明显,氮、磷之间主要为协同变化关系。氮净化与水量供给的正向主导驱动力为城镇建设用地和农村居民点密度,而植被覆盖度和水网密度具有显著负作用。氮、磷关系的主要影响因素为植被覆盖度,其次是耕地及林地比例,且均起到正向促进作用。主导驱动力识别有助于明确生态系统服务间关系的作用机制,为区域环境管理及生态保护规划提供科学依据。     

黄土高原生态工程对关键生态系统服务时空变化的影响——以延河流域为例

黄土高原地区生态工程的实施,使其生态环境得到显著改善,提高了生态系统服务能力。现有的研究成果中,生态工程对生态系统服务影响的定量评价比较匮乏。以延河流域为例,采用RUSLE模型和InVEST模型对生态系统服务进行评价,通过构建模型识别出不同时期的生态工程区,探究生态工程对生态系统服务的定量影响,为下一步生态工程的实施提供科学的指导和依据。结果表明：(1)2000—2018年,延河流域土壤保持服务和产水服务均呈波动增长趋势,但二者的变化并不同步。(2)4个时期内,能产生生态效应的生态工程区面积呈波动增长趋势,从854 km²增加到1343 km²,在整个流域均有分布,不同时期的重点分布区不同。(3)生态工程增强了区域的土壤保持能力和产水能力,土壤保持服务增强区面积从477.5 km²增加到1140.6 km²,保持的土壤总量从2.1×10⁷ t增加到5.6×10⁷ t;产水服务增强区面积从139.1 km²增加到485.5 km     

南京市景观时空动态变化及其驱动力

选择南京市1995、2000和2008年的土地利用/土地覆盖数据,对其土地利用斑块动态变化与景观斑块的空间稳定性进行探讨。结果表明,全市总体上以耕地、林地和草地的减少,以及建设用地和水域的大幅度增加为典型特征。1995-2000年间耕地、林地和草地面积净减少了12017.2 hm²,2000-2008年间净减少41029.13 hm²;全市的建设用地面积1995-2000年间仅增加了11981.3 hm², 2000-2008年间面积增加了41027.06 hm²。从斑块类型的面积稳定性来看,各种土地利用斑块的空间稳定性都非常高,除了2000-2008年间未利用土地的自身空间保持率为50.66%外,其他土地景观类型的保持率都在90%以上。从不稳定斑块的空间分布看,1995-2000年间,主要集中在南京市的栖霞区、浦口区和江宁区,呈离散分布特征; 2000-2008年间,不稳定斑块一方面呈现出集中连片的分布态势,这在江宁区和栖霞区最为典型,另一方面这种集中分布趋势也扩大到了六合区与高淳县。从不稳定斑块的转化方向看,均以建设用地所占比例最大,两个时段内转为建设用地的不稳定斑块分别占到了不稳定斑块总面积的87.09%和68.46%,其次的转化方向为水域,两个时间段内分别达到了10.61%和23.86%。从景观斑块变化的原因来看,对当地景观斑块稳定性影响最大的主要是人口与经济发展、城市规划的引导和重大的生态建设工程等人为因素。     

不同视角下基于InVEST模型的流域生态补偿标准核算——以渭河甘肃段为例

补偿标准是生态补偿项目研究的核心内容。为提高补偿标准的可操作性,以渭河干流甘肃段为研究区,基于流域生态补偿标准范围界定,在考虑耕地向多种林地转化情况下采用In VEST模型和加权法核算流域新增生态服务量,并以新增生态服务价值为补偿上限,农户机会成本为补偿下限,先从国家新一轮退耕还林政策导向视角,探讨15°—25°和25°以上坡耕地退耕还林的差异化补偿标准。再从利益公平分配视角,耦合机会成本投入与生态系统服务产出,探讨生态补偿净收益不同地区的差异化补偿标准。结果表明:(1)渭河干流甘肃段退耕还林,生态补偿标准范围为146.39—481.98万元km~(-2)a~(-1)。(2)重要水源地一级保护区和保留区迎水面15°—25°坡耕地全部退耕还林需补偿金额13.95—33.55亿元。25°以上坡耕地全部退耕还林需要补偿金额12.96—28.35亿元。(3)各县/区比较,退耕还林净收益麦积区最大,陇西县最小。其生态补偿范围分别为202.23—727.92、96.99—136.53万元km~(-2)a~(-1)。     

Island urbanization and its ecological consequences: A case study in the Zhoushan Island,East China

Islands, which provide multiple ecosystem services, are subject to increasing urbanization pressure due to the ongoing marine development, especially in developing countries. Insights into the island urbanization mechanism and its ecological consequences are essential to sustainable development. In the present paper, the satellite images, nighttime lights, and topographic data were integrated to characterize the spatially explicit urbanization process and mechanism during 1995–2011 in the Zhoushan Island, East China. Furthermore, the corresponding spatially explicit changes in ecosystem services, including net primary productivity (NPP), carbon sequestration and oxygen production (CSOP), nutrient cycling, crop production, and habitat quality, were quantified based on the Carnegie–Ames–Stanford Approach (CASA) and Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) models. The results showed that the Zhoushan Island had experienced a rapid urbanization over the years, with significant urban encroachment on the farmland and tidal flat. Moreover, the urban land expansion was positively correlated with that of the nighttime lights and negatively correlated with the elevation, slope, and the distance to shoreline. These indicated that the urban expansion was resulted from the enhancement of socioeconomic activities, and concentrated in the near-shore areas with low altitude and gentle slope. The urban encroachment on other land use types resulted in a decrease of 3.4 Gg C a⁻¹ NPP, 8.7 Gg a⁻¹ CSOP, 13.2 Gg a⁻¹ nutrient cycling, and 12.3 t a⁻¹ crop production, respectively. In addition, the habitat quality in 11% area of this island degraded substantially. Therefore, to achieve sustainable development of islands, it is urgent to implement more stringent policies, such as island spatial regulation, environmental impact assessment, intensive land use, and urban greening, etc.