全 文 ：
摇 摇 摇 摇 摇 生 态 学 报
摇 摇 摇 摇 摇 摇 摇 (SHENGTAI XUEBAO)
摇 摇 第 32 卷 第 24 期摇 摇 2012 年 12 月摇 (半月刊)
目摇 摇 次
从文献计量角度分析中国生物多样性研究现状 刘爱原,郭玉清,李世颖,等 (7635)……………………………
CO2 浓度升高和模拟氮沉降对青川箭竹叶营养质量的影响 周先容,汪建华,张摇 红,等 (7644)………………
陕西中部黄土高原地区空气花粉组成及其与气候因子的关系———以洛川县下黑木沟村为例
吕素青,李月从,许清海,等 (7654)
…………………
……………………………………………………………………………
长三角地区 1995—2007 年生态资产时空变化 徐昔保,陈摇 爽,杨桂山 (7667)…………………………………
基于智能体模型的青岛市林地生态格局评价与优化 傅摇 强,毛摇 锋,王天青,等 (7676)………………………
青藏高原高寒草地生态系统服务功能的互作机制 刘兴元,龙瑞军,尚占环 (7688)……………………………
北京城市绿地的蒸腾降温功能及其经济价值评估 张摇 彪,高吉喜,谢高地,等 (7698)…………………………
武汉市造纸行业资源代谢分析 施晓清,李笑诺,赵吝加,等 (7706)………………………………………………
丽江市家庭能耗碳排放特征及影响因素 王丹寅,唐明方,任摇 引,等 (7716)……………………………………
基于分布式水文模型和福利成本法的生态补偿空间选择研究 宋晓谕,刘玉卿,邓晓红,等 (7722)……………
设施塑料大棚风洞试验及风压分布规律 杨再强,张摇 波,薛晓萍,等 (7730)……………………………………
湖南珍稀濒危植物———珙桐种群数量动态 刘海洋,金晓玲,沈守云,等 (7738)…………………………………
云南岩陀及其近缘种质资源群体表型多样性 李萍萍,孟衡玲,陈军文,等 (7747)………………………………
沙埋和种子大小对柠条锦鸡儿种子萌发、出苗和幼苗生长的影响 杨慧玲,梁振雷,朱选伟,等 (7757)………
栗山天牛天敌花绒寄甲在栎林中的种群保持机制 杨忠岐,唐艳龙,姜摇 静,等 (7764)…………………………
基于相邻木排列关系的混交度研究 娄明华,汤孟平,仇建习,等 (7774)…………………………………………
三种回归分析方法在 Hyperion影像 LAI反演中的比较 孙摇 华,鞠洪波,张怀清,等 (7781)…………………
红松和蒙古栎种子萌发及幼苗生长对升温与降水综合作用的响应 赵摇 娟,宋摇 媛,孙摇 涛,等 (7791)………
新疆杨边材贮存水分对单株液流通量的影响 党宏忠,李摇 卫,张友焱,等 (7801)………………………………
火干扰对小兴安岭毛赤杨沼泽温室气体排放动态影响及其影响因素 顾摇 韩,牟长城,张博文 (7808)………
不同潮汐和盐度下红树植物幼苗秋茄的化学计量特征 刘滨尔,廖宝文,方展强 (7818)………………………
腾格里沙漠东南缘沙质草地灌丛化对地表径流及氮流失的影响 李小军,高永平 (7828)………………………
西双版纳人工雨林群落结构及其林下降雨侵蚀力特征 邓摇 云,唐炎林 ,曹摇 敏,等 (7836)…………………
西南高山地区净生态系统生产力时空动态 庞摇 瑞,顾峰雪,张远东,等 (7844)…………………………………
南北样带温带区栎属树种种子化学组成与气候因子的关系 李东胜,史作民,刘世荣,等 (7857)………………
模拟酸雨对龙眼叶片 PS域反应中心和自由基代谢的影响 李永裕,潘腾飞,余摇 东,等 (7866)………………
沈阳市城郊表层土壤有机污染评价 崔摇 健,都基众,马宏伟,等 (7874)…………………………………………
降雨对旱作春玉米农田土壤呼吸动态的影响 高摇 翔,郝卫平,顾峰雪,等 (7883)………………………………
冬季作物种植对双季稻根系酶活性及形态指标的影响 于天一,逄焕成,任天志,等 (7894)……………………
施氮量对小麦 /玉米带田土壤水分及硝态氮的影响 杨蕊菊,柴守玺,马忠明 (7905)……………………………
微山湖鸟类多样性特征及其影响因子 杨月伟,李久恩 (7913)……………………………………………………
新疆北部棉区作物景观多样性对棉铃虫种群的影响 吕昭智,潘卫林,张摇 鑫,等 (7925)………………………
杭州西湖北里湖沉积物氮磷内源静态释放的季节变化及通量估算 刘静静,董春颖,宋英琦,等 (7932)………
基于实码遗传算法的湖泊水质模型参数优化 郭摇 静,陈求稳,张晓晴,等 (7940)………………………………
气候环境因子和捕捞压力对南海北部带鱼渔获量变动的影响 王跃中,孙典荣,陈作志,等 (7948)……………
象山港南沙岛不同养殖类型沉积物酸可挥发性硫化物的时空分布 颜婷茹,焦海峰,毛玉泽,等 (7958)………
专论与综述
提高植物抗寒性的机理研究进展 徐呈祥 (7966)…………………………………………………………………
植被对多年冻土的影响研究进展 常晓丽,金会军,王永平,等 (7981)……………………………………………
凋落物分解主场效应及其土壤生物驱动 査同刚,张志强,孙摇 阁,等 (7991)……………………………………
街尘与城市降雨径流污染的关系综述 赵洪涛,李叙勇,尹澄清 (8001)…………………………………………
期刊基本参数:CN 11鄄2031 / Q*1981*m*16*374*zh*P* ￥ 70郾 00*1510*40*
室室室室室室室室室室室室室室
2012鄄12
封面图说: 永兴岛海滩植被———永兴岛是中国西沙群岛的主岛,也是西沙群岛及南海诸岛中最大的岛屿。 国务院 2012 年 6 月
批准设立的地级三沙市,管辖西沙群岛、中沙群岛、南沙群岛的岛礁及其海域,三沙市人民政府就驻西沙永兴岛。 永
兴岛岛上自然植被密布,野生植物有 148 种,占西沙野生植物总数的 89% ,主要树种有草海桐(羊角树)、麻枫桐、野
枇杷、海棠树和椰树等。 其中草海桐也称为羊角树,是多年生常绿亚灌木植物,它们总是喜欢倚在珊瑚礁岸或是与
其他滨海植物聚生于海岸沙滩,为典型的滨海植物。
彩图提供: 陈建伟教授摇 北京林业大学摇 E鄄mail: cites. chenjw@ 163. com
第 32 卷第 24 期
2012 年 12 月
生 态 学 报
ACTA ECOLOGICA SINICA
Vol. 32,No. 24
Dec. ,2012
http: / / www. ecologica. cn
基金项目:国家科技支撑计划课题(2012BAD22B01);林业公益性行业科研专项(201104006,200804001)资助
收稿日期:2012鄄04鄄09; 摇 摇 修订日期:2012鄄10鄄10
*通讯作者 Corresponding author. E鄄mail: zyd@ caf. ac. cn
DOI: 10. 5846 / stxb201204090496
庞瑞,顾峰雪,张远东,侯振宏,刘世荣.西南高山地区净生态系统生产力时空动态.生态学报,2012,32(24):7844鄄7856.
Pang R,Gu F X,Zhang Y D, Hou Z H, Liu S R. Temporal鄄spatial variations of net ecosystem productivity in alpine area of southwestern China. Acta
Ecologica Sinica,2012,32(24):7844鄄7856.
西南高山地区净生态系统生产力时空动态
庞摇 瑞1,顾峰雪2,张远东1,*,侯振宏1,刘世荣1
(1. 中国林业科学研究院森林生态环境与保护研究所,国家林业局森林生态环境重点实验室,北京摇 100091;
2. 中国农业科学院农业环境与可持续发展研究所,农业部旱作节水农业重点实验室,北京摇 100081)
摘要:西南高山地区生态系统类型丰富、地形复杂, 是响应全球气候变化的重点区域,对全球气候变化具有重要的指示作用。
应用生态系统模型(Carbon Exchange between Vegetation, Soil, and the Atmosphere,CEVSA)模型估算了 1954—2010 年西南高山
地区净生态系统生产力(NEP)的时空变化,分析了其对气候变化的响应。 结果表明:(1) 1954—2010 年西南高山地区 NEP 平
均为 29. 7 g C·m-2·a-1,其中低海拔地区常绿针叶林和常绿阔叶林 NEP 较高,而高海拔地区的草地覆盖类型 NEP 较低。 (2)西
南高山地区 NEP 总量的变动范围为-8. 36—29. 4Tg C / a,平均每年吸收碳 15. 4Tg C;NEP 年际下降趋势显著(P<0. 05),平均每
年减少 0. 187Tg C,下降显著的区域占研究地区总面积的 35. 2% (P<0. 05),其中草地(-0. 526 g C·m-2·a-2, P<0. 01)和常绿针
叶林(-0. 691 g C·m-2·a-2, P<0. 01)下降趋势极为显著。 (3) 年 NEP 总量的年际变化与年平均温度呈负相关( r = -0. 454,P<
0. 01),与年降水量呈正相关( r=0. 708,P<0. 01),与温度显著负相关的区域占 60. 3% (P<0. 05),与降水显著正相关的区域占
52. 1% (P<0. 05),其中草地和常绿针叶林均与温度极显著负相关( r= -0. 603, P<0. 01; r= -0. 485, P<0. 01),而与降水量极显
著正相关( r=0. 554, P<0. 01; r = 0. 749, P<0. 01)。 (4)西南高山地区是明显的碳汇区,但是由于土壤异养呼吸(HR,
heterotrophic respiration)的增长速度大于净初级生产力(NPP, net primary production)的增长速度,最近 20a 有部分地区开始由
碳汇转为碳源。
关键词:NEP;时空动态;气候变化;西南高山地区
Temporal鄄spatial variations of net ecosystem productivity in alpine area of
southwestern China
PANG Rui1,GU Fengxue2,ZHANG Yuandong1,*, HOU Zhenhong1, LIU Shirong1
1 Key Laboratory of Forest Ecology and Environment, State Forestry Administration, Institute of Forest Ecology, Environment and Protection, Chinese
Academy of Forestry, Beijing 100091,China
2 Key Laboratory of Dryland Agriculture, Ministry of Agriculture, Institute of Environment and Sustainable Development in Agriculture, Chinese Academy of
Agricultural Sciences, Beijing 100081, China
Abstract: There are various types of ecosystems and complex landform in alpine area of southwestern China, which make it
very sensitive to global climate change. Quantifying the spatial and temporal variation of terrestrial ecosystem carbon sink in
this place requires a detailed understanding of carbon exchange between vegetation, soil, and the atmosphere. A process鄄
based biogeochemical model CEVSA ( Carbon Exchange between Vegetation, Soil, and the Atmosphere) was used to
estimate temporal and spatial variations of net ecosystem production (NEP) in alpine area of southwestern China from 1954
to 2010. In this study, the set of independent environmental variables to run the model consisted of climate ( ten days mean
value of temperature, precipitation, relative humidity, cloudiness), atmospheric CO2 concentration, soil and vegetation
http: / / www. ecologica. cn
types at a spatial resolution of 0. 1毅. We first ran CEVSA model with average climate data during the period 1954—2010
until an ecological equilibrium was reached, then made dynamic simulations from 1954 to 2010 with transient changes in
climate at a time鄄step of ten days. To achieve the results, various kinds of computer software were applied, such as
ANUSPLIN4. 1, Arcgis9. 3, SPSS18. 0, Fortran 90. The results showed: (1)The mean value of NEP was 29. 7 g C·m-2·
a-1 during the period 1954—2010 in alpine area of southwestern China, with higher NEP in the evergreen needle鄄leaved and
evergreen broadleaved tree cover at low altitude and lower NEP in herbaceous cover at high altitude of study area. (2)
Annual total NEP ranged from -8. 36 to 29. 4 Tg C / a with a mean value of 15. 4 Tg C / a. NEP showed a statistically
significant decreasing trend with a reduction rate of 0. 187 Tg C / a (P <0. 05) during the period 1954—2010 and the
significant decreasing area of NEP accounted for 35. 2% (P < 0. 05). The decreasing trend was highly significant in
herbaceous cover(-0. 526 g C·m-2·a-2, P<0. 01) and evergreen needle鄄leaved tree cover ( -0. 691 g C·m-2·a-2, P<
0郾 01). (3)Annual total NEP was significantly negatively correlated with annual mean temperature (AMT) ( r = -0. 454,
P<0. 01) and positively correlated with annual total precipitation (ATP)( r = 0. 708, P<0. 01). The negative correlation
(P<0. 05) of annual NEP and AMT existed in 60. 3% of the study region while the positive correlation (P<0. 05) of
annual NEP with ATP existed in 52. 1% of the study region, respectively. Herbaceous cover and evergreen needle鄄leaved
tree cover as the two main vegetation types in study area were both correlated negatively with AMT( r = -0. 603, P<0. 01;
r= -0. 485, P<0. 01) and positively with ATP ( r = 0. 554, P < 0. 01; r = 0. 749, P < 0. 01), reaching an extremely
remarkable level. (4)Alpine area of southwestern China acted as a carbon sink, however, a small fraction of terrestrial
ecosystems here shifted from carbon sinks into carbon sources in recent 20 years, for that the NPP (net primary production)
growth rate was less than that of the soil HR ( heterotrophic respiration ) . These estimates indicate that the carbon
sequestration capacity in alpine area of southwestern China was weakened by the ongoing climate change.
Key Words: NEP;temporal鄄spatial variations;climate change;alpine area of southwestern China
陆地生态系统在全球碳循环中起着极其重要的作用[1鄄2],是全球大气 CO2浓度年际变异的主要控制因
素[3鄄4]。 陆地生态系统碳循环及其对气候变化的响应早期研究多集中于植被净初级生产力(NPP) [5鄄8],即单
位时间单位面积上植被所能生产的有机干物质总量,是光合作用所吸收的碳和自养呼吸所释放的碳之间的差
值[9],它是碳循环的原动力[10]。 后来,对陆地碳循环的研究逐渐开始关注净生态系统生产力(NEP) [11鄄13],即
生态系统碳的净吸收或净排放,代表着陆地生态系统和大气生态系统的净 CO2 交换量。 NEP 这一概念最早
是Woodwell等[14]在分析陆地生物圈源、汇的问题时提出的,用来表示较大尺度生态系统中碳的净贮存量,是
生态系统植被净初级生产力和土壤异氧呼吸(HR)之差[15]。 NEP 直接定性定量描述了陆地生态系统碳源汇
的能力,且通过 CO2 对气候系统产生直接影响[16]。
由于《京都议定书》将陆地生态系统碳吸收纳入到控制温室气体排放的框架之中,陆地碳汇空间分布及
其变化受到了空前的重视[17鄄18]。 虽然国内外对陆地生态系统碳收支进行了大量研究,普遍认为陆地生态系
统起着大气 CO2 碳汇的作用,但是对碳汇及其空间格局变化的估算仍存在很大的不确定性[19鄄22],所以有必要
对特定区域的碳吸收的分布及其变化进行研究。 我国位于北半球中高纬度地区,是全球陆地碳汇研究的重点
地区。 Cao利用 CEVSA模型得出 1981—2000 年间全国 NEP 年均值为 0. 07 Pg C / a (-0. 32—0. 25 Pg C / a),
碳吸收主要发生在东北地区、藏东南地区及黄淮海平原等地区,且气候变化降低了中国陆地生态系统的碳吸
收能力[18]。 陶波应用 CEVSA模型进一步研究表明东北平原中部和黄土高原的 NEP 有明显的下降趋势,而
在华中大部分地区有增加趋势[23]。 方精云利用资源清查资料和卫星遥感数据得出 1981—2000 年中国陆地
生态系统的总碳汇(植被和土壤)为 0. 136—0. 176 Pg C / a[24]。 Mu 运用 Biome鄄BGC 模型模拟结果表明中国
区域 1961—2000 年间 NEP 均值为 0. 16 Pg C / a (0. 07—0. 23 Pg C / a),平均每年下降-0. 02 T g C[25]。 中国
地区 NEP 总体为碳汇区,但碳汇分布格局及变化趋势存在空间差异。
5487摇 24 期 摇 摇 摇 庞瑞摇 等:西南高山地区净生态系统生产力时空动态 摇
http: / / www. ecologica. cn
西南高山地区位于“世界第三极冶青藏高原东南边缘,高山峡谷地貌,其显著的气候垂直分异、特殊的地
质地貌、不同发育阶段的土壤和植被,致使植物物候、生长、生殖与生理对温度和 CO2 浓度的变化都给予极强
的反馈,是响应全球气候变化的重点区域,对全球气候变化具有重要的调节和指示作用[26鄄29]。 关于西南高山
地区碳循环对气候变化的响应以往的研究主要集中在主要森林类型碳储量、碳密度、生物量及净初级生产力
的研究[30鄄35],对亚高山森林土壤碳排放也有少量研究[36],但缺少对于西南高山地区 NEP 时空动态的深入研
究,以及气候对 NEP 影响的分析。 本文应用高分辨率的气候数据库和生态系统过程模型估算了 1954—2010
年期间西南高山地区 NEP,分析了该地区碳源 /汇的分布格局及其对气候变化的响应。
1摇 研究地区和研究方法
1. 1摇 研究地区概况
本文采用吴中伦《中国森林》的林区划分方法,西南高山地区为我国第二大天然林区,地理范围北部西接
青藏高原,北侧及东侧与黄土高原相接壤,东侧南段与四川盆地分界,南界东段与云贵高原相接,西南段为国
境线,西南高山地区作为一个玉立玉的级区,向下又分为 4 个域级区,从北到南依次是洮河白龙江云杉冷杉林
区、岷江冷杉林区,大渡河雅砻江金沙江云杉冷杉林区,藏东南云杉冷杉林区[37](图 1)。 本研究中,统计的西
南高山地区总面积为 5. 19伊105km2,总的地势是西北高、向东南倾斜,主要为高山峡谷地貌,年平均气温约为
5—10益,年降水量约为 600—1000mm,且各地段差别较大。 辽阔的面积复杂的地理地貌孕育了本区复杂的
植被类型,主要森林类型为云杉林及冷杉林,低海拔处有阔叶林,位置偏北降水量较少的洮河林区和白龙江林
区主要是落叶阔叶林,大渡河向西到金沙江流域及雅鲁藏布江流域海拔较低处有常绿阔叶林分布。 林下灌丛
丰富,且森林线以上有高海拔灌丛带,灌木带以上为高山草甸或高山草原。
图 1摇 西南高山地区地理位置
Fig. 1摇 The location of alpine area of southwestern China
6487 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
http: / / www. ecologica. cn
1. 2摇 数据来源
本研究中,模型所需的气象数据(10d 平均温度,降水量,相对湿度和云量)来自国家气象信息中心,包括
1954—2010年全国约 756个台站的旬观测数据,然后应用 ANUSPLIN软件进行内插获得分辨率 0. 1毅的栅格数据
作为模型气象输入数据。 大气 CO2 浓度资料来自于美国夏威夷 Mauna Loa 观测所(http: / / cdiac. esd. ornl. gov /
ftp / maunaloa鄄co2 / maunaloa. co2) 。 土壤类型和质地资料取自 1颐100 万土壤类型图和第二次土壤普查数据,并使
之匹配于气候数据的空间分辨率。 植被数据来源于 http: / / www鄄gem. jrc. it / glc2000,是基于 2000 年的 spot 影像
制作的 1km分辨率的土地覆盖分类数据,在此基础上进行了插值重采样,提取了 0. 1毅分辨率的植被数据。
1. 3摇 模型介绍
CEVSA(Carbon Exchange between Vegetation, Soil, and Atmosphere)模型是一个基于生理生态过程模拟植
物鄄土壤鄄大气系统能量交换和水碳氮耦合循环的生物地球化学循环模型[38鄄40]。 CEVSA 模型通过 3 个子模型
描述这些过程及其与各种环境变量的相互作用:淤估算植物鄄土壤鄄大气之间水热交换、土壤含水量和气孔传
导等过程的生物物理子模型;于计算植物光合作用、呼吸作用、氮吸收速率、叶面积以及碳氮在植物各器官之
间分配、积累、周转和凋落物产生的植物生理生长子模型;盂估算土壤有机质分解与转化和有机氮矿化等过程
的土壤碳氮转化子模型[41]。 CEVSA模型已经在全球和中国区域应用于研究陆地生态系统碳循环对气候变
化的响应[18,23,38,40鄄44]。
1. 4摇 模型运行与验证
在模拟过程中,首先应用 1954—2010 平均气候数据将模型运行至生态系统平衡态,即各个状态变量如植
被、土壤碳贮量以及土壤含水量等年际变化量小于 0. 1% ,且 NPP、LT(凋落物量)分别与 HR 相等,然后再用
1954—2010 每 10d气候数据驱动模型来进行动态模拟,并反复运行模型以消除假定的生态系统状态变量的
初始值(即平衡态假设)对模拟结果的影响,之后即得到最终模拟结果。
摇 图 2摇 研究区域 1954—2010 年均净生态系统生产力(NEP)空间
分布
Fig. 2 摇 Spatial distribution of mean annual net ecosystem
production(NEP)during the period 1954—2010 in alpine area of
southwestern China
CEVSA模型应用的生物学和生态学原理、计算方程及参数均取自大量的实验室试验和野外试验观测。
该模型曾应用叶片或植株生理试验、样点初级生产力观测和遥感反演数据进行了验证,模型估算的 NPP、叶
面积指数及植被和土壤中的碳贮量与实地调查和测定所得的实际值有很好的一致性[38鄄39]。 从 20 世纪 90 年
代中期以后,大型环境控制试验和涡度相关技术的应用为模型的检验测试提供了连续、长期的生态系统尺度
碳通量观测数据[45鄄46]。 顾峰雪等[47]应用涡度相关通量观测数据对 CEVSA 模型的验证和评价表明,CEVSA
模型较好地模拟了不同类型生态系统水碳通量的动态特征及其对气候变化的响应。
1. 5摇 数据分析
相关性分析采用 Pearson相关系数法,对 CEVSA模
型估算出的 1954—2010 年净生态系统生产力与分别温
度、降水数据逐点计算相关系数,得到 NEP 与温度、降
水的相关系数图,然后设置 0. 05、0. 01 两个显著性水平
对相关系数进行分类。 NEP 年际变化趋势的分析采用
线性趋势倾向率方法[48鄄49],在研究区域内逐点进行线
性回归,最小二乘法拟合得到的一元线性方程斜率即为
趋势倾向率,并用相关系数检验法来确定变化是否显
著[50鄄53],同样设置 0. 05、0. 01 两个显著性水平对相关
系数进行分类,查相关系数临界值表得 r0. 05 = 0. 261,
r0. 01 =0. 339。
2摇 结果
2. 1摇 NEP 空间分布
图 2 为 1954—2010 时间段西南高山地区年平均
7487摇 24 期 摇 摇 摇 庞瑞摇 等:西南高山地区净生态系统生产力时空动态 摇
http: / / www. ecologica. cn
NEP 分布图,NEP 平均为 29. 7 g C·m-2·a-1。 西南高山地区为明显的碳汇区,NEP 为正值的区域占研究区域
总面积的 95. 6% 。 NEP 较高的区域为主要为洮河、白龙江及岷江流域的常绿针叶林、落叶阔叶林及定期水
淹灌丛 /草地地带(>45 g C·m-2·a-1)和大渡河流域的常绿阔叶林覆盖地区(>30 g C·m-2·a-1),以及雅鲁藏布
江流域大拐弯以南的高山深谷区(>45 g C·m-2·a-1),而负值主要零星分布在大渡河雅砻江金沙江林区及藏东
南地区的草地覆盖地区。
西南高山地区 NEP 最高为落叶阔叶林地带(57. 5 g C·m-2·a-1),但面积只占林区面积的 1. 04% ,其次为
定期水淹灌丛 /草地(53. 7 g C·m-2·a-1),最低为草地覆盖类型(17. 5 g C·m-2·a-1),且草地为西南高山地区面
积最大的植被类型(35. 5% )。 常绿针叶林(31. 2% )和常绿阔叶林(15. 3% )面积仅次于草地地区,NEP 分别
为 40郾 5 g C·m-2·a-1、41. 6 g C·m-2·a-1。
格点的 NEP 多年平均值与温度的多年平均值正相关( r = 0. 368,P<0. 01)而与降水量多年平均值和海拔
均负相关,相关系数分别为-0. 142(P<0. 01)、-0. 330(P<0. 01)。
表 1摇 西南亚高山地区不同植被类型 57 年 NEP均值、NEP趋势倾向率及其与温度、降水相关系数
Table 1摇 57鄄year annual mean NEP, NEP tendency rate and its correlation coefficient with mean annual temperature and annual precipitation
of different vegetation types during the period 1954—2010
植被类型
Vegetation type
面积百分比
Area ratio
/ %
单位面积 NEP
NEP in unit area
/ (g C·m-2·a-1)
NEP 总量
Total NEP
/ (Tg C / a)
趋势倾向率
Tendency rate
/ (g C·m-2·a-2)
NEP 与温度
相关系数
NEP 降水量
相关系数
常绿阔叶林 15. 3 41. 6 3. 31 0. 185 -0. 163 0. 515**
落叶阔叶林(郁闭) 1. 04 57. 5 0. 312 -1. 86* -0. 013 0. 454**
常绿针叶林 31. 2 40. 5 6. 56 -0. 691** -0. 485** 0. 749**
林木与其它天然植被嵌合区 1. 95 24. 2 0. 245 0. 670** -0. 005 0. 267*
常绿灌丛 3. 30 19. 7 0. 339 0. 762* -0. 077 0. 311*
落叶灌丛 0. 02 41. 7 0. 004 4. 27** 0. 281 -0. 058
草地 35. 5 17. 5 3. 23 -0. 526** -0. 603** 0. 554**
稀疏草地 /稀疏灌丛 0. 14 20. 0 0. 014 -0. 425* -0. 387** 0. 245
定期水淹灌丛 /草地 0. 99 53. 7 0. 276 0. 135 -0. 055 0. 113
耕作和管理区 4. 99 31. 8 0. 824 -0. 358 -0. 130 0. 659**
农田、林区与其它天然植被
嵌合区 2. 01 32. 3 0. 338 0. 587*
-0. 027 0. 167
裸地 0. 21 0 0
水体 1. 62 0 0
雪和冰 1. 74 0 0
摇 摇 常绿阔叶林 Tree Cover, broadleaved, evergreen;落叶阔叶林(郁闭)Tree Cover, broadleaved, deciduous, closed;常绿针叶林 Tree Cover, needle鄄
leaved, evergreen;林木与其它天然植被嵌合区 Mosaic: Tree cover / Other natural vegetation;常绿灌丛 Shrub Cover, closed鄄open, evergreen;落叶灌
丛 Shrub Cover, closed鄄open, deciduous;草地 Herbaceous Cover, closed鄄open;稀疏草地 /稀疏灌丛 Sparse Herbaceous or sparse Shrub Cover;定期水淹
灌丛 /草地 Regularly flooded Shrub and / or Herbaceous Cover;耕作和管理区 Cultivated and managed areas;农田、林区与其它天然植被嵌合区 Mosaic:
Cropland / Tree Cover / Other natural vegetation;裸地 Bare Areas;水体 Water Bodies;雪和冰 Snow and Ice;NEP 与温度相关系数,Correlation
coefficient between NEP and mean annual temperature;NEP 降水量相关系数 Correlation coefficient between NEP and annual precipitation; *, P<0. 05;
**, P<0. 01
2. 2摇 NEP 年际变化趋势
西南高山地区 NEP 总量各年平均为 15. 4Tg C / a(图 3),其中碳吸收量最大为常绿针叶林(6. 56Tg C / a),
其次为常绿阔叶林(3. 31Tg C / a)和草地(3. 23Tg C / a)植被类型。 1954—2010 时间段,年 NEP 总量最大值为
2010 年的 29. 4Tg C / a,最小值为 2009 年的-8. 36Tg C / a,下降趋势显著(P<0. 05),平均每年减少 0. 187Tg C。
研究时段内,西南高山地区温度上升趋势显著(P<0. 05),趋势倾向率为 0. 08益 / 10a,降水量下降趋势显
著(P<0. 01),趋势倾向率为 18. 3mm / 10a。 NEP 的年总量与年平均温度呈负相关( r= -0. 454,P<0. 01),与年
8487 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
http: / / www. ecologica. cn
摇 图 3摇 西南高山地区 1954—2010 净生态系统生产力总量(NEP)
时间动态
Fig. 3 摇 Interannual variations in annual total net ecosystem
production (NEP)during the period 1954—2010 in alpine area of
southwestern China
降水量呈正相关( r = 0. 708,P<0. 01)。 降水量下降,有
76. 9%的年份 NEP 下降,降水量升高,则 87. 1%的年份
NEP 也是升高的。 因此,西南高山地区 NEP 年际变化
主要受降水量影响。 NEP 最高的年份 2010 年,虽然该
年温度比 57a 平均水平高出 7. 78% ,但降水量充足
(867mm,1. 5% ),致使碳吸收量较高。 2009 年为 NEP
最低的年份,温度虽然和 2010 相差不大,但是降水量为
研究时段内最低值(717mm,-16. 1% )(图 4)。
2. 3摇 NEP 年际变化趋势空间分布
如图 5 所示,西南亚高山地区趋势倾向率以负值为
主,占研究区域总面积的 61. 0% ,显著降低的面积占
35. 2% (P<0. 05),主要分布在研究地区北部的洮河、白
龙江流域及岷江流域的常绿针叶林和落叶阔叶林(郁
闭)和川西大渡河、雅砻江、金沙江流域北部的草地地
带。 NEP 显著增加的区域仅占 12. 1% ,且主要发生在
藏东南雅鲁藏布江流域,主要植被覆盖为常绿灌丛、常
图 4摇 研究区域 1954—2010 年平均温度和年降水量距平变化
Fig. 4摇 Changes in mean annual temperature departure and annual precipitation departure during the period 1954—2010 in alpine area of
southwestern China
绿阔叶林及林木与其它天然植被嵌合区。
研究时段内,西南高山地区 NEP 下降极为显著的为常绿针叶林和草地覆盖地带(P<0. 01),平均每年分
别减少 0. 691 g C·m-2·a-1、0. 526 g C·m-2·a-1,而林木与其它天然植被嵌合区上升趋势最为显著(P<0. 01),年
增加量为 0. 670 g C·m-2·a- 1(表 1)。
近百年来全球主要有两个增温时段,第 1 个时段为 20 世纪 10—40 年代,第 2 个时段为从 70 年代后期开
始到现在;中国地区温度 50—60 年代出现了的相对冷期,从 1970 年代开始增温, 20 世纪 90 年代后期以后明
显偏暖[35,54]。 西南高山地区温度变化与全国范围温度变化趋势基本一致,1954—1970 温度降低趋势显著
(P<0. 01),1971—1990 为增加趋势,但不显著,1991—2010 增温趋势显著(P<0. 01),且 3 个时段变化趋势差
9487摇 24 期 摇 摇 摇 庞瑞摇 等:西南高山地区净生态系统生产力时空动态 摇
http: / / www. ecologica. cn
图 5摇 西南亚高山地区 1954—2010 年净生态系统生产力(NEP)变化趋势及相关系数法显著性检验
Fig. 5摇 Trend of annual net ecosystem production(NEP) changes, and its significance test with the correlation coefficient method during the
period 1954—2010 in alpine area of southwestern China
异显著(P<0. 01)。 图 6 为 3 个不同温度变化时段 NEP 空间分布图,结果表明 1954—1970 时间段,西南高山
地区除藏东南的雅鲁藏布江流域外均为明显的碳汇区,且 NEP 西北高于东南,总量为 19. 5Tg C / a。 1971—
1990 时间段,NEP 总量下降至 15. 6Tg C / a,且西北比东南下降趋势明显,然而雅鲁藏布江流域 NEP 有所增
加,由碳源转变成碳汇;1991—2010 时段 NEP 总量继续下降至 11. 8Tg C / a,这个时期碳源区占整个林区的
25. 4% ,主要为大渡河、雅砻江、金沙江及藏东南地区的草地覆盖地区,而藏东南雅鲁藏布江流域成为明显的
碳汇区。
图 6摇 研究区域 1954—1970、1971—1990、1991—2010 年净生态系统生产力(NEP)平均空间分布
Fig. 6摇 Spatial distribution of annual net ecosystem production(NEP)during the period 1954—1970,1971—1990,1991—2000 and 2001—
2010 in alpine area of Southwest China
0587 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
http: / / www. ecologica. cn
西南高山地区 NEP 与温度相关系数为负值的区域占 85. 4% ,通过显著性检验的面积占 60. 3% (P<
0郾 05),主要位于草地覆盖地区及雅鲁藏布江流域南部林区,且西北部较东南部与温度负相关性明显(图 7)。
NEP 与降水量主要为正相关,正值区占林区面积的量 76. 8% ,通过显著性检验为 52. 1% (P<0. 05),且主要发
生在洮河、白龙江及岷江流域南部的常绿针叶林和落叶阔叶林,大渡河流域的常绿阔叶林,雅砻江金沙江流域
的常绿针叶林及雅鲁藏布江流域南部(图 7)。
图 7摇 研究区域 1954—2010 年净生态系统生产力(NEP)与年平均温度和年降水量相关系数
Fig. 7摇 Correlation coefficient between annual net ecosystem production( NEP) and mean annual temperature and annual precipitation
during the period 1954—2010in alpine area of southwestern China
研究区域内,NEP 与温度极显著负相关的植被类型为常绿针叶林、草地和稀疏草地 /稀疏灌丛,相关系数
分别为-0. 485(P<0. 01)、-0. 603(P<0. 01)、-0. 387(P<0. 01),而与降水量极显著正相关的植被类型有常绿
阔叶林( r=0. 515,P<0. 01)、落叶阔叶林(郁闭)( r = 0. 454,P<0. 01)、常绿针叶林( r = 0. 749,P<0. 01)、草地
( r=0. 554,P<0. 01)及耕作和管理区( r=0. 659,P<0. 01)(表 1)。
4摇 讨论
西南高山地区平均每年碳吸收量为 15. 4Tg C / a,根据全国尺度碳吸收的估算数据[18,25,55], 占全国碳吸收
的 9%—22% 。 西南高山地区 NEP 均值为 29. 7 g C·m-2·a-1,略高于基于 Biome鄄BGC 模型和 BEPS 模型对东
北地区的估算结果 25. 8 g C·m-2·a-1 [56]。 虽然应用的土地覆盖分类数据不同,但植被类型 NEP 大小顺序基本
与应用 Biome鄄BGC模型对黄淮海地区研究结果基本一致,落叶阔叶林>常绿针叶林>草地[57]。 林区内落叶阔
叶林 NEP 较高,是因为其 NPP 较高(549 g C·m-2·a-1),而常绿针叶林和草地的 NPP 虽然相差不大,分别为
442 g C·m-2·a-1、430 g C·m-2·a-1,但是由于草地 HR(412 g C·m-2·a-1)与其 NPP 基本持平,所以导致草地 NEP
较低。 西南高山地区主要植被覆盖类型为草地、常绿针叶林,分别占林区面积的 35. 5% 、31. 2% ,碳吸收量则
分别占林区碳吸收总量的 20. 9% 、42. 4% 。 两种植被类型 NEP 均与温度极显著负相关,与降水量极显著正相
关,由于西南高山地区气候朝暖干化趋势发展,所以草地、常绿针叶林 NEP 均呈极显著下降趋势(P<0. 01),
导致大部分地区 NEP 趋势倾向率为负值。
NEP 时空动态与气候及其所致的 NPP 和 HR密切相关。
西南高山地区 NEP 的空间分布与 NPP( r=0. 253, P<0. 01)、HR( r=0. 130, P<0. 01)都为正相关,这与中
国南部地区 NEP 空间分布与 NPP 正相关、HR 负相关的结果不同[23],导致这一结果的原因源于西南高山地
区特殊的地理环境。 西南高山地区总的地势是西北高,向东南倾斜,温度降水自东南向西北逐次降低,相应的
植被类型则由落叶、常绿阔叶林和常绿针叶林过渡到高山、亚高山草甸,这一特殊的地理地貌,导致 NPP 与
HR均由东南向西北递减(图 8),且有很高的相关性,相关系数达 0. 992(P<0. 01)。
NPP 空间分布上与温度、降水均显著正相关,相关系数分别为 0. 750(P<0. 01),0. 532(P<0. 01),而与海
1587摇 24 期 摇 摇 摇 庞瑞摇 等:西南高山地区净生态系统生产力时空动态 摇
http: / / www. ecologica. cn
图 8摇 研究区域 1954—2010 年均净初级生产力(NPP)和土壤异氧呼吸(HR)空间分布图
Fig. 8摇 Spatial distribution of mean annual net primary production and soil heterotrophic respiration during the period 1954—2010 in
alpine area of southwestern China
拔高度呈显著负相关( r= -0. 706,P<0. 01),HR的空间分布与 NPP 相似,与温度、降水显著正相关,与海拔负
相关,相关系数分别为 0. 721(P<0. 01)、0. 564(P<0. 01)、-0. 681(P<0. 01)。 在 NPP、HR 的综合作用下,
NEP 的空间分布情况则比较复杂,同温度正相关( r=0. 368,P<0. 01),与降水量虽然负相关,但相关系数较低
( r= -0. 142, P<0. 01),而与海拔则负相关,相关系数为-0. 330(P<0. 01)。 相对于 NPP、HR由东南向西北递
减的规律性变化,NEP 的空间分布受植被类型的影响较为明显,高值区主要为低海拔处的常绿针叶林和常绿
阔叶林分布地区,而高海拔地区的草地覆盖类型 NEP 则较低。
图 9摇 1954—2010 年西南高山地区净初级生产力(NPP)和土壤异养呼吸(HR)时间动态及净初级生产力(NPP)和土壤异养呼吸(HR)对温
度敏感程度
Fig. 9摇 Interannual variations in net primary production ( NPP) and heterotrophic respiration ( HR) sensitivity to temperature of net
primary production(NPP)and heterotrophic respiration (HR) in alpine area of southwestern China
NEP 的年际变化主要由 HR决定,与 HR显著负相关( r= -0. 553,P<0. 01),而与 NPP 则不相关,HR上升
有 73. 5%的年份 NEP 是下降的,HR下降则有 65. 2%的年份 NEP 是上升的。 在年际尺度上 NPP、HR均增加
趋势显著(P<0. 01),平均每年分别增加 0. 523Tg C、0. 710Tg C(图 9),这与全国尺度上 NPP、HR的增长趋势
一致[18,55]。 西南高山地区 NPP、HR年际变化均与温度正相关(P<0. 01),相关系数分别为 0. 443(P<0. 01)、
0. 636(P<0. 01),其中 HR 较 NPP 对温度上升更为敏感(图 9)。 尽管有研究认为 NPP 对水分胁迫很敏
感[58],然而该地区 NPP 年际变化与降水量并不显著,这可能是因为西南高山地区降水丰富,水分并不是限制
植物生长的主要因子,这与应用 CASA 模型对降水量同样充沛的长江流域研究结果一致[59]。 土壤呼吸一般
来说是随降水量增加而升高的,但在降水量丰富的地区,增加的降水量可能会导致呼吸速率的降低[60],在雨
2587 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
http: / / www. ecologica. cn
水充沛的西南高山地区 HR年际变化与降水量呈负相关( r = -0. 369,P<0. 01)。 西南高山地区温度上升导致
HR的增加趋势大于 NPP 的增加趋势,而降水量的下降对 NPP 影响不大,但对 HR 有促进作用,二者综合作
用导致呈 NEP 呈下降趋势。 这与在 CEVSA、Biome鄄BGC 模型在全国尺度及 FORCCHN 模型在东北地区的研
究结果相一致[15,18,25]。 西南高山地区在 1991—2010 时段,已有 25. 4%的地区有碳汇转变成碳源,若是未来
气候继续朝暖干化发展,西南高山林区的碳汇能力可能会进一步减弱,将会有更多的地区由碳汇转变为碳源。
作为 NPP 和土壤呼吸的差值,NEP 估算结果的不确定性要远远大于 NPP [16],从模型本身来看,CEVSA
仅仅考虑了气候及 CO2 浓度变化对碳收支的影响,而并未考虑土地利用 / 土地覆被变化(LUCC)这一重要驱
动因子,但是对于 LUCC数据难以获得且过程机制尚不清楚,如在模型中考虑 LUCC 的因素,估算结果将会有
很大的变化[23]。 另外,西南高山地区复杂的地理环境导致温度和降水量地区间差异较大,虽然气象数据插值
时,考虑了海拔因素,但气候要素内插精度仍有待进一步提高。 复杂的地理环境孕育了西南高山地区复杂的
土地覆盖类型,植被类型不仅水平差异较大,且具有明显的垂直分布带,研究过程中虽然采用了 GLC 土地覆
盖分类数据,把研究地区土地覆盖类型细分到 14 种,但是鉴于研究地区复杂的植被类型分布,土地覆盖分类
数据有待于进一步完善。
5摇 结论
西南高山地区位于青藏高原东南边缘,独特地质地貌、水热条件分布,使本地区 NEP 的时空动态有其自
身特点。
(1)在空间分布上,NEP 与海拔负相关且与土地覆盖类型关系密切,NEP 高值区主要为低海拔处的常绿
针叶林和常绿阔叶林分布地区,而高海拔地区的草地覆盖类型 NEP 较低。
(2)西南高山地区大部分区域趋势倾向率为负值,其中主要植被类型草地和常绿针叶林在研究时间段内
下降趋势极为显著,且均与温度极显著负相关而与降水极显著正相关。
(3)西南高山地区是我国重要的碳汇区,碳汇量占已有文献估算量的 9%—22% 。。 在年际尺度上,由于
气候变化导致西南高山地区 HR增加速度大于 NPP 增加速度, NEP 在 1954—2010 时段下降显著。 若是未来
气候继续朝暖干化发展,西南高山林区的碳汇能力可能会进一步减弱,将会有更多的地区由碳汇转变为碳源。
References:
[ 1 ]摇 Grant R F, Black T A, Gaumont鄄Guay D, Klujn N, Barr A G, Morgenstern K, Nesic Z. Net ecosystem productivity of boreal aspen forests under
drought and climate change: mathematical modelling with ecosys. Agricultural and forest meteorology, 2006, 140(1 / 4): 152鄄170.
[ 2 ] 摇 BARR A G, Black T A, Hogg E H, Griffis T J, Morgenstern K, Kljun N, Theede A, Nesic Z. Climatic controls on the carbon and water balances
of a boreal aspen forest, 1994—2003. Global Change Biology, 2007, 13(3): 561鄄576.
[ 3 ] 摇 Kljun N, Black T A, Griffis T J, Barr A G, Gaumont鄄Guay D, Morgenstern K, McCaughey J H, Nesic Z. Response of net ecosystem productivity
of three boreal forest stands to drought. Ecosystems, 2006, 9(7): 1128鄄144.
[ 4 ] 摇 Bergeron O, Margolis H A, Black T A, Coursolle C, DUNN A L, BARR A G, WOFSY S C. Comparison of carbon dioxide fluxes over three boreal
black spruce forests in Canada. Global Change Biology, 2007, 13(1): 89鄄107.
[ 5 ] 摇 Raich J W, Rastetter E B, Melillo J M, Kicklighter D W, Steudler P A, Peterson B J, Grace A L, Moore Iii B, Vorosmarty C J. Potential net
primary productivity in South America: application of a global model. Ecological Applications, 1991, 1(4): 399鄄429.
[ 6 ] 摇 Melillo J M, McGuire A D, Kicklighter D W, Moore B, Vorosmarty C J, Schloss A L. Global climate change and terrestrial net primary
production. Nature, 1993, 363(6426): 234鄄240.
[ 7 ] 摇 Nemani R R, Keeling C D, Hashimoto H, Jolly W M, Piper S C, Tucker C J, Myneni R B, Running S W. Climate鄄driven increases in global
terrestrial net primary production from 1982 to 1999. Science, 2003, 300(5625): 1560鄄1563.
[ 8 ] 摇 Cramer W, Kicklighter D W, Bondeau A, Iii B M, Churkina G, Nemry B, Ruimy A, Schloss A L. Comparing global models of terrestrial net
primary productivity (NPP): overview and key results. Global Change Biology, 1999, 5(S1): 1鄄15.
[ 9 ] 摇 Chen B, Wang S Q, Liu R G, Song T. Study on Modeling and Spatial Pattern of Net Primary Production in China忆s Terrestrial Ecosystem.
Resources Science, 2007,29(6): 45鄄53.
[10] 摇 Field C B, Behrenfeld M J, Randerson J T, Falkowski P. Primary production of the biosphere: integrating terrestrial and oceanic components.
3587摇 24 期 摇 摇 摇 庞瑞摇 等:西南高山地区净生态系统生产力时空动态 摇
http: / / www. ecologica. cn
Science, 1998, 281(5374): 237鄄240.
[11] 摇 Falge E, Baldocchi D, Olson R, Anthoni P, Aubinet M, Bernhofer C, Burba G, Ceulemans R, Clement R, Dolman H. Gap filling strategies for
defensible annual sums of net ecosystem exchange. Agricultural and forest meteorology, 2001, 107(1): 43鄄69.
[12] 摇 Randerson J T, Chapin III F S, Harden J W, Neff J C, Harmon M E. Net ecosystem production: a comprehensive measure of net carbon
accumulation by ecosystems. Ecological applications, 2002, 12(4): 937鄄947.
[13] 摇 O忆Connell K E B, Gower S T, Norman J M. Net ecosystem production of two contrasting boreal black spruce forest communities. Ecosystems,
2003, 6(3): 248鄄260.
[14] 摇 Woodwell G M, Whitaker R H, Reiners W A, Likens G E, Delwich C C, Botkin D B. Biota and the world carbon budget. Science;(United
States), 1978, 199(4325):144鄄146.
[15] 摇 Zhao J F, Yan X D, Jia G S. Changes in carbon budget of Northeast China forest ecosystems under future climatic scenario. Chinese Journal of
Ecology, 2009,28(5): 781鄄787.
[16] 摇 Chang S L, Yang H X, Ge J P. Advance and questions in net ecosystem production. Journal of Beijing Normal University (Natural Science),
2005,41(5): 517鄄521.
[17] 摇 Schimel D S. Terrestrial ecosystems and the carbon cycle. Global Change Biology, 1995, 1(1): 77鄄91.
[18] 摇 Cao M K, Prince S D, Li K, Tao B, Small J, Shao X M. Response of terrestrial carbon uptake to climate interannual variability in China. Global
Change Biology, 2003, 9(4): 536鄄546.
[19] 摇 Tans P P, Fung I Y, Takahashi T. Observational contrains on the global atmospheric CO2 budget. Science, 1990, 247(4949): 1431鄄1438.
[20] 摇 Dai A, Fung I Y. Can Climate Variability Contribute to the" Missing" CO2 Sink? Global Biogeochemical Cycles, 1993, 7(3): 599鄄609.
[21] 摇 Goodale C L, Apps M J, Birdsey R A, Field C B, Heath L S, Houghton R A, Jenkins J C, Kohlmaier G H, Kurz W, Liu S. Forest carbon sinks
in the Northern Hemisphere. Ecological Applications, 2002, 12(3): 891鄄899.
[22] 摇 Schimel D S, House J I, Hibbard K A, Bousquet P, Ciais P, Peylin P, Braswell B H, Apps M J, Baker D, Bondeau A. Recent patterns and
mechanisms of carbon exchange by terrestrial ecosystems. Macmillan Magazines Ltd, 2001,414(6860): 169鄄172.
[23] 摇 Tao B, Cao M K, Li K R, Gu F X, Ji J J, Huang M, Zhang L M. Spatial patterns and changes of terrestrial net ecosystem productivity during
1981—2000 in China. Science in China (Series D: Earth Sciences), 2006, 36 (12): 1131鄄1139.
[24] 摇 Fang J Y, Guo Z D, Piao S L, Chen A P. Terrestrial vegetation carbon sinks in China from 1981 to 2000. Science in China (Series D ), 2007,37
(6): 804鄄812.
[25] 摇 Mu Q Z, Zhao M S, Running S W, Liu M L, Tian H Q. Contribution of increasing CO2 and climate change to the carbon cycle in China忆s
ecosystems. Journal of Geophysical Research, 2008, 113:G01018, doi:10. 1029 / 2006JG000316.
[26] 摇 Zhang Y D, Liu S R, Ma J M, Shi Z M, Liu X L. Wood land hydrological effects of birch forests in subalpine region of western Sichuan, China.
Acta Ecologica Sinica, 2005, 25(11): 2939鄄2946.
[27] 摇 Gu X P, Huang M, Jing J J, Wu Z P. The influence of climate change on vegetation net primary productivity in Southwestern China during Recent
20 Years Period. Journal of Natural Resources, 2007,22(2): 251鄄259.
[28] 摇 Tang X. Estimation of forest carbon storage and its spatial distribution in Sichuan province[D]. Yaan: Sichuan agricultural university, 2007.
[29] 摇 Liu B, Yang W Q, Wu F Z. Advances in the subalpine forest ecosystem processes. Acta Ecologica Sinica, 2010, 30(16): 4476鄄4483.
[30] 摇 Cheng G W, Luo J. The carbon accumulation and dissipation features of subalpine woodland in Mt. Gongga. Acta Geographica Sinica, 2003, 58
(2):179鄄185.
[31] 摇 Ma M D, Jiang H, Luo C D, Liu Y J. Preliminary study of carbon density, net production and carbon stock in natural spruce forests of northwest
subalpine Sichuan, china. Chinese Journal of Plant Ecology, 2007, 31(2): 305鄄312.
[32] 摇 Li Q. The characteristics of soil carbon and nitrogen and nutrient stock of different forest types in northwest of Sichuan[ D]. Yaan:Sichuan
agricultural university, 2008.
[33] 摇 Zhang G B, Liu S R, Zhang Y D, Miao N, Wang H. Dynamics of aboveground biomass of subalpine old鄄growth forest in the upper Minjiang River.
Acta Ecologica Sinica, 2008, 28(7): 3176鄄3184.
[34] 摇 Liu Q, Zhao C Z, Lin B, Cheng X Y, Yin H J, Chen J S, Zhou F F, Jiang F Y, Pan X L, Yin C Y Soil C Storages and their Responses to climate
warming in subalpine spruce plantation, Western Sichuan, China,2010. Mode shift of economy advance and independent innovation鄄China
Association for Science and Technology Annual Conference (Vol. I), 2010.
[35] 摇 Chen L X, Zhou X J, Li W L. Characteristics of the climate change and its formation mechanism in China in last 80 years. Acta Meteorologica
Sinica, 2004, 62(5): 634鄄646.
[36] 摇 Lu X Y, Huo C F, Fang J H, Cheng G W. Effects of climate change on soil greenhouse gas emissions in subalpine dark coniferous forest. Ecology
and Environmental Sciences, 2009, 18(6): 2194鄄2199.
4587 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
http: / / www. ecologica. cn
[37]摇 Wu Z L. China Forest. Beijing: China Forest Press, 1997: 486鄄498.
[38] 摇 Cao M K, Woodward F I. Net primary and ecosystem production and carbon stocks of terrestrial ecosystems and their responses to climate change.
Global Change Biology, 1998, 4(2): 185鄄198.
[39] 摇 Woodward F I, Smith T M, Emanuel W R. A global land primary productivity and phytogeography model. Global biogeochemical cycles, 1995, 9
(4): 471鄄490.
[40] 摇 Cao M K, Woodward F I. Dynamic responses of terrestrial ecosystem carbon cycling to global climate change. Nature, 1998, 393 (6682):
249鄄252.
[41] 摇 Li K R, Wang S Q, Cao M K. Carbon storage of vegetation and soil in china. Science in China (Series D), 2003, 33(1):72鄄80.
[42] 摇 Cao M K, Zhang Q, Shugart H H. Dynamic responses of African ecosystem carbon cycling to climate change. Climate Research, 2001, 17(2):
183鄄193.
[43] 摇 Cao M K, Tao B, Li K R, Shao X M, Prience S D. Interannual variation in terrestrial ecosystem carbon fluxes in China from 1981 to 1998. Acta
Botanica Sinica, 2003, 45(5): 552鄄560.
[44] 摇 Tao B, Li K R, Shao X M, Cao M K. Temporal and spatial pattern of net primary production of terrestrial ecosystems in China. Acta Geographica
Sinica, 2003,58(3): 372鄄380.
[45] 摇 Running S W, Baldocchi D D, Turner D P, Gower S T, Bakwin P S, Hibbard K A. A global terrestrial monitoring network integrating tower
fluxes, flask sampling, ecosystem modeling and EOS satellite data. Remote Sens Environ, 1999, 70(1): 108鄄127.
[46] 摇 Baldocchi D D, Wilson K B. Modeling CO2 and water vapor exchange of a temperate broadleaved forest across hourly to decadal time scales. Ecol
Model, 2001, 142(1): 155鄄184.
[47] 摇 Gu F X, Cao M K, Yu G R, Tao B, Wen X F,Liu Y F, Zhang L M. Modeling carbon exchange in different forest ecosystems by CEVSA model:
comparison with eddy covariance measurements. Advances in Earth Science, 2007,22(3): 313鄄321.
[48] 摇 Cai F, Zhang S J, Yu G R, Zhu Q L, Liu X A. Research of spatial鄄temporal evolvement characters of mean air temperature in China in recent 50
years based on spatialization technique. Plateau Meteorology, 2006, 25(6): 1168鄄1175.
[49] 摇 Yao Y B, Yang J H, Wang R Y, Lu D R. Change feature of net primary productivity of natural vegetation and its impact factors over the source
region of the Yangtze River in recent 50 years. Ecology and Environmental Sciences, 2010, 19(11): 2521鄄2528.
[50] 摇 He X Q, Liu W Q. Analysis of Applied Regression. Beijing: Chinese People忆s University Press, 2001: 18鄄40.
[51] 摇 Yang Y L. Temporal and spatial analysis of NDVI variation and correlation analysis between NDVI and temperature and rainfall in the natural grass
land of Qinghai Province [D]. Nanjing: Nanjing University of Information Science &Technology, 2008.
[52] 摇 Du J Q, Shu J M, Zhang L B, Guo Y. Responses of vegetation to climate change in the headwaters of China忆s Yellow River Basin based on zoning
of dry and wet climate. Chinese Journal of Plant Ecology, 2011, 35(11): 1192鄄1201
[53]摇 Piao S L, Fang J Y. Terrestrial net primary production and its spatio鄄temporal patterns in Qinghai鄄Xizang Plateau, China during 1982—1999.
Journal of Natural Resources, 2002, 17(3): 373鄄380.
[54] 摇 Zuo H C, Lu S H, Hu Y Q. Variations trend of yearly mean air temperature and precipitation in China in the last 50 years. Plateau Meteorology,
2004, 23(2): 238鄄244.
[55] 摇 Sun G D. Simulation of potential vegetation distribution and estimation of carbon flux in china from 1981 to 1998 with LPJ Dynamic Global
Vegetation Model. Climatic and Environmental Research, 2009, 14(4): 341鄄351.
[56] 摇 Wang J B. Chinese terrestrial net ecosystem productive model applied remote sensing [D]. Hangzhou: Zhejiang University, 2004.
[57] 摇 Hu B, Sun R, Chen Y J, Feng L C, Sun Liang. Estimation of the net Ecosystem Productivity in Huang鄄Huai Hai Region combining with Biome鄄
BGC Model and remote sensing data. Journal of Natural Resources, 2011, 26 (12):2062鄄2070.
[58] 摇 Churkina G, Running S W, Schloss A L. Comparing global models of terrestrial net primary productivity ( NPP): The importance of water
availability. Global Change Biology, 1999. 5(S1): 46鄄55.
[59] 摇 Ke J H, Piao S L, Fang J Y. NPP and its spatio鄄temporal patterns in the yangtze river watershed. Acta Phytoecologica Sinica, 2003, 27(6):
764鄄770.
[60] 摇 Wang W, Fang J Y. Soil respiration and human effects on global grasslands. Global and Planetary Change, 2009, 67(1 / 2): 20鄄28.
参考文献:
[ 9 ]摇 陈斌, 王绍强, 刘荣高, 宋婷. 中国陆地生态系统 NPP 模拟及空间格局分析. 资源科学, 2007,29(6): 45鄄53.
[15] 摇 赵俊芳, 延晓冬, 贾根锁. 未来气候情景下中国东北森林生态系统碳收支变化. 生态学杂志, 2009,28(5): 781鄄787.
[16] 摇 常顺利, 杨洪晓, 葛剑平. 净生态系统生产力研究进展与问题. 北京师范大学学报(自然科学版), 2005,41(5): 517鄄521.
[23] 摇 陶波,曹明奎 李克让 顾峰雪 季劲钧. 1981—2000 年中国陆地净生态系统生产力空间格局及其变化. 中国科学(D辑:地球科学), 2006,
5587摇 24 期 摇 摇 摇 庞瑞摇 等:西南高山地区净生态系统生产力时空动态 摇
http: / / www. ecologica. cn
36 (12): 131鄄1139.
[24] 摇 方精云, 郭兆迪, 朴世龙, 陈安平. 1981—2000 年中国陆地植被碳汇的估算. 中国科学(D辑:地球科学), 2007,37(6): 804鄄812.
[26] 摇 张远东, 刘世荣, 马姜明, 史作民, 刘兴良. 川西亚高山桦木林的林地水文效应. 生态学报, 2005, 25(11): 2939鄄2946.
[27] 摇 谷晓平, 黄玫, 季劲钧, 吴战平. 近 20 年气候变化对西南地区植被净初级生产力的影响. 自然资源学报, 2007,22(2): 251鄄259.
[28] 摇 唐宵. 四川森林植被碳储量估算及其空间分布特征[D]. 雅安:四川农业大学, 2007.
[29] 摇 刘彬, 杨万勤, 吴福忠. 亚高山森林生态系统过程研究进展.生态学报, 2010, 30(16): 4476鄄4483.
[30] 摇 程根伟, 罗辑. 贡嘎山亚高山林地碳的积累与耗散特征. 地理学报, 2003, 58(002): 179鄄185.
[31] 摇 马明东, 江洪, 罗承德, 刘跃建. 四川西北部亚高山云杉天然林生态系统碳密度, 净生产量和碳贮量的初步研究. 植物生态学报, 2007,
31(2): 305鄄312.
[32] 摇 李强. 四川西北部森林植被类型土壤碳氮特征和养分库研究[D]. 雅安:四川农业大学, 2008.
[33] 摇 张国斌, 刘世荣, 张远东, 缪宁, 王晖. 岷江上游亚高山林区老龄林地上生物量动态变化. 生态学报, 2008, 28(7): 3176鄄3184.
[34] 摇 刘庆, 赵春章, 林波, 程新颖, 尹华军, 陈劲松, 周非飞, 姜发艳, 潘新丽, 尹春英. 川西亚高山人工云杉林土壤碳库及其对气候变暖的
响应. 经济发展方式转变与自主创新———第十二届中国科学技术协会年会 (第一卷), 2010.
[35] 摇 陈隆勋, 周秀骥, 李维亮, 罗云峰, 朱文琴. 中国近 80 年来气候变化特征及其形成机制. 气象学报, 2004, 62(005): 634鄄646.
[36] 摇 鲁旭阳, 霍常富, 范继辉, 程根伟. 气候变化对亚高山暗针叶林土壤温室气体排放的影响. 生态环境学报, 2009, 18(6): 2194鄄2199.
[37] 摇 吴中伦. 中国森林. 北京: 中国林业出版社, 1997:486鄄498.
[41] 摇 李克让, 王绍强, 曹明奎. 中国植被和土壤碳贮量. 中国科学 D辑, 2003, 33(1): 72鄄80.
[44] 摇 陶波, 李克让, 邵雪梅, 曹明奎. 中国陆地净初级生产力时空特征模拟. 地理学报, 2003,58(3): 372鄄380.
[47] 摇 顾峰雪, 曹明奎, 于贵瑞, 陶波, 温学发, 刘允芬, 张雷明. 典型森林生态系统碳交换的机理模拟及其与观测的比较研究. 地球科学进
展, 2007,22(3): 313鄄321.
[48] 摇 蔡福, 张淑杰, 于贵瑞, 祝青林, 刘新安. 基于空间化技术对中国近 50 年平均气温时空演变特征的研究. 高原气象, 2006, 25(6):
1168鄄1175.
[49] 摇 姚玉璧, 杨金虎, 王润元, 陆登荣. 50 年长江源区域植被净初级生产力及其影响因素变化特征. 生态环境学报, 2010, 19 (11):
2521鄄2528.
[50] 摇 何晓群, 刘文卿. 应用回归分析. 北京: 中国人民大学出版社, 2001: 18鄄40.
[51] 摇 杨英莲. 青海省天然草地 NDVI 的时空化与气温和降水的关系分析[D].南京:南京信息工程大学, 2008.
[52] 摇 杜加强, 舒俭民, 张林波, 郭杨. 黄河上游不同干湿气候区植被对气候变化的响应. 植物生态学报, 2011, 35(11): 1192鄄1201.
[53] 摇 朴世龙, 方精云. 1982—1999 年青藏高原植被净第一性生产力及其时空变化. 自然资源学报, 2002, 17(3): 373鄄380.
[54] 摇 左洪超, 吕世华, 胡隐樵. 中国近 50 年气温及降水量的变化趋势分析. 高原气象, 2004, 23(2): 238鄄244.
[55] 摇 孙国栋. LPJ 模型对 1981—1998 年中国区域潜在植被分布和碳通量的模拟. 气候与环境研究, 2009, 14(4): 341鄄351.
[56] 摇 王军邦. 中国陆地净生态系统生产力遥感模型研究 [D]. 杭州: 浙江大学, 2004.
[57] 摇 胡波, 孙睿, 陈永俊, 冯丽超, 孙亮. 遥感数据结合 Biome鄄BGC 模型估算黄淮海地区生态系统生产力. 自然资源学报, 2011, 26(12):
2062鄄2070.
[59] 摇 柯金虎, 朴世龙, 方精云, 长江流域植被净第一性生产力及其时空格局研究. 植物生态学报, 2003, 27(6): 764鄄770.
6587 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
ACTA ECOLOGICA SINICA Vol. 32,No. 24 December,2012(Semimonthly)
CONTENTS
A bibliometric study of biodiversity research in China LIU Aiyuan, GUO Yuqing, LI Shiying,et al (7635)…………………………
Effects of elevated CO2 and nitrogen deposition on leaf nutrient quality of Fargesia rufa Yi
ZHOU Xianrong, WANG Jianhua, ZHANG Hong,et al (7644)
……………………………………………
……………………………………………………………………
Airborne pollen assemblages and their relationships with climate factors in the central Shaanxi Province of the Loess Plateau:
a case in Xiaheimugou, Luochuan County L譈 Suqing, LI Yuecong, XU Qinghai,et al (7654)…………………………………
Spatial and temporal change in ecological assets in the Yangtze River Delta of China 1995—2007
XU Xibao, CHEN Shuang, YANG Guishan (7667)
……………………………………
…………………………………………………………………………………
Evaluation and optimization of woodland ecological patterns for Qingdao based on the agent鄄based model
FU Qiang, MAO Feng, WANG Tianqing,et al (7676)
……………………………
……………………………………………………………………………
Interactive mechanism of service function of alpine rangeland ecosystems in Qinghai鄄Tibetan Plateau
LIU Xingyuan, LONG Ruijun, SHANG Zhanhuan (7688)
…………………………………
…………………………………………………………………………
Preliminary evaluation of air temperature reduction of urban green spaces in Beijing
ZHANG Biao, GAO Jixi, XIE Gaodi,et al (7698)
…………………………………………………
…………………………………………………………………………………
Resources metabolism analysis for the pulp and paper industry in Wuhan, China
SHI Xiaoqing,LI Xiaonuo,ZHAO Linjia,et al (7706)
………………………………………………………
………………………………………………………………………………
The characteristics and influential factors of direct carbon emissions from residential energy consumption: a case study of Lijiang
City, China WANG Danyin, TANG Mingfang, REN Yin, et al (7716)…………………………………………………………
Spatial targeting of payments for ecosystem services Based on SWAT Model and cost鄄benefit analysis
SONG Xiaoyu,LIU Yuqing,DENG Xiaohong,et al (7722)
…………………………………
…………………………………………………………………………
The wind tunnel test of plastic greenhouse and its surface wind pressure patterns
YANG Zaiqiang,ZHANG Bo,XUE Xiaoping,et al (7730)
………………………………………………………
…………………………………………………………………………
Population quantitative characteristics and dynamics of rare and endangered plant Davidia involucrata in Hunan Province
LIU Haiyang, JIN Xiaoling, SHEN Shouyun,et al (7738)
……………
…………………………………………………………………………
Phenotypic diversity in populations of germplasm resources of Rodgersia sambucifolia and related species
LI Pingping, MENG Hengling, CHEN Junwen,et al (7747)
……………………………
………………………………………………………………………
Effects of sand burial and seed size on seed germination, seedling emergence and growth of Caragana korshinskii Kom. (Fabaceae)
YANG Huiling, LIANG Zhenlei,ZHU Xuanwei,et al (7757)
…
………………………………………………………………………
Population鄄keeping mechanism of the parasitoid Dastarcus helophoroides (Coleoptera: Bothrideridae) of Massicus raddei
(Coleoptera: Cerambycidae) in oak forest YANG Zhongqi, TANG Yanlong, JIANG Jing,et al (7764)…………………………
Study of mingling based on neighborhood spatial permutation LOU Minghua, TANG Mengping, QIU Jianxi,et al (7774)……………
Comparison of three regression analysis methods for application to LAI inversion using Hyperion data
SUN Hua, JU Hongbo, ZHANG Huaiqing,et al (7781)
…………………………………
……………………………………………………………………………
Response of seed germination and seedling growth of Pinus koraiensis and Quercus mongolica to comprehensive action of warming
and precipitation ZHAO Juan, SONG Yuan, SUN Tao, et al (7791)……………………………………………………………
Impacts of water stored in sapwood Populus bolleana on its sap flux DANG Hongzhong, LI Wei,ZHANG Youyan,et al (7801)………
Dynamics of greenhouse gases emission and its impact factors by fire disturbance from Alnus sibirica forested wetland in
Xiaoxing忆an Mountains, Northeast China GU Han,MU Changcheng, ZHANG Bowen (7808)……………………………………
Different tide status and salinity alter stoichiometry characteristics of mangrove Kandelia candel seedlings
LIU Biner, LIAO Baowen, FANG Zhanqiang (7818)
……………………………
………………………………………………………………………………
Effects of shrub encroachment in desert grassland on runoff and the induced nitrogen loss in southeast fringe of Tengger Desert
LI Xiaojun, GAO Yongping (7828)
……
…………………………………………………………………………………………………
Community structure and throughfall erosivity characters of artificial rainforest in Xishuangbanna
DENG Yun, TANG Yanlin , CAO Min, et al (7836)
………………………………………
………………………………………………………………………………
Temporal鄄spatial variations of net ecosystem productivity in alpine area of southwestern China
PANG Rui,GU Fengxue,ZHANG Yuandong, et al (7844)
…………………………………………
…………………………………………………………………………
Relationships between chemical compositions of Quercus species seeds and climatic factors in temperate zone of NSTEC
LI Dongsheng, SHI Zuomin, LIU Shirong, et al (7857)
……………
……………………………………………………………………………
Effects of simulated acid rain stress on the PS域 reaction center and free radical metabolism in leaves of longan
LI Yongyu, PAN Tengfei, YU Dong, et al (7866)
……………………
…………………………………………………………………………………
Assessment of organic pollution for surface soil in Shenyang suburbs CUI Jian,DU Jizhong,MA Hongwei,et al (7874)………………
The impact of rainfall on soil respiration in a rain鄄fed maize cropland GAO Xiang, HAO Weiping, GU Fengxue, et al (7883)………
Effects of winter crops on enzyme activity and morphological characteristics of root in subsequent rice crops
YU Tianyi, PANG Huancheng,REN Tianzhi,et al (7894)
…………………………
…………………………………………………………………………
Dynamic changes of soil moisture and nitrate nitrogen in wheat and maize intercropping field under different nitrogen supply
YANG Ruiju, CHAI Shouxi, MA Zhongming (7905)
…………
………………………………………………………………………………
Characteristics of the bird diversity and the impact factors in Weishan Lake YANG Yuewei, LI Jiuen (7913)………………………
The effect of cropping landscapes on the population dynamics of the cotton bollworm Helicoverpa armigera (Lepidoptera,
Noctuidae) in the northern Xinjiang LU Zhaozhi, PAN Weilin, ZHANG Xin, et al (7925)……………………………………
The seasonal variations of nitrogen and phosphorus release and its fluxes from the sediments of the Beili Lake in the Hangzhou
West Lake LIU Jingjing,DONG Chunying,SONG Yingqi,et al (7932)……………………………………………………………
Optimization of lake model salmo based on real鄄coded genetic algorithm
GUO Jing, CHEN Qiuwen, ZHANG Xiaoqing, et al (7940)
………………………………………………………………
………………………………………………………………………
The influence of climatic environmental factors and fishing pressure on changes of hairtail catches in the northern South China
Sea WANG Yuezhong, SUN Dianrong, CHEN Zuozhi, et al (7948)………………………………………………………………
Seasonal and spatial distribution of acid volatile sulfide in sediment under different mariculture types in Nansha Bay, China
YAN Tingru, JIAO Haifeng, MAO Yuze, et al (7958)
…………
……………………………………………………………………………
Review and Monograph
Research progress on the mechanism of improving plant cold hardiness XU Chengxiang (7966)………………………………………
Influences of vegetation on permafrost: a review CHANG Xiaoli,JIN Huijun,WANG Yongping,et al (7981)…………………………
Home鄄field advantage of litter decomposition and its soil biological driving mechanism: a review
ZHA Tonggang, ZHANG Zhiqiang, SUN Ge, et al (7991)
………………………………………
…………………………………………………………………………
Research progress on the relationship of pollutants between road鄄deposited sediments and its washoff
ZHAO Hongtao, LI Xuyong, YIN Chengqing (8001)
…………………………………
………………………………………………………………………………
8008 摇 生摇 态摇 学摇 报摇 摇 摇 32 卷摇
《生态学报》2013 年征订启事
《生态学报》是中国生态学学会主办的生态学专业性高级学术期刊,创刊于 1981 年。 主要报道生态学研
究原始创新性科研成果,特别欢迎能反映现代生态学发展方向的优秀综述性文章;研究简报;生态学新理论、
新方法、新技术介绍;新书评介和学术、科研动态及开放实验室介绍等。
《生态学报》为半月刊,大 16 开本,300 页,国内定价 90 元 /册,全年定价 2160 元。
国内邮发代号:82鄄7,国外邮发代号:M670
标准刊号:ISSN 1000鄄0933摇 摇 CN 11鄄2031 / Q
全国各地邮局均可订阅,也可直接与编辑部联系购买。 欢迎广大科技工作者、科研单位、高等院校、图书
馆等订阅。
通讯地址: 100085 北京海淀区双清路 18 号摇 电摇 摇 话: (010)62941099; 62843362
E鄄mail: shengtaixuebao@ rcees. ac. cn摇 网摇 摇 址: www. ecologica. cn
摇 摇 编辑部主任摇 孔红梅摇 摇 摇 执行编辑摇 刘天星摇 段摇 靖
生摇 态摇 学摇 报
(SHENGTAI摇 XUEBAO)
(半月刊摇 1981 年 3 月创刊)
第 32 卷摇 第 24 期摇 (2012 年 12 月)
ACTA ECOLOGICA SINICA
摇
(Semimonthly,Started in 1981)
摇
Vol郾 32摇 No郾 24 (December, 2012)
编摇 摇 辑摇 《生态学报》编辑部
地址:北京海淀区双清路 18 号
邮政编码:100085
电话:(010)62941099
www. ecologica. cn
shengtaixuebao@ rcees. ac. cn
主摇 摇 编摇 冯宗炜
主摇 摇 管摇 中国科学技术协会
主摇 摇 办摇 中国生态学学会
中国科学院生态环境研究中心
地址:北京海淀区双清路 18 号
邮政编码:100085
出摇 摇 版摇
摇 摇 摇 摇 摇 地址:北京东黄城根北街 16 号
邮政编码:100717
印摇 摇 刷摇 北京北林印刷厂
发 行摇
地址:东黄城根北街 16 号
邮政编码:100717
电话:(010)64034563
E鄄mail:journal@ cspg. net
订摇 摇 购摇 全国各地邮局
国外发行摇 中国国际图书贸易总公司
地址:北京 399 信箱
邮政编码:100044
广告经营
许 可 证摇 京海工商广字第 8013 号
Edited by摇 Editorial board of
ACTA ECOLOGICA SINICA
Add:18,Shuangqing Street,Haidian,Beijing 100085,China
Tel:(010)62941099
www. ecologica. cn
Shengtaixuebao@ rcees. ac. cn
Editor鄄in鄄chief摇 FENG Zong鄄Wei
Supervised by摇 China Association for Science and Technology
Sponsored by摇 Ecological Society of China
Research Center for Eco鄄environmental Sciences, CAS
Add:18,Shuangqing Street,Haidian,Beijing 100085,China
Published by摇 Science Press
Add:16 Donghuangchenggen North Street,
Beijing摇 100717,China
Printed by摇 Beijing Bei Lin Printing House,
Beijing 100083,China
Distributed by摇 Science Press
Add:16 Donghuangchenggen North
Street,Beijing 100717,China
Tel:(010)64034563
E鄄mail:journal@ cspg. net
Domestic 摇 摇 All Local Post Offices in China
Foreign 摇 摇 China International Book Trading
Corporation
Add:P. O. Box 399 Beijing 100044,China
摇 ISSN 1000鄄0933
CN 11鄄2031 / Q
国内外公开发行 国内邮发代号 82鄄7 国外发行代号 M670 定价 70郾 00 元摇