基于InVEST模型的黔中喀斯特地区建设用地扩张碳储量变化研究.pdf

西北林学院学报2 0 2 4,3 9(1):1 8 5-1 9 2J o u r n a l o f N o r t h w e s t F o r e s t r y U n i v e r s i t y d o i:1 0.3 9 6 9/j.i s s n.1 0 0 1-7 4 6 1.2 0 2 4.0 1.2 5基于I n V E S T模型的黔中喀斯特地区建设用地扩张碳储量变化研究 收稿日期:2 0 2 2-1 1-1 4 修回日期:2 0 2 3-0 2-0 9 基金项目:贵州省教育厅青年科技人才成长项目(黔教合KY字2 0 2 22 0 2号);贵州财经大学在校学生科学研究项目(2 0 2 2 Z X S Y 1 1 0)。第一作者:李 月,博士,副教授。研究方向:土地资源管理与生态修复。E-m a i l:9 7 6 4 3 8 5 2 8q q.c o m李 月,罗红芬(贵州财经大学 公共管理学院,贵州 贵阳 5 5 0 0 2 5)摘 要:选择典型生态脆弱区黔中典型喀斯特地区普定县为研究对象,采用土地利用动态度和土地利用转移分析方法,应用I n V E S T模型测算2 0 0 3-2 0 2 0年的陆地生态系统碳储量及其空间分布格局,探讨建设用地扩张对碳储量变化的影响。结果表明:1)2 0 0 3-2 0 2 0年建设用地面积一直保持增加趋势,扩张速率表现为减缓态势,空间变化在中南部区域最为显著,主要侵占旱地、水田和灌木林地。2)普定县碳储量整体呈上升趋势,累计增加1 5.0 71 05t,主要因碳密度较高的灌木林地大幅挤占旱地。3)建设用地变化致使碳储量损失了2 2.3 71 04t,以中南部中心城区碳储量下降最为明显。因此,建议研究区在城市扩张过程中优先考虑碳密度较低的未利用地,同时限制建设用地对旱地、水田以及灌木林地的占用,以减少建设用地扩张对区域碳储存能力的影响。关键词:贵州;喀斯特;I n V E S T模型;建设用地;碳储量中图分类号:S 7 7 1.8 文献标志码:A 文章编号:1 0 0 1-7 4 6 1(2 0 2 4)0 1-0 1 8 5-0 8C a r b o n S t o r a g e C h a n g e o f C o n s t r u c t i o n L a n d E x p a n s i o n i n C e n t r a l G u i z h o u s K a r s t R e g i o n B a s e d o n I n V E S T M o d e lL I Y u e,L U O H o n g-f e n(S c h o o l o f P u b l i c A dm i n i s t r a t i o n,G u i z h o u U n i v e r s i t y o f F i n a n c e a n d E c o n o m i c s,G u i y a n g 5 5 0 0 2 5,G u i z h o u,C h i n a)A b s t r a c t:T o e l u c i d a t e t h e e f f e c t o f c o n s t r u c t i o n l a n d e x p a n s i o n o n c a r b o n s t o r a g e c a p a c i t y o f t e r r e s t r i a l e-c o s y s t e m i n t y p i c a l e c o l o g i c a l l y v u l n e r a b l e a r e a s i n C h i n a,P u d i n g C o u n t y,a t y p i c a l K a r s t a r e a i n c e n t r a l G u i z h o u,w a s s e l e c t e d a s t h e r e s e a r c h o b j e c t.T h e d y n a m i c a t t i t u d e i n d e x a n d l a n d u s e t r a n s f e r a n a l y s i s m e t h o d s w e r e u s e d t o e s t i m a t e a n d a n a l y z e t h e c a r b o n s t o r a g e a n d i t s s p a t i a l d i s t r i b u t i o n p a t t e r n o f d i f f e r-e n t l a n d t y p e s i n t h e c o u n t y f r o m 2 0 0 3 t o 2 0 2 0,c o m b i n e d w i t h t h e c a r b o n m o d e l I n V E S T,a n d t h e i n f l u e n c e o f c o n s t r u c t i o n l a n d e x p a n s i o n o n t h e c h a n g e o f c a r b o n s t o r a g e w a s d i s c u s s e d.T h e r e s u l t s s h o w e d t h a t 1)f r o m 2 0 0 3 t o 2 0 2 0,t h e c o n s t r u c t i o n l a n d m a i n t a i n e d a n e x p a n s i o n t r e n d,a n d t h e e x p a n s i o n r a t e s h o w e d a s l o w i n g t r e n d.S p a t i a l v a r i a t i o n w a s c o n c e n t r a t e d i n t h e c e n t r a l a n d s o u t h e r n r e g i o n s a n d e n c r o a c h m e n t m a i n l y o c c u r r e d o n d r y l a n d,p a d d y l a n d,a n d s h r u b l a n d.2)T h e e c o s y s t e m c a r b o n s t o r a g e i n P u d i n g C o u n t y s h o w e d a n o v e r a l l u p w a r d t r e n d f r o m 2 0 0 3 t o 2 0 2 0,w i t h a c u m u l a t i v e i n c r e a s e o f 1 5.0 71 05 t.T h e i n-c r e a s e i n s h r u b l a n d w a s t h e m a i n r e a s o n f o r t h e i n c r e a s e i n c a r b o n s t o r a g e.3)T h e c h a n g e o f c o n s t r u c t i o n l a n d r e s u l t e d i n a l o s s o f 2 2.3 71 04 t o f c a r b o n s t o r a g e,w h i c h w a s m a i n l y r e f l e c t e d i n t h e s i g n i f i c a n t d e-c r e a s e o f c a r b o n s t o r a g e i n t h e c e n t r a l u r b a n a r e a a n d t h e s o u t h e r n r e g i o n,w h i c h w a s d u e t o t h e c o n v e r s i o n o f s h r u b l a n d a n d p a d d y l a n d t o c o n s t r u c t i o n l a n d.T h e r e f o r e,i t i s s u g g e s t e d t h a t u n u s e d l a n d s h o u l d b e g i v e n p r i o r i t y i n t h e p r o c e s s o f u r b a n e x p a n s i o n i n t h e s t u d y a r e a,a n d t h e o c c u p a t i o n o f c o n s t r u c t i o n l a n d t o d r y l a n d,p a d d y f i e l d a n d s h r u b l a n d s h o u l d b e l i m i t e d t o r e d u c e t h e i m p a c t o f c o n s t r u c t i o n l a n d e x p a n-s i o n o n r e g i o n a l c a r b o n s t o r a g e c a p a c i t y.K e y w o r d s:G u i z h o u;K a r s t;I n V E S T m o d e l;c o n s t r u c t i o n l a n d;c a r b o n s t o r a g e 近百年的气候变化已经给全球自然生态系统和社会经济系统带来了深刻影响1。陆地生态系统中的碳储量能够通过吸收和释放大气中的温室气体,发挥巨大碳汇作用,对大气温室气体排放量以及碳循环的平衡具有重要意义2。土地利用覆被变化作为影响区域碳储量时空变化的关键原因之一3-7,土地利用变化的碳储量影响已成为全球气候变化及陆地生态系统碳汇研究的重点与热点。目前,国内外学者已经从土地利用的角度揭示了不同尺度下的碳储量空间变化2-1 1。B o l i n8通过对全球各地不同用地的碳储量进行分析,发现土地利用变化是导致大气环境C O2含量波动的主要因素;刘洋等9探究并预测了疏勒河流域1 9 9 0-2 0 1 5年及2 0 1 5-2 0 4 0年流域生态系统碳储量时空变化特征及其与土地利用方式之间的关系,发现未利用地向耕地和草地转化有利于碳储量增加,而草地向耕地和未利用地的转化则导致碳储量减少;侯瑞萍等1 2估算了2 0 2 0年长江经济带林地和其他生物质碳储量和碳汇量,并发现乔木林地的碳储量和碳汇量所占比例最大,以上研究均表明森林砍伐、林地转出会导致大量温室气体从陆地碳库排放至大气环境中。I n V E S T模型作为研究土地利用变化与生态系统服务功能相互关系较为成熟的生态模型之一,长期以来得到了广泛应用9-1 7。刘晓娟等1 6通过耦合F L U S-I n V E S T模型,从土地利用视角出发对中国2 1 0 0年的陆地生态系统碳储量进行模拟,并探讨其空间分异;邓喆等1 7基于I n V E S T模型C a r b o n模块计算因土地利用变化所导致的生态系统碳储量变化。大量研究表明,土地利用变化过程直接影响区域内植被、土壤等物质的固碳能力,进而影响整个区域碳储量的时空变化。纵观现有研究成果,内容主要聚焦在森林、农田、湿地等生态用地的碳储量时空变化规律以及驱动因素研究,而基于建设用地规模扩张的区域碳储量变化相对匮乏。在社会经济发展进程中,建设用地主要通过挤占林地、耕地等生态用地迅速扩张,致使陆地生态系统碳储量大幅流失,造成生态服务功能严重退化1 8。近年来,喀斯特地区受自身的脆弱性以及人类活动的影响,生态环境日益恶化,因此,探究喀斯特地区建设用地扩张引起的碳储量空间变化对区域碳平衡管理具有重要意义。贵州省作为国家生态文明试验区,具有人口容量小、生态系统敏感且脆弱等特征1 9,拥有巨大的固碳潜力。普定县地处黔中,岩溶分布广泛,土层浅薄,水土流失严重1 9,在快速城镇化和人口剧增的影响下,建设用地需求不断上升,区域土地利用方式的剧烈变化给敏感脆弱的生态环境带来巨大压力。本研究以黔中喀斯特区 普定县为对象,基于2 0 0 3、2 0 1 0、2 0 2 0年3期遥感影像数据,选取土地利用动态变化指数和土地利用转移矩阵方法,结合I n-V E S T模型C a r b o n模块,估算2 0 0 3-2 0 2 0年不同用地类型生态系统碳储量,探究区域生态系统碳储量空间分布格局和建设用地扩张对生态系统碳储量变化的影响,以期为区域生态可持续发展和典型生态脆弱区未来用地结构优化提供决策参考。1 材料与方法1.1 研究区概况普定县(1 0 5 2 7 4 9 -1 0 5 5 8 5 1 E,1 6 2 6 3 6 -2 6 3 1 4 2 N)位于贵州省西部,隶属安顺市,位于素有“黔之腹,滇之喉”之称的黔中腹地(图1),属亚热带季风湿润气候,年均降水量1 3 7 8.2 mm。年均气温1 5.1。地势总体呈南、北部高,中间低,全县总面积1 0 9 1 k m2,喀斯特地貌广泛发育,土层发育较薄,土地利用强度大。1.2 数据来源与处理所需基础数据主要包括2 0 0 3、2 0 1 0、2 0 2 0年3期遥感影像,碳库数据包括地上植被碳库Ca b o v e、地下植被碳库Cb e l o w、土壤碳有机碳库Cs o i l、死亡有机碳库Cd e a d 4大碳库(表1)。以3期遥感影像为基础数据,3期影像均采用近红外、红光与绿光波段进行标准假彩色合成,首先,以11 0 0 0 0地形图对3期影像进行几何精校正,并采用较为成熟的F L AA S H模型进行大气校正。其次,通过监督分类、人机交互式解译法进行目视解译提取地类,最后,结合普定县土地利用现状图、G o o g l e E a r t h高分辨率历史影像及全国国土调查主要数据公报等,对土地利用矢量图进行图斑修正。用地类型分为有林地、灌木林地、其他林地、草地、水田、旱地、水域、建设用地和未利用地等9种用地类型。1.3 土地利用动态度不同土地利用类型受自然和人为因素的干扰不同,其变化速率也不相同。而土地利用变化速率的区域差异可以用土地利用动态度模型来加以表述,单一土地利用动态度可以表示为2 9K=Ub-UaUa1T1 0 0%(1)681西北林学院学报3 9卷 式中:T为研究时段;K为研究时段某类用地类型的变化率;Ua和Ub分别为研究期开始和结束时段某土地利用类型的面积。1.4 土地利用转移矩阵转移矩阵是土地利用类型间相互转化的数量及方向定量研究的主要方法,能够具体反映土地利用变化的结构特征和各类型之间的转移方向。转移矩阵的数学形式可以表示为2 9Pi j=P1 1P1 2P1nP2 1P2 2P2nPn1Pn2Pn n (2)式中:P表示面积;n表示土地利用类型数;i、j分别表示研究初期与末期的土地利用类型;Pi j表示研究期间用地类型i向用地类型j的转移面积。图1 研究区地理位置F i g.1 O v e r v i e w o f t h e g e o g r a p h i c a l l o c a t i o n o f t h e s t u d y a r e a表1 研究区数据来源T a b l e 1 D a t a s o u r c e s f o r t h e s t u d y a r e a基础数据数据来源数据描述2 0 0 3、2 0 1 0、2 0 2 0年遥感影像地理空间数据云(h t t p:/w w w.g s c l o u d.c n/)L a n d s a t 7 E TM影像,分辨率3 0 m3 0 m生态系统碳库邹文涛等2 2、李默然等2 3、罗怀良等2 4、黄从红2 5、李义平2 6、解宪丽等2 7、李克让等2 8Ca b o v e、Cb e l o w、Cs o i l、Cd e a d1.5 I n V E S T模型法通过文献检索得到碳密度参数,土地利用碳密度数据主要收集研究区以及相邻区域范围内的研究成果,其中有林地、其他林地参考邹文涛等2 2的研究;灌木林地碳密度值参考李默然等2 3的研究成果;旱地、水田、草地、未利用地碳密度值主要参考罗怀良等2 4和黄从红2 5的研究成果;建设用地参考李义平2 6的研究;其他地类主要参考解宪丽等2 7和李克让等2 8的研究成果(表2)。土地利用碳密度数据结合土地利用数据(表2),基于I n V E S T模型中C a r b o n模块评估普定县碳储量的变化,碳储量计算公式为C=Cb o v e+Cb e l o w+Cs o i l+Cd e a d(3)式中:C表示生态系统碳储量;Ca b o v e表示地上植被碳库;Cb e l o w表示地下植被碳库;Cs o i l表示土壤有机碳库;Cd e a d表示死亡有机碳库。2 结果与分析2.1 土地利用变化分析2.1.1 土地利用动态度分析 利用土地利用动态度方法(式(1)对普定县各用地类型变化速率进行分析。由表3、图2可知,2 0 0 3-2 0 2 0年建设用地动态变化最剧烈,2 0 0 3-2 0 1 0年动态度高达1 4.1 3%,而2 0 1 0-2 0 2 0年变化速率出现下降,动态度降为1 1.1 3%;灌木林地面积在2 0 0 3-2 0 1 0年不断增加,且 动 态 变 化 程 度 仅 次 于 建 设 用 地,动 态 度 为1 0.0 9%,而水田、水域、其他林地、旱地在2 0 0 3-2 0 1 0年动态变化程度相对较小,且用地规模均呈减少的变化态势;2 0 1 0-2 0 2 0年水域变化速率上升至5.3 2%,其 他 林 地 变 化 速 率 大 幅 增 加,增 幅 达1 0 7.0 6%,表明其他林地在该时期以更快的速度缩减,水田和旱地的减少速率逐渐减弱,动态度分别为3.6 4%、0.8 4%。781第1期李 月 等:基于I n V E S T模型的黔中喀斯特地区建设用地扩张碳储量变化研究2.1.2 建设用地时空变化特征分析 对各个时期的土地利用类型转移情况进行分析(图3)。由图3(a)可知,2 0 0 3-2 0 1 0年普定县旱地转出最多,共流出2 2 3.1 9 k m2,其中有9.6 4%转入建设用地,水域转化面积最小,仅有1 2.2 3 k m2。由图3(b)得知,2 0 1 0-2 0 2 0年仍主要以旱地转出为主,但转出面积降低至1 4 4.7 8 k m2,其中有2 0.6 7%转入建设用地。对2 0 0 3-2 0 1 0、2 0 1 0-2 0 2 0年普定县建设用地转移情况进行分析(表4),发现研究区在2 0 0 3-2 0 2 0年建设用地经历了“快速扩张-平缓扩张”的变化过 程。其 中2 0 0 3-2 0 1 0年,建 设 用 地 净 扩 张2 3.0 5 k m2,旱地和水田的贡献率分别为7 3.8 8%和2 5.9 9%,说明占用旱地和水田是研究区该时期城镇表2 普定县土地利用类型碳密度T a b l e 2 C a r b o n d e n s i t y o f l a n d u s e t y p e s i n P u d i n g C o u n t yth m-2土地利用类型碳密度(以C计)地上地下土壤死亡有机物有林地3 1.9 26.3 81 4 6.8 23.5 0其他林地3 5.0 37.0 19 1.7 03.5 0灌木林地2 6.6 06 7.5 09 4.0 03.5 0草地0.8 20.8 78 9.2 01.0 0水田5.4 21.9 61 4 6.2 41.0 0旱地3.2 90.5 91 0 9.6 81.0 0水域0.0 00.0 00.0 00.0 0建设用地0.0 00.0 01 1 1.2 60.0 0未利用地0.7 40.1 36 9.9 20.0 0表3 2 0 0 3-2 0 2 0年土地利用动态度T a b l e 3 L a n d u s e d y n a m i c a t t i t u d e f r o m 2 0 0 3 t o 2 0 2 0%研究时段有林地其他林地灌木林地草地水田旱地水域建设用地未利用地2 0 0 3-2 0 1 0年8.5 4-3.4 21 0.0 92.3 9-4.8 6-3.2 0-3.6 51 4.1 30.7 22 0 1 0-2 0 2 0年3.5 4-7.0 83.9 2-6.0 8-3.6 4-0.8 45.3 21 1.1 3-8.6 9图2 2 0 0 3-2 0 2 0年普定县土地利用现状F i g.2 L a n d u s e s t a t u s m a p o f P u d i n g C o u n t y f r o m 2 0 0 3 t o 2 0 2 0图3 2 0 0 3-2 0 1 0年普定县土地利用转移矩阵弦图F i g.3 C h o r d d i a g r a m o f l a n d u s e t r a n s f e r m a t r i x i n P u d i n g C o u n t y f r o m 2 0 0 3 t o 2 0 1 0881西北林学院学报3 9卷 表4 2 0 0 3-2 0 2 0年普定县建设用地变化统计T a b l e 4 S t a t i s t i c s o n t h e c h a n g e o f c o n s t r u c t i o n l a n d i n P u d i n g C o u n t y f r o m 2 0 0 3 t o 2 0 2 0k m2土地利用类型2 0 0 3-2 0 1 0年转出转入总计2 0 1 0-2 0 2 0年转出转入总计有林地0.3 30.4 30.0 90.5 12.2 91.7 8其他林地0.2 50.8 30.5 90.0 32.0 21.9 9灌木林地2.3 21.3 7-0.9 51.3 51 0.6 79.3 2草地0.1 70.3 20.1 60.0 33.4 13.3 7水田5.6 01 1.5 85.9 90.0 62 1.7 32 1.6 7旱地5.4 42 2.4 71 7.0 31 7.2 62 9.9 31 2.6 7水域0.0 60.1 20.0 70.0 80.2 70.1 9未利用地0.0 70.1 50.0 80.0 00.5 70.5 6总计1 4.2 33 7.2 82 3.0 51 9.3 27 0.8 95 1.5 6扩张的主要方式。2 0 1 0-2 0 2 0年,建设用地共转入7 0.8 9 k m2,其中旱地、水田与灌木林地的贡献率最高,分别达到4 2.2 2%、3 0.6 5%、1 5.0 5%;共转出面积1 9.3 2 k m2,年均转出面积有所减少,该时期灌木林地也成为城镇扩张占用的主要用地类型之一,究其原因主要与建设用地需求上升有关。2 0 0 3-2 0 2 0年普定县建设用地扩张情况见图4,建设用地主要分布在地势相对平缓的中部和南部区域。2 0 1 0年建设用地面积相比2 0 0 3年扩张了约2倍,主要以南部区域最为显著,向东边呈带状蔓延扩张;其次为中部区域,由中心城区向南部以及东南部方向呈带状延伸发展,东部、北部则以零星点状零散分布,主要占用水田和旱地。2 0 2 0年建设用地面积上升至9 7.9 1 k m2,主要沿中心城区向东、南部区域呈带状发展,通过大规模占用水田向南部扩张,向东边扩张主要侵占旱地;南部伴有少量建设用地占用旱地。总的来说,普定县建设用地呈现南部率先大幅扩张后缩减、中部后期显著扩张、东部后期显著扩张的地域特征。图4 2 0 0 3-2 0 2 0年普定县建设用地扩张情况F i g.4 E x p a n s i o n o f c o n s t r u c t i o n l a n d i n P u d i n g C o u n t y f r o m 2 0 0 3 t o 2 0 2 02.2 碳储量时空变化规律特征分析由图5可见,研究区不同历史时期各土地利用类型碳 储 量 发 生 显 著 变 化。基 于I n V E S T模 型C a r b o n模块,得出2 0 0 3-2 0 2 0年普定县生态系统碳储 量 分 别 为1 5 0.2 21 05、1 5 8.3 41 05t和1 6 5.2 91 05t,整 体 呈 逐 期 递 增 趋 势,累 计 增 加1 5.0 71 05t。用地类型中灌木林地碳储量较多且增值最大,2 0 0 3-2 0 1 0年灌木林地碳储量年均增量达2.7 2 1 05t,未利用地碳储量变化最为稳定,2 0 0 3-2 0 2 0年碳储量仅下降1.2 11 05t。2 0 0 3-2 0 2 0年建设用地碳储量由2.2 91 05t增至1 0.8 71 05t,各期分别上升2.4 7 1 05、5.8 6 1 05t。由图6可知,普定县2 0 0 3-2 0 2 0年灌木林地面积增加1 9 1.3 5 k m2,碳储量上升3 5.1 7 1 05t,是总增量的2.3 3倍,而建设用地面积净增加7 4.6 1 k m2,碳储量仅增加8.8 51 05t,建设用地单位面积碳储量增值相对较少。水田碳储量损失最多,达1 8.9 11 05t,是用地规模萎缩最大、碳储量损失最严重的地类。因此,灌木林地是普定县碳储量上升的主要因素,2 0 2 0年普定县固碳能力较强的地类为灌木林981第1期李 月 等:基于I n V E S T模型的黔中喀斯特地区建设用地扩张碳储量变化研究地,而水田的固碳能力逐渐减弱。由图7可知,普定县碳储量较高的区域集中分布在东北部和西南部,碳储量低值区分布在中南部区域,这可能与东北、西南部建设用地扩张程度有关。进一步分析碳储量在空间上的变化特征发现(图8),2 0 0 3-2 0 2 0年,普定县碳储量逐期升高且空间演变极其明显,整体表现为东北部和西南部地区碳储量增加,中南部区域碳储量显著降低。通过叠加土地 利 用 现 状 图(图2)和 碳 储 量 空 间 分 布 图(图7)可知,普定县东北部和西南部区域的碳储量显著增加区域主要表现为大量灌木林地、有林地转入,而中北部地区城市扩张剧烈,建设用地广泛分布,该区域碳储量减少地区主要以灌木林地、水田转为建设用地为主。图5 2 0 0 3-2 0 2 0年普定县土地利用类型碳储量F i g.5 C a r b o n s t o r a g e o f l a n d u s e t y p e s i n P u d i n g C o u n t y f r o m 2 0 0 3 t o 2 0 2 0图6 2 0 0 3-2 0 2 0年普定县土地利用变化及其碳储量变化F i g.6 L a n d u s e c h a n g e a n d i t s c a r b o n s t o r a g e c h a n g e i n P u d i n g C o u n t y f r o m 2 0 0 3 t o 2 0 2 0图7 2 0 0 3-2 0 2 0年普定县碳储量空间分布F i g.7 S p a t i a l d i s t r i b u t i o n o f c a r b o n s t o r a g e i n P u d i n g C o u n t y f r o m 2 0 0 3 t o 2 0 2 0图8 2 0 0 3-2 0 2 0年普定县碳储量空间变化F i g.8 S p a t i a l v a r i a t i o n o f c a r b o n s t o r a g e i n P u d i n g C o u n t y f r o m 2 0 0 3 t o 2 0 2 02.3 建设用地扩张对碳储量的影响由表5可知,2 0 0 3-2 0 2 0年普定县建设用地变化致使碳储量损失了2 2.3 71 04t,说明建设用地在转出过程中的碳储增量远远小于其扩张引起的碳损091西北林学院学报3 9卷 失量。其中,2 0 0 3-2 0 1 0年建设用地变化碳储量流失2.5 31 04t,其转出过程中,建设用地向水田和灌木林地转移,使得水田、灌木林地碳库碳储量分别增加2.1 41 04、1.7 11 04t;在扩张进程中,建设用地仅占用草地、水域、未利用地,有利于生态系统碳储量增加,占用其他用地碳储量均有不同程度的流失,其中占用水田使得水田碳库减少碳储量4.6 21 04t。2 0 1 0-2 0 2 0年,建设用地变化致使碳储量净损失达1 9.8 4 1 04t,损失量约为上一阶段的8倍。主要因为该时期建设用地转出仅增加碳储量2.7 61 04t,且6 0.2 9%均来源于建设用地转灌木林地,而建设用地扩张碳储量损失了2 2.6 01 04t,是其转出碳储量增加的8.1 8倍,其中灌木林地和水田转为建设用地分别引起碳储量损失8.7 2 1 04、9.5 4 1 04t。表5 2 0 0 3-2 0 2 0年普定县建设用地转化的碳储量变化T a b l e 5 C a r b o n s t o r a g e c h a n g e s o f c o n s t r u c t i o n l a n d c o n v e r s i o n i n P u d i n g C o u n t y f r o m 2 0 0 3 t o 2 0 2 01 04t土地利用类型2 0 0 3-2 0 1 0年建设用地转出建设用地转入总计2 0 1 0-2 0 2 0年建设用地转出建设用地转入总计有林地0.1 8-0.1 80.0 00.5 7-2.1 8-1.6 0其他林地0.0 5-0.2 4-0.1 90.0 1-0.6 2-0.6 1灌木林地1.7 1-1.0 90.6 21.6 7-1 0.3 9-8.7 2草地-0.0 30.0 30.0 0-0.0 20.7 00.6 8水田2.1 4-4.6 2-2.4 80.0 7-9.6 0-9.5 4旱地0.1 5-0.6 9-0.5 40.5 8-1.0 1-0.4 4水域-0.0 60.0 90.0 3-0.1 10.3 40.2 3未利用地-0.0 20.0 50.0 30.0 00.1 60.1 6总计4.1 2-6.6 5-2.5 32.7 6-2 2.6 0-1 9.8 43 讨论3.1 建设用地扩张特征随着我国西部大开发的实施,黔中经济区普定县城镇化水平显著提升,社会经济条件的作用逐渐凸显,建设用地扩张剧烈,在地势条件相对平坦的中南部中心城区表现尤为明显。2 0 0 3-2 0 2 0年,建设用地一直保持扩张趋势,扩张速率表现为减缓态势,空间变化主要集中在中南部区域,集中连片分布明显,与城市发展范围保持一致;呈零点状分布的建设用地在2 0 1 0-2 0 2 0年大幅扩张,主要原因可能是由于各乡镇的城镇化水平显著提升所导致,建设用地面积增加主要由水田、灌木林地2种地类转入。虽然建设用地面积持续增加,但值得注意的是,2 0 1 0年南部呈带状分布的建设用地以及分布在中心城区西部的建设用地,于2 0 2 0年逐步恢复为旱地,这可能是由于在2 0 1 6年普定县实施城乡建设用地增减挂钩,拆旧区复垦耕地引起的。3.2 建设用地扩张的碳储量变化普定县在2 0 0 3-2 0 2 0年因建设用地扩张使得中南部中心城区碳储量大幅下降,用地类型变化具体表现为灌木林地、水田转为建设用地,可能由于该区域灌木林地、水田在转为建设用地的过程中,进入土壤的腐殖质减少,土壤腐殖质的矿化作用增加,进而使得土壤碳储量减少2 9,这与向书江等1对重庆主城区的研究结果较为一致,但不同的是,普定县建设用地扩张占用旱地几乎不对碳储量增加产生负面影响,而水田和灌木林地因建设用地侵占碳储量分别损失了1 4.2 21 04、1 1.4 81 04t,不利于区域碳汇。因此,为了提升区域碳汇能力,应注意在建设用地扩张过程中减少水田与灌木林地的占用,优先考虑固碳能力较弱的未利用地,促进土地布局优化,减少碳排放,提高区域生态效益,以达到生态与经济协调发展的目的。3.3 不确定分析本研究有效探究了普定县建设用地扩张的碳储量变化,但仍存在以下不确定性:首先,由于碳密度数据主要参考喀斯特地区研究成果,部分难以获取碳密度数据的地类通常参考非喀斯特地区,但其运行结果能够清晰地反映普定县近2 0 a来建设用地扩张的碳储量时空变化规律,为更好地反映区域的生态碳储存能力,在今后喀斯特地区碳储量研究中应加强喀斯特地区用地类型碳密度实测;其次,本研究土地利用数据经国土调查资料检验,且与秦罗义等3 1对普定县土地利用结构研究相似,精度达到一定要求,但土地利用数据是通过遥感影像解译获得的,土地利用数据精度会对研究结果有所影响,因此在今后的研究中可采用更高分辨率遥感影像,减少主观误差,以提升用地类型碳储量估算精度。4 结论2 0 0 3-2 0 2 0年建设用地一直保持扩张趋势,扩张面积达7 4.6 1 k m2,扩张速率表现为减缓态势,空间变化主要集中在中南部中心城区,通过侵占旱地、191第1期李 月 等:基于I n V E S T模型的黔中喀斯特地区建设用地扩张碳储量变化研究水田和灌木林地得以发展。2 0 0 3-2 0 2 0年普定县生态系统碳储量整体呈上升趋势且空间演变极其明显,累计增加1 5.0 71 05t,灌木林地增加是碳储量上升的主要原因,空间变化整体表现为东北部和西南部地区碳储量增加,中南部区域碳储量降低。建设用地转化致使碳储量损失了2 2.3 71 04t,其中2 0 1 0-2 0 2 0年损失碳储量占8 8.6 9%,主要表现为中南部区域碳储量大幅下降,源于灌木林地、水田转为建设用地。参考文献:1 向书江,张骞,王丹,等.近2 0年重庆市主城区碳储量对土地利用/覆被变化的响应及脆弱性分析J.自然资源学报,2 0 2 2,3 7(5):1 1 9 8-1 2 1 3.X I NA G S J,Z HANG Q,WAN G D,e t a l.R e s p o n s e a n d v u l n e r-a b i l i t y a n a l y s i s o f c a r b o n s t o r a g e t o L U C C i n t h e m a i n u r b a n a r e a o f C h o n g q i n g d u r i n g 2 0 0 0-2 0 2 0 y e a r sJ.J o u r n a l o f N a t u r a l R e s o u r c e s,2 0 2 2,3 7(5):1 1 9 8-1 2 1 3.(i n C h i n e s e)2 唐睿,彭开丽.土地利用变化对区域陆地碳储量的影响研究综述J.江苏农业科学,2 0 1 8,4 6(1 9):5-1 1.3 李姣,汪杰,李朗,等.洞庭湖生