天山北麓新人工绿洲扩张过程中的碳库格局变化与碳库迁移

张 芳, 熊黑钢, 冯 娟, 许仲林

1 新疆大学资源与环境科学学院,乌鲁木齐 830046 2 北京联合大学应用文理学院,北京 100083 3 教育部新疆绿洲生态重点实验室,乌鲁木齐 830046

天山北麓新人工绿洲扩张过程中的碳库格局变化与碳库迁移

张 芳1,3, 熊黑钢2,3,*, 冯 娟1,3, 许仲林1,3

1 新疆大学资源与环境科学学院,乌鲁木齐 830046 2 北京联合大学应用文理学院,北京 100083 3 教育部新疆绿洲生态重点实验室,乌鲁木齐 830046

基于地面调查数据、气象数据和遥感数据分析了40年来天山北麓新人工绿洲迅速扩张对区域植被和土壤的碳库格局、碳库迁移及碳储量的影响,结果如下：1976—2016年,研究区农田面积占比从3.25%增加到40%以上,沼泽消失,水库干涸,灌丛大幅减少,土壤盐碱化过程停止,裸盐碱地面积40年里减少近70%,新生草地在裸盐碱地上形成；2016年6月植被有效碳储量约为0.122Tg,比1976年下降了15.7%,成为一个弱的碳源；土壤碳库在干旱地区陆地生态系统碳库中占绝对主导地位,1976年研究区(0—60cm)的土壤有机碳储量为7.814Tg,其中耕地仅占4.2%,1996年土壤有机碳储量比1976年减少12.4%,呈明显碳源特征,到2016年,耕地土壤有机碳储量占比超过50%,总碳储量与1996年相比微降1.2%,土壤碳库重新趋于稳定。其中,裸盐碱地土壤碳库在向耕地、灌丛地、草地土壤碳库的迁移过程中贡献了1.265Tg的碳储量增量,草地土壤碳库在向耕地、灌丛地碳库的迁移过程中贡献了0.894Tg的碳储量增量；若维持现有耕地规模不再进行新的垦荒活动,40a垦荒造成的土壤碳库损失可以逐渐恢复并重新形成碳汇效应。

新人工绿洲扩张；碳库格局变化；碳库迁移；天山北麓

人类活动对碳循环的影响是气候变化研究的重要核心问题之一。大量研究表明,开垦会在一定程度上产生碳源效应[1- 4]。例如,三江平原沼泽湿地开垦初期土壤有机碳损失较快,经过15a以上稳定而长期的耕作后,有机碳最终趋于一个相对的稳定值[5- 6]；利用历史土地覆盖数据和簿记模型对东北地区过去300年来开垦导致碳排放量的估算研究表明,东北全境约38%的草地、20%的森林和灌丛被开垦,开垦导致的碳排放量中土壤碳库是较大的排放源[7]；开垦在新疆干旱荒漠地区也可能产生有利于碳汇的情况,自然土壤垦殖后最初5年有机碳含量增加明显,年均增加0.65g/kg以上,5年后维持在相对平衡的水平[8]。

新疆对于我国国家安全、经济发展及“一带一路”建设均具有极其重要的战略地位,是我国21世纪经济增长的重要支点和资源接替区,也是保障国家粮食安全的重要耕地储备基地[9]。1976—2016年,由于新人工绿洲快速扩张,面积仅占全疆8.22%的天山北麓经济带,人口和耕地面积超过全疆的25%[10]。开垦是对绿洲生态系统影响最大的人类活动,本研究尝试采用地面土壤调查数据、遥感数据和气象数据,结合光能利用率模型CASA(Carnegie-Ames-Stanford-Approach),研究天山北麓新人工绿洲迅速扩张对区域植被、土壤碳库格局和碳库转移的影响,以增进对人类活动影响下干旱区绿洲生态系统碳库反馈效应的理解和认识。

1 数据与方法

1.1 研究区概况

20世纪中叶以后形成的绿洲称为新绿洲[11]。从天山北麓绿洲分布和扩张轨迹来看,新绿洲多分布在冲积平原中下部[12]。近几十年来,奇台绿洲的高强度农业开发在天山北麓具有代表性,因此研究区选择在奇台平原绿洲下部与古尔班通古特沙漠东部毗邻的交错带区域,范围：43°57′—44°13′N,89°24′—90°09′E,其中非沙漠区域面积为1582.6km2,年均气温5℃左右,年均降水量约176mm,年平均蒸发势达2141mm,属温带大陆性干旱/半干旱气候。

