东北东部森林碳循环过程的集水区尺度模拟

【摘要】大气中C02浓度的升高引发了一系列严重的环境问题,使得碳循环研究受到了空前的重视。中高纬度的东北东部森林生态系统处于全球变化的敏感区域,在维护区域生态环境中起到重要作用,是碳循环研究的重点区域。利用模型模拟碳循环过程是碳循环研究的重要方法,其中IBIS(IntegratedBiosphereSimulator)模型是森林碳循环模拟的有利工具,代表着全球碳循环模型的发展方向。利用IBIS模型模拟东北东部森林生态系统碳循环对陆地碳循环和全球变化研究具有重要意义。为使IBIS模型更好地体现碳循环过程的空间异质性,适用于东北东部集水区尺度的碳循环模拟,本研究融合GIS、VisualFoxPro、VisualC#等建模技术,改进IBIS模型。在气象、地形、植被和土壤等数据的支持下,利用改造后的IBIS模型模拟张家沟集水区5种森林类型的碳循环过程,以实测数据对IBIS模拟结果进行验证,分析集水区碳通量的时空格局及其随地形因子的变化。主要结果和结论如下：(1)IBIS模型改进。由于IBIS模型只能模拟离散的样点(plot),无法针对连续的植被空间开展模拟,没有考虑地形对土壤水分流动、太阳辐射差异的影响,忽略了森林碳动态的空间异质性,造成森林碳循环模拟机理上欠缺。因此,本研究进一步完善IBIS模型机理,应用GIS技术,以一个集水区作为模拟尺度,以一个栅格单元作为最小模拟单位,增加了IBIS模型的地形分析模块,改进了土壤水分再分配和地表太阳辐射计算模块,筛选了适合中国东北地区的8种植被功能型(PFTs),细化了模型物候算法和模型参数,优化了模型编程语言和数据存储格式。使得改进后的IBIS模型能够利用实际地貌进行模拟,机理性更强,实现了集水区尺度的连续植被空间碳水过程的细致模拟,充分体现不同地形条件下森林碳循环过程的空间异质性,为集水区尺度的森林生态系统碳循环模拟提供了有利工具,为今后碳循环研究的尺度上推奠定了坚实基础。(2)IBIS模型运行系统的实现。应用VisualC#编程语言和VisualFoxPro数据库技术,实现了IBIS模型图形界面运行系统。该系统可以推广到全国其他地区模拟运算,主要功能有：模型输入参数的读取、编辑与管理,模型模拟运算,模型输出结果的设置、读取与导出,以及完善的系统出错提示功能等。总体来看,该系统具有结构清晰、界面美观、使用方便等特点。(3)IBIS模型验证。本研究从土壤呼吸、净初级生产力(NPP)两个方面,采用IBIS模型模拟数据与实测数据对比的方法验证模型模拟精度。结果表明：改造后的IBIS模型模拟效果较好,模拟结果基本可靠,适用于东北东部集水区森林碳循环的估算。(4)土壤呼吸模拟。研究区土壤呼吸年平均值为571gC/(m2·a)。年土壤呼吸的空间格局与生长季土壤呼吸空间格局相似,均表现为高值区主要分布在北部、西南和东南区域,低值区主要分布在沟谷附近,该格局与集水区的地形、植被及其组合等因素有关。研究区土壤呼吸的季节性变化呈单峰曲线形式,土壤呼吸峰值均出现在7月,其中落叶松林峰值最低,为85.5gC/m2,杂木林峰值最高,为146.3gC/m2。研究区年总土壤呼吸为865.44tC/a,杂木林分布面积最大,年平均土壤呼吸值最高(628gC/(m2·a)),年总土壤呼吸量为486.39tC/a；落叶松林分布面积最小,年平均土壤呼吸值最低(356gC/(m2·a)),年总土壤呼吸量为0.29tC/a。5种森林类型的土壤呼吸值与5cm深土壤温度存在极显著的指数关系,与土壤湿度的相关性较低,土壤温度的变化可以解释土壤呼吸约70%的季节变化。从地形因素看,坡度和坡向对土壤呼吸的影响大于海拔对土壤呼吸的影响,这与研究区的整体海拔不高且高差较小有直接关系。(5)净初级生产力(NPP)和净生态系统生产力(NEP)模拟。研究区森林NPP年平均值为375gC/(m2·a),森林NPP在空间上大体呈现为东北和西北两区域高,然后逐渐向中心沟谷地区减少；生长季中,不同月份的NPP空间格局变化显著。森林NPP季节变化与气温、降水量和太阳辐射等密切相关,NPP月均值变化曲线呈双峰形式,峰值出现在5月和8月。研究区森林NPP年总量为541.72tC/a,杂木林NPP年总量为260.00tC/a,占研究区NPP总量的48%。各植被类型中,杨桦林的年平均NPP最高(405gC/(m2·a)),杂木林最低(336gC/(m2·a))。从NPP与气候因子的关系看,NPP与降水量的相关性(R2=0.6657,P<0.01)高于与温度的相关性(R2=0.5383,P<0.05),说明降水量是决定集水区NPP值的主要限制因子。从地形因素看,海拔和坡度对森林NPP影响不大,而坡向对森林NPP影响最明显,阳坡的森林NPP较阴坡高41%。集水区NEP年平均值为93gC/(m2-a),东南区域是NEP最低的区域,中心区是NEP相对较高的区域。集水区年总NEP为81.19tC/a,杂木林和杨桦林所占面积最大,它们的NEP约占研究区NEP总量的58%。5种森林类型中,落叶松林的NEP值最高,碳汇能力最强,杂木林的NEP值最低,碳汇能力最弱。
【作者】郭丽娟；
【导师】国庆喜；
【作者基本信息】东北林业大学，生态学，2013，博士
【关键词】IBIS模型；集水区；土壤呼吸；净初级生产力；净生态系统生产力；时空格局；

