留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于系统动力学的校园碳收支测算模型及应用

胡玥昕 江洪 王颖

胡玥昕, 江洪, 王颖. 基于系统动力学的校园碳收支测算模型及应用[J]. 浙江农林大学学报, 2014, 31(6): 850-859. doi: 10.11833/j.issn.2095-0756.2014.06.005
引用本文: 胡玥昕, 江洪, 王颖. 基于系统动力学的校园碳收支测算模型及应用[J]. 浙江农林大学学报, 2014, 31(6): 850-859. doi: 10.11833/j.issn.2095-0756.2014.06.005
HU Yuexin, JIANG Hong, WANG Ying. System dynamics modeling for carbon budget at a university campus[J]. Journal of Zhejiang A&F University, 2014, 31(6): 850-859. doi: 10.11833/j.issn.2095-0756.2014.06.005
Citation: HU Yuexin, JIANG Hong, WANG Ying. System dynamics modeling for carbon budget at a university campus[J]. Journal of Zhejiang A&F University, 2014, 31(6): 850-859. doi: 10.11833/j.issn.2095-0756.2014.06.005

基于系统动力学的校园碳收支测算模型及应用

doi: 10.11833/j.issn.2095-0756.2014.06.005
基金项目: 

国家自然科学重大基金资助项目 61190114

国家重点基础研究发展规划("973"计划)项目 2011CB302705

国家重点基础研究发展规划("973"计划)项目 2011CB302705

国家重点基础研究发展规划("973"计划)项目 2010CB950702

国家重点基础研究发展规划("973"计划)项目 2010CB428503

国家自然科学基金资助项目 41171324

上海市战略性新兴产业重大项目 重大2013-14号

详细信息
    作者简介: 胡玥昕, 从事生态系统碳循环研究。E-mail:wendy-hoo@hotmail.com
    通信作者: 江洪, 教授, 博士, 博士生导师, 从事生态系统碳循环研究。E-mail:jianghong_china@hotmail.com
  • 中图分类号: S718.5