1.2 数据资料

①遥感数据：根据数据质量(云量、条带等)和时相匹配(研究区植物生长最旺盛季节为6月中旬前后)进行综合考虑,选用1976年6月6日的Landsat MSS影像、1996年6月28日和2016年6月19日的Landsat TM影像共3期卫星数据,时间间隔为20a。遥感数据使用前进行严格的预处理,包括大气辐射校正和几何精校正,几何校正误差不超过0.2个像元。

②地面土壤实测数据：考虑到与影像过境时间相对应,土壤数据选用于2009年6月9—23日采集的野外自然土壤和农田土壤样本所测定的结果：采样时在研究区未受开垦干扰的荒地定位布点60个,样点尽可能规则遍及所有荒地类型,用四分法采集剖面上0—20、20—40cm和40—60cm的土壤样本；从西向东在研究区选择耕种10年左右和耕种>20年的农田各挖取土壤剖面25个,同样采集0—20、20—40、40—60cm的土壤样本。所有土壤样本在实验室内自然晾干、磨碎,过2mm筛后,用于土壤理化分析。

③其他数据：1976—2016气象数据、2004年以来历年定点观测、入村走访记录和景观照片库等。

1.3 研究方法

①土地覆被变化分类：利用历年野外定点观测、走访记录和实地景观照片库,采用监督分类方法对各期影像进行地物覆盖类型分类与验证,计算通过遥感专业处理软件ENVI 5.0完成。

②植被地上生物量与地下生物量估算：用植被净初级生产力(Net Primary Productivity, NPP)估算植被地上生物量。NPP指单位面积的绿色植物在单位时间内累积的有机物量,是判定生态系统碳汇的主要因子[13]。目前,基于光能利用率的遥感模型被广泛应用于估算区域及全球尺度的NPP[14- 15],本研究采用CASA模型进行NPP估算,公式如下：

NPP=APAR×Tε1×Tε2×Wε×ε*

(1)

式中,APAR是光合有效辐射,Tε1和Tε2为温度胁迫系数,Wε为水分胁迫影响系数,最大光能转化率ε*取值0.389gC/MJ,具体计算过程和相关气象参数详见文献[16],需要说明的是,新疆绿洲为灌溉农业,灌溉过程由生产者根据土壤墒情决定,不受降水条件限制,因此本研究中,农田与沼泽的Wε取值为1。

不同植被类型具有不同的地下与地上部分生物量的比例系数,通过比例系数可以计算植被的地下生物量。本文采用7.89、4.42、15.68作为温性荒漠草原、灌丛、沼泽类的地下与地上部分生物量的比例系数[17]；农作物根茎比采用0.1[18]。

③土壤有机碳密度计算：土壤有机碳用重铬酸钾-外加热法测定,土壤容重通过环刀法测定,土壤有机碳密度计算公式如下[19]：

(2)

式中,SOCD为有机碳密度(kg/m2)；i为土层；Di为土层厚度(cm)；BDi为第i层土壤容重(g/cm3)；SOCCi为第i层土壤有机碳百分含量(%)。

2 研究结果

2.1 新人工绿洲扩张与地表覆盖格局演变

根据研究区的历年定点观测和入户访查资料,利用监督分类方法对3期影像进行分类(图1),分类精度分别为：96.38%、92.75%和98.45%。研究区目前主要的覆被类型为农田、草地、灌丛与裸盐碱地。种植面积最大的农作物有小麦、玉米、油葵及打瓜等；研究区天然植被主要为草本、灌木或半灌木的盐(碱)生植物,草地以叉毛蓬、猪毛菜、碱蓬、骆驼蓬、骆驼刺等低矮草本植物为主；灌丛由株形较大的灌木/半灌木与草本混合群落构成,例如柽柳+芨芨草+苦豆子+白刺等。研究区土壤呈碱性,土壤pH值>9.5的地方主要是光板地,几乎无植物生长,因此将其划分为裸盐碱地[20]。部分田边、路边有防护林成线状分布,但未能形成超过一个像元宽度(30m)的可提取面积信息,因此未将林地单独分类。

二是教学内容更新不足.当前教学内容仍以传统教材为主体，着重讲授地图学的基础理论和传统地图制作方法.随着大数据的到来和GIS技术的迅速发展，传统地图学教学内容已无法适应信息时代下地图制图的要求[10]，如传统纸质地图制图方法及工艺过程已逐渐被数字制图新理论、新技术和新方法所取代.因此在教学内容上应及时做出更新和调整，引入和补充数字地图的相关理论和制图流程，并将最新的地图制印和工艺流程纳入课堂教学中，为后续实习做好准备.