【参考文献】
[1]李宏平.南高丛越橘品种“V_3”（Vaccinium austral）四倍体诱导及鉴定[D].西南大学,园林植物与观赏园艺,2013,硕士.
[2]李鹏.数字媒体内容管理平台设计与实现[D].电子科技大学,软件工程（专业学位）,2012,硕士.
[3]倪志伟,蔡庆生.范例推理中的知识发现技术[J].小型微型计算机系统,2002,02:159-162.
[4]汤学亮.适用于低功耗广覆盖网络下的同步方案设计[D].北京交通大学,2015.
[5]王玲莉.区域经济视野下高等职业教育专业建设[D].湖南农业大学,高等教育学,2012,硕士.
[6]程鸿.数理统计在商业保险发展中的应用[J].保险职业学院学报,2005,05:29-33.
[7]郭斐.建立第三方C-C电子商务信用评级体系：关于必要性与可行性的研究[D].北京邮电大学,产业经济学,2013,硕士.
[8]李伟荣.AZ31B镁合金搅拌摩擦焊焊缝微弧氧化处理的研究[D].兰州理工大学,材料加工工程,2012,硕士.
[9]孙良,刘海龙,潘德惠.证券投资决策的自由边界问题[J].控制与决策,1999,04:57-61.
[10]王磊,王文云,王乔.基于软件构件技术的KDD算法的设计与实现[J].计算机工程与设计,2005,09:2467-2469.
[11]张娟.小学随班就读外籍学生管理现状及对策研究[D].辽宁师范大学,教育管理（专业学位）,2012,硕士.
[12]陈建军,杨振山.基于智能体和数据仓库的群决策支持系统[J].微机发展,2004,01:1-3+65.
[13]娄佳佳.基于安装误差敏感性的直齿锥齿轮齿面设计与实验分析[D].河南科技大学,机械设计及理论,2014,硕士.
[14]胡飞飞.江苏省典型铅酸蓄电池企业职业危害现状及对周边环境影响调查[D].南京医科大学,劳动卫生与环境卫生学,2013,硕士.
[15]王冰瑞.压电复合材料界面特性理论建模及扩展有限元模拟实现研究[D].西南交通大学,机械设计及理论,2014,硕士.
[16]汪强.基于uClinux的嵌入式SIP服务器的研发[D].浙江工业大学,计算机应用技术,2004,硕士.
[17]吴兴源.拉萨地体晚古生代早期岩浆岩成因及其对松多特提斯洋开启的启示[D].中国地质大学（北京）,矿物学、岩石学、矿床学,2014,硕士.
[18]刘传虎,刘福贵,王超明.巷道水平井的地震设计方法[J].石油地球物理勘探,1993,01:105-111+126.
[19]邝栗山.AZ31B和AZ91D镁合金低温力学性能研究[D].哈尔滨工业大学,空间材料与加工,2013,硕士.
[20]葛恬.文化话语的视域下基于中美对达尔富尔问题言论的考察[D].浙江大学,法语语言文学,2014,硕士.
[21]张海琴.南通市崇川区小学“体育、艺术2+1项目”中体育项目开展的现状研究[D].苏州大学,体育教育训练学（专业学位）,2014,硕士.
[22]贾萌君.移动电商消费者参与度在购买决策中的影响[D].首都经济贸易大学,工商管理（专业学位）,2014,硕士.
[23]葛少博.基于压缩感知理论的智能机器人视觉图像的重构算法研究[D].兰州理工大学,测试计量技术及仪器,2014,硕士.
[24]李秋嘉.再创造数学教学模式的实践[D].东北师范大学,教育,2003,硕士.
[25]胡志恒.多速率多射频无线Mesh网降低组播延时技术研究[D].哈尔滨工业大学,计算机技术,2013,硕士.
[26]田野.废旧塑料浮选过程中气泡粘附行为的可视化研究[D].长安大学,环境工程,2013,硕士.
[27]张驰.300MW核电主泵维护装置及工艺关键技术研究[D].浙江工业大学,机械工程,2011,硕士.
[28]周杰.青藏高原高寒草甸冻融退化过程中群落及化学计量特征的研究[D].兰州大学,农业推广（专业学位）,2013,硕士.
[29]徐赢光.新兴外来词反映下的日本变化[D].吉林大学,外国语言学及应用语言学,2014,硕士.
[30]陈金川.我国应高起点建设物流信息系统[N].中国交通报,2003-01-14.
[31]唐秋旻,黄淑伟,邓甦.基于知识发现的高校科研决策支持系统[J].辽宁大学学报(自然科学版),2004,02:153-155+179.
[32]陈雪波,李树生.多变量系统稳定性分析[J].自动化学报,1987,05:388-392.
[33]孟迪.多孔铜/碳纳米相复合负极的制备及其电化学性能的研究[D].天津大学,2012.
[34]张杲辉.亚波长金属阵列结构的太赫兹透射特性模拟和实验研究[D].首都师范大学,光学,2013,硕士.
[35]田孝武.VDSL系统的设计与实现[D].天津大学,通信与信息系统,2004,硕士.
[36]杨学兰.氨甲喋呤联合中药保守治疗异位妊娠的临床研究[D].山西医科大学,中西医结合妇科,2003,硕士.
[37]颜雷.农村人力资源开发中的政府作用探析[D].苏州大学,行政管理,2012,硕士.
[38]方美娥.插值问题的统一求解算法及复杂曲面B样条重建[D].湖南师范大学,计算数学,2004,硕士.
[39]刘建涛.半导体量子环动力学生长分析及电子结构特性研究[D].北京邮电大学,物理电子学,2013,硕士.
[40]贾瑾.基于空气质量数据解析大气复合污染时空特征及过程序列[D].浙江大学,2014.
[41]廖雪义.固定化微生物捷径生物脱氮的研究[D].广西大学,微生物学,2004,硕士.
[42]李晓辉.声乐套曲《天鹅之歌》演唱研究[D].天津音乐学院,音乐表演,2014,硕士.
[43]刘丹.船用导航雷达视频信号处理系统研究[D].西安科技大学,通信与信息系统,2014,硕士.
[44]岳宁.失肌皮神经损伤胸大肌部分移位重建屈肘术后供受肌内MuRF1和MAFbx的表达变化[D].遵义医学院,人体解剖与组织胚胎学,2013,硕士.
[45]周艳利.Nrf2在溃疡性结肠炎患者中的表达及其与氧化应激的关系[D].河北医科大学,内科学,2013,硕士.
[46]邱佳慧,陶成,刘留,谈振辉.U型槽无线信道多径传播特性测量与建模方法的研究[J].铁道学报,2014,01:40-48.
[47]张宁宁.文化安全背景下美剧的意识形态问题研究[D].安徽大学,马克思主义基本原理,2014,硕士.
[48]张爱英.~(99m)Tc-MIBI、~(18)F-FDG及~(18)F-FLT摄取率变化早期评价耐药人乳腺癌细胞化疗疗效的实验研究[D].苏州大学,影像医学与核医学,2013,硕士.
[49]王金涛.长江流域棉区红铃虫抗性基因频率检测[D].华中农业大学,农业昆虫与害虫防治,2014,硕士.
[50]张鹏.改良Woodward术治疗不同年龄先天性高肩胛症[D].山西医科大学,骨外科学（专业学位）,2013,硕士.