System dynamics modeling for carbon budget at a university campus

  • 摘要: 以系统动力学VensimPLE软件为技术平台, 结合Li COR-6400与Li COR-8100观测数据, 以浙江农林大学东湖校区为研究对象, 建立校园碳收支测算模型, 对校园中人类活动系统碳排放与自然生态系统碳收支分别进行了测算, 并提出低碳建设参考路径。模型测算结果表明:校园全年整体表现为碳源, 年碳排放量为4 147.48 t·a-1, 平均单位面积碳排放量为3.24 t·m-2·a-1; 人类活动碳排放量5 555.45 t·a-1, 人均碳排放强度为24.0 kg·人-1·月-1, 其中, 生活办公用电排碳量最大。自然生态系统年固碳量1 407.96 t·a-1, 季节变化呈现双峰特征, 春秋两季较高, 冬季最低。单位绿化面积碳吸收量为187.0 g·m-2·月-1。同时, 低碳校园情景模拟中通过绿地碳汇、人均用能与交通用能等因素的共同控制, 可比原始状态再实现29%的碳排放削减, 说明校园低碳建设仍有一定发展空间。
  • 图  1  浙江农林大学校园航拍遥感影像与影像分类结果

    Figure  1  Aerial remote sensing image and classification results

    图  2  校园碳收支系统动力学模型详细流图

    Figure  2  The stock-flow diagram of campus carbon budget

    图  3  人类活动系统碳排放总量及组成成分

    Figure  3  Carbon emission and composition of human activity system

    图  4  人类活动系统人均碳排放量

    Figure  4  Per capital carbon emission in human activity system

    图  5  校园草地碳收支季节变化

    Figure  5  Seasonal changes of carbon budget of grassland in campus

    图  6  校园林地碳收支季节变化

    Figure  6  Seasonal changes of carbon budget of woodland in campus

    图  7  校园自然生态系统碳汇季节变化特征

    Figure  7  Monthly carbon budget of natural ecosystem in campus

    图  8  植物光合速率季节变化

    Figure  8  Monthly average plant photosynthetic rate of natural ecosystem in campus

    图  9  植物与土壤呼吸速率季节变化

    Figure  9  Monthly average vergetable and soil respiration rate of natural ecosystem in campus

    图  10  校园碳收支季节变化

    Figure  10  Monthly carbon budget of whole campus

    表  1  浙江农林大学校园土地覆盖分类结果

    Table  1.   Land cover classification results

    土地覆盖类型 面积/m2 百分比/%
    林地 154 755.0 12.70
    草地 471 856.0 38.72
    水域 52 107.8 4.28
    建筑用地 539 769.0 44.30
    合计 1218 487.8
    下载: 导出CSV

    表  2  自然生态系统参数拟合结果

    Table  2.   Parameterized equations in natural ecosystem

    内容 公式 精度R2
    林地植物光合速率 P-wood=∫11(-0.000 036×Par2 +0.023 0×Par+2.16)dt 0.809
    草地植物光合速率 P-wood=∫11(-0.000 056×Par2 +0.001 6×Par+8.46)dt 0.759
    林地植物呼吸速率 VR-wood=∫11 (0.42Swe0.051×f(T, h, PAR, σ) dt 0.937
    草地植物呼吸速率 VR-grass=∫11 (1.99Sge0.041×f(T, h, PAR, σ) dt 0.951
    林地土壤呼吸速率 SR-wood=∫11 (0.87Swe0.055×f(T, h, PAR, σ) dt 0.958
    草地土壤呼吸速率 SR-grass=∫11 (1.52Sge0.06×f(T, h, PAR, σ) dt 0.958
    说明:T为温度,PAR为光合有效辐射,σ为二氧化碳浓度,t为时常,Sw为林地面积,Sg为草地面积。
    下载: 导出CSV

    表  3  碳排放组成成分贡献度

    Table  3.   Contribution of each kind of carbon emission's composition

    时间阶段 Cp /% Cfood /% Ctra /% Clife /%
    工作日 32.67 6.48 1.13 59.72
    寒暑假 26.00 11.92 5.81 56.27
    下载: 导出CSV

    表  4  低碳情景与原状态碳收支结果比较

    Table  4.   Comparation between low-carbon scenario and original state

    状态 基础数据 计算结果
    绿地率/% 单位绿化面积碳汇/(g·m-2 人均生活用能碳排放/(g·人-1 ·d-1 日均交通碳排放/(kg·d-1 绿地碳汇/(t·a-1 生活用能碳排放(包括锅炉)/(t·a-1 交通碳排放/(t·a-1 碳收支/(t·a-1 人均碳排放量/(g·人-1 ·d-1
    原状态 49 187 535.69 203.8 -1 408 5 555.48 73.35 4 147.48 439.27
    低碳情景 60 224 428.55 163.0 -2 066 3 138.95 58.68 2 929.44 310.27
    下载: 导出CSV
  • [1] 何华.华南居住区绿地碳汇作用研究及其在全生命周期碳收支评价中的应用[D]:重庆:重庆大学, 2010.

    HE Hua. Research on Green Space Carbon Sink of Residential Area in South China and its Application in Carbon Budget of Residential Area for Whole Life Circle[D]. Chongqing:Chongqing University, 2010.
    [2] 马巾英, 尹锴, 吝涛.城市复合生态系统碳氧平衡分析[J].环境科学学报, 2011, 31(8):1808-1816.

    MA Jinying, YI Kai, LIN Tao. Analysis of the carbon and oxygen balance of acomplex urban ecosystem:a case study in the coastal city of Xiamen[J]. Acta Sci Circumst, 2011, 31(8):1808-1816.
    [3] FENG Yueyi, CHEN Shaoqing, ZHANG Lixiao. System dynamics modeling for urban energy consumption and CO2 emis-sions:a case study of Beijing, China[J]. Ecol Model, 2012, 252(1):44-52.
    [4] 周志.基于系统动力学的广东省低碳经济发展路径选择[D].广州:华南理工大学, 2011.

    ZHOU Zhi. The Choosing of Low-carbon Economy Development Path of Guangdong Province Based by System Dynamics Theory[D]. Guangzhou:South China University of Technology, 2011.
    [5] ZHAO Min, KONG Zhenghong, ESCOBEDO F J, et al. Impacts of urban forests on offsetting carbon emissions from industrial energy use in Hangzhou, China[J]. J Environ Manage, 2010, 91(4):807-813.
    [6] 李霞, 孙睿, 李远, 等.北京海淀公园绿地二氧化碳通量[J].生态学报, 2010, 30(24):6715-6725.

    LI Xia, SUN Rui, LI Yuan, et al. Carbon dioxide fluxes on green space in Haidian Park, Beijing[J]. Acta Ecol Sin, 2010, 30(24):6715-6725.
    [7] 吴飞, 吴泽民, 吴文友.马鞍山市城市森林碳储量及分布格局研究[J].安徽农业大学学报, 2012, 39(4):519-526.

    WU Fei, WU Zemin, WU Wenyou. Appraisal of carbon storage in urban forest patches and its distribution pattern in Maanshan city[J]. J Anhui Agric Univ, 2012, 39(4):519-526.
    [8] 王其藩.系统动力学[M].北京:清华大学出版社, 1994.
    [9] 赵荣钦, 黄贤金.城市系统碳循环:特征, 机理与理论框架[J].生态学报, 2013, 33(2):358-366.

    ZHAO Rongqin, HUANG Xianjin. Carboncycle of urbansystem:characteristics, mechanism and theoretical framework[J]. Acta Ecol Sin, 2013, 33(2):358-366.
    [10] 李雪玉, 姚林如, 杨锐琦, 等.城市碳循环清单及相关对策研究:以南昌市为例[J].全国商情·理论研究, 2012(20):12-14.

    LI Xueyu, YAO Linru, YANG Ruiqi, et al. Research on urban carbon cycle list and relevant countermeasures[J]. Business Inf Theor Res, 2012(20):12-14.
    [11] 何华, 龙天渝, 鄢涛, 等.居住区使用阶段碳收支研究[J].煤气与热力, 2010, 30(4):40-44.

    HE Hua, LONG Tianyu, YAN Tao, et al. Research on carbon budget during usestage of residential areas[J]. Gas & Heat, 2010, 30(4):40-43.
    [12] 蒋运生, 柴胜丰, 唐辉, 等.光照强度对广西地不容光合特性和生长的影响[J].广西植物, 2009, 29(6):792-796.

    JIANG Yunsheng, CHAI Shengfeng, TANG Hui, et al. Effect of different light intensities on photosynthetic characteristics and growth of Steohania kwangsiensis[J]. Guihaia, 2009, 29(6):792-796.
    [13] 张光飞, 翟书华, 苏文华, 等.光照和温度对鸟巢蕨光合速率的影响研究[J].昆明学院学报, 2008, 30(4):62-63.

    ZHANG Guangfei, ZHAI Shuhua, SU Wenhua, et al. Effects of light intensity and air temperature on the photosynthesis of Neottopters nidus[J]. J Kunming Univ, 2008, 30(4):62-63.
    [14] 陈柯, 王小德.常春油麻藤、中华常春藤和扶芳藤的光合特性比较分析[J].安徽农业大学学报, 2008, 35(2):196-199.

    CHEN Ke, WANG Xiaode. Comparative analysis on photosynthetic characteristics and shade tolerance between Mucuna sempervirens, Hedera nepalensis var. sinensi and Euonymus fortunei[J]. J Anhui Agric Univ, 2008, 35(2):196-199.
    [15] 王红兵, 王励勉, 秦俊, 等.上海地区九种木兰科植物光合特性的研究[J].北方园艺, 2008(8):54.

    WANG Hongbing, WANG Limian, QIN Jun, et al. Research on photosynthesis characteristic of nine species of Magnoliaceae plants in Shanghai[J]. Northern Hortic, 2008(8):54.
    [16] 穆大刚, 付翔, 许华林. 5种红树植物光合产量模型研究[J].安徽农业科学, 2008, 36(35):1535-1535.

    MU Dagang, FU Xiang, XU Hualin. Study on phtotsynthetic yield model for five Mangrove species[J]. J Anhui Agric Sci, 2008, 36(35):1535-1535.
    [17] 陈磊, 潘青华, 金洪.温湿度对紫叶黄栌光合特性变化的影响[J].中国农学通报, 2008, 24(6):124-128.

    CHEN Lei, PAN Qinghua, JIN Hong. Research on influence of relative humidity and air temperature on photosynthetic characteristics of Cotinus coggygria prupureus[J]. Chin Agric Sci Bull, 2008, 24(6):124-128.
    [18] 郭新想.居住区绿化植物固碳能力评价方法研究[D].重庆:重庆大学, 2010.

    GUO Xinxiang. Study on Green Plants CO2-fixing Capacity Evaluation Methods in Residential Area[D]. Chongqing:Chongqing University, 2010.
    [19] 陈志银.红麻光合、呼吸强度与气象因子关系的研究[J].中国农业气象, 1990, 11(4):29-32.

    CHEN Zhiyin. Research on the relationship between photosynthesis, respiration intensity of Kenaf and meteorological factors[J]. Chin J Agrometeorol, 1990, 11(4):29-32.
    [20] 王淼, 刘亚琴, 郝占庆, 等.长白山阔叶红松林生态系统的呼吸速率[J]. 2006, 17(10):1789-1795.

    WANG Miao, LIU Yaqin, HAO Zhanqing, et al. Respiration rate of broad leaved Korean pine fo rest ecosystem in Changbai Mountains[J]. Chin J Appl Ecol, 2006, 17(10):1789-1795.
    [21] 崔骁勇, 陈佐忠, 陈四清.草地土壤呼吸研究进展[J].生态学报, 2001, 21(2):315-325.

    CUI Xiaoyong, CHEN Zuozhong, CHEN Siqing. Progress in research on soil respiration of grasslands[J]. Acta Ecol Sin, 2001, 21(2):315-325.
    [22] 杜丽君.红壤几种典型利用方式CO2的排放规律及其影响因素[D].武汉:华中农业大学, 2006.

    DU Lijun. CO2 Emission Rule and its Influencing Factors of Red Soil in Different Use Types[D]. Wuhan:Huazhong Agricultural University, 2006.
    [23] 刘绍辉, 方精云, 清田信.北京山地温带森林的土壤呼吸[J].植物生态学报, 1998, 22(2):119-126.

    LIU Shaohui, FANG Jingyun, QING Tianxin. Soil respritation of mountainous temperate forests in Beijing, China[J]. Acta Phytoecol Sin, 1998, 22(2):119-126.
    [24] 张宪权.东北地区落叶松人工林土壤呼吸的时空异质性研究[D].上海:华东师范大学, 2005.

    ZHANG Xianquan. Larch Plantation's Spatial and Temporal Heterogeneity of Soil Respiration in Northeast China[D]. Shanghai:East China Normal University, 2005.
    [25] KUCERA C L, KIRKHAM D R. Soil respiration studies in tallgrass prairie in Missouri[J]. Ecology, 1971, 52(5):912-915.
    [26] 宋永昌, 由文辉, 王祥荣.城市生态学[M].上海:华东师范大学出版社, 2000.
    [27] 李志鹏.基于系统动力学的城市交通能源消耗与碳排放预测[D].天津:天津大学, 2012.

    LI Zhipeng. The Forecast of the Energy Consumption and Carbon Emission of Urban Traffic Based on System Dynamics-An Empirical Analysis on Tianjin[D]. Tianjin:Tianjin University, 2012.
    [28] DHAKAL S. Urban energy use and carbon emissions from cities in China and policy implications[J]. Energ Polic, 2009, 37(11):4208-4219.
    [29] 毛军需.太阳辐射度与光照度换算浅析[J].河南农业科学, 1995(1):11-12.

    MAO Junxu. The analysis of the conversion between solar radiation and light inensity[J]. J Henan Agric Sci, 1995(1):11-12.
    [30] 陈云飞, 江洪, 周国模, 等.人工高效经营雷竹林CO2通量估算及季节变化特征[J].生态学报, 2012, 33(11):3434-3444.

    CHEN Yunfei, JIANG Hong, ZHOU Guomo, et al. Estimation of CO2 fluxes and its seasonal variations from the effective management lei bamboo (Phyllostachys violascens)[J]. Acta Ecol Sin, 2012, 33(11):3434-3444.
  • [1] 裴建川, 张书廷, 杨金艳, 张进.  立体生态模块处理杭州市玉皇山南基金小镇水体氮的效果 . 浙江农林大学学报, 2018, 35(6): 987-996. doi: 10.11833/j.issn.2095-0756.2018.06.001
    [2] 汤旭, 郑洁, 冯彦, 李燕坤, 王时军, 张大红.  云南省县域森林生态安全评价与空间分析 . 浙江农林大学学报, 2018, 35(4): 684-694. doi: 10.11833/j.issn.2095-0756.2018.04.014
    [3] 彭琛琛, 蔡玉梅, 边振兴.  烟台市生态红线划定研究 . 浙江农林大学学报, 2018, 35(5): 818-828. doi: 10.11833/j.issn.2095-0756.2018.05.005
    [4] 金文奖, 侯平, 张伟, 梁立成, 俞飞.  温州鳌江流域表层底泥及河岸土壤重金属空间分布与生态风险评价 . 浙江农林大学学报, 2017, 34(6): 963-971. doi: 10.11833/j.issn.2095-0756.2017.06.001
    [5] 曾琪, 刘健, 余坤勇, 张今朝, 郑文英, 陈樟昊.  福建金丝湾森林公园功能区划与生态脆弱性耦合研究 . 浙江农林大学学报, 2017, 34(5): 833-840. doi: 10.11833/j.issn.2095-0756.2017.05.009
    [6] 高常军, 魏龙, 贾朋, 田惠玲, 李树光.  基于去重复性分析的广东省滨海湿地生态系统服务价值估算 . 浙江农林大学学报, 2017, 34(1): 152-160. doi: 10.11833/j.issn.2095-0756.2017.01.021
    [7] 王帆, 江洪, 牛晓栋.  大气水汽稳定同位素组成在生态系统水循环中的应用 . 浙江农林大学学报, 2016, 33(1): 156-165. doi: 10.11833/j.issn.2095-0756.2016.01.021
    [8] 陈永根, 周传斌, 朱慧芳, 王如松.  发达地区农村固体废弃物管理与资源化策略 . 浙江农林大学学报, 2015, 32(6): 940-946. doi: 10.11833/j.issn.2095-0756.2015.06.018
    [9] 王琳, 景元书, 张悦.  基于MODIS的长江中下游地区植被净第一性生产力时空变化规律 . 浙江农林大学学报, 2015, 32(6): 829-836. doi: 10.11833/j.issn.2095-0756.2015.06.002
    [10] 晏闻博, 柳丹, 彭丹莉, 李松, 陈俊任, 叶正钱, 吴家森, 王海龙.  重金属矿山生态治理与环境修复技术进展 . 浙江农林大学学报, 2015, 32(3): 467-477. doi: 10.11833/j.issn.2095-0756.2015.03.021
    [11] 易武英, 苏维词, 周文龙, 唐金刚, 张凤太.  基于元胞自动机模型的贵阳市花溪区生态安全预警模拟研究 . 浙江农林大学学报, 2015, 32(3): 369-375. doi: 10.11833/j.issn.2095-0756.2015.03.006
    [12] 王适, 吴永波, 陈杰, 丁晓叶.  上海崇明生态农业园区花菜田的碳平衡 . 浙江农林大学学报, 2014, 31(2): 190-195. doi: 10.11833/j.issn.2095-0756.2014.02.005
    [13] 李德会, 李贤伟, 王巧, 荣丽, 杨渺, 刘朔.  林木根系呼吸影响因素及根系呼吸对全球变化的响应 . 浙江农林大学学报, 2007, 24(2): 231-238.
    [14] 李正才, 徐德应, 杨校生, 傅懋毅, 孙雪忠, 奚金荣.  7 种不同林农土地利用类型残体的有机碳储量 . 浙江农林大学学报, 2007, 24(5): 581-586.
    [15] 戴建兵, 俞益武, 曹群.  湿地保护与管理研究综述 . 浙江农林大学学报, 2006, 23(3): 328-333.
    [16] 贺佳飞, 周伟, 李明会.  多疣狭口蛙繁殖生态初步观察 . 浙江农林大学学报, 2006, 23(3): 311-315.
    [17] 付晓萍, 田大伦, 黄智勇.  模拟酸雨对植物形态学效应的影响 . 浙江农林大学学报, 2006, 23(5): 521-526.
    [18] 蒋文伟, 刘彤, 丁丽霞, 温国胜, 张万荣, 钟泰林.  景观生态空间异质性的研究进展 . 浙江农林大学学报, 2003, 20(3): 311-314.
    [19] 张万荣, 方红, 冯维维.  高等学校校园绿化探讨 . 浙江农林大学学报, 2003, 20(2): 177-181.
    [20] 田信桥.  论市场经济与校园文化优化 . 浙江农林大学学报, 1996, 13(4): 467-472.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2014.06.005

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2014/6/850

图(10) / 表(4)
计量
  • 文章访问数:  4093
  • HTML全文浏览量:  517
  • PDF下载量:  478
  • 被引次数: 0
出版历程
  • 收稿日期:  2013-12-03
  • 修回日期:  2014-02-18
  • 刊出日期:  2014-12-20

基于系统动力学的校园碳收支测算模型及应用

doi: 10.11833/j.issn.2095-0756.2014.06.005
    基金项目:

    国家自然科学重大基金资助项目 61190114

    国家重点基础研究发展规划("973"计划)项目 2011CB302705

    国家重点基础研究发展规划("973"计划)项目 2011CB302705

    国家重点基础研究发展规划("973"计划)项目 2010CB950702

    国家重点基础研究发展规划("973"计划)项目 2010CB428503

    国家自然科学基金资助项目 41171324

    上海市战略性新兴产业重大项目 重大2013-14号

    作者简介:

    胡玥昕, 从事生态系统碳循环研究。E-mail:wendy-hoo@hotmail.com

    通信作者: 江洪, 教授, 博士, 博士生导师, 从事生态系统碳循环研究。E-mail:jianghong_china@hotmail.com
  • 中图分类号: S718.5

摘要: 以系统动力学VensimPLE软件为技术平台, 结合Li COR-6400与Li COR-8100观测数据, 以浙江农林大学东湖校区为研究对象, 建立校园碳收支测算模型, 对校园中人类活动系统碳排放与自然生态系统碳收支分别进行了测算, 并提出低碳建设参考路径。模型测算结果表明:校园全年整体表现为碳源, 年碳排放量为4 147.48 t·a-1, 平均单位面积碳排放量为3.24 t·m-2·a-1; 人类活动碳排放量5 555.45 t·a-1, 人均碳排放强度为24.0 kg·人-1·月-1, 其中, 生活办公用电排碳量最大。自然生态系统年固碳量1 407.96 t·a-1, 季节变化呈现双峰特征, 春秋两季较高, 冬季最低。单位绿化面积碳吸收量为187.0 g·m-2·月-1。同时, 低碳校园情景模拟中通过绿地碳汇、人均用能与交通用能等因素的共同控制, 可比原始状态再实现29%的碳排放削减, 说明校园低碳建设仍有一定发展空间。

English Abstract

胡玥昕, 江洪, 王颖. 基于系统动力学的校园碳收支测算模型及应用[J]. 浙江农林大学学报, 2014, 31(6): 850-859. doi: 10.11833/j.issn.2095-0756.2014.06.005
引用本文: 胡玥昕, 江洪, 王颖. 基于系统动力学的校园碳收支测算模型及应用[J]. 浙江农林大学学报, 2014, 31(6): 850-859. doi: 10.11833/j.issn.2095-0756.2014.06.005
HU Yuexin, JIANG Hong, WANG Ying. System dynamics modeling for carbon budget at a university campus[J]. Journal of Zhejiang A&F University, 2014, 31(6): 850-859. doi: 10.11833/j.issn.2095-0756.2014.06.005
Citation: HU Yuexin, JIANG Hong, WANG Ying. System dynamics modeling for carbon budget at a university campus[J]. Journal of Zhejiang A&F University, 2014, 31(6): 850-859. doi: 10.11833/j.issn.2095-0756.2014.06.005
  • 城市生态系统碳收支的测算是全球碳循环研究中的重要问题,也是正在完成的一项挑战,同时,低碳建设也已成为世界热点问题之一[1]。通过构建数学模型测算不同区域碳收支的研究结果表明:碳源远大于碳汇,且呈现增长趋势[2]。现有城市区域碳收支研究结果认为:工业与交通是碳排放最大的2个部门[3-4];城市绿地具有一定碳汇能力[5-6],且碳汇强度表现为山体林>行道树>公园[7]。目前,针对城市生态系统碳收支的测算研究主要为省、市及以上尺度,社区等小尺度研究相对较少,而生物量法、涡度相关法以及遥感反演法在小尺度碳汇研究上均有一定局限。如何测算并动态反映社区尺度碳收支情况是需要进一步突破的问题。系统动力学(system dynamics)认为系统是由单元、单元的运动和信息组成,能定性与定量地分析研究系统,以“结构—功能”模拟为突出特点,从系统的微观结构入手建模,构造系统的基本结构,进而模拟与分析系统的动态行为[8]。本研究以浙江农林大学东湖校区为研究对象,以系统动力学VensimPLE软件为技术平台,提出校园碳收支计算的方法并反映其碳收支变化的特征,为“零碳”和“低碳”社区的建设评价提供了量化工具,又为未来校园以及居住区等的低碳建设提供参考。

    • 模型以浙江农林大学东湖校区为系统边界,覆盖时间为2012年,时间步长以月为单位。研究对象坐落于浙江省临安市,校舍建筑面积近60万m2,在校师生约为27 000人。研究对象具有校园与植物园两重属性,是一座集教学、科研、生物多样性保护等功能于一体的现代化生态校园。校园碳收支系统主要包括校园范围内各项人类活动产生的碳排放及自然生态系统的碳收支等2个方面。

      人类活动子系统以学校的作息规律造成的人员数量变动为基础,具体包括4个模块:人类呼吸模块主要考虑校园内在校师生日常呼吸产生的碳排放;交通模块主要考虑进入校园的车辆行驶所产生的碳排放;生活办公耗能模块主要考虑学生住宿及教学办公消耗电力产生的碳排放;锅炉能耗模块主要考虑学校食堂化石燃料燃烧产生的碳排放。

      自然生态系统以校园内植被、土壤的碳收支活动为基础,具体包括3个模块:植被光合碳吸收估算模块,包含校园内所有植被光合作用净固碳量;植被呼吸碳排放估算模块,包含植被呼吸作用所排碳量;土壤呼吸碳排放估算模块,包含所有土壤呼吸作用所排碳量。

    • 人类作为城市生态系统内最重要的活动对象,其自身的呼吸作用作为城市生态系统自然垂直碳输出的一部分,以二氧化碳(CO2)的形式将碳排入大气[9]。通过人均呼吸碳排放量与人口总量共同计算人类呼吸碳排放量的方法已被广泛使用[2, 9-10]Cp=∫(Cpmon)dt; Cpmon=Cp0×P×t。其中:Cp为人类呼吸碳排量(t);Cpmon为月呼吸碳排量(t);P为在校总人数;t为时长(月平均天数);Cp0为人均呼吸碳排量每天(t·d-1)。交通系统的碳排放与交通量、交通类型、行驶里程数都密切相关。本研究采用何华等[11]对居住区交通碳排放量的测算方法:Ctra=∫(Ctra-m)dtCtra-m=Tra×Ct0×t×ntra×Dtra。其中:Ctra为交通碳排量(t);Ctra-m为月交通碳排量(t);Tra为汽车总数;t为时长(月平均天数);Cto为单位里程(1 km)汽车平均碳排量(kg·km-1);ntra为汽车日均通行校园的次数;Dtra为汽车次均校园内行驶千米数(km)。生活办公能耗主要以电力为主。模型利用电力碳排放系数与用电总量共同计算生活办公模块的碳排放:Clife=∫C(life-m)dt; Clife-m=Elife×Ielc; Elife=Estu+Esta。其中:Clife为生活办公碳排量(t);Clife-m为月生活办公碳排量(t);Elife为耗电量(kW·h);Ielc为电碳排放系数[t·(kW·h)-1];Estu为生活区耗电量(kW·h);Esta为教学办公区耗电量(kW·h)。食堂锅炉采用化石燃料燃烧与电力相组合的能源消耗方式,具体为:Cfood=∫(Cfood-m)dt; Cfood-mEi×Ii。其中:Cfood为锅炉燃烧碳排量(t);Cfood-m为月锅炉燃烧碳排放量(t);Ei为锅炉燃烧各类能耗量;Ii为各类能源碳排放系数。校园内人类活动系统碳排放总量为:Chuman=Cp+Ctra+Clife+Cfood。其中:Chuman为人类活动系统碳排总量。

    • 自然生态系统主要指研究范围内植被及土壤共同组成的绿地系统。本研究建立了逐时植物呼吸碳排放方程、逐时土壤呼吸碳排放方程和逐时植物光合碳吸收方程。植物光合碳吸收测算。物种差异、光强、温度和湿度是影响植物光合作用的主要因素,其中光强、温度和湿度均表现为最适情况下促进光合[12-14],过高或过低均抑制光合[15-17]。植物光合碳吸收量具体计算方式如下:Pij=a×Pari2+b×Par+c; ${C_P} = \sum\limits_i^n {\left( {{P_{ij}} \times {A_{ij}}} \right)} $ 。其中:Pij为第j种植被类型第i时的光合速率(μmol·m-2·s-1);Pari为第i时的光合有效辐射(μmol·m-2·s-1);a, b, c为系数项。Aij为第j种植被类型第i时的叶面积(m2);Cp为植被光合碳吸收总量(t)。叶面积是指该种植物的所有叶片的单面叶面积之和。假设常绿植物全年总叶面积不变,估算方法参考居住区绿地碳汇能力评价中叶面积指数的计算方法[1, 18]。植物呼吸碳排放测算:Rik=ak∫ebktikdt[19-20]; ${C_v} = \sum\limits_i^n {\left( {{R_{ij}} \times {A_{ij}}} \right)} $ 。其中Rij为第j种植被类型第i时的呼吸速率(μmol·m-2·s-1);ti为第i时的温度(℃);ab为系数项;Aij为第j种植被类型第i时的叶面积(m2);Cv为植被光合碳吸收总量(t)。土壤呼吸碳排放测算。土壤呼吸一般表现为土壤向大气净释放二氧化碳[21]。物种差异[22]、温度[23-24]是影响土壤呼吸速率的主要因素。具体计算方式如下:RSik=akebktik[25]; ${C_S} = \sum\limits_i^m {\left( {{R_{{\rm{Sik}}}} \times {A_{{\rm{Sik}}}}} \right)} $ 。其中RSik为第j种植被类型第k时的呼吸速率(μmol·m-2·s-1);ti为第i时的温度(℃);ab为系数项;ASik为第k种植被类型第i时的叶面积(m2);CS为土壤碳排放总量(t)。

    • 人类活动系统中,在校人口总数、车辆通行数量与行驶里数、生活办公用电量、锅炉能源消耗量均由学校相关部门提供。人数调节系数取值0.2。人均呼吸碳排放量为900 g·人-1·d-1[26]。车辆单位里程碳排放量为0.06 kg·km-1[11, 27]。在生活办公与锅炉燃烧模块中,电力与柴油碳排放因子分别为7.72×104 t·(kW·h)-1与0.862 6 t·t-1(二氧化碳/标准煤)[28]

    • ①观测数据。选取校园林地与草地的主流植物种类,以及不同疏密程度与不同环境的绿地,利用Li COR-6400与Li COR-8100分别观测植物光合速率、植物夜间呼吸速率与土壤呼吸速率。同时,记录Li COR-6400与Li COR-8100所测得的温度、湿度与光照强度数据。②校园不同土地利用类型的面积。各种植被类型的面积通过航拍遥感影像分类获得。影像空间分辨率为1 m,选择非监督分类方法ISOData进行分类(图 1),获取林地、草地、水体、建筑用地4种土地利用类型的面积数据(表 1)。③植被叶面积指数。根据乔木与草地的叶面积计算公式,通过校园植被调查数据获取乔木的株数、年龄、冠幅及草地的高度,来计算植被的总叶面积。乔木和草地的叶面积指数平均值分别采用3.41与2.46 [18]

      图  1  浙江农林大学校园航拍遥感影像与影像分类结果

      Figure 1.  Aerial remote sensing image and classification results

      表 1  浙江农林大学校园土地覆盖分类结果

      Table 1.  Land cover classification results

      土地覆盖类型 面积/m2 百分比/%
      林地 154 755.0 12.70
      草地 471 856.0 38.72
      水域 52 107.8 4.28
      建筑用地 539 769.0 44.30
      合计 1218 487.8
    • 校园碳收支估算模型详细流图见图 2。人类活动系统中,人均用电量、汽车通行量等作为辅助常数变量输入模型,模型按照核算公式运行后可直接计算得到人类活动系统碳排放结果。自然生态系统碳收支的计算相对复杂。在运行模型前,需要在整理观测数据的基础上,取所有植物种类(林地与草地分开)和土壤同一时刻光合作用速率、土壤呼吸速率的均值来分别代表校园区域植被光合速率与土壤呼吸速率。利用SPSS统计软件,选取温度、湿度、光合有效辐射(由光照强度数据转换所得[29])等环境因素数据与各类速率数据拟合。拟合结果(表 2)通过统计学检验后作为计算公式输入模型,温度、湿度等环境因素数据作为辅助变量输入模型,光合与土壤呼吸速率为模型运行后输出的中间结果,自然生态系统碳吸收量为最终输出结果,进行自然生态系统碳收支的核算。

      图  2  校园碳收支系统动力学模型详细流图

      Figure 2.  The stock-flow diagram of campus carbon budget

      表 2  自然生态系统参数拟合结果

      Table 2.  Parameterized equations in natural ecosystem

      内容 公式 精度R2
      林地植物光合速率 P-wood=∫11(-0.000 036×Par2 +0.023 0×Par+2.16)dt 0.809
      草地植物光合速率 P-wood=∫11(-0.000 056×Par2 +0.001 6×Par+8.46)dt 0.759
      林地植物呼吸速率 VR-wood=∫11 (0.42Swe0.051×f(T, h, PAR, σ) dt 0.937
      草地植物呼吸速率 VR-grass=∫11 (1.99Sge0.041×f(T, h, PAR, σ) dt 0.951
      林地土壤呼吸速率 SR-wood=∫11 (0.87Swe0.055×f(T, h, PAR, σ) dt 0.958
      草地土壤呼吸速率 SR-grass=∫11 (1.52Sge0.06×f(T, h, PAR, σ) dt 0.958
      说明:T为温度,PAR为光合有效辐射,σ为二氧化碳浓度,t为时常,Sw为林地面积,Sg为草地面积。
    • ①校园人类活动系统2012年碳排放总量约为5 555.48 t图 3)。2,7,8月随着在校人数减少出现骤减现象。在正常工作日,4-6月与10-12月碳排放量均呈现上升趋势。这可能是因为夏天与冬天校园内空调等电器电力消耗增多。②校园人类活动系统中生活办公能耗产生的碳排放量最大(3 312.49 t·a-1);其次为人类呼吸排放(1 797.81 t·a-1);锅炉燃烧与交通产生的碳排放量较小,分别为371.80 t·a-1与73.35 t·a-1表 3)。另外,正常工作时段与寒暑假时段校园碳排放的构成成分的贡献度有轻微差别。在寒暑假时段,人类呼吸碳排放贡献度较正常工作时段下降较多,而交通碳排放与锅炉燃烧碳排放的贡献度较正常工作时段有所上升。这是由于寒暑假时段在校人数骤减程度大于交通与锅炉燃烧碳排放下降程度造成。③校园人均碳排放量约为24 kg·人-1·月-1图 4),其在8月达到最大(31 kg·人-1·月-1,4月最小(21 kg·人-1·月-1)。与碳排放量月变化特征相反,人均碳排放量在寒暑假时段(2,7,8月)较高。4-6月和10-12月期间,人均碳排放量有缓慢上升趋势。

      图  3  人类活动系统碳排放总量及组成成分

      Figure 3.  Carbon emission and composition of human activity system

      表 3  碳排放组成成分贡献度

      Table 3.  Contribution of each kind of carbon emission's composition

      时间阶段 Cp /% Cfood /% Ctra /% Clife /%
      工作日 32.67 6.48 1.13 59.72
      寒暑假 26.00 11.92 5.81 56.27

      图  4  人类活动系统人均碳排放量

      Figure 4.  Per capital carbon emission in human activity system

    • 自然生态系统碳收支由林地、草地2种植被类型决定,林地与草地碳收支中各种作用季节变化特征见图 5图 6。①校园自然生态系统林地碳吸收年累计量达到-258.14 t·a-1,草地碳吸收年累计量为-1 149.82 t·a-1。林地与草地碳吸收量在季节变化上均呈现出双峰特征,春、秋两季较高,林地最高达到-31.82 t·月-1,草地最高达到-168.21 t·月-1。林地与草地碳吸收量均在1月降至最低,分别为0.87 t·月-1与49.32 t·月-1,成为微弱的碳源。校园内林地与草地的光合作用在季节变化上与碳吸收量季节变化特征一致,呈现双峰,且均在春、秋两季表现出较强的光合作用。夏季由于高温环境,光合作用产生“午休”现象,林地与草地光合作用均有所下降。冬季受低温、低光照的环境影响,林地与草地的光合作用均在1月达到全年最低。校园内植物呼吸与土壤呼吸作用在季节变化上呈现单峰特征,夏季7月出现顶峰。受冬季低温、积雪等影响,校园内植物呼吸与土壤呼吸碳排量均在1月达到全年最小值。②校园自然生态系统总碳汇年累计量为-1 407.96 t·a-1,呈现双峰(图 7),4月和10月碳汇量达到2个峰值,分别为-200.02 t·月-1与-190.71 t·月-1;1月和7月碳汇量最低,处于波谷,分别为50.18 t·月-1与9.16 t·月-1。陈云飞等[30]对临安太湖源雷竹Phyllostachys violascens林的二氧化碳通量季节变化进行了研究,表明雷竹生态系统净生态系统交换量(NEE)呈现双峰特征,冬季1月为最大排放源,春、夏、秋季为碳汇,但夏季相对较弱。这与校园自然生态系统固碳特征一致。李霞等[6]针对北京海淀公园二氧化碳通量进行了研究,认为其二氧化碳通量季节变化呈现单峰特征,3-10月为碳汇,11月-翌年2月为碳源,于8-9月碳吸收量达到峰值。这可能与南北气候差异有关,临安夏季7-8月酷暑高温,二氧化碳通量于10:00左右达到峰值后下降[30]。而北京海淀公园夏季二氧化碳通量于12:00左右达到峰值[6]。在碳汇构成成分上,草地碳汇量高于林地碳汇量,这可能由两方面原因造成:一是草地月均光合、呼吸速率均比林地大(图 8图 9);二是校园内草地面积是林地面积的3.5倍(表 1)。③单位绿化面积碳吸收量月平均值为-187 g·m-2,最小值出现在1月,为80 g·m-2;最大值出现在4月,为-319 g·m-2

      图  5  校园草地碳收支季节变化

      Figure 5.  Seasonal changes of carbon budget of grassland in campus

      图  6  校园林地碳收支季节变化

      Figure 6.  Seasonal changes of carbon budget of woodland in campus

      图  7  校园自然生态系统碳汇季节变化特征

      Figure 7.  Monthly carbon budget of natural ecosystem in campus

      图  8  植物光合速率季节变化

      Figure 8.  Monthly average plant photosynthetic rate of natural ecosystem in campus

      图  9  植物与土壤呼吸速率季节变化

      Figure 9.  Monthly average vergetable and soil respiration rate of natural ecosystem in campus

    • 由于人类活动密集,浙江农林大学东湖校区全年表现为碳源,年累计碳排放量为4 147.48 t·a-1。植被的碳吸收量相对校园人类活动所造成的碳排放量较小,不能弥补其排碳量,碳平衡主要受到校园活动人口及燃料燃烧的影响形成与人类活动系统碳排放特征季节性变化。在2月和8月,校园碳排放强度较弱;1月,碳排放强度最强,碳排放量达到649.04 t·月-1图 10)。

      图  10  校园碳收支季节变化

      Figure 10.  Monthly carbon budget of whole campus

    • 如果要建设低碳校园,则需要提出低碳校园的标准。这个标准实质上就是校园碳源与碳汇的基准,它是对校园碳排放与碳吸收的量化。由于研究对象为已建成使用的校园,不考虑建筑碳源,在明确校园规模与人均生命活动碳排放基准的前提下,低碳校园建设的影响因素主要为绿地碳汇、人均生活用能与交通用能[1]

      基于上述影响因素,假设低碳校园情景的指标取值如下:①校园绿地率提升至60%,单位绿化面积碳汇能力提升20%;②人均生活用能碳排放量削减20%,日均交通碳排放削减20%。在系统动力学模型中,任何相关变量的变化都会带来整个系统相应变化,为低碳校园政策效果评估提供了直观工具。按照上述指标取值调整模型中的校园绿化面积、人均生活用能、交通通行量、平均交通碳排放量等变量后运行模型,得到低碳校园情境下碳收支情况如表 4

      表 4  低碳情景与原状态碳收支结果比较

      Table 4.  Comparation between low-carbon scenario and original state

      状态 基础数据 计算结果
      绿地率/% 单位绿化面积碳汇/(g·m-2 人均生活用能碳排放/(g·人-1 ·d-1 日均交通碳排放/(kg·d-1 绿地碳汇/(t·a-1 生活用能碳排放(包括锅炉)/(t·a-1 交通碳排放/(t·a-1 碳收支/(t·a-1 人均碳排放量/(g·人-1 ·d-1
      原状态 49 187 535.69 203.8 -1 408 5 555.48 73.35 4 147.48 439.27
      低碳情景 60 224 428.55 163.0 -2 066 3 138.95 58.68 2 929.44 310.27

      在低碳校园情境下,校园碳排放量下降至2 929.44 t·a-1,比原状态减少了29%以上。这说明校园的低碳建设仍有余地,但建设空间已有局限。这可能是由于东湖校区人口密集、现状植被覆盖率相对较高、植被体系已相对成熟等原因造成的。

    • 浙江农林大学东湖校区全年碳收支表现为碳源,年累计碳排放约为4 147.48 t·a-1,其中,人类活动系统产生碳排放约为5 555.45 t·a-1,自然生态系统碳吸收量约为-1 407.96 t·a-1。校园全年碳收支季节变化规律与人类活动子系统变化规律相似:寒暑假骤减,冬、夏季排放量较春、秋季高的特征。校园整体碳收支最大值出现在1月,为649.04 t·月-1;最小值出现在2月和8月,约为22.00 t·月-1。在校园人类活动子系统中,生活办公耗能是最主要的碳排放来源,贡献率达到55%以上。人均碳排放量平均为24 kg·人-1·月-1。在自然生态子系统中,林地年碳汇量为258.14 t·a-1,草地年碳汇量为1 149.82 t·a-1。碳汇呈现双峰的季节变化特征:春、秋两季较高,冬季达到最低。在4月达到-200.02 t·月-1,7月降低至9.16 t·月-1。单位绿化面积碳吸收量月平均值为-187.00 g·m-2。在完成校园碳收支测算的基础上,模拟低碳校园情景,从绿地碳汇能力、人均生活用能与交通用能3个方面控制校园碳收支,达到比原始状态削减29%以上碳排放量的结果。

      本研究虽利用Vensim模型得到浙江农林大学东湖校区全年碳收支的动态结果,但未能很好地反映碳汇系统与碳源系统间相互作用与反馈机制,未能研究确定碳汇能力最大化的植物配置模式。今后有必要进一步进行该方面的研究。另外,本研究虽给出了低碳校园建设与评价基准的参考,但未深入考究该标准的操作可行性和区域的适应性,这是今后低碳城市建设有待研究和解决的问题。

参考文献 (30)

目录

    /

    返回文章
    返回