研究区东南部海拔较高,向西北处逐渐降低[21]。由图1可知,1976年时研究区地表水资源较为丰富,西北角建有一平原水库,泉水溢出带由西南向东北弧状分布,形成长短不一、流向向北的泉水河,泉水溢出带在影像上表现为沼泽和灌丛特征,经实地走访,沼泽及泉水河周围生长有茂密的芦苇+芨芨草混合灌丛；1976年分类图上的大面积白色区域为高反射率的裸盐碱地；1976年研究区农田稀少,呈零星斑块在中、西部小范围散布,说明当时人类活动强度很弱,研究区主要表现为天然草地+灌丛+沼泽+裸盐碱地的荒漠草原景观。1996年,新垦耕地在原有田块周围迅速扩张,泉水溢出带消失且大部分被开垦,水库干涸,灌丛明显减少；2016年,地表覆盖主要由农田+草地+裸盐碱地构成,适合大型农业机械作业的条田连接成片,原来以荒漠草原景观为主的天然绿洲被新人工农业绿洲所替代。

图1 研究区地表覆盖分类图Fig.1 The classified images of study area from 1976 to 2016

提取各年份不同地物类别的覆盖面积(表1),1976年有农田5.15×103hm2,仅占研究区面积的3.25%,1996年,农田增加至34.49×103hm2,是1976年的6.7倍,2016年农田占比已超过研究区面积的40%。由于高密度打井开采地下水来灌溉新垦耕地,到1996年时研究区地下水位降幅超过10m[22],土壤中的毛管悬着水与毛管上升水不再连接,土壤盐碱化进程停止,活性盐碱土变成残余盐碱土,在研究区近180mm年降水量的淋洗下,表层土壤盐碱化程度逐渐减弱,植物开始生长,裸盐碱地逐渐向新生草地演变甚至被开垦、改良成为农田。40a里,裸盐碱地面积减少41.06×103hm2,减少近70%,同期,灌丛则从17.87×103hm2减少到5.20×103hm2,主要散布于田间地头,已很少成片分布。

表1 1976—2016年地表覆盖类型的面积变化

通过迁移矩阵来观察40年来地表覆被分布格局的变化(表2),1976年沼泽面积中的75.58%、水体面积的77.96%、灌丛面积的66.93%、草地面积中的41.43%和裸盐碱地中的23.86%被开垦,构成了2016年农田；1976年草地中的43.62%在40年后仍旧是草地,这部分草地被新垦农田分割呈碎片化分布,与59.38%裸盐碱地演替而来的新生草地重新连接成片,加上来自21.44%的灌丛、16.16%的水体、16.57%的沼泽和7.40%的农田退化形成的草地,共同构成2016年的草地景观。1976年的裸盐碱地中仅有15.01%到2016年时依然存在,其余的裸盐碱地来自少量的草地、沼泽和灌丛等的退化。新人工绿洲扩张使农田大幅增加,区域水、土环境及覆被格局均发生巨大改变,因此,区域碳库平衡和碳库构成也必然随之发生改变。

表2 2016年相比1976年的地物覆被类型迁移/%

2.2 新人工绿洲扩张过程中的植被碳库反馈效应

植被碳库由植物的地上生物量和地下生物量两部分构成。地上生物量用研究区6月NPP来计算,地下生物量采用地下与地上生物量比例系数来计算,植被碳储量按照植物有机干物质中碳占的比重(国际上常用转换率为 0.45)转换为碳量[23- 24]。

2.2.1 地上与地下生物量变化

分别计算除水体和裸盐碱地外各地物类别中各个像元的NPP值(图2),介于120—70g·C/m2之间的主要为沼泽、芨芨草+苦豆子+柽柳+樟味藜+盐爪爪等混生灌丛群落,≤70g·C/m2的主要为荒漠草原,草地的单位面积NPP约在30—40g·C/m2范围内波动,相对较稳定。农田单位面积NPP最高,多高于120g·C/m2。

图2 研究区不同年份NPP分布Fig.2 The spatial patterns of above-ground NPP indifferent years

计算不同时期植被的地上与地下生物量总量(表3),1976年6月,研究区地上生物量总计为41.02×103t,到2016年6月达到90.60×103t,增加了120.88%,其中,农田扩张产生的地上生物量增量达55.56×103t。1976年农作物对研究区地上生物量的贡献率为8.62%,1996年急增至49.09%,2016年则达到了65.23%。由于农作物地下生物量较小,2016年农田地下生物量在整个研究区地下生物量中占比不足2.5%。研究区1976年的地下生物量为284.08×103t,1996年迅速下降到211.73×103t,2016年有所回升为245.30×103t,比1976年仍减少13.65%。40年里,灌丛和沼泽净损失了11.59×103t地上生物量和88.57×103t地下生物量,而草地的地上、地下生物量均呈增加态势,分别增加了5.61×103t和44.24×103t。从植被总生物量来看,1976年研究区的地上与地下生物量之和为325.10×103t,1996年为268.26×103t,2016年为335.90×103t,比1976年略增3.3%,说明研究区的新人工绿洲扩张过程中,植被总生物量响应特征为迅速下降—逐渐恢复—相对稳定,植被生物量分布格局显著改变,地上生物量大幅增加,地下生物量明显减少。

表3 1976—2016研究区6月份地上与地下NPP/103t

2.2.2 植被碳库迁移与碳储量变化

表4 研究区不同时期植被碳库构成/GgC

方精云等认为,过去20年里中国农作物生物量碳库和固碳效率增加,但这些增加的生物量碳绝大部分在短期内通过分解释放回大气中,因此将农作物生物量碳汇设定为零[23]。从这个角度出发,研究区的农作物碳储量在40年里增加了约15倍,均属于无效碳汇,而灌丛、沼泽等天然植被碳库向农作物碳库迁移,实际上是一个碳汇无效化过程。表5计算了2016年相较于1976年研究区植被碳库的迁移状况。1976年的农田中有1.55Gg作物碳在40年后依然是作物碳,只是固碳量增加了0.51Gg,属无效碳汇；1976年有17.82Gg灌丛碳迁移到了2016年的作物碳库中成为无效碳汇,迁移过程中固碳量损失12.41Gg,固碳能力下降了近70%；沼泽和草地在向2016年作物碳库的迁移过程中产生的碳排放分别为16.48Gg和25.56Gg,并各自将18.83Gg和38.53Gg的有效碳汇变成了无效碳汇；裸盐碱地在向草地和灌丛地演化过程中为植被碳库提供了58.36Gg有效碳汇补偿,占2016年有效碳汇的近一半。

综上所述,农田扩张是植被碳汇无效化过程,天然植被减少或消失使区域植被固碳潜力降低,而裸盐碱地向新生草地演替可产生植被有效碳汇补偿效应。从植被有效碳汇储量来看,1976年为144.55Gg,1996年为106.98Gg,下降了37.57Gg,2016年为121.91Gg,较1996年增加13.96%,说明,在新人工绿洲扩张初期,研究区植被碳库损耗最快,之后,由于裸盐碱地的碳汇补偿作用使碳损耗速度趋缓。2016年植被有效碳库较1976年下降了22.64Gg,说明人工绿洲扩张使研究区植被固碳潜力降低,呈弱碳源特征。

表5 2016年相较于1976年的植物碳库迁移/GgC

列数据为1976年不同植被碳库向2016年各类植被碳库的迁移,行数据为2016年植被碳库对来自1976年各类植被碳的接受情况； “+”、“-”符号表示碳库迁移过程中碳储量的增加或减少

2.3 新人工绿洲扩张过程中的土壤碳库反馈效应

研究区地带性土壤为荒漠灰钙土,发育在黄土状物质覆盖的山前倾斜平原和古老冲积平原上,矿质营养丰富,是新疆优质的土地资源之一[25]。

2.3.1 研究区土壤碳库构成

表6为实测的不同地类土壤有机碳密度分布情况。灌丛地土壤有机碳表聚性最高,其次是草地。灌丛地群落高度大、结构相对稳定,每年输入土壤的凋落物较多,因此土壤有机碳密度最大,0—60cm厚度土层中平均可达9.067kg/m2,其固碳量是草地的近2.4倍、裸盐碱地4.3倍。耕作10年左右的土壤有机碳密度与草地相近,比裸盐碱地高61.64%,但是仅为灌丛地的37.51%,耕作时间超过20年后,土壤有机碳密度明显提高,达6.404kg/m2,是裸盐碱地的3.04倍。

表6 土壤有机碳密度均值分布/(kg/m2)

研究区沼泽在1996年以前就全部消失,因此无法获得沼泽土壤有机碳的实测数据。对已发表文献中相关研究进行分析,全球湿地土壤(土层厚度0—100cm)的有机碳密度在42.5—66.1kg/m2之间[26-27]。牟长城等于5 月下旬在大兴安岭实测的天然灌丛沼泽土壤有机碳密度为33.06kg/m2(0—50cm)[28],霍莉莉对黑龙江省东南部未垦殖天然沼泽湿地土壤实测的有机碳密度为 29.1kg/m2(土层厚度0—60cm)[29]。考虑到淹水条件下,沼泽湿地有机碳积累主要受温度影响,霍莉莉的采样区纬度与本文研究区较接近且土层厚度一致[29],因此本文采用有机碳密度29.1kg/m2来估算研究区沼泽消失造成的土壤有机碳损失。

计算不同时期各地类土壤碳储量(表7),1996年和1976年耕地用耕种>20年的碳密度均值计算,1996年之后开垦的耕地用耕种10年的碳密度计算。1976年研究区土壤碳总储量为7814.1Gg,耕地、草地、灌丛地、裸盐碱地和沼泽地分别占4.2%、33.6%、20.7%、15.9%和25.6%,草地为最大土壤碳库。1996年由于沼泽消失,灌丛地减少等,研究区的土壤有机碳储量下降到6843.2Gg,呈碳源特征。2016年,由于耕地面积的大幅增加,耕地土壤有机碳储量占比近50%,超过草地,总碳储量与1996年相比仅微降1.2%,土壤碳库重新趋于稳定。

2.3.2 土壤碳库迁移与碳储量变化

新人工绿洲的大规模扩张导致地下水被大量超采,整个天山北麓尤以奇台县地下水超采最为严重[30],近几年来,研究区的垦荒行为开始被严格控制。

表7 研究区不同时期土壤碳库构成/GgC

假设研究区维持现有生态景观格局,不再进行新的开垦活动,已垦新农田在持续耕种后土壤碳汇效率趋于稳定和平衡(表8,其中耕地碳库均采用耕种>20年土壤有机碳密度计算),耕地土壤碳储量可达4126.8Gg,占土壤总碳库的53.9%,土壤总碳储量可达7653.9Gg,仅比1976年略低160.2Gg,与1996年相比回升11.8%,40年中开垦造成的土壤碳库损失可基本得到恢复。其中,裸盐碱地土壤碳库在向耕地、灌丛地、草地土壤碳库的迁移过程中可分别为总土壤碳储量贡献603.9、68.9、592.6Gg的碳汇增量,对土壤碳库的补偿贡献达16.53%。草地土壤碳库在向耕地、灌丛地碳库的迁移过程中则分别贡献了748Gg和145.5Gg的碳储量增量。沼泽地是干旱区最为脆弱的土壤碳库,在碳库迁移过程中造成的碳汇损失也最大。

表8 2016年覆被格局稳定条件下相较于1976年的土壤碳库迁移/GgC

列数据为1976年不同类型土壤碳库向2016年各类土壤碳库的迁移,行数据为2016年土壤碳库对来自1976年各类土壤碳的接受情况；“+”、“-”符号表示碳库迁移过程中碳储量的增加或减少

3 讨论

3.1 气候因素对研究区碳库的可能影响

40年来,天山北麓的气候呈现一定的暖湿化趋势,但升温主要在冬季,影响植物生长的春季温度仅上升了0.14℃/10a、夏季气温还略有下降,为-0.14℃/10a[31],说明温度对研究区植被生物量和植被碳库的影响较小。研究区近40年来的年降水量升率为15.8mm/10a[31],不足2mm/a的降水升率对干旱区植被碳库的影响也相对有限。

温度相对稳定情况下,水分和土壤结构的变化是影响土壤有机碳的分解或积累的直接因素。研究区原来碳密度较高的沼泽地、灌丛地、草地垦为耕地后,淹水条件和地表水供给消失,土壤透气性增加,将促进土壤中原有有机碳的分解和矿化形成碳排放；耕作土壤通过灌溉和种植可以保障水分和有机质的定期输入,长期耕作后的土壤会维持一个相对稳定的碳密度水平。因此,人类开垦活动是驱动研究区生态系统碳库格局改变和碳库迁移的决定性因素。

3.2 研究区土壤脱盐碱化与碳失汇

科学家们在全球碳循环研究中发现大气中CO2存在收支不平衡,认为可能存在一个未知的碳汇(失汇Missing sink)[32]。许多研究指出这个较大的未知碳汇可能存在于北半球中纬度地区的陆地生态系统中[33-34]。Houghton认为土地耕作制度的变化增加了全球碳循环研究的不确定性,促使 “碳失汇”的形成[35]。我国学者也认为土地利用方式改变对 “碳失汇”成因具有重要作用,揭示“碳失汇”将有助于理解碳的迁移和转化过程[36]。但是由于碳失汇问题的复杂性,目前许多研究尚有待于通过增加实验观测获取土地利用、气候以及植被、土壤碳氮动态详细资料来进行验证[37]。

本研究中,天山北麓强烈的人类开垦活动改变了地表水分的供应和分布,盐碱土发育过程中断,使原来大面积分布在荒漠-绿洲交错带的盐碱土进入脱盐碱化的逆向演变,盐碱土碳库逐渐向新生草地碳库、耕作土壤碳库等具有较高固碳潜力的碳库迁移,为开垦造成的碳排放提供了有效碳汇补偿,这种补偿能在多大程度上为 “碳失汇”做出贡献尚有待今后在更大区域尺度上进行进一步的综合观测和研究。

4 结论

1976—2016年,研究区新人工绿洲快速扩张,植被总生物量经历了从迅速下降到逐渐恢复的过程,在这个过程中,整个区域的生物量分布格局发生显著改变,植被地上生物量大幅增加,地下生物量大幅减少。开垦导致的天然植被碳库向作物碳库迁移是植被碳汇无效化的过程,植被无效碳汇增加使研究区植被固碳潜力降低,成为一个弱的碳源。

1976年研究区(0—60cm)的土壤碳总储量为7.814Tg,其中耕地、草地、灌丛地、裸盐碱地和沼泽地分别占4.2%、33.6%、20.7%、15.9%和25.6%,到2016年,耕地土壤有机碳储量占比近50%。裸盐碱地土壤碳库在向耕地、灌丛地、草地土壤碳库的迁移过程中为总土壤碳库增加了1.265Tg的碳储量,草地土壤碳库在向耕地、灌丛地碳库的迁移过程中增加了0.894Tg的碳储量。碳排放效应减弱,土壤碳库开始重新趋于稳定。在研究区水资源的制约下,若维持现有耕地规模不再进行新的垦荒活动,40a来开垦造成的土壤碳库损失可以逐渐得到恢复。

[1] 张利, 周广胜, 汲玉河, 白永飞. 中国草地碳储量时空动态模拟研究. 中国科学: 地球科学, 2016, 46(10): 1392- 1405.

[2] Poeplau C, Don A. Sensitivity of soil organic carbons tocks and fractions to different land-use changes across Europe. Geoderma, 2013, 192(1): 189- 201.

[3] Houghton R A, Hackler J L. Sources and sinks of carbon from land-use change in China. Global Biogeochemical Cycles, 2003, 17(2): 1034.

[4] Sellers P J, Dickinson R E, Randall D A, Betts A K, Hall F G, Berry J A, Collatz G J, Denning A S, Mooney H A, Nobre C A, Sato N, Field C B, Henderson-Sellers A. Modeling the exchanges of energy, water, and carbon between continents and the atmosphere. Science, 1997, 275(5299): 502- 509.

[5] 张金波, 宋长春, 杨文燕. 沼泽湿地垦殖对土壤碳动态的影响. 地理科学, 2006, 26(3): 340- 344.

[6] 江长胜, 王跃思, 郝庆菊, 王定勇. 土地利用对沼泽湿地土壤碳影响的研究. 水土保持学报, 2009, 23(5): 248- 252.

[7] 李蓓蓓, 方修琦, 叶瑜, 张学珍. 中国东北地区过去300年耕地开垦导致的碳收支. 中国科学: 地球科学, 2014, 44(9): 1987- 1996.

[8] 徐万里, 唐光木, 盛建东, 梁智, 周勃, 朱敏. 垦殖对新疆绿洲农田土壤有机碳组分及团聚体稳定性的影响. 生态学报, 2010, 30(7): 1773- 1779.

[9] 刘新平. 新疆绿洲土地资源可持续利用的经济学分析. 北京: 中国文史出版社, 2005: 1- 4.

[10] 程维明, 周成虎, 李建新. 天山北麓经济发展与绿洲扩张. 地理学报, 2002, 57(5): 561- 568.

[11] 樊自立, 马英杰, 艾力西尔·库尔班, 沈玉玲. 试论中国荒漠区人工绿洲生态系统的形成演变和可持续发展. 中国沙漠, 2004, 24(1): 10- 16.

[12] 韩德麟. 关于绿洲若干问题的认识. 干旱区资源与环境, 1995, 9(3): 13- 31.

[13] 杨延征, 马元丹, 江洪, 朱求安, 刘金勋, 彭长辉. 基于IBIS模型的1960- 2006年中国陆地生态系统碳收支格局研究. 生态学报, 2016, 36(13): 3911- 3922.

[14] Parton W J, Scurlock J M O, Ojima D S, Gilmanov T G, Scholes R J, Schimel D S, Kirchner T, Menaut J C, Seastedt T, Garcia MoyaE, Kamnalrut A, Kinyamario J I. Observations and modeling of biomass and soil organic matter dynamics for the grassland biome worldwide. Global Biogeochemical Cycles, 1993, 7(4): 785- 809.

[15] Fang J Y, Piao S L, Field C B, Pan Y D, Guo Q H, Zhou L M, Peng C H, Tao S. Increasing net primary production in China from 1982 to 1999. Frontiers in Ecology and the Environment,2003, 1(6): 293- 297.

[16] 朴世龙, 方精云, 郭庆华. 利用CASA模型估算我国植被净第一性生产力. 植物生态学报, 2001, 25(5): 603- 608.

[17] 朴世龙, 方精云, 贺金生, 肖玉. 中国草地植被生物量及其空间分布格局. 植物生态学报, 2004, 28(4): 491- 498.

[18] Jackson R B, Canadell J, Ehleringer J R, Mooney H A, Sala O E, Schulze E D. A global analysis of root distributions for terrestrial biomes. Oecologia, 1996, 108(3): 389- 411.

[19] 徐敏云, 李培广, 谢帆, 黄顶, 王芳, 晏宗明, 王堃. 土地利用和管理方式对农牧交错带土壤碳密度的影响. 农业工程学报, 2011, 27(7): 320- 325.

[20] 张芳, 熊黑钢, 安放舟, 夏倩柔. 基于盐(碱)生植被盖度的土壤碱化分级. 土壤学报, 2012, 49(4): 665- 672.

[21] 张芳, 熊黑钢, 田源, 栾福明. 区域尺度地形因素对奇台绿洲土壤盐渍化空间分布的影响. 环境科学研究, 2011, 24(7): 731- 739.

[22] 夏倩柔, 熊黑钢, 张芳. 绿洲不同尺度地下水时空动态变化特征及成因分析. 新疆农业科学, 2013, 50(6): 1137- 1144.

[23] 方精云, 郭兆迪, 朴世龙, 陈安平. 1981—2000年中国陆地植被碳汇的估算. 中国科学 D辑: 地球科学, 2007, 37(6): 804- 812.

[24] 柴曦, 梁存柱, 梁茂伟, 韩伟华, 李智勇, 苗百岭, 王炜, 王立新. 内蒙古草甸草原与典型草原地下生物量与生产力季节动态及其碳库潜力. 生态学报, 2014, 34(19): 5530- 5540.

[25] 中国科学院新疆综合考察队, 中国科学院土壤研究所. 新疆土壤地理. 北京: 科学出版社, 1965: 55- 56.

[26] Mitra S, Wassmann R, Vlek P L G. An appraisal of global wetland area and its organic carbon stock. Current Science, 2005, 88(1): 25- 35.

[27] 郑姚闽, 牛振国, 宫鹏, 戴永久, 上官微. 湿地碳计量方法及中国湿地有机碳库初步估计. 科学通报, 2013, 58(2): 170- 180.

[28] 牟长城, 王彪, 卢慧翠, 包旭, 崔巍. 大兴安岭天然沼泽湿地生态系统碳储量. 生态学报, 2013, 33(16): 4956- 4965.

[29] 霍莉莉. 沼泽湿地垦殖前后土壤有机碳垂直分布及其稳定性特征研究[D]. 长春: 中国科学院研究生院(东北地理与农业生态研究所), 2013.

[30] 王涛, 徐明. 天山北麓东段小河流域平原区浅层地下水超采区评价. 地下水, 2016, 38(5): 37- 39.

[31] 宋文娟, 熊黑钢, 戚冉. 新疆奇台县近42a气候变化特征分析. 干旱区资源与环境, 2009, 23(4): 100- 106.

[32] Houghton R A, Davidson E A,Woodwell G M. Missing sinks, feedbacks, and understanding the role of terrestrial ecosystems in the global carbon balance. Global Biogeochemical Cycles, 1998, 12(1): 25- 34.

[33] Siegenthaler U, Sarmiento J L. Atmospheric carbon dioxide and the ocean. Nature, 1993, 365(6442): 119- 125.

[34] Tans P P, Fung I Y, Takahashi T. Observational constraints on the global atmospheric CO2budget.Science, 1990, 247(4949): 1431- 1438.

[35] Houghton R A. Changes in the storage of terrestrial carbon since 1850//Lal R, Kimble J, Levine E, Stewart B A,eds. Soils and Global Change. BocaRaton:Lewis Publishers, 1995:45- 65.

[36] 徐小锋, 宋长春. 全球碳循环研究中“碳失汇”研究进展. 中国科学院研究生院学报, 2004, 21(2): 145- 152.

[37] 王效科, 白艳莹, 欧阳志云, 苗鸿. 全球碳循环中的失汇及其形成原因. 生态学报, 2002, 22(1): 94- 103.

ChangesinthecarbonpoolpatternandcarbonpoolmigrationduringtheexpansionofnewartificialoasesinthenorthernfoothillsofTianshanMountainswithinthelast40years

ZHANG Fang1,3, XIONG Heigang2,3,*, FENG Juan1,3, XU Zhonglin1,3

1CollegeofResources&EnvironmentScience,XinjiangUniversity,Urumqi830046,China2CollegeofArtandScience,BeijingUnionUniversity,Beijing100083,China3KeyLaboratoryofOasisEcology(XinjiangUniversity)MinistryofEducation,Urumqi830046,China

Based on ground survey data, meteorological data, and remote sensing data, this study analyzed the impact of regional vegetation and soil from the rapid expansion of new artificial oases in the northern foothills of Tianshan Mountains in the past 40 years on carbon pool pattern and migration and carbon sequestration.The results are as follows: from 1976 to 2016, the proportion of farmland in the study area increased from 3.25% to more than 40%. As the marsh disappeared, the reservoir dried up, the shrub cover was drastically reduced, the soil salinization process ceased, and the degree of salinization gradually decreased. Finally, new grassland has formed over bare saline-alkali soil, and there has been an approximate 70% reduction in the area of naked saline alkali over the last 40 years. The effective carbon sink of vegetation was decreased by 15.7% compared with that in 1976. The soil(0—60cm) total carbon storage in the study area in 1976 was 7.814TgC, of which, farmland only accounted for 4.2%. In 2016, the farmland soil organic carbon storage accounted for more than 50%, and the total soil carbon storage decreased slightly by 1.2% compared with that in 1996. The soil carbon pool appears to have stopped declining and has achieved a new balance. The soil carbon pools of bare soil salinity increased by 1.265TgC during the transfer of soil carbon to farmland, shrub land, and grassland soil carbon pools. The carbon stocks of grassland soils increased by 0.8934TgC during migration to the farmland and shrub land carbon pool. If the existing scale of cultivated land is well controlled, the loss of the soil carbon pool caused by the 40-year reclamation can be gradually recovered and the carbon sink effect can be re-formed. Human activities are decisive factors that influence the pattern of NPP distribution in the study area. Reclamation activities are generally carried out in arid oasis, especially in the desert intertidal zone around the old oasis during 1976 to 2016. Once oases are reclaimed and farmland extends to desert, a new artificial oasis is dramatically expanded, with a drastic change in water consumption patterns and temporal and spatial distributions, causing significant changes in the ecological environment. The overdraft of groundwater leads to the disappearance of spring overflow and the drying of reservoirs, which has negative effects on the ecological environment in the research area. Therefore, the sustainable development of water resources and future development trend of ecosystem carbon pool feedback in the process of new artificial oasis expansion should be taken seriously and demands further studies.

expansion of new artificial oases; carbon pool pattern; carbon pool migration; northern foothills of Tianshan Mountains

国家自然科学基金(41261049,41761041,41671198,41361098)；中国博士后基金面上资助(2013M532100)

2017- 03- 31;

2017- 07- 11

*通讯作者Corresponding author.E-mail: xhg1956@sohu.com

10.5846/stxb201703310565

张芳, 熊黑钢, 冯娟, 许仲林.天山北麓新人工绿洲扩张过程中的碳库格局变化与碳库迁移.生态学报,2017,37(23):7862- 7871.

Zhang F, Xiong H G, Feng J, Xu Z L.Changes in the carbon pool pattern and carbon pool migration during the expansion of new artificial oases in the northern foothills of Tianshan Mountains within the last 40 years.Acta Ecologica Sinica,2017,37(23):7862- 7871.